Gibbs state

.gitignore

@@ -3,6 +3,8 @@ tests/data/temp/
 tests/data/performance_results/
 tests/data/plots/
 *.hdf5
+*.pkl
+*.pdf
 
 # local development
 local_dev/

examples/system-correlations.py

@@ -6,7 +6,7 @@
 
 import numpy as np
 import matplotlib.pyplot as plt
-from IPython import embed
+#from IPython import embed
 
 import oqupy
 import oqupy.operators as op
@@ -33,6 +33,8 @@ def gaussian_shape(t, area = 1.0, tau = 1.0, t_0 = 0.0):
                                 cutoff=omega_cutoff,
                                 cutoff_type='gaussian',
                                 temperature=temperature)
+print('type=',type(correlations))
+print(issubclass(type(correlations), oqupy.correlations.BaseCorrelations))
 bath = oqupy.Bath(op.sigma("z")/2.0, correlations)
 
 

oqupy/__init__.py

@@ -31,6 +31,9 @@
     'Control',
     'Dynamics',
     'FileProcessTensor',
+    'GibbsParameters',
+    'GibbsTempo',
+    'gibbs_tempo_compute',
     'guess_tempo_parameters',
     'helpers',
     'import_process_tensor',
@@ -111,6 +114,9 @@
 
 from oqupy.tempo import Tempo
 from oqupy.tempo import TempoParameters
+from oqupy.tempo import GibbsTempo
+from oqupy.tempo import GibbsParameters
 from oqupy.tempo import MeanFieldTempo
 from oqupy.tempo import guess_tempo_parameters
 from oqupy.tempo import tempo_compute
+from oqupy.tempo import gibbs_tempo_compute