10: does not converge

This commit is contained in:
Valentin Boettcher 2022-09-05 15:55:13 +02:00
parent cc07e05e2f
commit 449b48c8e8

View file

@ -39,7 +39,7 @@ Here we try to reproduce the anti zeno engine from the paper.
ε=1,#.1,
ω_c=1,
ε_couple=.68,
n=60,
n=20,
detune=.5,
ω_0=20,
T_c=5e3,
@ -48,7 +48,7 @@ Here we try to reproduce the anti zeno engine from the paper.
γ=.2,
switch_cycles=1,
therm_initial_state=False,
#ε_init=.001/2,
ε_init=.001/2,
terms=7,
dt=0.01/2,
sp_tol=1e-4,
@ -77,11 +77,13 @@ Here we try to reproduce the anti zeno engine from the paper.
#+RESULTS:
#+begin_src jupyter-python
%matplotlib tk
flow, interaction, int_pow, system, sys_pow = model.all_energies_online_from_cache()
fig, ax = plt.subplots()
print(sys_pow.N)
pu.plot_with_σ(model.t, (sys_pow+int_pow.sum_baths()).integrate(model.t), ax=ax)
pu.plot_with_σ(model.t, (-1 * flow).sum_baths().integrate(model.t) + interaction.sum_baths() + system + 10, ax=ax)
# pu.plot_with_σ(model.t, (sys_pow+int_pow.sum_baths()).integrate(model.t), ax=ax)
pu.plot_with_σ(model.t, (sys_pow+int_pow.sum_baths()), ax=ax)
# pu.plot_with_σ(model.t, (-1 * flow).sum_baths().integrate(model.t) + interaction.sum_baths() + system + 10, ax=ax)
# with aux.get_data(model) as data:
# #pu.plot_with_σ(model.t, model.total_energy_from_power(data), ax=ax)
# pu.plot_with_σ(model.t, model.interaction_energy(data).sum_baths() - interaction.sum_baths(), ax=ax)
@ -93,9 +95,8 @@ Here we try to reproduce the anti zeno engine from the paper.
#+RESULTS:
:RESULTS:
: 10
| hline | <AxesSubplot:> | ((<matplotlib.lines.Line2D at 0x7fec3f75eb80>) <matplotlib.collections.PolyCollection at 0x7fec3f75e9d0>) |
[[file:./.ob-jupyter/72099bcb22465077395cd8acccbecc8f6e768f84.svg]]
: 66590
| hline | <AxesSubplot:> | ((<matplotlib.lines.Line2D at 0x7fe28cb48250>) <matplotlib.collections.PolyCollection at 0x7fe28cb484c0>) |
:END:
@ -104,20 +105,14 @@ Here we try to reproduce the anti zeno engine from the paper.
#+end_src
#+RESULTS:
:RESULTS:
| <Figure | size | 520x320 | with | 1 | Axes> | <AxesSubplot:xlabel= | $\tau$ | > |
[[file:./.ob-jupyter/d1345d964ba6362247d8f34e3ea98e5b5f91a7db.svg]]
:END:
#+begin_src jupyter-python
plot_az_sd_overview(model, params)
#+end_src
#+RESULTS:
:RESULTS:
| <Figure | size | 520x320 | with | 1 | Axes> | <AxesSubplot:xlabel= | $\omega$ | ylabel= | Spectral Density | > |
[[file:./.ob-jupyter/afa180fb667b63a127dbecc86f43980e65eabb6a.svg]]
:END:
#+begin_src jupyter-python
@ -136,8 +131,502 @@ Here we try to reproduce the anti zeno engine from the paper.
#+RESULTS:
:RESULTS:
| <Figure | size | 520x320 | with | 1 | Axes> | <AxesSubplot:> |
[[file:./.ob-jupyter/5ab38efefb7c03476efb3dea22c64715e4c7131f.svg]]
#+begin_example
Process ForkPoolWorker-27:
Process ForkPoolWorker-26:
Process ForkPoolWorker-25:
Traceback (most recent call last):
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hiro_models/model_auxiliary.py", line 148, in get_data
return HIData(path, read_only=read_only, robust=False, **kwargs)
Process ForkPoolWorker-28:
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/core/hierarchy_data.py", line 252, in __init__
with h5py.File(hdf5_name, "r", libver="latest") as h5File:
Traceback (most recent call last):
Traceback (most recent call last):
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/h5py/_hl/files.py", line 533, in __init__
fid = make_fid(name, mode, userblock_size, fapl, fcpl, swmr=swmr)
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/h5py/_hl/files.py", line 226, in make_fid
fid = h5f.open(name, flags, fapl=fapl)
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hiro_models/model_auxiliary.py", line 148, in get_data
return HIData(path, read_only=read_only, robust=False, **kwargs)
Traceback (most recent call last):
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hiro_models/model_auxiliary.py", line 148, in get_data
return HIData(path, read_only=read_only, robust=False, **kwargs)
File "h5py/_objects.pyx", line 54, in h5py._objects.with_phil.wrapper
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/core/hierarchy_data.py", line 252, in __init__
with h5py.File(hdf5_name, "r", libver="latest") as h5File:
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hiro_models/model_auxiliary.py", line 148, in get_data
return HIData(path, read_only=read_only, robust=False, **kwargs)
File "h5py/_objects.pyx", line 55, in h5py._objects.with_phil.wrapper
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/core/hierarchy_data.py", line 252, in __init__
with h5py.File(hdf5_name, "r", libver="latest") as h5File:
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/core/hierarchy_data.py", line 252, in __init__
with h5py.File(hdf5_name, "r", libver="latest") as h5File:
File "h5py/h5f.pyx", line 106, in h5py.h5f.open
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/h5py/_hl/files.py", line 533, in __init__
fid = make_fid(name, mode, userblock_size, fapl, fcpl, swmr=swmr)
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/h5py/_hl/files.py", line 226, in make_fid
fid = h5f.open(name, flags, fapl=fapl)
OSError: Unable to open file (file is already open for write (may use <h5clear file> to clear file consistency flags))
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/h5py/_hl/files.py", line 533, in __init__
fid = make_fid(name, mode, userblock_size, fapl, fcpl, swmr=swmr)
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/h5py/_hl/files.py", line 533, in __init__
fid = make_fid(name, mode, userblock_size, fapl, fcpl, swmr=swmr)
File "h5py/_objects.pyx", line 54, in h5py._objects.with_phil.wrapper
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/h5py/_hl/files.py", line 226, in make_fid
fid = h5f.open(name, flags, fapl=fapl)
During handling of the above exception, another exception occurred:
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/h5py/_hl/files.py", line 226, in make_fid
fid = h5f.open(name, flags, fapl=fapl)
File "h5py/_objects.pyx", line 55, in h5py._objects.with_phil.wrapper
File "h5py/_objects.pyx", line 54, in h5py._objects.with_phil.wrapper
File "h5py/_objects.pyx", line 54, in h5py._objects.with_phil.wrapper
File "h5py/_objects.pyx", line 55, in h5py._objects.with_phil.wrapper
File "h5py/h5f.pyx", line 106, in h5py.h5f.open
File "h5py/h5f.pyx", line 106, in h5py.h5f.open
Traceback (most recent call last):
File "h5py/_objects.pyx", line 55, in h5py._objects.with_phil.wrapper
OSError: Unable to open file (file is already open for write (may use <h5clear file> to clear file consistency flags))
File "h5py/h5f.pyx", line 106, in h5py.h5f.open
OSError: Unable to open file (file is already open for write (may use <h5clear file> to clear file consistency flags))
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
OSError: Unable to open file (file is already open for write (may use <h5clear file> to clear file consistency flags))
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
During handling of the above exception, another exception occurred:
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 125, in worker
result = (True, func(*args, **kwds))
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 125, in worker
result = (True, func(*args, **kwds))
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 125, in worker
result = (True, func(*args, **kwds))
Traceback (most recent call last):
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 48, in mapstar
return list(map(*args))
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 48, in mapstar
return list(map(*args))
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 48, in mapstar
return list(map(*args))
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/stocproc.py", line 1064, in alpha_times_pi
return alpha(tau) * np.pi
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/stocproc.py", line 1064, in alpha_times_pi
return alpha(tau) * np.pi
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/stocproc.py", line 1064, in alpha_times_pi
return alpha(tau) * np.pi
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 377, in __call__
,* zeta_func(
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 377, in __call__
,* zeta_func(
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 18, in zeta_func
return lerch_phi(shift, s, a)
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 125, in worker
result = (True, func(*args, **kwds))
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 18, in zeta_func
return lerch_phi(shift, s, a)
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 377, in __call__
,* zeta_func(
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/zeta_function.py", line 35, in _arb_lerch_phi
if _arb_c_lerch_phi(ctypes.byref(cy), *args, 0):
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/zeta_function.py", line 35, in _arb_lerch_phi
if _arb_c_lerch_phi(ctypes.byref(cy), *args, 0):
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 18, in zeta_func
return lerch_phi(shift, s, a)
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 48, in mapstar
return list(map(*args))
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/zeta_function.py", line 35, in _arb_lerch_phi
if _arb_c_lerch_phi(ctypes.byref(cy), *args, 0):
KeyboardInterrupt
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/stocproc.py", line 1064, in alpha_times_pi
return alpha(tau) * np.pi
KeyboardInterrupt
KeyboardInterrupt
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 377, in __call__
,* zeta_func(
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 18, in zeta_func
return lerch_phi(shift, s, a)
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/zeta_function.py", line 35, in _arb_lerch_phi
if _arb_c_lerch_phi(ctypes.byref(cy), *args, 0):
KeyboardInterrupt
Process ForkPoolWorker-32:
Traceback (most recent call last):
Process ForkPoolWorker-29:
Process ForkPoolWorker-30:
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
Traceback (most recent call last):
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
Traceback (most recent call last):
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
Process ForkPoolWorker-31:
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 125, in worker
result = (True, func(*args, **kwds))
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 125, in worker
result = (True, func(*args, **kwds))
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 48, in mapstar
return list(map(*args))
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 48, in mapstar
return list(map(*args))
Traceback (most recent call last):
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 377, in __call__
,* zeta_func(
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/stocproc.py", line 1064, in alpha_times_pi
return alpha(tau) * np.pi
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/stocproc.py", line 1064, in alpha_times_pi
return alpha(tau) * np.pi
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 125, in worker
result = (True, func(*args, **kwds))
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 377, in __call__
,* zeta_func(
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 377, in __call__
,* zeta_func(
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 18, in zeta_func
return lerch_phi(shift, s, a)
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 18, in zeta_func
return lerch_phi(shift, s, a)
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/zeta_function.py", line 35, in _arb_lerch_phi
if _arb_c_lerch_phi(ctypes.byref(cy), *args, 0):
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 125, in worker
result = (True, func(*args, **kwds))
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 48, in mapstar
return list(map(*args))
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/zeta_function.py", line 35, in _arb_lerch_phi
if _arb_c_lerch_phi(ctypes.byref(cy), *args, 0):
KeyboardInterrupt
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 48, in mapstar
return list(map(*args))
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/stocproc.py", line 1064, in alpha_times_pi
return alpha(tau) * np.pi
KeyboardInterrupt
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 18, in zeta_func
return lerch_phi(shift, s, a)
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/stocproc.py", line 1064, in alpha_times_pi
return alpha(tau) * np.pi
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 377, in __call__
,* zeta_func(
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/zeta_function.py", line 35, in _arb_lerch_phi
if _arb_c_lerch_phi(ctypes.byref(cy), *args, 0):
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 18, in zeta_func
return lerch_phi(shift, s, a)
KeyboardInterrupt
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/zeta_function.py", line 35, in _arb_lerch_phi
if _arb_c_lerch_phi(ctypes.byref(cy), *args, 0):
KeyboardInterrupt
#+end_example
# [goto error]
#+begin_example
---------------------------------------------------------------------------
OSError Traceback (most recent call last)
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hiro_models/model_auxiliary.py:148, in get_data(model, data_path, read_only, **kwargs)
 147 try:
--> 148 return HIData(path, read_only=read_only, robust=False, **kwargs)
 149 except:
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/core/hierarchy_data.py:252, in HIData.__init__(self, hdf5_name, read_only, hi_key, hi_key_bin, hi_key_bin_hash, check_consistency, overwrite_key, robust, stream_file)
 251 if not hi_key:
--> 252 with h5py.File(hdf5_name, "r", libver="latest") as h5File:
 253 try:
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/h5py/_hl/files.py:533, in File.__init__(self, name, mode, driver, libver, userblock_size, swmr, rdcc_nslots, rdcc_nbytes, rdcc_w0, track_order, fs_strategy, fs_persist, fs_threshold, fs_page_size, page_buf_size, min_meta_keep, min_raw_keep, locking, alignment_threshold, alignment_interval, **kwds)
 530 fcpl = make_fcpl(track_order=track_order, fs_strategy=fs_strategy,
 531 fs_persist=fs_persist, fs_threshold=fs_threshold,
 532 fs_page_size=fs_page_size)
--> 533 fid = make_fid(name, mode, userblock_size, fapl, fcpl, swmr=swmr)
 535 if isinstance(libver, tuple):
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/h5py/_hl/files.py:226, in make_fid(name, mode, userblock_size, fapl, fcpl, swmr)
 225 flags |= h5f.ACC_SWMR_READ
--> 226 fid = h5f.open(name, flags, fapl=fapl)
 227 elif mode == 'r+':
File h5py/_objects.pyx:54, in h5py._objects.with_phil.wrapper()
File h5py/_objects.pyx:55, in h5py._objects.with_phil.wrapper()
File h5py/h5f.pyx:106, in h5py.h5f.open()
OSError: Unable to open file (file is already open for write (may use <h5clear file> to clear file consistency flags))
During handling of the above exception, another exception occurred:
KeyboardInterrupt Traceback (most recent call last)
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/method_ft.py:558, in get_dt_for_accurate_interpolation(t_max, tol, ft_ref, diff_method)
 556 ft_ref_n[::2] = ft_ref_n_old
--> 558 ft_ref_n_new = np.asarray(pool.map(ft_ref, tau[1::2]))
 560 ft_ref_n[1::2] = np.array(ft_ref_n_new)
File /nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py:364, in Pool.map(self, func, iterable, chunksize)
 360 '''
 361 Apply `func` to each element in `iterable`, collecting the results
 362 in a list that is returned.
 363 '''
--> 364 return self._map_async(func, iterable, mapstar, chunksize).get()
File /nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py:765, in ApplyResult.get(self, timeout)
 764 def get(self, timeout=None):
--> 765 self.wait(timeout)
 766 if not self.ready():
File /nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py:762, in ApplyResult.wait(self, timeout)
 761 def wait(self, timeout=None):
--> 762 self._event.wait(timeout)
File /nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/threading.py:581, in Event.wait(self, timeout)
 580 if not signaled:
--> 581 signaled = self._cond.wait(timeout)
 582 return signaled
File /nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/threading.py:312, in Condition.wait(self, timeout)
 311 if timeout is None:
--> 312 waiter.acquire()
 313 gotit = True
KeyboardInterrupt:
During handling of the above exception, another exception occurred:
KeyboardInterrupt Traceback (most recent call last)
Input In [15], in <cell line: 1>()
----> 1 plot_excited_state(model, params, with_init=True)
File ~/Documents/Projects/UNI/master/masterarb/python/energy_flow_proper/10_antizeno_engine/plot_utils.py:40, in wrap_plot.<locals>.wrapped(ax, setup_function, *args, **kwargs)
 37 if not ax:
 38 fig, ax = setup_function()
---> 40 ret_val = f(*args, ax=ax, **kwargs)
 41 return (fig, ax, ret_val) if ret_val else (fig, ax)
File ~/Documents/Projects/UNI/master/masterarb/python/energy_flow_proper/10_antizeno_engine/anti_zeno_engine.py:286, in plot_excited_state(model, params, ax, with_init)
 284 @wrap_plot
 285 def plot_excited_state(model, params, ax=None, with_init=False):
--> 286 with aux.get_data(model) as data:
 287 mask = model.t >= 0 if with_init else model.t > params.τ_init
 289 plot_with_σ(
 290 model.t[mask],
 291 abs(
 (...)
 299 ax=ax,
 300 )
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hiro_models/model_auxiliary.py:152, in get_data(model, data_path, read_only, **kwargs)
 148 return HIData(path, read_only=read_only, robust=False, **kwargs)
 149 except:
 150 return HIData(
 151 path,
--> 152 hi_key=model.hops_config,
 153 read_only=False,
 154 check_consistency=False,
 155 overwrite_key=True,
 156 robust=False,
 157 **kwargs,
 158 )
 160 else:
 161 raise RuntimeError(f"No data found for model with hash '{hexhash}'.")
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hiro_models/one_qubit_model.py:842, in QubitModelMutliBath.hops_config(self)
 833 default_solver_args.update(self.solver_args)
 835 integration = params.IntP(t=self.t, **default_solver_args)
 837 return params.HIParams(
 838 SysP=system,
 839 IntP=integration,
 840 HiP=hierarchy,
 841 Eta=[self.driving_process(i) for i in range(self.num_baths)],
--> 842 EtaTherm=self.thermal_processes,
 843 )
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hiro_models/one_qubit_model.py:726, in QubitModelMutliBath.thermal_processes(self)
 719 @property
 720 def thermal_processes(self) -> list[Optional[hopsflow.hopsflow.StocProc]]:
 721 """
 722  The thermal noise stochastic processes for each bath.
 723  :any:`None` means zero temperature.
 724  """
--> 726 return [self.thermal_process(i) for i in range(self.num_baths)]
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hiro_models/one_qubit_model.py:726, in <listcomp>(.0)
 719 @property
 720 def thermal_processes(self) -> list[Optional[hopsflow.hopsflow.StocProc]]:
 721 """
 722  The thermal noise stochastic processes for each bath.
 723  :any:`None` means zero temperature.
 724  """
--> 726 return [self.thermal_process(i) for i in range(self.num_baths)]
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hiro_models/one_qubit_model.py:706, in QubitModelMutliBath.thermal_process(self, i)
 703 if self.T[i] == 0:
 704 return None
--> 706 return (
 707  sp.StocProc_TanhSinh
 708  if self.therm_methods[i] == "tanhsinh"
 709  else sp.StocProc_FFT
 710 )(
 711  spectral_density=self.thermal_spectral_density(i),
 712  alpha=self.thermal_correlations(i),
 713  t_max=self.t.max(),
 714  intgr_tol=self.thermal_process_tolerances[i].integration,
 715  intpl_tol=self.thermal_process_tolerances[i].interpolation,
 716  negative_frequencies=False,
 717 )
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/stocproc.py:111, in StocProc.__init__(self, *args, **kwargs)
 108 if "invalidate_cache" in kwargs:
 109 del kwargs["invalidate_cache"]
--> 111 self.__post_init__(*args, **kwargs)
 113 if not hasattr(self, "_from_cache") and hasattr(self, "_cache_path"):
 114 with open(self._cache_path, "wb") as f:
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/stocproc.py:568, in StocProc_FFT.__post_init__(self, spectral_density, t_max, alpha, intgr_tol, intpl_tol, seed, negative_frequencies, scale, calc_deriv)
 566 if not negative_frequencies:
 567 log.info("non neg freq only")
--> 568 a, b, N, dx, dt = method_ft.calc_ab_N_dx_dt(
 569  integrand=spectral_density,
 570  intgr_tol=intgr_tol,
 571  intpl_tol=intpl_tol,
 572  t_max=t_max,
 573  ft_ref=ft_ref,
 574  opt_b_only=True,
 575  )
 576 else:
 577 log.info("use neg freq")
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/method_ft.py:633, in calc_ab_N_dx_dt(integrand, intgr_tol, intpl_tol, t_max, ft_ref, opt_b_only, diff_method)
 630 c = t_max
 632 c = min(c, t_max)
--> 633 dt_tol = get_dt_for_accurate_interpolation(
 634  t_max=c, tol=intpl_tol, ft_ref=ft_ref, diff_method=diff_method
 635 )
 637 log.debug("requires dt < {:.3e}".format(dt_tol))
 639 log.debug("get_N_a_b_for_accurate_fourier_integral, please wait ...")
File /nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/method_ft.py:580, in get_dt_for_accurate_interpolation(t_max, tol, ft_ref, diff_method)
 578 finally:
 579 pool.close()
--> 580 pool.join()
File /nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py:662, in Pool.join(self)
 660 elif self._state not in (CLOSE, TERMINATE):
 661 raise ValueError("In unknown state")
--> 662 self._worker_handler.join()
 663 self._task_handler.join()
 664 self._result_handler.join()
File /nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/threading.py:1060, in Thread.join(self, timeout)
 1057 raise RuntimeError("cannot join current thread")
 1059 if timeout is None:
-> 1060 self._wait_for_tstate_lock()
 1061 else:
 1062 # the behavior of a negative timeout isn't documented, but
 1063 # historically .join(timeout=x) for x<0 has acted as if timeout=0
 1064 self._wait_for_tstate_lock(timeout=max(timeout, 0))
File /nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/threading.py:1080, in Thread._wait_for_tstate_lock(self, block, timeout)
 1077 return
 1079 try:
-> 1080 if lock.acquire(block, timeout):
 1081 lock.release()
 1082 self._stop()
KeyboardInterrupt:
Process ForkPoolWorker-35:
Process ForkPoolWorker-36:
Process ForkPoolWorker-33:
Traceback (most recent call last):
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
Traceback (most recent call last):
Traceback (most recent call last):
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 114, in worker
task = get()
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/queues.py", line 365, in get
res = self._reader.recv_bytes()
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 125, in worker
result = (True, func(*args, **kwds))
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/connection.py", line 221, in recv_bytes
buf = self._recv_bytes(maxlength)
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 114, in worker
task = get()
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 48, in mapstar
return list(map(*args))
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/queues.py", line 364, in get
with self._rlock:
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/connection.py", line 419, in _recv_bytes
buf = self._recv(4)
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/stocproc.py", line 1064, in alpha_times_pi
return alpha(tau) * np.pi
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/connection.py", line 384, in _recv
chunk = read(handle, remaining)
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 377, in __call__
* zeta_func(
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/synchronize.py", line 95, in __enter__
return self._semlock.__enter__()
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 18, in zeta_func
return lerch_phi(shift, s, a)
KeyboardInterrupt
KeyboardInterrupt
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/zeta_function.py", line 35, in _arb_lerch_phi
if _arb_c_lerch_phi(ctypes.byref(cy), *args, 0):
KeyboardInterrupt
Process ForkPoolWorker-34:
Traceback (most recent call last):
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 125, in worker
result = (True, func(*args, **kwds))
File "/nix/store/m3c6f5apfwnlnc5ms0ijk3k045ma0pqf-python3-3.9.13/lib/python3.9/multiprocessing/pool.py", line 48, in mapstar
return list(map(*args))
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/stocproc/stocproc.py", line 1064, in alpha_times_pi
return alpha(tau) * np.pi
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 377, in __call__
* zeta_func(
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/bcf.py", line 18, in zeta_func
return lerch_phi(shift, s, a)
File "/nix/store/jhaz6wa06sid3pc6c0g3wmid421vldwj-python3-3.9.13-env/lib/python3.9/site-packages/hops/util/zeta_function.py", line 35, in _arb_lerch_phi
if _arb_c_lerch_phi(ctypes.byref(cy), *args, 0):
KeyboardInterrupt
#+end_example
:END:
#+begin_src jupyter-python