diff --git a/python/energy_flow_proper/08_dynamic_one_bath/mod_freq.py b/python/energy_flow_proper/08_dynamic_one_bath/mod_freq.py index 7e02b74..080c480 100644 --- a/python/energy_flow_proper/08_dynamic_one_bath/mod_freq.py +++ b/python/energy_flow_proper/08_dynamic_one_bath/mod_freq.py @@ -18,11 +18,11 @@ logging_setup(logging.INFO, show_stocproc=False) from hops.util.dynamic_matrix import SmoothStep, Periodic, Harmonic, ConstantMatrix -Δs = np.linspace(1,10,20)#np.sort(np.concatenate((np.linspace(1, 5, 20), np.linspace(5, 7, int(20/5 * 2))))) +Δs = np.linspace(1,10,10)#np.sort(np.concatenate((np.linspace(1, 5, 20), np.linspace(5, 7, int(20/5 * 2))))) Δ_models = [] strobe_ts = [] strobe_indices_s = [] -δs = [2, 1,.5, .1] +δs = np.linspace(.1, 2, 10) for Δ in Δs: for δ in δs: proto, strobe_t, strobe_indices = es.energy_shovel(Δ, periods=20, k_max=5, modulate_system=False) @@ -62,10 +62,10 @@ for (model, data), Δ, strobe_t, strobe_indices in zip( # ax.plot(model.t, model.L(model.t)[:,0,1]) # print(strobe_t[1]) -plt.ylabel(r"$\frac{\Delta E}{T}$") -ax.set_xlabel(r"$\Delta$") -ax.set_ylabel(r"$P_\mathrm{max}$") -for i, energy, σ in zip(itertools.count(0), final_e, final_e_error): - ax.errorbar(Δs, energy, σ, label=rf"$\alpha(0)$ = {δs[i]}") -ax.legend() -fs.export_fig("delta_dependence") +# plt.ylabel(r"$\frac{\Delta E}{T}$") +# ax.set_xlabel(r"$\Delta$") +# ax.set_ylabel(r"$P_\mathrm{max}$") +# for i, energy, σ in zip(itertools.count(0), final_e, final_e_error): +# ax.errorbar(Δs, energy, σ, label=rf"$\alpha(0)$ = {δs[i]}") +# ax.legend() +# fs.export_fig("delta_dependence") diff --git a/python/energy_flow_proper/08_dynamic_one_bath/shovel_experiments.org b/python/energy_flow_proper/08_dynamic_one_bath/shovel_experiments.org index 65e5f37..09472f2 100644 --- a/python/energy_flow_proper/08_dynamic_one_bath/shovel_experiments.org +++ b/python/energy_flow_proper/08_dynamic_one_bath/shovel_experiments.org @@ -33,6 +33,7 @@ #+end_src #+RESULTS: +: 2022-08-24 11:35:56,913 INFO worker.py:956 -- Connecting to existing Ray cluster at address: 141.30.17.16:6379 #+begin_src jupyter-python Δ = 5 @@ -70,7 +71,23 @@ : \(λ=0.1, Δ=5, T=5\) #+begin_src jupyter-python - aux.integrate_multi(models, 2000) + print(es.models_table(models)) +#+end_src + +#+RESULTS: +: \begin{tabular}{llll} +: \hline +: $ω_c$ & $2$ & $2$ & $2$ \\ +: $α(0)$ & $0.7$ & $0.7$ & $0.7$ \\ +: $T$ & $5$ & $5$ & $5$ \\ +: $N$ & $10000$ & $10000$ & $10000$ \\ +: $k_{\mathrm{max}}$ & $5$ & $5$ & $5$ \\ +: \hline +: \end{tabular} + + +#+begin_src jupyter-python + aux.integrate_multi(models, 10_000) #+end_src #+RESULTS: @@ -105,7 +122,11 @@ #+end_src #+RESULTS: -[[file:./.ob-jupyter/14ef222d0be01acd82256fc3625773767dc4d381.svg]] +:RESULTS: +: /nix/store/9qpa0hd0d43m1cdv97chak1hmls63kmv-python3-3.9.13-env/lib/python3.9/site-packages/tikzplotlib/_cleanfigure.py:166: UserWarning: Cleaning Line Collections (scatter plot) is not supported yet. +: warnings.warn( +[[file:./.ob-jupyter/ee12898f4f9fbb35ce5766124d28698de214dc25.svg]] +:END: #+begin_src jupyter-python pu.plot_energy_overview(without_friction, strobe_frequency=Δ) @@ -128,8 +149,22 @@ #+RESULTS: :RESULTS: +#+begin_example + /nix/store/9qpa0hd0d43m1cdv97chak1hmls63kmv-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/9qpa0hd0d43m1cdv97chak1hmls63kmv-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[0] = start + /nix/store/9qpa0hd0d43m1cdv97chak1hmls63kmv-python3-3.9.13-env/lib/python3.9/site-packages/matplotlib/axes/_axes.py:5225: ComplexWarning: Casting complex values to real discards the imaginary part + pts[N + 1] = end + /nix/store/9qpa0hd0d43m1cdv97chak1hmls63kmv-python3-3.9.13-env/lib/python3.9/site-packages/matplotlib/axes/_axes.py:5228: ComplexWarning: Casting complex values to real discards the imaginary part + pts[1:N+1, 1] = dep1slice + /nix/store/9qpa0hd0d43m1cdv97chak1hmls63kmv-python3-3.9.13-env/lib/python3.9/site-packages/matplotlib/axes/_axes.py:5230: ComplexWarning: Casting complex values to real discards the imaginary part + pts[N+2:, 1] = dep2slice[::-1] + /nix/store/9qpa0hd0d43m1cdv97chak1hmls63kmv-python3-3.9.13-env/lib/python3.9/site-packages/numpy/ma/core.py:3377: ComplexWarning: Casting complex values to real discards the imaginary part + _data[indx] = dval +#+end_example | hline | | (None ) | -[[file:./.ob-jupyter/878176ffdb8b483e2661b8b0c45ccd184123096f.svg]] +[[file:./.ob-jupyter/d0b99487451a16c0dd865c4ab5f79ea373b11011.svg]] :END: ** No System mod vs System Mod @@ -286,14 +321,14 @@ import scipy import pathlib - ωs = np.linspace(.5, 2.5, 6) + ωs = np.linspace(.5, 2.1, 15) ω_models = [] Δ = 5 for ω_c in ωs: proto, strobe_t, strobe_indices = es.energy_shovel( Δ, - periods=3, + periods=4, modulate_system=False, δ=0.1, ω_c = ω_c, @@ -307,11 +342,12 @@ def cost(δ): proto.δ = δ - aux.integrate(proto, 200) + aux.integrate(proto, 500) with aux.get_data(proto) as data: inter = proto.interaction_energy(data) return (abs(inter.value).max() - goal) ** 2 + cache_path = pathlib.Path(f"./.cache/{proto.hexhash}_delta_{goal}.npy") if cache_path.exists(): with cache_path.open("rb") as cache: @@ -324,7 +360,7 @@ final, strobe_t, strobe_indices = es.energy_shovel( Δ, - periods=40, + periods=50, modulate_system=False, δ=δ, ω_c=ω_c @@ -338,413 +374,39 @@ #+RESULTS: #+begin_example \(λ=0.1, Δ=5, T=5\) - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 111 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. - 100% 111/111 [00:10<00:00, 10.61it/s] - Loading: 100% 7/7 [00:00<00:00, 9.73it/s] - [INFO root 2175806] Writing cache to: results/interaction_3d6b814a4a40a7ec39b3fe256bc5e8b963b5f3a92c40ad33648103cc75a935e5_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - - Func-count x f(x) Procedure - 1 1.2077 0.0654011 initial - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:16<00:00, 12.09it/s] -Loading: 100% 7/7 [00:00<00:00, 8.85it/s] - [INFO root 2175806] Writing cache to: results/interaction_94b5a6aaee35e3476809d36916bfc614188b213a6cf4ebc3a5390e9b6febbe2f_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 2 1.8923 0.0308913 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:19<00:00, 10.43it/s] -Loading: 100% 7/7 [00:00<00:00, 8.75it/s] - [INFO root 2175806] Writing cache to: results/interaction_f2ee401c1340d1af5103ad410180519b98d9beac58bace10bacd851b03b2bc2a_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 3 2.3154 0.0161309 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:17<00:00, 11.41it/s] -Loading: 100% 7/7 [00:00<00:00, 8.64it/s] - [INFO root 2175806] Writing cache to: results/interaction_e1539c4fd1ede415ebaf31771d48e8643ae6cd0b4ea27faef191e439246c4085_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 4 2.5769 0.00942849 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:21<00:00, 9.14it/s] -Loading: 100% 7/7 [00:00<00:00, 9.00it/s] - [INFO root 2175806] Writing cache to: results/interaction_f7111f79d5a4ac86ad4d3b7dbec54e89042dc3ed40b193bde10a81ecc0bb1208_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 5 2.73851 0.00618262 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:17<00:00, 11.59it/s] -Loading: 100% 7/7 [00:00<00:00, 8.80it/s] - [INFO root 2175806] Writing cache to: results/interaction_d6f31be0c61bea2aae7919360a558441cb07bbbf6588eac9738b9e4111a968e4_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 6 2.83839 0.00452204 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:18<00:00, 11.09it/s] -Loading: 100% 7/7 [00:00<00:00, 7.59it/s] - [INFO root 2175806] Writing cache to: results/interaction_8b91035fbb9300c17ee6ad3e2185472a42ee80ddbf197879eedfbd42652e8cb0_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 7 2.90012 0.0036268 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:17<00:00, 11.51it/s] -Loading: 100% 7/7 [00:00<00:00, 8.38it/s] - [INFO root 2175806] Writing cache to: results/interaction_1bf959d0d9c03c59501e33d8841ae6e91fa501c0d1fca6e6fa2443d13f99887e_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 8 2.93827 0.00312336 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:17<00:00, 11.60it/s] -Loading: 100% 7/7 [00:00<00:00, 8.92it/s] - [INFO root 2175806] Writing cache to: results/interaction_c06ad5f302e53a3ba4d69c50f122b54c41bf053fff0b550b2135f61a56130b6a_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 9 2.96185 0.00283121 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:18<00:00, 10.81it/s] -Loading: 100% 7/7 [00:00<00:00, 8.81it/s] - [INFO root 2175806] Writing cache to: results/interaction_194c65d76110d7b5a6c9095c3078259f408b729e8fb7762409cfa223d31bd860_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 10 2.97642 0.00265789 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:17<00:00, 11.35it/s] -Loading: 100% 7/7 [00:00<00:00, 8.48it/s] - [INFO root 2175806] Writing cache to: results/interaction_2af5601f9bf69c6a7c6cba45db9ef0a9b568e6228fff64c9536cb24bf53083f3_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 11 2.98543 0.00255353 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:22<00:00, 8.96it/s] -Loading: 100% 7/7 [00:00<00:00, 8.71it/s] - [INFO root 2175806] Writing cache to: results/interaction_24e183371857095ced0bfbcea07f31b8dfff9eb864de4aae0dd17f95585d1f0c_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 12 2.99099 0.00249009 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:13<00:00, 15.05it/s] -Loading: 100% 7/7 [00:00<00:00, 8.90it/s] - [INFO root 2175806] Writing cache to: results/interaction_ec56f2bc38a257bc406d3bbf95397d5d731e173d64625b6434a2ddd9e0758441_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 13 2.99443 0.00245129 golden - - Optimization terminated successfully; - The returned value satisfies the termination criteria - (using xtol = 0.01 ) \(λ=0.1, Δ=5, T=5\) \(λ=0.1, Δ=5, T=5\) - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:11<00:00, 17.27it/s] -Loading: 100% 7/7 [00:00<00:00, 8.90it/s] - [INFO root 2175806] Writing cache to: results/interaction_fc89ddac14b1ee58802eda98ac31c051ad91495af8cdd919bd2fd8cc037c3f7a_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + [INFO hops.core.integration 205608] Choosing the nonlinear integrator. + [INFO hops.core.integration 205608] Using 21 integrators. + [INFO hops.core.integration 205608] Some 500 trajectories have to be integrated. + [INFO hops.core.integration 205608] Using 792 hierarchy states. + 100% 500/500 [00:50<00:00, 9.88it/s] + Loading: 100% 21/21 [00:02<00:00, 8.37it/s] + [INFO root 205608] Writing cache to: results/interaction_0d64c622542a492e583af382f6ccff07bcb12c2bcc8cefd7d617e4746e4269d6_interaction_energy_task_500_None_f76dc3d000fede3c15900b379894a79f5537ebdb64427f8b0299801907ab4181.npy Func-count x f(x) Procedure - 1 1.2077 0.0181862 initial - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:14<00:00, 13.41it/s] -Loading: 100% 7/7 [00:00<00:00, 9.02it/s] - [INFO root 2175806] Writing cache to: results/interaction_c383c6db04f79a2b79e442cc1093fe15f3d9b951943ed5306f4b06714fd1bd65_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 2 1.8923 0.000113858 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:12<00:00, 16.28it/s] -Loading: 100% 7/7 [00:00<00:00, 8.84it/s] - [INFO root 2175806] Writing cache to: results/interaction_1962605cecb7d5d4adc2c3f9f2df5383968ae7115126d9249fdc19cb55c7bbc4_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 3 2.3154 0.00986562 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:18<00:00, 11.06it/s] -Loading: 100% 7/7 [00:00<00:00, 9.08it/s] - [INFO root 2175806] Writing cache to: results/interaction_639b1e85650272ff0b1fd8ac5ce398e3d9470ba2b01e75cc7b853f33e8786d69_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 4 1.84569 7.24876e-07 parabolic - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:15<00:00, 12.66it/s] -Loading: 100% 7/7 [00:00<00:00, 8.87it/s] - [INFO root 2175806] Writing cache to: results/interaction_844a7bf0877c61239058b28081d76775298b2cb586b4948cd90340b1eb2c250e_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 5 1.84135 4.02743e-09 parabolic - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:12<00:00, 16.22it/s] -Loading: 100% 7/7 [00:00<00:00, 8.34it/s] - [INFO root 2175806] Writing cache to: results/interaction_b8451128d3d3052da1f36c92a6e6cc66abbb3975502883ca2ef8400ddc423229_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 6 1.83802 5.87104e-07 parabolic - - Optimization terminated successfully; - The returned value satisfies the termination criteria - (using xtol = 0.01 ) - \(λ=0.1, Δ=5, T=5\) - \(λ=0.1, Δ=5, T=5\) - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:11<00:00, 17.99it/s] -Loading: 100% 7/7 [00:00<00:00, 8.83it/s] - [INFO root 2175806] Writing cache to: results/interaction_6bd9835068b25b2c102a75a18461f638e93747868b77715824c5c98985faaf93_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - - Func-count x f(x) Procedure - 1 1.2077 0.00113026 initial - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:17<00:00, 11.34it/s] -Loading: 100% 7/7 [00:00<00:00, 9.12it/s] - [INFO root 2175806] Writing cache to: results/interaction_7df854ae50a57fa03e7fcca4c48f0dd66187944217190ebaceb43fa3a8e2763e_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 2 1.8923 0.0287167 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:17<00:00, 11.54it/s] -Loading: 100% 7/7 [00:00<00:00, 8.77it/s] - [INFO root 2175806] Writing cache to: results/interaction_0faf4ca4557dc1a61dad8255619b8a3f3a78087194b844d1607a587eb7fbaec8_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 3 0.784597 0.0256398 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:14<00:00, 13.72it/s] -Loading: 100% 7/7 [00:00<00:00, 9.12it/s] - [INFO root 2175806] Writing cache to: results/interaction_01c475ecd620d978eddc7ebbbfebe7b03bb5c8746c2503f8efb0f81784f90e49_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 4 1.32279 4.11683e-07 parabolic - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:10<00:00, 18.63it/s] -Loading: 100% 7/7 [00:00<00:00, 8.99it/s] - [INFO root 2175806] Writing cache to: results/interaction_3ebac6753d288b31249005f09a5a31126585f7bad68164d70d53251edb3cc8a1_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 5 1.31945 1.22156e-07 parabolic - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:18<00:00, 10.61it/s] -Loading: 100% 7/7 [00:00<00:00, 8.37it/s] - [INFO root 2175806] Writing cache to: results/interaction_d0ed5cf72e278b10cd19099ff86d6247f03135344961ae240fbddd5e3875e650_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 6 1.31612 1.79751e-06 parabolic - - Optimization terminated successfully; - The returned value satisfies the termination criteria - (using xtol = 0.01 ) - \(λ=0.1, Δ=5, T=5\) - \(λ=0.1, Δ=5, T=5\) - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:10<00:00, 18.86it/s] -Loading: 100% 7/7 [00:00<00:00, 8.63it/s] - [INFO root 2175806] Writing cache to: results/interaction_1362a2eaf290edd955fd6194376db340caefe3be0892210770d60e37cc7ff09a_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - - Func-count x f(x) Procedure - 1 1.2077 0.0100288 initial - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:10<00:00, 18.44it/s] -Loading: 100% 7/7 [00:00<00:00, 9.20it/s] - [INFO root 2175806] Writing cache to: results/interaction_1e4e7a01f41028445031de0ad7a1158bbac5d5e7694a7f7de2b6b39f9bcb4641_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 2 1.8923 0.126393 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:17<00:00, 11.41it/s] -Loading: 100% 7/7 [00:00<00:00, 8.99it/s] - [INFO root 2175806] Writing cache to: results/interaction_3c18810832478fd1f709a56e3e7f2998c69230747e7f61d8c3fc1e2d43239f22_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 3 0.784597 0.00420412 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:14<00:00, 13.76it/s] -Loading: 100% 7/7 [00:00<00:00, 9.16it/s] - [INFO root 2175806] Writing cache to: results/interaction_ed769efd49a9e4b11cbcb308a4aea1d898ab0624fce070c9ce353626612eed55_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 4 0.947339 2.26549e-07 parabolic - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:18<00:00, 10.99it/s] -Loading: 100% 7/7 [00:00<00:00, 8.30it/s] - [INFO root 2175806] Writing cache to: results/interaction_44b0aba0e7df72cdc9885474420f085b17a43de542b4bcdc4ef7d8c2c7d50a36_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 5 0.950891 8.40922e-07 parabolic - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:18<00:00, 11.06it/s] -Loading: 100% 7/7 [00:01<00:00, 6.77it/s] - [INFO root 2175806] Writing cache to: results/interaction_b3bb2ba793b1c1200bbaa39fcc4b412ba8048feb26907908964904ccc3b9f3d2_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 6 0.944006 3.18144e-06 parabolic - - Optimization terminated successfully; - The returned value satisfies the termination criteria - (using xtol = 0.01 ) - \(λ=0.1, Δ=5, T=5\) - \(λ=0.1, Δ=5, T=5\) - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:10<00:00, 19.11it/s] -Loading: 100% 7/7 [00:00<00:00, 8.98it/s] - [INFO root 2175806] Writing cache to: results/interaction_51062b294c91b4302a81e75710c316e053b18138996ceae34eccea46c036b65b_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - - Func-count x f(x) Procedure - 1 1.2077 0.056141 initial - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:15<00:00, 12.56it/s] -Loading: 100% 7/7 [00:00<00:00, 9.16it/s] - [INFO root 2175806] Writing cache to: results/interaction_e35c886dee85b36a5a9d21230b0b567950bd4dbf962fb8be7afd9edc5aeb8b76_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 2 1.8923 0.309729 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:17<00:00, 11.44it/s] -Loading: 100% 7/7 [00:00<00:00, 9.15it/s] - [INFO root 2175806] Writing cache to: results/interaction_8e4dffe69139a934526015d4e0816e9d3d13b2dc9c6348a1a237f3370bacd8a7_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 3 0.784597 0.00076149 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:10<00:00, 19.26it/s] -Loading: 100% 7/7 [00:00<00:00, 9.58it/s] - [INFO root 2175806] Writing cache to: results/interaction_b33af14e2b7d088ceb178c89ea5039b5ba892253fc66fb5ca8d509c408934b71_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 4 0.693505 0.000368138 parabolic - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:16<00:00, 12.29it/s] -Loading: 100% 7/7 [00:01<00:00, 6.31it/s] - [INFO root 2175806] Writing cache to: results/interaction_e3acb0a0eb86c95d221f44d6cc3c1f1bfc29176fe579f04444861e49668956c8_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 5 0.73028 4.85496e-08 parabolic - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:18<00:00, 10.88it/s] -Loading: 100% 7/7 [00:00<00:00, 8.41it/s] - [INFO root 2175806] Writing cache to: results/interaction_a87ad84893a778283058b1ae985a923922d4640fb39b3c4f2d9aa0cc0b83ce21_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 6 0.733613 2.2303e-06 parabolic - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:15<00:00, 13.11it/s] -Loading: 100% 7/7 [00:00<00:00, 8.39it/s] - [INFO root 2175806] Writing cache to: results/interaction_38b796ef49ae7155eebac7efe7333657c4cfbf1a57866afbb6260d12e0fe7ee4_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 7 0.726947 3.7442e-06 parabolic - - Optimization terminated successfully; - The returned value satisfies the termination criteria - (using xtol = 0.01 ) - \(λ=0.1, Δ=5, T=5\) - \(λ=0.1, Δ=5, T=5\) - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:11<00:00, 17.25it/s] -Loading: 100% 7/7 [00:01<00:00, 6.50it/s] - [INFO root 2175806] Writing cache to: results/interaction_8c1d054ba36e00ae216460bec970cf00b44f3078e283cf4bae4b3c8a3c9f81bc_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - - Func-count x f(x) Procedure - 1 1.2077 0.0869517 initial - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:11<00:00, 16.91it/s] -Loading: 100% 7/7 [00:00<00:00, 8.17it/s] - [INFO root 2175806] Writing cache to: results/interaction_6f192fd0cf946684c8d564f7833e28457086f75b2b4322741c6a959c7dc727f6_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 2 1.8923 0.416013 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:14<00:00, 13.52it/s] -Loading: 100% 7/7 [00:00<00:00, 9.03it/s] - [INFO root 2175806] Writing cache to: results/interaction_5702c098b5eec0490e8b5f6a6614ef5dbea3e4813b39d4fe9ff8330fd740b923_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 3 0.784597 0.00411023 golden - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:18<00:00, 10.94it/s] -Loading: 100% 7/7 [00:00<00:00, 9.25it/s] - [INFO root 2175806] Writing cache to: results/interaction_3c96ae3a6d958e60a27ca7011d7df9150316911bf2ee294a05f61c4f4e3bbc5b_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 4 0.615481 0.00101504 parabolic - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:15<00:00, 12.97it/s] -Loading: 100% 7/7 [00:00<00:00, 8.93it/s] - [INFO root 2175806] Writing cache to: results/interaction_bf9f0303b85efeae264ef318c82372e1ce942062fd45d2dc0e2252cc3ae9af94_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 5 0.669505 9.06649e-07 parabolic - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:13<00:00, 14.53it/s] -Loading: 100% 7/7 [00:00<00:00, 8.53it/s] - [INFO root 2175806] Writing cache to: results/interaction_b580450f0947a7a5349a6afc3fa2b2f590f7bd75962009baca9a8ac81ce552a4_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 6 0.672839 8.97062e-07 parabolic - [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. - [INFO hops.core.integration 2175806] Using 21 integrators. - [INFO hops.core.integration 2175806] Some 200 trajectories have to be integrated. - [INFO hops.core.integration 2175806] Using 792 hierarchy states. -100% 200/200 [00:11<00:00, 16.86it/s] -Loading: 100% 7/7 [00:00<00:00, 8.92it/s] - [INFO root 2175806] Writing cache to: results/interaction_63ffd1b4c435def7bb61d3f2bdec62c397d380e17129b96908e308ce659cc7cf_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy - 7 0.676172 8.09719e-06 parabolic - - Optimization terminated successfully; - The returned value satisfies the termination criteria - (using xtol = 0.01 ) - \(λ=0.1, Δ=5, T=5\) + 1 1.2077 0.0482932 initial + [INFO hops.core.integration 205608] Choosing the nonlinear integrator. + [INFO hops.core.integration 205608] Using 21 integrators. + [INFO hops.core.integration 205608] Some 500 trajectories have to be integrated. + [INFO hops.core.integration 205608] Using 792 hierarchy states. +100% 500/500 [00:47<00:00, 10.55it/s] +Loading: 100% 21/21 [00:02<00:00, 7.56it/s] + [INFO root 205608] Writing cache to: results/interaction_403926748d4ac79cfaafc734fdc421404835cb84e1156f980e394a6860dbbb59_interaction_energy_task_500_None_f76dc3d000fede3c15900b379894a79f5537ebdb64427f8b0299801907ab4181.npy + 2 1.8923 0.0143301 golden + [INFO hops.core.integration 205608] Choosing the nonlinear integrator. + [INFO hops.core.integration 205608] Using 21 integrators. + [INFO hops.core.integration 205608] Some 500 trajectories have to be integrated. + [INFO hops.core.integration 205608] Using 792 hierarchy states. +100% 500/500 [00:50<00:00, 9.88it/s] +Loading: 100% 21/21 [00:02<00:00, 7.49it/s] + [INFO root 205608] Writing cache to: results/interaction_38d5b05e454d6ee954c83a2924fb52122f2858b4578bece15efc5845e6fcc946_interaction_energy_task_500_None_f76dc3d000fede3c15900b379894a79f5537ebdb64427f8b0299801907ab4181.npy + 3 2.3154 0.00343145 golden + [INFO hops.core.integration 205608] Choosing the nonlinear integrator. + [INFO hops.core.integration 205608] Using 21 integrators. + [INFO hops.core.integration 205608] Some 500 trajectories have to be integrated. + [INFO hops.core.integration 205608] Using 792 hierarchy states. +100% 500/500 [00:48<00:00, 10.28it/s] #+end_example #+begin_src jupyter-python :tangle night.py @@ -752,11 +414,14 @@ Loading: 100% 7/7 [00:00<00:00, 8.92it/s] #+end_src #+RESULTS: -: [INFO hops.core.integration 2175806] Choosing the nonlinear integrator. -: [INFO hops.core.integration 2175806] Using 21 integrators. -: [INFO hops.core.integration 2175806] Some 10000 trajectories have to be integrated. -: [INFO hops.core.integration 2175806] Using 792 hierarchy states. -: 26% 2564/10000 [26:23<1:16:32, 1.62it/s] + + +#+begin_src jupyter-python + print(es.models_table(ω_models)) +#+end_src + +#+RESULTS: +: 712fb5f5-49f4-4383-9e5c-13eca64a73ca #+begin_src jupyter-python :tangle night.py from pprint import pprint @@ -784,32 +449,33 @@ Loading: 100% 7/7 [00:00<00:00, 8.92it/s] #+RESULTS: :RESULTS: -: -[[file:./.ob-jupyter/f4892164e6434a0713ec3095b99b9bce035d0ad7.svg]] +: +[[file:./.ob-jupyter/30920399a6c7f58b871c27cdd51448731ea68388.svg]] :END: #+begin_src jupyter-python for model in ω_models: - vals = abs(model.full_thermal_bcf(model.t)) - plt.plot(model.t, vals / vals.max(), linewidth=1) + vals = abs(model.bcf(model.t)) + plt.plot(model.t, model.bcf_scale * vals, linewidth=1) plt.xlim(-.1, 6) #+end_src #+RESULTS: -[[file:./.ob-jupyter/9ee1e4c90f02ba70df29cd4a6986f8910bc9d559.svg]] +[[file:./.ob-jupyter/e5ef6781d602578992dff67446244729599a952b.svg]] #+begin_src jupyter-python ωs = np.linspace(.01, 30, 1000) for model in ω_models: - plt.plot(ωs, model.full_thermal_spectral_density(ωs), linewidth=1, label=fr"$\omega_c={model.ω_c:.2f}$") + plt.plot(ωs, model.bcf_scale * model.spectral_density(ωs), linewidth=1, label=fr"$\omega_c={model.ω_c:.2f}$") plt.legend() + plt.xlabel(r"$\omega$") + plt.ylabel(r"$J(\omega)$") + + fs.export_fig("omega_sd", tikz=False, x_scaling=.49, y_scaling=.3) #+end_src #+RESULTS: -:RESULTS: -: -[[file:./.ob-jupyter/61dc5eced82b18a3290de486f000e10bb5e30063.svg]] -:END: +[[file:./.ob-jupyter/c6f8a23bf778858fc9c1a4469706922a5ec5106d.svg]] The time scale of the bath and the modulation is also interesting. @@ -830,6 +496,7 @@ The time scale of the bath and the modulation is also interesting. *** Energy stuff #+begin_src jupyter-python fig, ax = plt.subplots() + ax.axhline(0, color="gray") for model, data in aux.model_data_iterator(ω_models): _, _, bar = pu.plot_with_σ( model.t, @@ -837,19 +504,27 @@ The time scale of the bath and the modulation is also interesting. ax=ax, label=fr"$\omega_c={model.ω_c:.2f}$", ) + ax.set_xlim(0, model.t[strobe_indices][3]) + ax.set_ylabel(r"$\langle H_\mathrm{I}\rangle$") + ax.set_xlabel(r"$\tau$") + ax.legend() + #fs.export_fig("omega_interactions", fig, tikz=False, x_scaling=.49, y_scaling=.3) #+end_src #+RESULTS: -: 3a394dbb-cf59-4f13-9cbf-066c05b681b1 +: 3ca75dfb-222a-4c19-8741-0a0200bc2f32 +Total Energy Extraction #+begin_src jupyter-python fig, ax = plt.subplots() import scipy + for model, data in aux.model_data_iterator(ω_models): + energy = model.total_energy_from_power(data) * (1 / es.ergo(model.T)) _, _, bar = pu.plot_with_σ( model.t, - model.total_energy_from_power(data) * (1/es.ergo(model.T)), + energy, ax=ax, label=fr"$\omega_c={model.ω_c:.2f}$", strobe_frequency=Δ, @@ -858,7 +533,7 @@ The time scale of the bath and the modulation is also interesting. pu.plot_with_σ( model.t, - (model.total_energy_from_power(data) * (1/es.ergo(model.T))), + energy, ax=ax, # label=fr"$ω_c={model.ω_c:.2f}$", # strobe_frequency=Δ, @@ -867,32 +542,121 @@ The time scale of the bath and the modulation is also interesting. color=bar.lines[0].get_color(), ) - τ_off = 1 / (model.bcf_coefficients()[1][0].real.min()) - a0 = abs(model.bcf(0)) - τ_off = scipy.optimize.fsolve(lambda τ: abs(model.bcf(τ)) - (a0/300), 10) + # τ_off = 1 / (model.bcf_coefficients()[1][0].real.min()) + # # a0 = abs(model.bcf(0)) + # # τ_off = scipy.optimize.fsolve(lambda τ: abs(model.bcf(τ)) - (a0/300), 10) + # ax.axvline( + # τ_off, + # color=bar.lines[0].get_color(), + # linestyle="dotted", + # linewidth=1, + # zorder=-1, + # ) + + ts = model.t[strobe_indices] + e_s = energy.value[strobe_indices] + + + ind_min = e_s.argmin() + t_min = ts[ind_min] + e_min = e_s[ind_min] + ax.axvline( - τ_off, + t_min, color=bar.lines[0].get_color(), + ymin=0, + ymax=abs((e_min + 0.41) / (0.45 + 0.1)), linestyle="dotted", - linewidth=1, - zorder=-1, ) + ax.axhline( + e_min, + color=bar.lines[0].get_color(), + xmin=0, + xmax=t_min / model.t.max(), + linestyle="dotted", + ) + ax.set_xlabel(r"$\tau$") ax.set_ylabel(r"$(\langle H\rangle_\tau -\langle H\rangle_0)/\mathcal{W}_\mathrm{max}$") # ax.axhline(, color="grey", linewidth=1, linestyle="dashed") # ax.plot(model.t, model.L.operator_norm(model.t), linewidth=1) + ax.set_xlim(0, model.t.max()) + ax.set_ylim(-0.41, 0.1) ax.legend() - #fs.export_fig("omegas_total") + fs.export_fig("omegas_total", tikz=False, y_scaling=0.4) #+end_src #+RESULTS: -: fd39097d-e2f5-466b-86a2-4979c5fb9726 +: 914dca52-374e-4b3d-ae08-e6f330d76c5e - #+begin_src jupyter-python +Now we plot max power vs memory time. +#+begin_src jupyter-python + import scipy + + opt_powers = [] + opt_energies = [] + memories = [] + for model, data in aux.model_data_iterator(ω_models): + energy = model.total_energy_from_power(data) + ts = model.t[strobe_indices] + e_s = energy.value[strobe_indices] + + powers = np.divide(e_s, ts, where=ts > 0) + p_min = powers.min() + ind_min = powers.argmin() + # t_min = ts[ind_min] + # e_min = e_s[ind_min] + + opt_powers.append(abs(p_min)) + opt_energies.append(abs(e_s.min()) * (1 / es.ergo(model.T))) + + # bcf = model.full_thermal_bcf + # α0 = bcf(0).real + # obj = lambda τ: (abs(bcf(τ))/α0 - 1e-1) ** 2 + # τ_off = scipy.optimize.basinhopping( + # obj, 1 + # ).x + memories.append(3 / model.ω_c) + # memories.append(τ_off) + + fig, ax = plt.subplots() + lines = ax.plot(memories, opt_powers, marker="o", markersize=3) + ax.set_xlabel(r"$\tau_\mathrm{B}=3 / \omega_c$") + ax.set_ylabel(r"$P_\mathrm{max}$", color=lines[0].get_color()) + + ax2 = ax.twinx() + + ax2.plot([], []) + lines = ax2.plot(memories, opt_energies, marker="o", markersize=3) + ax2.set_ylabel( + r"$|\Delta \langle H\rangle|_\mathrm{max} / \mathcal{W}_\mathrm{max}$", + color=lines[0].get_color(), + ) + + # def f(t, s, h): + # t = t-s + # out = np.zeros_like(t) + # np.log(t, out=out, where=t>1) + + # return out * h + + # opt,_ = scipy.optimize.curve_fit(f, memories, opt_powers) + # smooth_mem = np.linspace(memories[0], memories[-1], 1000) + # ax.plot(smooth_mem, f(smooth_mem, *opt)) # + plt.axvline(1 * np.pi * 2/Δ, color="gray", linestyle="dotted") + fs.export_fig("omega_energies_and_powers", fig, tikz=False, x_scaling=.5, y_scaling=.35) +#+end_src + +#+RESULTS: +: 839b9591-0955-409b-9e52-ab3f0656fed3 + + + +#+begin_src jupyter-python fig, ax = plt.subplots() for model, data in aux.model_data_iterator(ω_models): - fs.plot_with_σ( + pu.plot_with_σ( model.t, model.system_energy(data), ax=ax, @@ -902,29 +666,6 @@ The time scale of the bath and the modulation is also interesting. ax.legend() #+end_src -#+RESULTS: -:RESULTS: -# [goto error] -#+begin_example - --------------------------------------------------------------------------- - AttributeError Traceback (most recent call last) - Input In [46], in () -  1 fig, ax = plt.subplots() -  2 for model, data in aux.model_data_iterator(ω_models): - ----> 3 fs.plot_with_σ( -  4 model.t, -  5 model.system_energy(data), -  6 ax=ax, -  7 label=fr"$ω_c={model.ω_c:.2f}$", -  8 strobe_frequency=Δ, -  9 ) -  10 ax.legend() - - AttributeError: module 'figsaver' has no attribute 'plot_with_σ' -#+end_example -[[file:./.ob-jupyter/d617a8bf271a07559fa626d2691271db7dcf8ce4.svg]] -:END: - #+begin_src jupyter-python @@ -943,20 +684,23 @@ The time scale of the bath and the modulation is also interesting. #+begin_src jupyter-python fig, ax = plt.subplots() int_0 = ω_models[0].interaction_energy(aux.get_data(ω_models[0])).for_bath(0) - for model, data in aux.model_data_iterator(ω_models): - fs.plot_with_σ( + for model, data in aux.model_data_iterator(ω_models[1:]): + pu.plot_with_σ( model.t, - (model.interaction_energy(data) - int_0) * (1 / abs(int_0).max.value), + (abs(model.interaction_energy(data)) - abs(int_0)) * (1 / abs(int_0).max.value), ax=ax, bath=0, - label=fr"$ω_c={model.ω_c:.2f}$", + label=fr"$\omega_c={model.ω_c:.2f}$", linewidth=0.5, ) - #ax.legend() + ax.legend() #+end_src #+RESULTS: -[[file:./.ob-jupyter/dada909ebf98f84e3ab6ff83510d6e169fe46d1e.svg]] +:RESULTS: +: +[[file:./.ob-jupyter/a70df5465a7e39d80dd2ffbbc3b7863cbafe9e1c.svg]] +:END: - pretty similar - substantial :P @@ -1105,21 +849,55 @@ The time scale of the bath and the modulation is also interesting. :END: *** State stuff +#+begin_src jupyter-python + from hops.util.utilities import relative_entropy +#+end_src + +#+RESULTS: + #+begin_src jupyter-python fig, ax = plt.subplots() for model,data in aux.model_data_iterator(ω_models): - fs.plot_with_σ(model.t, EnsembleValue(relative_entropy(data, strobe_indices[-1])), ax=ax, label=model.ω_c, strobe_frequency=Δ) + ax.plot(model.t, relative_entropy(data, strobe_indices[-1])[0], ax=ax, label=model.ω_c, strobe_frequency=Δ) ax.legend() #+end_src #+RESULTS: :RESULTS: -: /nix/store/l77qg6h691cglkdhnl4bn06f7flrwkr0-python3-3.9.13-env/lib/python3.9/site-packages/scipy/linalg/_matfuncs_inv_ssq.py:827: LogmExactlySingularWarning: The logm input matrix is exactly singular. -: warnings.warn(exact_singularity_msg, LogmExactlySingularWarning) -: /home/hiro/src/hops/hops/util/utilities.py:114: RuntimeWarning: invalid value encountered in sqrt -: np.sqrt(np.trace(_square_mat(Δρ @ (logm2(ρ) + 1 - logm2(ρ_ref))).real)) -: -[[file:./.ob-jupyter/db58e138f71c5d6e18f644e3090e6f3fb6dd5863.svg]] +# [goto error] +#+begin_example + --------------------------------------------------------------------------- + RuntimeError Traceback (most recent call last) + Input In [10], in () +  1 fig, ax = plt.subplots() + ----> 2 for model,data in aux.model_data_iterator(ω_models): +  3 ax.plot(model.t, relative_entropy(data, strobe_indices[-1])[0], ax=ax, label=model.ω_c, strobe_frequency=Δ) +  4 ax.legend() + + File ~/src/two_qubit_model/hiro_models/model_auxiliary.py:160, in model_data_iterator(models, *args, **kwargs) +  152 """ +  153 Yields tuples of ``model, data``, where ``data`` is already opened +  154 and will be closed automatically. +  155 +  156 For the rest of the arguments see :any:`get_data`. +  157 """ +  159 for model in models: + --> 160 with get_data(model, *args, **kwargs) as data: +  161 yield model, data + + File ~/src/two_qubit_model/hiro_models/model_auxiliary.py:146, in get_data(model, data_path, read_only, **kwargs) +  135 return HIData( +  136 path, +  137 hi_key=model.hops_config, +  (...) +  142 **kwargs, +  143 ) +  145 else: + --> 146 raise RuntimeError(f"No data found for model with hash '{hexhash}'.") + + RuntimeError: No data found for model with hash 'b0c0e289e5e8222d53973e3d565a0dffb4d53dd4e82b04a09d8b5b36f1f12e9e'. +#+end_example +[[file:./.ob-jupyter/392bec586936678b5b4d8a88746cf9458c91ce03.svg]] :END: Whole different steady state -> phase transition? @@ -1173,6 +951,758 @@ Whole different steady state -> phase transition? *** TODO Find out what happens when we use the same peak heights +*** Optimization for weaker coupling +:PROPERTIES: +:header-args: :tangle cutoff_dep_weak.py :session 08_shovel_cutoff_dep_weak :async yes :noweb yes +:END: + +#+begin_src jupyter-python :noweb yes + <> +#+end_src + +#+RESULTS: +: 2022-08-24 11:56:14,755 INFO worker.py:956 -- Connecting to existing Ray cluster at address: 141.30.17.16:6379 + +- found that with low temperature only a small cutoff produces a net + energy gain + +#+begin_src jupyter-python :tangle night.py + import scipy + import pathlib + + ωs = np.linspace(.5, 2.5, 6) + ω_models = [] + Δ = 5 + + for ω_c in ωs: + proto, strobe_t, strobe_indices = es.energy_shovel( + Δ, + periods=3, + modulate_system=False, + δ=0.01, + ω_c = ω_c, + ) + + proto.k_max = 5 + proto.therm_method = "fft" + + goal = 0.01 + + + def cost(δ): + proto.δ = δ + aux.integrate(proto, 200) + with aux.get_data(proto) as data: + inter = proto.interaction_energy(data) + + return (abs(inter.value).max() - goal) ** 2 + cache_path = pathlib.Path(f"./.cache/{proto.hexhash}_delta_{goal}.npy") + if cache_path.exists(): + with cache_path.open("rb") as cache: + δ = np.load(cache, allow_pickle=True) + else: + δ = scipy.optimize.minimize_scalar(cost, bounds=(.1, 3), method="bounded", options=dict(xatol=1e-2, disp=3)).x + cache_path.parent.mkdir(parents=True, exist_ok=True) + with cache_path.open("wb") as cache: + np.save(cache, δ) + + final, strobe_t, strobe_indices = es.energy_shovel( + Δ, + periods=40, + modulate_system=False, + δ=δ, + ω_c=ω_c + ) + final.k_max = 5 + final.therm_method = "fft" + + ω_models.append(final.copy()) +#+end_src + +#+RESULTS: +#+begin_example + \(λ=0.1, Δ=5, T=5\) + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. + 0it [00:00, ?it/s] + + Func-count x f(x) Procedure + 1 1.2077 0.0180267 initial + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + 2 1.8923 0.0458991 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:18<00:00, 10.69it/s] +Loading: 100% 7/7 [00:00<00:00, 8.49it/s] + [INFO root 3001754] Writing cache to: results/interaction_63c415a41b49902825b954d812a2c2b8a30364ec6a55a3a1014c86b3e2157de4_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 3 0.784597 0.0070849 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:10<00:00, 18.90it/s] +Loading: 100% 7/7 [00:00<00:00, 9.54it/s] + [INFO root 3001754] Writing cache to: results/interaction_7ddc83088888eb108cabc3e5f049f5efadda793bfd4e5d50e1b33c9e8167838a_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 4 0.523104 0.00280658 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:10<00:00, 19.26it/s] +Loading: 100% 7/7 [00:00<00:00, 8.53it/s] + [INFO root 3001754] Writing cache to: results/interaction_5b8b089e7497660dcfe0cd06e6694994d487f6a8c34c53b7d8d51ffef58eba03_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 5 0.361493 0.00112972 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:10<00:00, 19.30it/s] +Loading: 100% 7/7 [00:00<00:00, 8.10it/s] + [INFO root 3001754] Writing cache to: results/interaction_99803887724d5c9d32c8437ea30bc4b2f9a4e39a09041c1df19aa14c42a3750c_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 6 0.261611 0.000466631 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:10<00:00, 19.91it/s] +Loading: 100% 7/7 [00:01<00:00, 6.13it/s] + [INFO root 3001754] Writing cache to: results/interaction_2d8cb0d378823d491ddf84cbf5fd602c4d59dd57df6f6d914067deb0c64f3505_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 7 0.199881 0.000200614 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:09<00:00, 20.42it/s] +Loading: 100% 7/7 [00:00<00:00, 8.31it/s] + [INFO root 3001754] Writing cache to: results/interaction_1e8178e38f3316b51cf94218644435d71b485725e7ab871c00787c0121d2940b_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 8 0.16173 9.14146e-05 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:17<00:00, 11.43it/s] +Loading: 100% 7/7 [00:00<00:00, 8.48it/s] + [INFO root 3001754] Writing cache to: results/interaction_13c06f8632cde3ecc3cc47fe4b9567285ef153d52529c35022c2aef7394c0a7e_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 9 0.138151 4.50809e-05 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:09<00:00, 20.29it/s] +Loading: 100% 7/7 [00:00<00:00, 9.13it/s] + [INFO root 3001754] Writing cache to: results/interaction_a7f515d0b7062fdc312c19acf47f958639b7c370f765a614eea88d7edddc5a8c_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 10 0.123579 2.45413e-05 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:15<00:00, 12.53it/s] +Loading: 100% 7/7 [00:01<00:00, 5.42it/s] + [INFO root 3001754] Writing cache to: results/interaction_d1cec2b0043f9ff291edd3436a537481168229ddf9e223528de2e2fa52f4f5f1_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 11 0.114572 1.49433e-05 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 13.84it/s] +Loading: 100% 7/7 [00:00<00:00, 8.64it/s] + [INFO root 3001754] Writing cache to: results/interaction_4d497d01ab97816b9d5a89545b6b7372949c41814dbfa75301835b537b4c279c_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 12 0.109006 1.01949e-05 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:09<00:00, 20.45it/s] +Loading: 100% 7/7 [00:00<00:00, 9.22it/s] + [INFO root 3001754] Writing cache to: results/interaction_ec2228808637db38145961fbe203067c8410205018909661e47379103f111983_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 13 0.105566 7.71256e-06 golden + + Optimization terminated successfully; + The returned value satisfies the termination criteria + (using xtol = 0.01 ) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + + Func-count x f(x) Procedure + 1 1.2077 0.0650984 initial + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + 2 1.8923 0.160537 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:18<00:00, 10.93it/s] +Loading: 100% 7/7 [00:00<00:00, 8.94it/s] + [INFO root 3001754] Writing cache to: results/interaction_265bd9e02f5dfb6f939e3fdfbe8b62a27630dee4a418304c1aa4bf9e1812842e_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 3 0.784597 0.0268188 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:17<00:00, 11.67it/s] +Loading: 100% 7/7 [00:00<00:00, 8.75it/s] + [INFO root 3001754] Writing cache to: results/interaction_5a41cbf4faa8a4291f8e521445c40c4c34fe7454b818a7c229e5ee2c2ba6ce13_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 4 0.523104 0.0113515 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:09<00:00, 21.85it/s] +Loading: 100% 7/7 [00:01<00:00, 6.60it/s] + [INFO root 3001754] Writing cache to: results/interaction_380ea4179d4729c6d96ccddf6485d3df5b84da5de6011c67d47d20f5e1d93251_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 5 0.361493 0.0050189 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:15<00:00, 12.56it/s] +Loading: 100% 7/7 [00:00<00:00, 8.59it/s] + [INFO root 3001754] Writing cache to: results/interaction_a3eca842f741a401bd0091a23e9a621177023c9539674f2b0298e9106f6291aa_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 6 0.261611 0.00236651 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 13.74it/s] +Loading: 100% 7/7 [00:00<00:00, 8.78it/s] + [INFO root 3001754] Writing cache to: results/interaction_0e593a98325c9fafb97d22158ad1facaad6fc5ca137ed7c0848d83df7bfa0441_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 7 0.199881 0.00121629 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:13<00:00, 14.39it/s] +Loading: 100% 7/7 [00:00<00:00, 9.00it/s] + [INFO root 3001754] Writing cache to: results/interaction_cc2fb0211e56e64b2e7bbb7f239b5b9ccd6afa9250f45915a4aabd3a25cd7f5a_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 8 0.16173 0.000693993 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:15<00:00, 13.19it/s] +Loading: 100% 7/7 [00:00<00:00, 8.75it/s] + [INFO root 3001754] Writing cache to: results/interaction_3cf56f8e85bdaf33518d91bbe629349725a77f5e29a61ba90da4f8a377402032_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 9 0.138151 0.000443652 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 14.04it/s] +Loading: 100% 7/7 [00:00<00:00, 8.82it/s] + [INFO root 3001754] Writing cache to: results/interaction_e1ad251d113ec55b3559794bb5ef4297ac074ee003a4187173ab86900b5c72d4_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 10 0.123579 0.000316711 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 13.69it/s] +Loading: 100% 7/7 [00:00<00:00, 9.06it/s] + [INFO root 3001754] Writing cache to: results/interaction_19421f97d2d66ab9056056c938fdb27e224fb0f4f84e87c28883f3acd3172c0c_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 11 0.114572 0.000248892 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:08<00:00, 23.61it/s] +Loading: 100% 7/7 [00:00<00:00, 9.37it/s] + [INFO root 3001754] Writing cache to: results/interaction_8f156a6e5ca8957b5f8ae62e9af16e70bef53a58566c0cb45a6b3f532ae6ef86_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 12 0.109006 0.000211045 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:08<00:00, 22.69it/s] +Loading: 100% 7/7 [00:00<00:00, 8.60it/s] + [INFO root 3001754] Writing cache to: results/interaction_d5227d8a7ad266c1341d94471437b05bb437516a2d4dde33b3daa9bdfdc2363b_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 13 0.105566 0.000189209 golden + + Optimization terminated successfully; + The returned value satisfies the termination criteria + (using xtol = 0.01 ) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + + Func-count x f(x) Procedure + 1 1.2077 0.127007 initial + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + 2 1.8923 0.312995 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + 3 0.784597 0.0528429 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:09<00:00, 20.62it/s] +Loading: 100% 7/7 [00:00<00:00, 9.59it/s] + [INFO root 3001754] Writing cache to: results/interaction_d67dc8d25a1df6616a4f143b79e61cbfb966a8e3e7d9370583b1fd02a39f00f9_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 4 0.523104 0.0227445 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:08<00:00, 22.75it/s] +Loading: 100% 7/7 [00:00<00:00, 9.39it/s] + [INFO root 3001754] Writing cache to: results/interaction_a0dce3990864557ede695e55fa65b86bcb62a6b853613145fe7fe3d4d953fa32_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 5 0.361493 0.0103056 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 13.42it/s] +Loading: 100% 7/7 [00:00<00:00, 9.39it/s] + [INFO root 3001754] Writing cache to: results/interaction_55c1d3d9479e427e7b49aeae88045ac0cc3ca0ae26ae8abcd6cddda6e3f832b2_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 6 0.261611 0.00502798 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:16<00:00, 12.42it/s] +Loading: 100% 7/7 [00:00<00:00, 9.46it/s] + [INFO root 3001754] Writing cache to: results/interaction_d00b26def8af27c64f6254303ec287ca90c509937793c65fa1bc5bc0dbaf0544_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 7 0.199881 0.00269536 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:11<00:00, 17.33it/s] +Loading: 100% 7/7 [00:00<00:00, 9.38it/s] + [INFO root 3001754] Writing cache to: results/interaction_7bcdae3069cd1b2eb3c0bcf0ce16bd0596bfe6ff351a36b44b7bbfd4e5f09572_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 8 0.16173 0.00161191 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 13.89it/s] +Loading: 100% 7/7 [00:00<00:00, 9.42it/s] + [INFO root 3001754] Writing cache to: results/interaction_28c820848f5aeb94e721e8e6a43784ce32d7eb5ea67dce1be64608e512ea0bcf_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 9 0.138151 0.00108002 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:08<00:00, 23.85it/s] +Loading: 100% 7/7 [00:00<00:00, 9.42it/s] + [INFO root 3001754] Writing cache to: results/interaction_f224bd5222f43a6f7b9fd819aed5cf63889d2362356afc41090124642ca937a8_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 10 0.123579 0.000804111 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:13<00:00, 14.65it/s] +Loading: 100% 7/7 [00:00<00:00, 9.06it/s] + [INFO root 3001754] Writing cache to: results/interaction_4dccca9091fa8dff53947f82545961ef119dba126fade69e0013ec53b10128e0_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 11 0.114572 0.000653813 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:15<00:00, 13.28it/s] +Loading: 100% 7/7 [00:00<00:00, 9.57it/s] + [INFO root 3001754] Writing cache to: results/interaction_e118ac1347755f8bd115eaf171a4293c537a93cea3f53f94fcc0b42cbdbca012_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 12 0.109006 0.000568659 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 14.26it/s] +Loading: 100% 7/7 [00:00<00:00, 9.56it/s] + [INFO root 3001754] Writing cache to: results/interaction_38e3ae1c50348944f7698e7c04d8aefa9ce8b6010052cba303c1e5ff4f556158_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 13 0.105566 0.00051899 golden + + Optimization terminated successfully; + The returned value satisfies the termination criteria + (using xtol = 0.01 ) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + + Func-count x f(x) Procedure + 1 1.2077 0.240241 initial + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + 2 1.8923 0.555797 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + 3 0.784597 0.10573 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:09<00:00, 20.91it/s] +Loading: 100% 7/7 [00:00<00:00, 9.25it/s] + [INFO root 3001754] Writing cache to: results/interaction_d7fc1b89ad19a829331fb718d93b397264b7eb0dd6fea3bd8e032a86993d4115_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 4 0.523104 0.0478651 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:09<00:00, 21.46it/s] +Loading: 100% 7/7 [00:00<00:00, 9.07it/s] + [INFO root 3001754] Writing cache to: results/interaction_059abd637a3764be06e30427fa1df20728fb9b832c0134ed0a201966f39abcbc_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 5 0.361493 0.022689 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:12<00:00, 16.09it/s] +Loading: 100% 7/7 [00:00<00:00, 9.20it/s] + [INFO root 3001754] Writing cache to: results/interaction_806d9ec974b346bab4c0fcf4a67b484546a9104fa95213c0a8b19d5650ee9dc6_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 6 0.261611 0.0115502 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:09<00:00, 22.13it/s] +Loading: 100% 7/7 [00:00<00:00, 9.47it/s] + [INFO root 3001754] Writing cache to: results/interaction_c74f5159e0314282bb5bb4f573318c55be52ae9d5dabb61945c6154c5b397bc4_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 7 0.199881 0.00646096 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 13.69it/s] +Loading: 100% 7/7 [00:00<00:00, 9.18it/s] + [INFO root 3001754] Writing cache to: results/interaction_5e603929d7d5ed3ce5c50eb819f4081be6a5d5237f9b0482a408d186c18447b7_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 8 0.16173 0.00402601 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:13<00:00, 14.71it/s] +Loading: 100% 7/7 [00:00<00:00, 9.32it/s] + [INFO root 3001754] Writing cache to: results/interaction_7f73962ba037b71013729bb478fcc6f0994de1f2c2c609072f46c95b10fcd6f1_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 9 0.138151 0.00279997 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 14.27it/s] +Loading: 100% 7/7 [00:00<00:00, 9.23it/s] + [INFO root 3001754] Writing cache to: results/interaction_af5cca670fb2d48cf123202f1a45003eefcdd4068d1db9c6519f14510a3395d5_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 10 0.123579 0.0021506 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:15<00:00, 12.74it/s] +Loading: 100% 7/7 [00:00<00:00, 8.06it/s] + [INFO root 3001754] Writing cache to: results/interaction_8395edd87d975c9d58da04f59dbf0b6c97988d07a56c5767b1953df0f660d70b_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 11 0.114572 0.00179111 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 13.88it/s] +Loading: 100% 7/7 [00:00<00:00, 8.89it/s] + [INFO root 3001754] Writing cache to: results/interaction_397f7776d460b54982c18f746ab0accddfcfec50d0c124b142fb185983693ce4_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 12 0.109006 0.00158503 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:17<00:00, 11.56it/s] +Loading: 100% 7/7 [00:00<00:00, 9.65it/s] + [INFO root 3001754] Writing cache to: results/interaction_d35682e768d6d151ed21a0a5698ecd045f4b258bbc67aa21168b171895dd6eba_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 13 0.105566 0.00146383 golden + + Optimization terminated successfully; + The returned value satisfies the termination criteria + (using xtol = 0.01 ) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + + Func-count x f(x) Procedure + 1 1.2077 0.393055 initial + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + 2 1.8923 0.895926 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + 3 0.784597 0.174386 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 13.67it/s] +Loading: 100% 7/7 [00:00<00:00, 9.60it/s] + [INFO root 3001754] Writing cache to: results/interaction_705469e542ccb61bd42c07b6eb73c995b9b1b4c7f55f7da915b5578be9e56aa9_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 4 0.523104 0.0792455 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:19<00:00, 10.49it/s] +Loading: 100% 7/7 [00:00<00:00, 9.28it/s] + [INFO root 3001754] Writing cache to: results/interaction_8761ae6608c67df925507d996f751cad27f3c24413da822ae6289228f39123f3_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 5 0.361493 0.0377972 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 13.83it/s] +Loading: 100% 7/7 [00:00<00:00, 9.25it/s] + [INFO root 3001754] Writing cache to: results/interaction_b170b42b40a7423836efc97a9a477a5f71da9b3dea9da3e964e2e3bab2211e2c_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 6 0.261611 0.0193978 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:10<00:00, 19.06it/s] +Loading: 100% 7/7 [00:00<00:00, 9.61it/s] + [INFO root 3001754] Writing cache to: results/interaction_c470a2cf287c41e3782289024fe7b456e0538f72e9ba01bcf8063cd9f401922e_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 7 0.199881 0.0109588 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:10<00:00, 18.21it/s] +Loading: 100% 7/7 [00:00<00:00, 9.55it/s] + [INFO root 3001754] Writing cache to: results/interaction_b1fbf2a1433d530a4d9e9db8547336a4f9880c36aa1b30a69fedbcba9c503c31_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 8 0.16173 0.00690704 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 14.01it/s] +Loading: 100% 7/7 [00:00<00:00, 9.33it/s] + [INFO root 3001754] Writing cache to: results/interaction_ea19bb55e554fe8f6d03c390eaba2abeb25728bbc435a11ad780e81165f9f758_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 9 0.138151 0.00485795 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:10<00:00, 18.23it/s] +Loading: 100% 7/7 [00:00<00:00, 9.68it/s] + [INFO root 3001754] Writing cache to: results/interaction_1084c48945b1179bdbf34a244fcf88e66d952b1819efa1ccb73d103776f46344_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 10 0.123579 0.00376785 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:15<00:00, 13.02it/s] +Loading: 100% 7/7 [00:00<00:00, 9.56it/s] + [INFO root 3001754] Writing cache to: results/interaction_0142776638fbea4aa8f96e8e45af2d52f1d9c20b4c23804bd1660e561f4f524d_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 11 0.114572 0.00316207 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:19<00:00, 10.10it/s] +Loading: 100% 7/7 [00:00<00:00, 9.54it/s] + [INFO root 3001754] Writing cache to: results/interaction_d6d4bbb458339a6a32484072cd9cff3d7b48891fb504e5e539b1f3ce6510a6b6_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 12 0.109006 0.00281377 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:10<00:00, 19.16it/s] +Loading: 100% 7/7 [00:00<00:00, 9.93it/s] + [INFO root 3001754] Writing cache to: results/interaction_c6cb44c07a77a3999db1a05c6d05648e8845f27ab6890c50eb2b21117be620b4_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 13 0.105566 0.0026085 golden + + Optimization terminated successfully; + The returned value satisfies the termination criteria + (using xtol = 0.01 ) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + + Func-count x f(x) Procedure + 1 1.2077 0.469055 initial + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + 2 1.8923 1.07121 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +0it [00:00, ?it/s] + 3 0.784597 0.206217 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:17<00:00, 11.32it/s] +Loading: 100% 7/7 [00:00<00:00, 9.22it/s] + [INFO root 3001754] Writing cache to: results/interaction_e00ac21086846ce0750dd1cb231c526f577baac66481748fea1e0eaba7c6b75f_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 4 0.523104 0.0928712 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:20<00:00, 9.98it/s] +Loading: 100% 7/7 [00:00<00:00, 9.44it/s] + [INFO root 3001754] Writing cache to: results/interaction_74b250b1757fc1d17cde1f7d9bc046f0fb410bfc62336e8befa3043533be96e9_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 5 0.361493 0.0439751 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:20<00:00, 9.63it/s] +Loading: 100% 7/7 [00:00<00:00, 9.20it/s] + [INFO root 3001754] Writing cache to: results/interaction_f54c823fb2732cb62fbee6c24faed1c8c0afd854bea126e43dd3a47ac85a7c1a_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 6 0.261611 0.0224839 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:14<00:00, 13.65it/s] +Loading: 100% 7/7 [00:00<00:00, 9.34it/s] + [INFO root 3001754] Writing cache to: results/interaction_bd26e2c29b70c3a54d4390e3a379cbe5cea315ea2f9aca6619cf1b985b629b80_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 7 0.199881 0.012688 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:18<00:00, 10.92it/s] +Loading: 100% 7/7 [00:00<00:00, 9.56it/s] + [INFO root 3001754] Writing cache to: results/interaction_fa5c4f672756e0a607a206a278a94331332b8a8f838ee3396220921762fcf4e8_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 8 0.16173 0.00800075 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:11<00:00, 17.78it/s] +Loading: 100% 7/7 [00:00<00:00, 10.04it/s] + [INFO root 3001754] Writing cache to: results/interaction_54c280072b90b96ad06197bb08948cdcb93eed6b5f6ee30682c99c8719b59869_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 9 0.138151 0.00563403 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:20<00:00, 10.00it/s] +Loading: 100% 7/7 [00:00<00:00, 9.28it/s] + [INFO root 3001754] Writing cache to: results/interaction_50876017eb4018d566091e44c80ef1f62a90d0fc6ba7c9d36041691e28a6781e_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 10 0.123579 0.00437572 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:17<00:00, 11.37it/s] +Loading: 100% 7/7 [00:00<00:00, 9.03it/s] + [INFO root 3001754] Writing cache to: results/interaction_cd1023121add760f208a4fb2ffd3e204190a5577a8303f20c9a699c010e76be5_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 11 0.114572 0.00367655 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:19<00:00, 10.14it/s] +Loading: 100% 7/7 [00:00<00:00, 9.60it/s] + [INFO root 3001754] Writing cache to: results/interaction_6971f1723ab7c35c5f54e2dc9078bfdbf41f2cfe57d0035a96aa35434cc32c78_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 12 0.109006 0.00327452 golden + [INFO hops.core.integration 3001754] Choosing the nonlinear integrator. + [INFO hops.core.integration 3001754] Using 21 integrators. + [INFO hops.core.integration 3001754] Some 200 trajectories have to be integrated. + [INFO hops.core.integration 3001754] Using 792 hierarchy states. +100% 200/200 [00:19<00:00, 10.15it/s] +Loading: 100% 7/7 [00:00<00:00, 9.35it/s] + [INFO root 3001754] Writing cache to: results/interaction_37387c15aeee6fadce9ec178949f80ba01136ecae0e379dfc5f07c02f819f7ed_interaction_energy_task_200_None_1e4380ee26552450e9766cdb6c79bb8f5e5765417c0885509798d121d4719bfc.npy + 13 0.105566 0.00303758 golden + + Optimization terminated successfully; + The returned value satisfies the termination criteria + (using xtol = 0.01 ) + \(λ=0.1, Δ=5, T=5\) +#+end_example + +#+begin_src jupyter-python + for model in ω_models: + print(model.δ) +#+end_src + +#+RESULTS: +: 0.10556619830394895 +: 0.10556619830394895 +: 0.10556619830394895 +: 0.10556619830394895 +: 0.10556619830394895 +: 0.10556619830394895 + + +#+begin_src jupyter-python + ωs = np.linspace(.01, 30, 1000) + for model in ω_models: + plt.plot(ωs, model.bcf_scale * model.spectral_density(ωs), linewidth=1, label=fr"$\omega_c={model.ω_c:.2f}$") + plt.legend() + plt.xlabel(r"$\omega$") + plt.ylabel(r"Spectral Density") + + fs.export_fig("omega_sd_weak", tikz=False, x_scaling=.49, y_scaling=.4) +#+end_src + +#+RESULTS: +[[file:./.ob-jupyter/994947171f7bb6b02104068f9f07e18897411c55.svg]] + +#+begin_src jupyter-python + for model in ω_models: + plt.plot(model.t, model.bcf_scale * model.full_thermal_bcf(model.t).real, label=model.ω_c) + plt.legend() +#+end_src + +#+RESULTS: +:RESULTS: +: +[[file:./.ob-jupyter/e739b1271ce8b62aed01552f5b0705b83682b536.svg]] +:END: + ** Modulation Frequency Dependence :PROPERTIES: :header-args: :tangle mod_freq.py :session 08_mod_freq :noweb yes :async yes @@ -1188,11 +1718,11 @@ Whole different steady state -> phase transition? from hops.util.dynamic_matrix import SmoothStep, Periodic, Harmonic, ConstantMatrix - Δs = np.linspace(1,10,20)#np.sort(np.concatenate((np.linspace(1, 5, 20), np.linspace(5, 7, int(20/5 * 2))))) + Δs = np.linspace(1,10,10)#np.sort(np.concatenate((np.linspace(1, 5, 20), np.linspace(5, 7, int(20/5 * 2))))) Δ_models = [] strobe_ts = [] strobe_indices_s = [] - δs = [2, 1,.5, .1] + δs = np.linspace(.1, 2, 10) for Δ in Δs: for δ in δs: proto, strobe_t, strobe_indices = es.energy_shovel(Δ, periods=20, k_max=5, modulate_system=False) @@ -1208,78 +1738,98 @@ Whole different steady state -> phase transition? \(λ=0.1, Δ=1.0, T=5\) \(λ=0.1, Δ=1.0, T=5\) \(λ=0.1, Δ=1.0, T=5\) - \(λ=0.1, Δ=1.4736842105263157, T=5\) - \(λ=0.1, Δ=1.4736842105263157, T=5\) - \(λ=0.1, Δ=1.4736842105263157, T=5\) - \(λ=0.1, Δ=1.4736842105263157, T=5\) - \(λ=0.1, Δ=1.9473684210526314, T=5\) - \(λ=0.1, Δ=1.9473684210526314, T=5\) - \(λ=0.1, Δ=1.9473684210526314, T=5\) - \(λ=0.1, Δ=1.9473684210526314, T=5\) - \(λ=0.1, Δ=2.4210526315789473, T=5\) - \(λ=0.1, Δ=2.4210526315789473, T=5\) - \(λ=0.1, Δ=2.4210526315789473, T=5\) - \(λ=0.1, Δ=2.4210526315789473, T=5\) - \(λ=0.1, Δ=2.894736842105263, T=5\) - \(λ=0.1, Δ=2.894736842105263, T=5\) - \(λ=0.1, Δ=2.894736842105263, T=5\) - \(λ=0.1, Δ=2.894736842105263, T=5\) - \(λ=0.1, Δ=3.3684210526315788, T=5\) - \(λ=0.1, Δ=3.3684210526315788, T=5\) - \(λ=0.1, Δ=3.3684210526315788, T=5\) - \(λ=0.1, Δ=3.3684210526315788, T=5\) - \(λ=0.1, Δ=3.8421052631578947, T=5\) - \(λ=0.1, Δ=3.8421052631578947, T=5\) - \(λ=0.1, Δ=3.8421052631578947, T=5\) - \(λ=0.1, Δ=3.8421052631578947, T=5\) - \(λ=0.1, Δ=4.315789473684211, T=5\) - \(λ=0.1, Δ=4.315789473684211, T=5\) - \(λ=0.1, Δ=4.315789473684211, T=5\) - \(λ=0.1, Δ=4.315789473684211, T=5\) - \(λ=0.1, Δ=4.789473684210526, T=5\) - \(λ=0.1, Δ=4.789473684210526, T=5\) - \(λ=0.1, Δ=4.789473684210526, T=5\) - \(λ=0.1, Δ=4.789473684210526, T=5\) - \(λ=0.1, Δ=5.263157894736842, T=5\) - \(λ=0.1, Δ=5.263157894736842, T=5\) - \(λ=0.1, Δ=5.263157894736842, T=5\) - \(λ=0.1, Δ=5.263157894736842, T=5\) - \(λ=0.1, Δ=5.7368421052631575, T=5\) - \(λ=0.1, Δ=5.7368421052631575, T=5\) - \(λ=0.1, Δ=5.7368421052631575, T=5\) - \(λ=0.1, Δ=5.7368421052631575, T=5\) - \(λ=0.1, Δ=6.2105263157894735, T=5\) - \(λ=0.1, Δ=6.2105263157894735, T=5\) - \(λ=0.1, Δ=6.2105263157894735, T=5\) - \(λ=0.1, Δ=6.2105263157894735, T=5\) - \(λ=0.1, Δ=6.684210526315789, T=5\) - \(λ=0.1, Δ=6.684210526315789, T=5\) - \(λ=0.1, Δ=6.684210526315789, T=5\) - \(λ=0.1, Δ=6.684210526315789, T=5\) - \(λ=0.1, Δ=7.157894736842105, T=5\) - \(λ=0.1, Δ=7.157894736842105, T=5\) - \(λ=0.1, Δ=7.157894736842105, T=5\) - \(λ=0.1, Δ=7.157894736842105, T=5\) - \(λ=0.1, Δ=7.63157894736842, T=5\) - \(λ=0.1, Δ=7.63157894736842, T=5\) - \(λ=0.1, Δ=7.63157894736842, T=5\) - \(λ=0.1, Δ=7.63157894736842, T=5\) - \(λ=0.1, Δ=8.105263157894736, T=5\) - \(λ=0.1, Δ=8.105263157894736, T=5\) - \(λ=0.1, Δ=8.105263157894736, T=5\) - \(λ=0.1, Δ=8.105263157894736, T=5\) - \(λ=0.1, Δ=8.578947368421051, T=5\) - \(λ=0.1, Δ=8.578947368421051, T=5\) - \(λ=0.1, Δ=8.578947368421051, T=5\) - \(λ=0.1, Δ=8.578947368421051, T=5\) - \(λ=0.1, Δ=9.052631578947368, T=5\) - \(λ=0.1, Δ=9.052631578947368, T=5\) - \(λ=0.1, Δ=9.052631578947368, T=5\) - \(λ=0.1, Δ=9.052631578947368, T=5\) - \(λ=0.1, Δ=9.526315789473683, T=5\) - \(λ=0.1, Δ=9.526315789473683, T=5\) - \(λ=0.1, Δ=9.526315789473683, T=5\) - \(λ=0.1, Δ=9.526315789473683, T=5\) + \(λ=0.1, Δ=1.0, T=5\) + \(λ=0.1, Δ=1.0, T=5\) + \(λ=0.1, Δ=1.0, T=5\) + \(λ=0.1, Δ=1.0, T=5\) + \(λ=0.1, Δ=1.0, T=5\) + \(λ=0.1, Δ=1.0, T=5\) + \(λ=0.1, Δ=2.0, T=5\) + \(λ=0.1, Δ=2.0, T=5\) + \(λ=0.1, Δ=2.0, T=5\) + \(λ=0.1, Δ=2.0, T=5\) + \(λ=0.1, Δ=2.0, T=5\) + \(λ=0.1, Δ=2.0, T=5\) + \(λ=0.1, Δ=2.0, T=5\) + \(λ=0.1, Δ=2.0, T=5\) + \(λ=0.1, Δ=2.0, T=5\) + \(λ=0.1, Δ=2.0, T=5\) + \(λ=0.1, Δ=3.0, T=5\) + \(λ=0.1, Δ=3.0, T=5\) + \(λ=0.1, Δ=3.0, T=5\) + \(λ=0.1, Δ=3.0, T=5\) + \(λ=0.1, Δ=3.0, T=5\) + \(λ=0.1, Δ=3.0, T=5\) + \(λ=0.1, Δ=3.0, T=5\) + \(λ=0.1, Δ=3.0, T=5\) + \(λ=0.1, Δ=3.0, T=5\) + \(λ=0.1, Δ=3.0, T=5\) + \(λ=0.1, Δ=4.0, T=5\) + \(λ=0.1, Δ=4.0, T=5\) + \(λ=0.1, Δ=4.0, T=5\) + \(λ=0.1, Δ=4.0, T=5\) + \(λ=0.1, Δ=4.0, T=5\) + \(λ=0.1, Δ=4.0, T=5\) + \(λ=0.1, Δ=4.0, T=5\) + \(λ=0.1, Δ=4.0, T=5\) + \(λ=0.1, Δ=4.0, T=5\) + \(λ=0.1, Δ=4.0, T=5\) + \(λ=0.1, Δ=5.0, T=5\) + \(λ=0.1, Δ=5.0, T=5\) + \(λ=0.1, Δ=5.0, T=5\) + \(λ=0.1, Δ=5.0, T=5\) + \(λ=0.1, Δ=5.0, T=5\) + \(λ=0.1, Δ=5.0, T=5\) + \(λ=0.1, Δ=5.0, T=5\) + \(λ=0.1, Δ=5.0, T=5\) + \(λ=0.1, Δ=5.0, T=5\) + \(λ=0.1, Δ=5.0, T=5\) + \(λ=0.1, Δ=6.0, T=5\) + \(λ=0.1, Δ=6.0, T=5\) + \(λ=0.1, Δ=6.0, T=5\) + \(λ=0.1, Δ=6.0, T=5\) + \(λ=0.1, Δ=6.0, T=5\) + \(λ=0.1, Δ=6.0, T=5\) + \(λ=0.1, Δ=6.0, T=5\) + \(λ=0.1, Δ=6.0, T=5\) + \(λ=0.1, Δ=6.0, T=5\) + \(λ=0.1, Δ=6.0, T=5\) + \(λ=0.1, Δ=7.0, T=5\) + \(λ=0.1, Δ=7.0, T=5\) + \(λ=0.1, Δ=7.0, T=5\) + \(λ=0.1, Δ=7.0, T=5\) + \(λ=0.1, Δ=7.0, T=5\) + \(λ=0.1, Δ=7.0, T=5\) + \(λ=0.1, Δ=7.0, T=5\) + \(λ=0.1, Δ=7.0, T=5\) + \(λ=0.1, Δ=7.0, T=5\) + \(λ=0.1, Δ=7.0, T=5\) + \(λ=0.1, Δ=8.0, T=5\) + \(λ=0.1, Δ=8.0, T=5\) + \(λ=0.1, Δ=8.0, T=5\) + \(λ=0.1, Δ=8.0, T=5\) + \(λ=0.1, Δ=8.0, T=5\) + \(λ=0.1, Δ=8.0, T=5\) + \(λ=0.1, Δ=8.0, T=5\) + \(λ=0.1, Δ=8.0, T=5\) + \(λ=0.1, Δ=8.0, T=5\) + \(λ=0.1, Δ=8.0, T=5\) + \(λ=0.1, Δ=9.0, T=5\) + \(λ=0.1, Δ=9.0, T=5\) + \(λ=0.1, Δ=9.0, T=5\) + \(λ=0.1, Δ=9.0, T=5\) + \(λ=0.1, Δ=9.0, T=5\) + \(λ=0.1, Δ=9.0, T=5\) + \(λ=0.1, Δ=9.0, T=5\) + \(λ=0.1, Δ=9.0, T=5\) + \(λ=0.1, Δ=9.0, T=5\) + \(λ=0.1, Δ=9.0, T=5\) + \(λ=0.1, Δ=10.0, T=5\) + \(λ=0.1, Δ=10.0, T=5\) + \(λ=0.1, Δ=10.0, T=5\) + \(λ=0.1, Δ=10.0, T=5\) + \(λ=0.1, Δ=10.0, T=5\) + \(λ=0.1, Δ=10.0, T=5\) \(λ=0.1, Δ=10.0, T=5\) \(λ=0.1, Δ=10.0, T=5\) \(λ=0.1, Δ=10.0, T=5\) @@ -1363,13 +1913,13 @@ Whole different steady state -> phase transition? # ax.plot(model.t, model.L(model.t)[:,0,1]) # print(strobe_t[1]) - plt.ylabel(r"$\frac{\Delta E}{T}$") - ax.set_xlabel(r"$\Delta$") - ax.set_ylabel(r"$P_\mathrm{max}$") - for i, energy, σ in zip(itertools.count(0), final_e, final_e_error): - ax.errorbar(Δs, energy, σ, label=rf"$\alpha(0)$ = {δs[i]}") - ax.legend() - fs.export_fig("delta_dependence") + # plt.ylabel(r"$\frac{\Delta E}{T}$") + # ax.set_xlabel(r"$\Delta$") + # ax.set_ylabel(r"$P_\mathrm{max}$") + # for i, energy, σ in zip(itertools.count(0), final_e, final_e_error): + # ax.errorbar(Δs, energy, σ, label=rf"$\alpha(0)$ = {δs[i]}") + # ax.legend() + # fs.export_fig("delta_dependence") #+end_src #+RESULTS: