From 449b48c8e855e118db6ead8330a046f1c5f140c4 Mon Sep 17 00:00:00 2001 From: Valentin Boettcher Date: Mon, 5 Sep 2022 15:55:13 +0200 Subject: [PATCH] 10: does not converge --- .../10_antizeno_engine/continously_couple.org | 519 +++++++++++++++++- 1 file changed, 504 insertions(+), 15 deletions(-) diff --git a/python/energy_flow_proper/10_antizeno_engine/continously_couple.org b/python/energy_flow_proper/10_antizeno_engine/continously_couple.org index 51d7c58..167264b 100644 --- a/python/energy_flow_proper/10_antizeno_engine/continously_couple.org +++ b/python/energy_flow_proper/10_antizeno_engine/continously_couple.org @@ -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 | | (() ) | -[[file:./.ob-jupyter/72099bcb22465077395cd8acccbecc8f6e768f84.svg]] +: 66590 +| hline | | (() ) | :END: @@ -104,20 +105,14 @@ Here we try to reproduce the anti zeno engine from the paper. #+end_src #+RESULTS: -:RESULTS: |
| | -[[file:./.ob-jupyter/d1345d964ba6362247d8f34e3ea98e5b5f91a7db.svg]] -:END: #+begin_src jupyter-python plot_az_sd_overview(model, params) #+end_src #+RESULTS: -:RESULTS: |
| | -[[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: -|
| | -[[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 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 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 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 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 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 () + ----> 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..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 (.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