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 d6e790c..0b4968d 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 @@ -148,12 +148,12 @@ model, Δ, (τ_mod, τ_c, τ_bath, cycles, model.ω_s, ω_0, τ_s, τ_off, n, Δ Δ=1.3, ε=.5, ω_c=0.2, - ε_couple=0.7, + ε_couple=0.8, n=3, - detune=-0.19, + detune=-0.10, ω_0=2, T_c=3, - T_h=11, + T_h=6, δ=[1] * 2, γ=0.5 / 10, switch_cycles=1, 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 0fcf668..a2937e2 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,7 @@ 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-06-28_09-55-52_051747_8907/sockets/plasma_store', 'raylet_socket_name': '/tmp/ray/session_2022-06-28_09-55-52_051747_8907/sockets/raylet', 'webui_url': '', 'session_dir': '/tmp/ray/session_2022-06-28_09-55-52_051747_8907', 'metrics_export_port': 61788, 'gcs_address': '141.30.17.225:42071', 'address': '141.30.17.225:42071', 'node_id': '5a08a875431de9ddeac7253b9760300dd31cc6a01c2245d351a5ce52'}) +: 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-06-29_11-08-32_984823_8907/sockets/plasma_store', 'raylet_socket_name': '/tmp/ray/session_2022-06-29_11-08-32_984823_8907/sockets/raylet', 'webui_url': '', 'session_dir': '/tmp/ray/session_2022-06-29_11-08-32_984823_8907', 'metrics_export_port': 53835, 'gcs_address': '141.30.17.225:62738', 'address': '141.30.17.225:62738', 'node_id': '06bd5dfb2bc3dd0bd249d2952cdeb83f270068f2fea685f1cd984eee'}) #+begin_src jupyter-python :results none from hops.util.logging_setup import logging_setup @@ -168,12 +168,12 @@ Init ray and silence stocproc. Δ=1.3, ε=.5, ω_c=0.2, - ε_couple=0.7, + ε_couple=0.8, n=3, - detune=-0.19, + detune=-0.10, ω_0=2, T_c=3, - T_h=11, + T_h=6, δ=[1] * 2, γ=0.5 / 10, switch_cycles=1, @@ -196,7 +196,7 @@ Let's test the assumptions of the paper. #+end_src #+RESULTS: -: 6 +: 3 ** BCFs and Modulations #+begin_src jupyter-python :tangle nil @@ -212,8 +212,8 @@ Let's test the assumptions of the paper. #+RESULTS: :RESULTS: -| | -[[file:./.ob-jupyter/9bbb032a9a9418bd7edc1a954321cdc566e9ca78.svg]] +| | +[[file:./.ob-jupyter/d3eb36a9ceab4bb3e3725c25708186507f69046d.svg]] :END: #+begin_src jupyter-python :tangle nil @@ -226,13 +226,13 @@ Let's test the assumptions of the paper. plt.plot(ωs, np.sinc((ωs - ω_0 - Δ) * τ_s * cycles)) plt.plot(ωs, np.sinc((ωs - ω_0 - Δ) * τ_s * cycles * 2), color="orange", linewidth=.5) plt.plot(ωs, np.sinc((ωs - ω_0 - Δ) * τ_s * cycles * 10), color="yellow", linewidth=.4) - #plt.xlim(3, 4) + plt.xlim(3, 4) #+end_src #+RESULTS: :RESULTS: -| | -[[file:./.ob-jupyter/1f608368d4d13aa49b8f0263e3f9cf29ad99c4b7.svg]] +| 3.0 | 4.0 | +[[file:./.ob-jupyter/49f9da04e6433af73f0213f44347dd5f57e2a285.svg]] :END: #+begin_src jupyter-python @@ -308,8 +308,8 @@ Let's test the assumptions of the paper. #+RESULTS: :RESULTS: -: -[[file:./.ob-jupyter/8429b60bde6a26d5e0d2746d8a9b762934de8b34.svg]] +: +[[file:./.ob-jupyter/2a2d5d18230fa86567e88bbb847b3f053ea84623.svg]] :END: - **too fast decoupling kills it** @@ -335,7 +335,19 @@ Let's test the assumptions of the paper. #+end_src #+RESULTS: -[[file:./.ob-jupyter/8268ecba9ee72d6cc0c5d7a8007dac5f4294c8d6.svg]] +:RESULTS: +: /nix/store/av32kn66py1m065lq75rmmlxrsqvgbaa-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/av32kn66py1m065lq75rmmlxrsqvgbaa-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/av32kn66py1m065lq75rmmlxrsqvgbaa-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/av32kn66py1m065lq75rmmlxrsqvgbaa-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/av32kn66py1m065lq75rmmlxrsqvgbaa-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/b26c30bb6a05d98389cced5435c826ad2882953b.svg]] +:END: - no steady state ... but we have to average... @@ -347,8 +359,8 @@ Let's test the assumptions of the paper. #+RESULTS: :RESULTS: -| | -[[file:./.ob-jupyter/bd601a9265e37d0819f52afc09917dfdf582e956.svg]] +| | +[[file:./.ob-jupyter/3c9e5433c4bac6944b1ac0393acd7c6bb9d39b90.svg]] :END: ** TODO Power and Efficiency @@ -373,7 +385,7 @@ We need the time points where we sample the total energy. #+end_src #+RESULTS: -[[file:./.ob-jupyter/c1bd776c9a91ec8149ba76b8314408be6d2fbf5f.svg]] +[[file:./.ob-jupyter/0df98d7d13813b6c85c31084ec99ea6b95941b01.svg]]