From 9908405f4b61d04ef73b30fbaa8c94cb46c0d708 Mon Sep 17 00:00:00 2001 From: Valentin Boettcher Date: Mon, 11 Jul 2022 18:23:17 +0200 Subject: [PATCH] aaand now, let's try anti-zeno --- .../10_antizeno_engine/10_first_anti_zeno.py | 5 +- .../10_antizeno_engine/anti_zeno_engine.org | 64 +++++++++++++++---- 2 files changed, 56 insertions(+), 13 deletions(-) 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 545a74e..3cd3236 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 @@ -179,7 +179,7 @@ def anti_zeno_engine( (τ_mod, τ_c, τ_bath, cycles, model.ω_s, ω_0, τ_s, τ_off, n, Δ_switch, τ_init), ) = anti_zeno_engine( Δ=11, - ε=1, + ε=.5, ω_c=2, ε_couple=0.2, n=3, @@ -250,8 +250,9 @@ with aux.get_data(model) as data: fs.plot_with_σ(model.t, model.bath_energy(data), bath=0, ax=ax) fs.plot_with_σ(model.t, model.bath_energy(data), bath=1, ax=ax) fs.plot_with_σ(model.t, model.bath_energy(data).sum_baths(), ax=ax) + fs.plot_with_σ(model.t, model.total_energy_from_power(data), ax=ax) - #fs.plot_with_σ(model.t, model.total_energy(data), ax=ax) + fs.plot_with_σ(model.t, model.total_energy(data), ax=ax) #fs.plot_with_σ(model.t, model.interaction_energy(data).for_bath(1), ax=ax) #fs.plot_with_σ(model.t, model.system_energy(data), ax=ax) 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 b8778f1..874cb94 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 @@ -25,7 +25,10 @@ Init ray and silence stocproc. #+end_src #+RESULTS: -: RayContext(dashboard_url='', python_version='3.9.13', ray_version='1.13.0', ray_commit='e4ce38d001dbbe09cd21c497fedd03d692b2be3e', address_info={'node_ip_address': '141.30.17.225', 'raylet_ip_address': '141.30.17.225', 'redis_address': None, 'object_store_address': '/tmp/ray/session_2022-07-11_17-56-08_601005_343893/sockets/plasma_store', 'raylet_socket_name': '/tmp/ray/session_2022-07-11_17-56-08_601005_343893/sockets/raylet', 'webui_url': '', 'session_dir': '/tmp/ray/session_2022-07-11_17-56-08_601005_343893', 'metrics_export_port': 33700, 'gcs_address': '141.30.17.225:52666', 'address': '141.30.17.225:52666', 'node_id': 'c868ef9cb4594a88f4d75d22b97e6e30132618843dcccc5031cf3e62'}) +:RESULTS: +: 2022-07-11 18:20:21,473 INFO worker.py:956 -- Connecting to existing Ray cluster at address: 141.30.17.16:6379 +: RayContext(dashboard_url='', python_version='3.9.13', ray_version='1.13.0', ray_commit='e4ce38d001dbbe09cd21c497fedd03d692b2be3e', address_info={'node_ip_address': '141.30.17.225', 'raylet_ip_address': '141.30.17.225', 'redis_address': None, 'object_store_address': '/tmp/ray/session_2022-07-07_19-11-40_930040_398742/sockets/plasma_store.1', 'raylet_socket_name': '/tmp/ray/session_2022-07-07_19-11-40_930040_398742/sockets/raylet', 'webui_url': '', 'session_dir': '/tmp/ray/session_2022-07-07_19-11-40_930040_398742', 'metrics_export_port': 63895, 'gcs_address': '141.30.17.16:6379', 'address': '141.30.17.16:6379', 'node_id': 'ca50ebda6b5629030bd1ad7d4047a1a1fde48b5b2041bf7976b7d105'}) +:END: #+begin_src jupyter-python :results none from hops.util.logging_setup import logging_setup @@ -199,7 +202,7 @@ Init ray and silence stocproc. (τ_mod, τ_c, τ_bath, cycles, model.ω_s, ω_0, τ_s, τ_off, n, Δ_switch, τ_init), ) = anti_zeno_engine( Δ=11, - ε=1, + ε=.5, ω_c=2, ε_couple=0.2, n=3, @@ -245,8 +248,8 @@ Let's test the assumptions of the paper. #+RESULTS: :RESULTS: -| | -[[file:./.ob-jupyter/f197c60a4dd01166b654064ca4d794d7645d96df.svg]] +| | +[[file:./.ob-jupyter/dd12f4cd01aa9f158355ddff8a1151fa1fb0107e.svg]] :END: #+begin_src jupyter-python :tangle nil @@ -264,8 +267,8 @@ Let's test the assumptions of the paper. #+RESULTS: :RESULTS: -| | -[[file:./.ob-jupyter/0b52df5f4c1b72a077f2a791adbb3c8c1bd2fe32.svg]] +| | +[[file:./.ob-jupyter/1c64ce4d6cddbe1bf47c0c9900650b58860bdbb7.svg]] :END: #+begin_src jupyter-python @@ -322,7 +325,33 @@ Let's test the assumptions of the paper. #+end_src #+RESULTS: -: 602a65a4-50ca-485c-84f9-0f0adcce05f1 +#+begin_example + [INFO hops.core.integration 537693] Choosing the nonlinear integrator. + [INFO hops.core.integration 537693] Using 21 integrators. + [INFO hops.core.integration 537693] Some 100 trajectories have to be integrated. + [INFO hops.core.integration 537693] Using 1820 hierarchy states. + 0% 0/100 [00:00