make stocproc happy

This commit is contained in:
Valentin Boettcher 2022-07-11 16:59:15 +02:00
parent 109903d184
commit b56f4e12ac
2 changed files with 26 additions and 34 deletions

View file

@ -146,8 +146,8 @@ def anti_zeno_engine(
k_max=4,
bcf_terms=[terms, terms],
truncation_scheme="simplex",
driving_process_tolerances=[StocProcTolerances(1e-3, 1e-3)] * 2,
thermal_process_tolerances=[StocProcTolerances(1e-3, 1e-3)] * 2,
driving_process_tolerances=[StocProcTolerances(1e-4, 1e-4)] * 2,
thermal_process_tolerances=[StocProcTolerances(1e-4, 1e-4)] * 2,
T=[T_c, T_h],
L=[L] * 2,
H=H,
@ -178,12 +178,12 @@ def anti_zeno_engine(
Δ,
(τ_mod, τ_c, τ_bath, cycles, model.ω_s, ω_0, τ_s, τ_off, n, Δ_switch, τ_init),
) = anti_zeno_engine(
Δ=10,
Δ=15,
ε=1,
ω_c=1,
ω_c=2,
ε_couple=0.2,
n=6,
detune=-0.9,
detune=-2,
ω_0=20,
T_c=8,
T_h=40,
@ -220,7 +220,7 @@ vs = np.linspace(0.1, 10, 100)
plt.plot(vs, chi(vs, ω_0))
plt.plot(vs, G_h(vs))
aux.integrate(model, 1000)
aux.integrate(model, 1)
# _, ax = fs.plot_energy_overview(model, markersize=1, ensemble_args=dict(gc_sleep=0.05))

View file

@ -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-07-11_13-58-20_276114_35102/sockets/plasma_store', 'raylet_socket_name': '/tmp/ray/session_2022-07-11_13-58-20_276114_35102/sockets/raylet', 'webui_url': '', 'session_dir': '/tmp/ray/session_2022-07-11_13-58-20_276114_35102', 'metrics_export_port': 63919, 'gcs_address': '141.30.17.225:60837', 'address': '141.30.17.225:60837', 'node_id': '282b42830d5ff26f04c6280c3cd20e672efbdb52867015623c35edde'})
: 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_16-54-20_590904_33936/sockets/plasma_store', 'raylet_socket_name': '/tmp/ray/session_2022-07-11_16-54-20_590904_33936/sockets/raylet', 'webui_url': '', 'session_dir': '/tmp/ray/session_2022-07-11_16-54-20_590904_33936', 'metrics_export_port': 61718, 'gcs_address': '141.30.17.225:55788', 'address': '141.30.17.225:55788', 'node_id': '714c1cafb7153ccef4b9cc0a9481d02e096e4cf1720657c269036390'})
#+begin_src jupyter-python :results none
from hops.util.logging_setup import logging_setup
@ -162,8 +162,8 @@ Init ray and silence stocproc.
k_max=4,
bcf_terms=[terms, terms],
truncation_scheme="simplex",
driving_process_tolerances=[StocProcTolerances(1e-3, 1e-3)] * 2,
thermal_process_tolerances=[StocProcTolerances(1e-3, 1e-3)] * 2,
driving_process_tolerances=[StocProcTolerances(1e-4, 1e-4)] * 2,
thermal_process_tolerances=[StocProcTolerances(1e-4, 1e-4)] * 2,
T=[T_c, T_h],
L=[L] * 2,
H=H,
@ -198,12 +198,12 @@ Init ray and silence stocproc.
Δ,
(τ_mod, τ_c, τ_bath, cycles, model.ω_s, ω_0, τ_s, τ_off, n, Δ_switch, τ_init),
) = anti_zeno_engine(
Δ=10,
Δ=15,
ε=1,
ω_c=1,
ω_c=2,
ε_couple=0.2,
n=6,
detune=-0.9,
detune=-2,
ω_0=20,
T_c=8,
T_h=40,
@ -221,6 +221,7 @@ Init ray and silence stocproc.
#+end_src
#+RESULTS:
: 2c065e45-b64f-4824-93e5-b627a6e8db14
Let's test the assumptions of the paper.
@ -229,7 +230,7 @@ Let's test the assumptions of the paper.
#+end_src
#+RESULTS:
: 39
: 74fc4e96-b52e-4719-ba95-715290a68745
** BCFs and Modulations
#+begin_src jupyter-python :tangle nil
@ -244,10 +245,7 @@ Let's test the assumptions of the paper.
#+end_src
#+RESULTS:
:RESULTS:
| <matplotlib.lines.Line2D | at | 0x7fc84443f850> |
[[file:./.ob-jupyter/99eef71824c4bba4e10908bdebb9f52cb2ec3cc9.svg]]
:END:
: 10a2c713-bc80-47a6-be8f-b154b9485736
#+begin_src jupyter-python :tangle nil
ωs = np.linspace(0.01, 4*ω_0, 10000)
@ -259,14 +257,11 @@ 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(20, 30)
plt.xlim(30, 40)
#+end_src
#+RESULTS:
:RESULTS:
| <matplotlib.lines.Line2D | at | 0x7fc880d29d30> |
[[file:./.ob-jupyter/aade01d9e60146b2c865140bab73ec1a58a82b0a.svg]]
:END:
: 788c9bc5-86c5-414c-a47b-457bd878c63f
#+begin_src jupyter-python
ts = np.linspace(0,50,1000)
@ -318,21 +313,18 @@ Let's test the assumptions of the paper.
** TODO Integration
#+begin_src jupyter-python
aux.integrate(model, 1000)
aux.integrate(model, 1)
#+end_src
#+RESULTS:
#+begin_example
[INFO hops.core.integration 11515] Using 4 integrators.
[INFO hops.core.integration 11515] Some 1 trajectories have to be integrated.
[INFO hops.core.integration 11515] Using 91 hierarchy states.
100% 1/1 [01:28<00:00, 88.94s/it]
[INFO hops.core.integration 11515] Choosing the nonlinear integrator.
[INFO hops.core.integration 11515] Using 4 integrators.
[INFO hops.core.integration 11515] Some 0 trajectories have to be integrated.
[INFO hops.core.integration 11515] Using 91 hierarchy states.
0it [00:00, ?it/s]
#+end_example
: [INFO hops.core.integration 33936] Choosing the nonlinear integrator.
: [INFO hops.core.integration 33936] Using 4 integrators.
: [INFO hops.core.integration 33936] Some 1 trajectories have to be integrated.
: [INFO hops.core.integration 33936] Using 1820 hierarchy states.
: 0% 0/1 [00:00<?, ?it/s][INFO hops.core.signal_delay 33936] caught sig 'SIGINT'
: [INFO hops.core.signal_delay 33936] caught sig 'SIGINT'
: [INFO hops.core.signal_delay 33936] caught sig 'SIGINT'
: [INFO hops.core.signal_delay 33936] caught sig 'SIGINT'
#+begin_src jupyter-python