diff --git a/python/energy_flow_proper/10_antizeno_engine/10_first_anti_zeno.py b/python/energy_flow_proper/10_antizeno_engine/10_first_anti_zeno.py index 4efccea..15c74bd 100644 --- a/python/energy_flow_proper/10_antizeno_engine/10_first_anti_zeno.py +++ b/python/energy_flow_proper/10_antizeno_engine/10_first_anti_zeno.py @@ -288,7 +288,7 @@ vs = np.linspace(0.1, 10, 100) plt.plot(vs, chi(vs, ω_0)) plt.plot(vs, G_h(vs)) -aux.integrate(model, 100) +aux.integrate(model, 10000) #_, ax = fs.plot_energy_overview(model, markersize=1, ensemble_args=dict(gc_sleep=0.05)) diff --git a/python/energy_flow_proper/10_antizeno_engine/anti_zeno_engine.org b/python/energy_flow_proper/10_antizeno_engine/anti_zeno_engine.org index b6704aa..d01bff7 100644 --- a/python/energy_flow_proper/10_antizeno_engine/anti_zeno_engine.org +++ b/python/energy_flow_proper/10_antizeno_engine/anti_zeno_engine.org @@ -316,8 +316,8 @@ Let's test the assumptions of the paper. #+RESULTS: :RESULTS: -: -[[file:./.ob-jupyter/52319af8927d7fbb9f93dd71384478ec06ef6faf.svg]] +: +[[file:./.ob-jupyter/8139a949bde14728a4e6d9f0d2640fc8edef8218.svg]] :END: #+begin_src jupyter-python :tangle nil @@ -369,8 +369,8 @@ Let's test the assumptions of the paper. #+RESULTS: :RESULTS: -[[file:./.ob-jupyter/0e7e96064672b6bc919ca89d52825c59e9543cfb.svg]] -[[file:./.ob-jupyter/00d9a4ef93705497b6722f5e6ab91c0542aea18e.svg]] +[[file:./.ob-jupyter/9b211b2cea3136d4aa15c794b0b0caa6ae2d4677.svg]] +[[file:./.ob-jupyter/8c422814e6d81aa1252f97e948775b58da037a44.svg]] :END: #+begin_src jupyter-python @@ -424,145 +424,38 @@ Let's test the assumptions of the paper. ** Integration #+begin_src jupyter-python - aux.integrate(model, 100) + aux.integrate(model, 10000) #+End_src #+RESULTS: -:RESULTS: #+begin_example - [INFO hops.core.integration 241007] Choosing the nonlinear integrator. - [INFO hops.core.integration 241007] Using 21 integrators. - [INFO hops.core.integration 241007] Some 90 trajectories have to be integrated. - [INFO hops.core.integration 241007] Using 165 hierarchy states. - 0% 0/90 [00:00::_M_invoke() - (raylet) std::_Function_handler<>::_M_invoke() - (raylet) std::_Function_handler<>::_M_invoke() - (raylet) ray::rpc::ClientCallImpl<>::OnReplyReceived() - (raylet) std::_Function_handler<>::_M_invoke() - (raylet) EventTracker::RecordExecution() - (raylet) std::_Function_handler<>::_M_invoke() - (raylet) boost::asio::detail::completion_handler<>::do_complete() - (raylet) boost::asio::detail::scheduler::do_run_one() - (raylet) boost::asio::detail::scheduler::run() - (raylet) boost::asio::io_context::run() - (raylet) main - (raylet) __libc_start_call_main - (raylet) - (raylet) [2022-07-13 10:41:27,801 C 578762 578762] (raylet) node_manager.cc:170: This node has beem marked as dead. - (raylet) *** StackTrace Information *** - (raylet) ray::SpdLogMessage::Flush() - (raylet) ray::RayLog::~RayLog() - (raylet) std::_Function_handler<>::_M_invoke() - (raylet) std::_Function_handler<>::_M_invoke() - (raylet) std::_Function_handler<>::_M_invoke() - (raylet) ray::rpc::ClientCallImpl<>::OnReplyReceived() - (raylet) std::_Function_handler<>::_M_invoke() - (raylet) EventTracker::RecordExecution() - (raylet) std::_Function_handler<>::_M_invoke() - (raylet) boost::asio::detail::completion_handler<>::do_complete() - (raylet) boost::asio::detail::scheduler::do_run_one() - (raylet) boost::asio::detail::scheduler::run() - (raylet) boost::asio::io_context::run() - (raylet) main - (raylet) __libc_start_call_main - (raylet) - (raylet) [2022-07-20 10:41:29,167 E 122306 122344] (raylet) agent_manager.cc:107: The raylet exited immediately because the Ray agent failed. The raylet fate shares with the agent. This can happen because the Ray agent was unexpectedly killed or failed. See `dashboard_agent.log` for the root cause. - 20% 18/90 [01:56<07:45, 6.47s/it][INFO hops.core.signal_delay 241007] caught sig 'SIGINT' - 21% 19/90 [01:59<07:26, 6.29s/it] - 2022-07-20 15:03:10,268 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,271 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,273 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,277 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,278 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,281 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,285 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - [INFO hops.core.signal_delay 241007] caught 1 signal(s) - 2022-07-20 15:03:10,286 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - [INFO hops.core.signal_delay 241007] emit signal 'SIGINT' - 2022-07-20 15:03:10,288 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,289 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,291 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,292 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,293 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,343 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,365 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,412 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,435 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,468 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,475 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,513 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,521 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,539 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,577 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. - 2022-07-20 15:03:10,583 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information. + 100% 1/1 [02:31<00:00, 151.78s/it] + [INFO hops.core.integration 241007] Choosing the nonlinear integrator. + [INFO hops.core.integration 241007] Using 21 integrators. + [INFO hops.core.integration 241007] Some 99 trajectories have to be integrated. + [INFO hops.core.integration 241007] Using 165 hierarchy states. + 0% 0/99 [00:00() - ----> 1 aux.integrate(model, 100) - - File ~/src/two_qubit_model/hiro_models/model_auxiliary.py:108, in integrate(model, n, data_path, clear_pd) -  98 # with model_db(data_path) as db: -  99 # if hash in db and "data" db[hash] -  101 supervisor = HOPSSupervisor( -  102 model.hops_config, -  103 n, -  104 data_path=data_path, -  105 data_name=hash, -  106 ) - --> 108 supervisor.integrate(clear_pd) -  110 with supervisor.get_data(True) as data: -  111 with model_db(data_path) as db: - - File ~/src/hops/hops/core/integration.py:1288, in HOPSSupervisor.integrate(self, clear_pd) -  1285 break -  1287 integration.update() - -> 1288 data.new_samples( -  1289 idx=index, -  1290 incomplete=incomplete, -  1291 psi0=psi0, -  1292 aux_states=aux_states, -  1293 stoc_proc=stoc_proc, -  1294 result_type=self.params.HiP.result_type, -  1295 normed=self._normed_average, -  1296 rng_seed=seed, -  1297 ) - - File ~/src/hops/hops/core/signal_delay.py:87, in sig_delay.__exit__(self, exc_type, exc_val, exc_tb) -  84 if len(self.sigh.sigs_caught) > 0 and self.handler is not None: -  85 self.handler(self.sigh.sigs_caught) - ---> 87 self._restore() - - File ~/src/hops/hops/core/signal_delay.py:68, in sig_delay._restore(self) -  66 for i, s in enumerate(self.sigs): -  67 signal.signal(s, self.old_handlers[i]) - ---> 68 self.sigh.emit() - - File ~/src/hops/hops/core/signal_delay.py:42, in SigHandler.emit(self) -  40 for s in self.sigs_caught: -  41 log.info("emit signal '{}'".format(SIG_MAP[s])) - ---> 42 os.kill(os.getpid(), s) - - KeyboardInterrupt: -#+end_example -:END: #+Begin_src jupyter-python @@ -578,21 +471,11 @@ Let's test the assumptions of the paper. #+RESULTS: :RESULTS: -#+begin_example - Loading: 90% 9/10 [00:00<00:00, 138.66it/s] - (_ensemble_remote_function pid=22270, ip=141.30.17.8) WARNING:root:Using mpmath for the zeta (lerch phi) function. Install a recent version of ``arb`` for more performance. - (_ensemble_remote_function pid=22267, ip=141.30.17.8) WARNING:root:Using mpmath for the zeta (lerch phi) function. Install a recent version of ``arb`` for more performance. - (_ensemble_remote_function pid=22265, ip=141.30.17.8) WARNING:root:Using mpmath for the zeta (lerch phi) function. Install a recent version of ``arb`` for more performance. - (_ensemble_remote_function pid=22268, ip=141.30.17.8) WARNING:root:Using mpmath for the zeta (lerch phi) function. Install a recent version of ``arb`` for more performance. - (_ensemble_remote_function pid=22264, ip=141.30.17.8) WARNING:root:Using mpmath for the zeta (lerch phi) function. Install a recent version of ``arb`` for more performance. - (_ensemble_remote_function pid=22266, ip=141.30.17.8) WARNING:root:Using mpmath for the zeta (lerch phi) function. Install a recent version of ``arb`` for more performance. - (_ensemble_remote_function pid=22271, ip=141.30.17.8) WARNING:root:Using mpmath for the zeta (lerch phi) function. Install a recent version of ``arb`` for more performance. - (_ensemble_remote_function pid=22269, ip=141.30.17.8) WARNING:root:Using mpmath for the zeta (lerch phi) function. Install a recent version of ``arb`` for more performance. - [INFO root 241007] Writing cache to: results/interaction_power_67cca8aef6bb8d5c44cd03d2e435820ea5a6ad6df4c336139cd691b49fb51540_interaction_energy_task_29_None_1132783dc3e4b63333ab34011167d84b9c35e6e83409bc4c77ac7f48df69c157.npy - Loading: 80% 4/5 [00:00<00:00, 171.20it/s] - [INFO root 241007] Writing cache to: results/f8ad5eb5a8a6d9fb98ca17682691e53506757d7244d54c67c53f9f3314c68a05_67cca8aef6bb8d5c44cd03d2e435820ea5a6ad6df4c336139cd691b49fb51540_op_exp_task_29_None_8eaaa0856f1a481703397da647054bbc3035c456b4fe49fb302a4fb76fa7c865.npy -#+end_example -[[file:./.ob-jupyter/b3e65afd3656a3af59f525458b4df39d5222d518.svg]] +: Loading: 98% 98/100 [00:13<00:00, 7.24it/s] +: [INFO root 241007] Writing cache to: results/interaction_power_7a2f8149e57b9baf77b692150756e3ffd42b69a9c58bf684d60dc7a29b750ae8_interaction_energy_task_100_None_a1ebae168ba6027d405a66e94ca963c84f8aedb16621da5954b5631bcd1636c4.npy +: Loading: 98% 49/50 [00:01<00:00, 37.90it/s] +: [INFO root 241007] Writing cache to: results/f8ad5eb5a8a6d9fb98ca17682691e53506757d7244d54c67c53f9f3314c68a05_7a2f8149e57b9baf77b692150756e3ffd42b69a9c58bf684d60dc7a29b750ae8_op_exp_task_100_None_1ff6d34b7de893995dfd0ae8efe2ed770a6773ef52b6e717b4b4ee9b9e5d285d.npy +[[file:./.ob-jupyter/8d4453ac0d9d598221ae19b35505623ee1e35dbe.svg]] :END: - **too fast decoupling kills it** @@ -619,7 +502,19 @@ Let's test the assumptions of the paper. #+end_src #+RESULTS: -[[file:./.ob-jupyter/be0e2c3ba2a64e4022395469cd31f4dc5d5e37f2.svg]] +:RESULTS: +: /nix/store/k12zz30f8mdh9da2r17c5aw87f7kagl9-python3-3.9.13-env/lib/python3.9/site-packages/matplotlib/cbook/__init__.py:1298: ComplexWarning: Casting complex values to real discards the imaginary part +: return np.asarray(x, float) +: /nix/store/k12zz30f8mdh9da2r17c5aw87f7kagl9-python3-3.9.13-env/lib/python3.9/site-packages/matplotlib/axes/_axes.py:5218: ComplexWarning: Casting complex values to real discards the imaginary part +: pts[0] = start +: /nix/store/k12zz30f8mdh9da2r17c5aw87f7kagl9-python3-3.9.13-env/lib/python3.9/site-packages/matplotlib/axes/_axes.py:5219: ComplexWarning: Casting complex values to real discards the imaginary part +: pts[N + 1] = end +: /nix/store/k12zz30f8mdh9da2r17c5aw87f7kagl9-python3-3.9.13-env/lib/python3.9/site-packages/matplotlib/axes/_axes.py:5222: ComplexWarning: Casting complex values to real discards the imaginary part +: pts[1:N+1, 1] = dep1slice +: /nix/store/k12zz30f8mdh9da2r17c5aw87f7kagl9-python3-3.9.13-env/lib/python3.9/site-packages/matplotlib/axes/_axes.py:5224: ComplexWarning: Casting complex values to real discards the imaginary part +: pts[N+2:, 1] = dep2slice[::-1] +[[file:./.ob-jupyter/0ed19e6213d611d35fa44ceba7ce2c0255bf0bce.svg]] +:END: - no steady state ... but we have to average... @@ -631,8 +526,8 @@ Let's test the assumptions of the paper. #+RESULTS: :RESULTS: -| | -[[file:./.ob-jupyter/8c2b46e2adbdadab36d7fdf6bc2d97f63c2217b3.svg]] +| | +[[file:./.ob-jupyter/2d5aee00e6f799251447c4d0206a214bbd5c18b0.svg]] :END: ** TODO Power and Efficiency @@ -715,8 +610,8 @@ We need the time points where we sample the total energy. #+RESULTS: :RESULTS: -: \(P=-0.12\pm 0.05\) -[[file:./.ob-jupyter/ee342b28ea9884126abbf873ede39644af783db9.svg]] +: \(P=-0.006\pm 0.012\) +[[file:./.ob-jupyter/6128cf745a54bf219e524c94c44516724ee55f45.svg]] :END: