try slower mod

This commit is contained in:
Valentin Boettcher 2022-06-21 15:55:15 +02:00
parent f82deb97a2
commit bce70ccd5c
2 changed files with 46 additions and 31 deletions

View file

@ -144,14 +144,14 @@ def anti_zeno_engine(
)
model, Δ, (τ_mod, τ_c, τ_bath, cycles, model.ω_s, ω_0, τ_s, τ_off, n, Δ_switch) = anti_zeno_engine(
Δ=1.5,
Δ=1,
ε=.5 ,
ω_c=0.2,
ε_couple=0.85,
n=10,
n=5,
detune=-0.19,
ω_0=2,
T_c=1.6,
T_c=4,
T_h=10,
δ=[1] * 2,
γ=0.5 / 10,
@ -184,7 +184,7 @@ vs = np.linspace(0.1, 10, 100)
plt.plot(vs, chi(vs, ω_0))
plt.plot(vs, G_h(vs))
aux.integrate(model, 2_000)
aux.integrate(model, 2000)
_, ax = fs.plot_energy_overview(model, markersize=1, ensemble_args=dict(gc_sleep=0))

View file

@ -25,7 +25,7 @@ Init ray and silence stocproc.
#+end_src
#+RESULTS:
: RayContext(dashboard_url='', python_version='3.9.12', ray_version='1.12.1', ray_commit='4863e33856b54ccf8add5cbe75e41558850a1b75', 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-21_13-08-14_458638_1925013/sockets/plasma_store', 'raylet_socket_name': '/tmp/ray/session_2022-06-21_13-08-14_458638_1925013/sockets/raylet', 'webui_url': '', 'session_dir': '/tmp/ray/session_2022-06-21_13-08-14_458638_1925013', 'metrics_export_port': 64729, 'gcs_address': '141.30.17.225:49099', 'address': '141.30.17.225:49099', 'node_id': 'ec4fc290c3ec934fac3c68412b4d14e7d7ee724ef9ced6b50a979a49'})
: RayContext(dashboard_url='', python_version='3.9.12', 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-21_15-53-28_658570_2169805/sockets/plasma_store', 'raylet_socket_name': '/tmp/ray/session_2022-06-21_15-53-28_658570_2169805/sockets/raylet', 'webui_url': '', 'session_dir': '/tmp/ray/session_2022-06-21_15-53-28_658570_2169805', 'metrics_export_port': 40469, 'gcs_address': '141.30.17.225:47762', 'address': '141.30.17.225:47762', 'node_id': '67301e784de4fea24ad52ad9a29741b018133f9b17ada3daf66a313c'})
#+begin_src jupyter-python :results none
from hops.util.logging_setup import logging_setup
@ -164,14 +164,14 @@ Init ray and silence stocproc.
* Model Definition
#+begin_src jupyter-python
model, Δ, (τ_mod, τ_c, τ_bath, cycles, model.ω_s, ω_0, τ_s, τ_off, n, Δ_switch) = anti_zeno_engine(
Δ=1.5,
Δ=1,
ε=.5 ,
ω_c=0.2,
ε_couple=0.85,
n=10,
n=5,
detune=-0.19,
ω_0=2,
T_c=1.6,
T_c=4,
T_h=10,
δ=[1] * 2,
γ=0.5 / 10,
@ -193,7 +193,7 @@ Let's test the assumptions of the paper.
#+end_src
#+RESULTS:
: 3
: 2
** BCFs and Modulations
#+begin_src jupyter-python :tangle nil
@ -209,8 +209,8 @@ Let's test the assumptions of the paper.
#+RESULTS:
:RESULTS:
| <matplotlib.lines.Line2D | at | 0x7ffa2bdb6b50> |
[[file:./.ob-jupyter/29d05104979912a12b959c1f097a852591c6054e.svg]]
| <matplotlib.lines.Line2D | at | 0x7f682cd0ad30> |
[[file:./.ob-jupyter/639832825a9dbdc531cad6e7fddeeccd8f62b97d.svg]]
:END:
#+begin_src jupyter-python :tangle nil
@ -223,13 +223,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:
| 3.0 | 4.0 |
[[file:./.ob-jupyter/ff58aca678dd6432231e8576520b612762b4c3f0.svg]]
| <matplotlib.lines.Line2D | at | 0x7f682cc57c40> |
[[file:./.ob-jupyter/131703f1c352eb45f976c10644c8072bee3c5fb0.svg]]
:END:
#+begin_src jupyter-python
@ -241,8 +241,8 @@ Let's test the assumptions of the paper.
#+RESULTS:
:RESULTS:
| <Figure | size | 432x288 | with | 1 | Axes> | <AxesSubplot:> |
[[file:./.ob-jupyter/33d008316e291c8d17151ead0418b51dfbd8b754.svg]]
[[file:./.ob-jupyter/3b942f66139250435ad213dc65f50c7309e1e663.svg]]
[[file:./.ob-jupyter/701494660db131c6c13d5c161e1812f1504ab5fd.svg]]
[[file:./.ob-jupyter/f6c628244bed62a53e1a5c355ed8af995d4a9219.svg]]
:END:
#+begin_src jupyter-python
@ -256,7 +256,7 @@ Let's test the assumptions of the paper.
#+RESULTS:
:RESULTS:
| <Figure | size | 432x288 | with | 1 | Axes> | <AxesSubplot:> |
[[file:./.ob-jupyter/ebe254a5be95214174933422d4682e3714882012.svg]]
[[file:./.ob-jupyter/052e3d18439fca2d3a800f3d94fdb168abe468ae.svg]]
:END:
@ -275,16 +275,31 @@ Let's test the assumptions of the paper.
#+RESULTS:
:RESULTS:
| <matplotlib.lines.Line2D | at | 0x7f807c0240a0> |
[[file:./.ob-jupyter/d30f7ce8409cc860d687e4ee003fa7d955d054fa.svg]]
| <matplotlib.lines.Line2D | at | 0x7f00fdbf7370> |
[[file:./.ob-jupyter/1f7aa2b79184210d6a4f95c9d65e0d0725319071.svg]]
:END:
** TODO Integration
#+begin_src jupyter-python
aux.integrate(model, 2_000)
aux.integrate(model, 2000)
#+end_src
#+RESULTS:
#+begin_example
[INFO hops.core.integration 2137276] Choosing the nonlinear integrator.
/home/hiro/src/hops/hops/core/hierarchy_data.py:497: UserWarning: test_file_version FAILED with exception Unable to open file (file is already open for write/SWMR write (may use <h5clear file> to clear file consistency flags))
warnings.warn(
/home/hiro/src/hops/hops/core/hierarchy_data.py:500: UserWarning: hdf5_name .data/5dbbf4ccb2b94cd79c9f1367996141a291047f11b5e36a39afc55d07a2ebc90b/_2/5dbbf4ccb2b94cd79c9f1367996141a291047f11b5e36a39afc55d07a2ebc90b_2535941ad73bf6f57abe33e03f90f548_1.h5
warnings.warn("hdf5_name {}".format(self.hdf5_name))
/home/hiro/src/hops/hops/core/hierarchy_data.py:503: UserWarning: Moving .data/5dbbf4ccb2b94cd79c9f1367996141a291047f11b5e36a39afc55d07a2ebc90b/_2/5dbbf4ccb2b94cd79c9f1367996141a291047f11b5e36a39afc55d07a2ebc90b_2535941ad73bf6f57abe33e03f90f548_1.h5 to .data/5dbbf4ccb2b94cd79c9f1367996141a291047f11b5e36a39afc55d07a2ebc90b/_2/5dbbf4ccb2b94cd79c9f1367996141a291047f11b5e36a39afc55d07a2ebc90b_2535941ad73bf6f57abe33e03f90f548_1.h5backup_1655818803.6054995 and starting fresh.
warnings.warn(
[INFO hops.core.integration 2137276] Using 4 integrators.
[INFO hops.core.integration 2137276] Some 1 trajectories have to be integrated.
[INFO hops.core.integration 2137276] Using 1001 hierarchy states.
0% 0/1 [00:00<?, ?it/s]
#+end_example
#+begin_src jupyter-python
_, ax = fs.plot_energy_overview(model, markersize=1, ensemble_args=dict(gc_sleep=0))
@ -297,8 +312,8 @@ Let's test the assumptions of the paper.
#+RESULTS:
:RESULTS:
: <matplotlib.legend.Legend at 0x7ffa308d2e50>
[[file:./.ob-jupyter/4aa25d0d35da92d90b13e4144ae4ead13d102566.svg]]
: <matplotlib.legend.Legend at 0x7f683b8c6f40>
[[file:./.ob-jupyter/51ea7ea8fa1c89b34891f4d3d5e792eeece1ce4b.svg]]
:END:
- **too fast decoupling kills it**
@ -321,17 +336,17 @@ Let's test the assumptions of the paper.
#+RESULTS:
:RESULTS:
: /nix/store/pwhaggpgvxhy410r6vlx448v9lp8n08s-python3-3.9.12-env/lib/python3.9/site-packages/matplotlib/cbook/__init__.py:1298: ComplexWarning: Casting complex values to real discards the imaginary part
: /nix/store/j9m9rb3yfvdx891s95j90alz8pazq5kw-python3-3.9.12-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/pwhaggpgvxhy410r6vlx448v9lp8n08s-python3-3.9.12-env/lib/python3.9/site-packages/matplotlib/axes/_axes.py:5218: ComplexWarning: Casting complex values to real discards the imaginary part
: /nix/store/j9m9rb3yfvdx891s95j90alz8pazq5kw-python3-3.9.12-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/pwhaggpgvxhy410r6vlx448v9lp8n08s-python3-3.9.12-env/lib/python3.9/site-packages/matplotlib/axes/_axes.py:5219: ComplexWarning: Casting complex values to real discards the imaginary part
: /nix/store/j9m9rb3yfvdx891s95j90alz8pazq5kw-python3-3.9.12-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/pwhaggpgvxhy410r6vlx448v9lp8n08s-python3-3.9.12-env/lib/python3.9/site-packages/matplotlib/axes/_axes.py:5222: ComplexWarning: Casting complex values to real discards the imaginary part
: /nix/store/j9m9rb3yfvdx891s95j90alz8pazq5kw-python3-3.9.12-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/pwhaggpgvxhy410r6vlx448v9lp8n08s-python3-3.9.12-env/lib/python3.9/site-packages/matplotlib/axes/_axes.py:5224: ComplexWarning: Casting complex values to real discards the imaginary part
: /nix/store/j9m9rb3yfvdx891s95j90alz8pazq5kw-python3-3.9.12-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/5984b33291841aa3665cf7d8749dd0d66370ed8a.svg]]
[[file:./.ob-jupyter/cce36b28d371016575407492deb9316748064d73.svg]]
:END:
@ -345,8 +360,8 @@ Let's test the assumptions of the paper.
#+RESULTS:
:RESULTS:
| <matplotlib.lines.Line2D | at | 0x7ffa3017d6a0> |
[[file:./.ob-jupyter/e5661de260b83a368f954a1307035a66735b5272.svg]]
| <matplotlib.lines.Line2D | at | 0x7f682c9e7880> |
[[file:./.ob-jupyter/ce2fdae5dbd2a2377604b4c283036fafc81240af.svg]]
:END:
** TODO Power and Efficiency
@ -371,7 +386,7 @@ We need the time points where we sample the total energy.
#+end_src
#+RESULTS:
[[file:./.ob-jupyter/bd5332cacf1a951d281eb6a02c13885fa26f0c51.svg]]
[[file:./.ob-jupyter/d833ee1e584b08c25336774111f2799ff4466b6d.svg]]