From 520adcb1a756858db3dc4621890cd7e4f704809d Mon Sep 17 00:00:00 2001 From: Valentin Boettcher Date: Fri, 26 Aug 2022 14:00:05 +0200 Subject: [PATCH] milked 08:mod_freq, have to do more detailed simulations now --- .../08_dynamic_one_bath/mod_freq.py | 89 +- .../shovel_experiments.org | 1189 +++++++++++++---- 2 files changed, 1025 insertions(+), 253 deletions(-) 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 080c480..1b49307 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 @@ -9,7 +9,7 @@ import plot_utils as pu import ray ray.shutdown() -ray.init(address="auto") +ray.init() from hops.util.logging_setup import logging_setup import logging @@ -17,39 +17,57 @@ logging_setup(logging.INFO, show_stocproc=False) from hops.util.dynamic_matrix import SmoothStep, Periodic, Harmonic, ConstantMatrix - -Δs = np.linspace(1,10,10)#np.sort(np.concatenate((np.linspace(1, 5, 20), np.linspace(5, 7, int(20/5 * 2))))) +Δs = [ + 1, + 2, + 3, + 4, + 5, + 6, + 7, + 8, + 9, + 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 = np.linspace(.1, 2, 10) +δs = np.linspace(0.1, 1.5, 10) for Δ in Δs: for δ in δs: - proto, strobe_t, strobe_indices = es.energy_shovel(Δ, periods=20, k_max=5, modulate_system=False) + proto, strobe_t, strobe_indices = es.energy_shovel( + Δ, + periods=int(Δ / (2 * np.pi) * 20), + k_max=3, + bcf_terms=4, + modulate_system=False, + bcf_norm_method="unit_therm", + ) proto.δ = δ strobe_ts.append(strobe_t) strobe_indices_s.append(strobe_indices) Δ_models.append(proto) -aux.integrate_multi(Δ_models, 10_000) +Δ_models[1].bcf(0) * Δ_models[-1].bcf_scale -final_e = [[], []] -final_e_error = [[], []] +aux.integrate_multi(Δ_models, 2000) + +dim_models = len(δs) +power = np.zeros((len(Δs), dim_models)) + +# final_e_error = [[], []] ensemble_arg = dict(overwrite_cache=False) fig, ax = plt.subplots() -for (model, data), Δ, strobe_t, strobe_indices in zip( +for (model, data), Δ, strobe_t, strobe_indices, idx in zip( aux.model_data_iterator(Δ_models), np.array([[Δ]*len(δs) for Δ in Δs]).flatten(), strobe_ts, strobe_indices_s, + range(len(Δ_models)) ): energies = model.total_energy_from_power(data, **ensemble_arg) - idx = energies.value[strobe_indices].argmin() - energies = energies * (1 / strobe_t[idx]) - - energy_idx = δs.index(model.δ) - final_e[energy_idx].append(energies.value[strobe_indices[idx]]) - final_e_error[energy_idx].append(energies.σ[strobe_indices[idx]]) + powers = energies.value[strobe_indices[1:]] * (1/strobe_t[1:]) + power[int(idx/dim_models), idx % dim_models] = abs(powers[powers <= 0]).max() # fs.plot_energy_overview(model, ensemble_args=ensemble_arg) # fig, ax = plt.subplots() @@ -69,3 +87,44 @@ for (model, data), Δ, strobe_t, strobe_indices in zip( # ax.errorbar(Δs, energy, σ, label=rf"$\alpha(0)$ = {δs[i]}") # ax.legend() # fs.export_fig("delta_dependence") + +normed_power = power.copy() +for i in range(dim_models): + normed_power[:, i] /= power[:, i].max() + +for i in range(dim_models): + plt.plot(Δs, power[:, i]) + +f, (a, a2) = plt.subplots(ncols=2, sharey=True) +im = a.imshow( + power, + interpolation="gaussian", + interpolation_stage="data", + aspect="auto", + origin="lower", + cmap="plasma", + extent=[δs.min(), δs.max(), min(Δs), max(Δs)], +) +f.colorbar(im, ax=a) + +im2 = a2.imshow( + normed_power, + interpolation="gaussian", + interpolation_stage="data", + aspect="auto", + origin="lower", + cmap="plasma", + extent=[δs.min(), δs.max(), min(Δs), max(Δs)], +) +f.colorbar(im2, ax=a2) +a.set_ylabel(r"$\Delta$") +for ax in (a, a2): + ax.set_xlabel(r"$\alpha_\beta(0)$") + +a.set_title(r"$\bar{P}$") +a2.set_title(r"$\bar{P}/\bar{P}_{\mathrm{\max},\alpha}$") +fs.export_fig("power_heatmap", f, tikz=False) + +for model, data in aux.model_data_iterator(Δ_models): + inter = model.interaction_energy(data) + print(abs(inter.value).max()) 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 09472f2..31d285f 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 @@ -13,7 +13,7 @@ import ray ray.shutdown() - ray.init(address="auto") + ray.init() from hops.util.logging_setup import logging_setup import logging @@ -303,7 +303,7 @@ :END: * Only Coupling Modulation -** Dependence on Cutoff +** TODO Dependence on Cutoff :PROPERTIES: :header-args: :tangle cutoff_dep.py :session 08_shovel_cutoff_dep :async yes :noweb yes :END: @@ -317,7 +317,7 @@ - found that with low temperature only a small cutoff produces a net energy gain -#+begin_src jupyter-python :tangle night.py +#+begin_src jupyter-python import scipy import pathlib @@ -376,44 +376,144 @@ \(λ=0.1, Δ=5, T=5\) \(λ=0.1, Δ=5, T=5\) \(λ=0.1, Δ=5, T=5\) - [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.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] + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) #+end_example -#+begin_src jupyter-python :tangle night.py +#+begin_src jupyter-python + ω_models[5].δ +#+end_src + +#+RESULTS: +: array(1.59883336) + +#+begin_src jupyter-python aux.integrate_multi(ω_models, 10_000) #+end_src #+RESULTS: +:RESULTS: +#+begin_example + [INFO hops.core.integration 254439] Choosing the nonlinear integrator. + [INFO hops.core.integration 254439] Using 21 integrators. + [INFO hops.core.integration 254439] Some 9966 trajectories have to be integrated. + [INFO hops.core.integration 254439] Using 792 hierarchy states. + 0% 0/9966 [00:00() + ----> 1 aux.integrate_multi(ω_models, 10_000) + + File ~/src/two_qubit_model/hiro_models/model_auxiliary.py:84, in integrate_multi(models, *args, **kwargs) +  77 """Integrate the hops equations for the ``models``. +  78 Like :any:`integrate` just for many models. +  79 +  80 A call to :any:`ray.init` may be required. +  81 """ +  83 for model in models: + ---> 84 integrate(model, *args, **kwargs) + + File ~/src/two_qubit_model/hiro_models/model_auxiliary.py:108, in integrate(model, n, data_path, clear_pd) +  98 # with model_db(data_path) as db: +  99 # if hash in db and "data" db[hash] +  101 supervisor = HOPSSupervisor( +  102 model.hops_config, +  103 n, +  104 data_path=data_path, +  105 data_name=hash, +  106 ) + --> 108 supervisor.integrate(clear_pd) +  110 with supervisor.get_data(True) as data: +  111 with model_db(data_path) as db: + + File ~/src/hops/hops/core/integration.py:1288, in HOPSSupervisor.integrate(self, clear_pd) +  1285 break +  1287 integration.update() + -> 1288 data.new_samples( +  1289 idx=index, +  1290 incomplete=incomplete, +  1291 psi0=psi0, +  1292 aux_states=aux_states, +  1293 stoc_proc=stoc_proc, +  1294 result_type=self.params.HiP.result_type, +  1295 normed=self._normed_average, +  1296 rng_seed=seed, +  1297 ) + + File ~/src/hops/hops/core/signal_delay.py:87, in sig_delay.__exit__(self, exc_type, exc_val, exc_tb) +  84 if len(self.sigh.sigs_caught) > 0 and self.handler is not None: +  85 self.handler(self.sigh.sigs_caught) + ---> 87 self._restore() + + File ~/src/hops/hops/core/signal_delay.py:68, in sig_delay._restore(self) +  66 for i, s in enumerate(self.sigs): +  67 signal.signal(s, self.old_handlers[i]) + ---> 68 self.sigh.emit() + + File ~/src/hops/hops/core/signal_delay.py:42, in SigHandler.emit(self) +  40 for s in self.sigs_caught: +  41 log.info("emit signal '{}'".format(SIG_MAP[s])) + ---> 42 os.kill(os.getpid(), s) + + KeyboardInterrupt: +#+end_example +:END: #+begin_src jupyter-python @@ -423,11 +523,6 @@ Loading: 100% 21/21 [00:02<00:00, 7.49it/s] #+RESULTS: : 712fb5f5-49f4-4383-9e5c-13eca64a73ca -#+begin_src jupyter-python :tangle night.py - from pprint import pprint - print_diff_tbl(ω_models) -#+end_src - *** BCFs #+begin_src jupyter-python @@ -512,7 +607,48 @@ The time scale of the bath and the modulation is also interesting. #+end_src #+RESULTS: -: 3ca75dfb-222a-4c19-8741-0a0200bc2f32 +:RESULTS: +# [goto error] +#+begin_example + --------------------------------------------------------------------------- + RuntimeError Traceback (most recent call last) + Input In [20], in () +  1 fig, ax = plt.subplots() +  2 ax.axhline(0, color="gray") + ----> 3 for model, data in aux.model_data_iterator(ω_models): +  4 _, _, bar = pu.plot_with_σ( +  5 model.t, +  6 model.interaction_energy(data).for_bath(0), +  7 ax=ax, +  8 label=fr"$\omega_c={model.ω_c:.2f}$", +  9 ) +  10 ax.set_xlim(0, model.t[strobe_indices][3]) + + 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 'fa41de4dce918a19bae348fe2636d0e3c092a90bf7a4c5f11e378c4b1f38d20f'. +#+end_example +[[file:./.ob-jupyter/6bb1ad6b59df8fa850749bc64d469393c08158bf.svg]] +:END: Total Energy Extraction @@ -588,7 +724,48 @@ Total Energy Extraction #+end_src #+RESULTS: -: 914dca52-374e-4b3d-ae08-e6f330d76c5e +:RESULTS: +# [goto error] +#+begin_example + --------------------------------------------------------------------------- + RuntimeError Traceback (most recent call last) + Input In [21], in () +  1 fig, ax = plt.subplots() +  2 import scipy + ----> 4 for model, data in aux.model_data_iterator(ω_models): +  5 energy = model.total_energy_from_power(data) * (1 / es.ergo(model.T)) +  6 _, _, bar = pu.plot_with_σ( +  7 model.t, +  8 energy, +  (...) +  12 markersize=3, +  13 ) + + 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 'fa41de4dce918a19bae348fe2636d0e3c092a90bf7a4c5f11e378c4b1f38d20f'. +#+end_example +[[file:./.ob-jupyter/e50959fe53d7ef60c6d7441a095c12433620df7a.svg]] +:END: Now we plot max power vs memory time. #+begin_src jupyter-python @@ -649,7 +826,42 @@ Now we plot max power vs memory time. #+end_src #+RESULTS: -: 839b9591-0955-409b-9e52-ab3f0656fed3 +:RESULTS: +# [goto error] +#+begin_example + --------------------------------------------------------------------------- + RuntimeError Traceback (most recent call last) + Input In [22], in () +  4 opt_energies = [] +  5 memories = [] + ----> 6 for model, data in aux.model_data_iterator(ω_models): +  7 energy = model.total_energy_from_power(data) +  8 ts = model.t[strobe_indices] + + 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 'fa41de4dce918a19bae348fe2636d0e3c092a90bf7a4c5f11e378c4b1f38d20f'. +#+end_example +:END: @@ -966,7 +1178,7 @@ Whole different steady state -> phase transition? - found that with low temperature only a small cutoff produces a net energy gain -#+begin_src jupyter-python :tangle night.py +#+begin_src jupyter-python import scipy import pathlib @@ -1717,15 +1929,32 @@ Loading: 100% 7/7 [00:00<00:00, 9.35it/s] #+begin_src jupyter-python from hops.util.dynamic_matrix import SmoothStep, Periodic, Harmonic, ConstantMatrix - - Δs = np.linspace(1,10,10)#np.sort(np.concatenate((np.linspace(1, 5, 20), np.linspace(5, 7, int(20/5 * 2))))) + Δs = [ + 1, + 2, + 3, + 4, + 5, + 6, + 7, + 8, + 9, + 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 = np.linspace(.1, 2, 10) + δs = np.linspace(0.1, 1.5, 10) for Δ in Δs: for δ in δs: - proto, strobe_t, strobe_indices = es.energy_shovel(Δ, periods=20, k_max=5, modulate_system=False) + proto, strobe_t, strobe_indices = es.energy_shovel( + Δ, + periods=int(Δ / (2 * np.pi) * 20), + k_max=3, + bcf_terms=4, + modulate_system=False, + bcf_norm_method="unit_therm", + ) proto.δ = δ strobe_ts.append(strobe_t) strobe_indices_s.append(strobe_indices) @@ -1734,173 +1963,640 @@ Loading: 100% 7/7 [00:00<00:00, 9.35it/s] #+RESULTS: #+begin_example - \(λ=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, Δ=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\) - \(λ=0.1, Δ=10.0, T=5\) + \(λ=0.1, Δ=1, T=5\) + \(λ=0.1, Δ=1, T=5\) + \(λ=0.1, Δ=1, T=5\) + \(λ=0.1, Δ=1, T=5\) + \(λ=0.1, Δ=1, T=5\) + \(λ=0.1, Δ=1, T=5\) + \(λ=0.1, Δ=1, T=5\) + \(λ=0.1, Δ=1, T=5\) + \(λ=0.1, Δ=1, T=5\) + \(λ=0.1, Δ=1, T=5\) + \(λ=0.1, Δ=2, T=5\) + \(λ=0.1, Δ=2, T=5\) + \(λ=0.1, Δ=2, T=5\) + \(λ=0.1, Δ=2, T=5\) + \(λ=0.1, Δ=2, T=5\) + \(λ=0.1, Δ=2, T=5\) + \(λ=0.1, Δ=2, T=5\) + \(λ=0.1, Δ=2, T=5\) + \(λ=0.1, Δ=2, T=5\) + \(λ=0.1, Δ=2, T=5\) + \(λ=0.1, Δ=3, T=5\) + \(λ=0.1, Δ=3, T=5\) + \(λ=0.1, Δ=3, T=5\) + \(λ=0.1, Δ=3, T=5\) + \(λ=0.1, Δ=3, T=5\) + \(λ=0.1, Δ=3, T=5\) + \(λ=0.1, Δ=3, T=5\) + \(λ=0.1, Δ=3, T=5\) + \(λ=0.1, Δ=3, T=5\) + \(λ=0.1, Δ=3, T=5\) + \(λ=0.1, Δ=4, T=5\) + \(λ=0.1, Δ=4, T=5\) + \(λ=0.1, Δ=4, T=5\) + \(λ=0.1, Δ=4, T=5\) + \(λ=0.1, Δ=4, T=5\) + \(λ=0.1, Δ=4, T=5\) + \(λ=0.1, Δ=4, T=5\) + \(λ=0.1, Δ=4, T=5\) + \(λ=0.1, Δ=4, T=5\) + \(λ=0.1, Δ=4, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=5, T=5\) + \(λ=0.1, Δ=6, T=5\) + \(λ=0.1, Δ=6, T=5\) + \(λ=0.1, Δ=6, T=5\) + \(λ=0.1, Δ=6, T=5\) + \(λ=0.1, Δ=6, T=5\) + \(λ=0.1, Δ=6, T=5\) + \(λ=0.1, Δ=6, T=5\) + \(λ=0.1, Δ=6, T=5\) + \(λ=0.1, Δ=6, T=5\) + \(λ=0.1, Δ=6, T=5\) + \(λ=0.1, Δ=7, T=5\) + \(λ=0.1, Δ=7, T=5\) + \(λ=0.1, Δ=7, T=5\) + \(λ=0.1, Δ=7, T=5\) + \(λ=0.1, Δ=7, T=5\) + \(λ=0.1, Δ=7, T=5\) + \(λ=0.1, Δ=7, T=5\) + \(λ=0.1, Δ=7, T=5\) + \(λ=0.1, Δ=7, T=5\) + \(λ=0.1, Δ=7, T=5\) + \(λ=0.1, Δ=8, T=5\) + \(λ=0.1, Δ=8, T=5\) + \(λ=0.1, Δ=8, T=5\) + \(λ=0.1, Δ=8, T=5\) + \(λ=0.1, Δ=8, T=5\) + \(λ=0.1, Δ=8, T=5\) + \(λ=0.1, Δ=8, T=5\) + \(λ=0.1, Δ=8, T=5\) + \(λ=0.1, Δ=8, T=5\) + \(λ=0.1, Δ=8, T=5\) + \(λ=0.1, Δ=9, T=5\) + \(λ=0.1, Δ=9, T=5\) + \(λ=0.1, Δ=9, T=5\) + \(λ=0.1, Δ=9, T=5\) + \(λ=0.1, Δ=9, T=5\) + \(λ=0.1, Δ=9, T=5\) + \(λ=0.1, Δ=9, T=5\) + \(λ=0.1, Δ=9, T=5\) + \(λ=0.1, Δ=9, T=5\) + \(λ=0.1, Δ=9, T=5\) + \(λ=0.1, Δ=10, T=5\) + \(λ=0.1, Δ=10, T=5\) + \(λ=0.1, Δ=10, T=5\) + \(λ=0.1, Δ=10, T=5\) + \(λ=0.1, Δ=10, T=5\) + \(λ=0.1, Δ=10, T=5\) + \(λ=0.1, Δ=10, T=5\) + \(λ=0.1, Δ=10, T=5\) + \(λ=0.1, Δ=10, T=5\) + \(λ=0.1, Δ=10, T=5\) #+end_example + #+begin_src jupyter-python - aux.integrate_multi(Δ_models, 10_000) + Δ_models[1].bcf(0) * Δ_models[-1].bcf_scale #+end_src #+RESULTS: -:RESULTS: -: [INFO hops.core.integration 2187665] Choosing the nonlinear integrator. -# [goto error] -#+begin_example - --------------------------------------------------------------------------- - RuntimeError Traceback (most recent call last) - Input In [6], in () - ----> 1 aux.integrate_multi(Δ_models, 10_000) - - File ~/src/two_qubit_model/hiro_models/model_auxiliary.py:84, in integrate_multi(models, *args, **kwargs) -  77 """Integrate the hops equations for the ``models``. -  78 Like :any:`integrate` just for many models. -  79 -  80 A call to :any:`ray.init` may be required. -  81 """ -  83 for model in models: - ---> 84 integrate(model, *args, **kwargs) - - File ~/src/two_qubit_model/hiro_models/model_auxiliary.py:108, in integrate(model, n, data_path, clear_pd) -  98 # with model_db(data_path) as db: -  99 # if hash in db and "data" db[hash] -  101 supervisor = HOPSSupervisor( -  102 model.hops_config, -  103 n, -  104 data_path=data_path, -  105 data_name=hash, -  106 ) - --> 108 supervisor.integrate(clear_pd) -  110 with supervisor.get_data(True) as data: -  111 with model_db(data_path) as db: - - File ~/src/hops/hops/core/integration.py:1241, in HOPSSupervisor.integrate(self, clear_pd) -  1238 num_integrators = int(ray.available_resources().get("CPU", 0)) -  1240 if num_integrators == 0: - -> 1241 raise RuntimeError("No cpu available for integration!") -  1243 log.info(f"Using {num_integrators} integrators.") -  1245 indices = self.get_job_args(data) - - RuntimeError: No cpu available for integration! -#+end_example -:END: +| 1.240875157568988+0j | #+begin_src jupyter-python - final_e = [[], []] - final_e_error = [[], []] + aux.integrate_multi(Δ_models, 2000) +#+end_src + +#+RESULTS: +#+begin_example + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. + 0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] + [INFO hops.core.integration 453705] Choosing the nonlinear integrator. + [INFO hops.core.integration 453705] Using 4 integrators. + [INFO hops.core.integration 453705] Some 0 trajectories have to be integrated. + [INFO hops.core.integration 453705] Using 35 hierarchy states. +0it [00:00, ?it/s] +#+end_example + +#+begin_src jupyter-python + dim_models = len(δs) + power = np.zeros((len(Δs), dim_models)) + + # final_e_error = [[], []] ensemble_arg = dict(overwrite_cache=False) fig, ax = plt.subplots() - for (model, data), Δ, strobe_t, strobe_indices in zip( + for (model, data), Δ, strobe_t, strobe_indices, idx in zip( aux.model_data_iterator(Δ_models), np.array([[Δ]*len(δs) for Δ in Δs]).flatten(), strobe_ts, strobe_indices_s, + range(len(Δ_models)) ): energies = model.total_energy_from_power(data, **ensemble_arg) - idx = energies.value[strobe_indices].argmin() - energies = energies * (1 / strobe_t[idx]) - - energy_idx = δs.index(model.δ) - final_e[energy_idx].append(energies.value[strobe_indices[idx]]) - final_e_error[energy_idx].append(energies.σ[strobe_indices[idx]]) + powers = energies.value[strobe_indices[1:]] * (1/strobe_t[1:]) + power[int(idx/dim_models), idx % dim_models] = abs(powers[powers <= 0]).max() # fs.plot_energy_overview(model, ensemble_args=ensemble_arg) # fig, ax = plt.subplots() @@ -1923,48 +2619,65 @@ Loading: 100% 7/7 [00:00<00:00, 9.35it/s] #+end_src #+RESULTS: -:RESULTS: -# [goto error] -#+begin_example - --------------------------------------------------------------------------- - RuntimeError Traceback (most recent call last) - Input In [7], in () -  3 ensemble_arg = dict(overwrite_cache=False) -  4 fig, ax = plt.subplots() - ----> 5 for (model, data), Δ, strobe_t, strobe_indices in zip( -  6 aux.model_data_iterator(Δ_models), -  7 np.array([[Δ]*len(δs) for Δ in Δs]).flatten(), -  8 strobe_ts, -  9 strobe_indices_s, -  10 ): -  11 energies = model.total_energy_from_power(data, **ensemble_arg) -  12 idx = energies.value[strobe_indices].argmin() +[[file:./.ob-jupyter/3678e7450bf95cb9f19cad5c334a5b8cb0fedf07.svg]] - 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 +#+begin_src jupyter-python :results none + normed_power = power.copy() + for i in range(dim_models): + normed_power[:, i] /= power[:, i].max() +#+end_src + + +#+begin_src jupyter-python + for i in range(dim_models): + plt.plot(Δs, power[:, i]) +#+end_src + +#+RESULTS: +[[file:./.ob-jupyter/314ad7fe53881fb46e1526f595a91fd3d82e8e90.svg]] + +#+begin_src jupyter-python + f, (a, a2) = plt.subplots(ncols=2, sharey=True) + im = a.imshow( + power, + interpolation="gaussian", + interpolation_stage="data", + aspect="auto", + origin="lower", + cmap="plasma", + extent=[δs.min(), δs.max(), min(Δs), max(Δs)], + ) + f.colorbar(im, ax=a) + + im2 = a2.imshow( + normed_power, + interpolation="gaussian", + interpolation_stage="data", + aspect="auto", + origin="lower", + cmap="plasma", + extent=[δs.min(), δs.max(), min(Δs), max(Δs)], + ) + f.colorbar(im2, ax=a2) + a.set_ylabel(r"$\Delta$") + for ax in (a, a2): + ax.set_xlabel(r"$\alpha_\beta(0)$") + + a.set_title(r"$\bar{P}$") + a2.set_title(r"$\bar{P}/\bar{P}_{\mathrm{\max},\alpha}$") + fs.export_fig("power_heatmap", f, tikz=False) +#+end_src + +#+RESULTS: +[[file:./.ob-jupyter/0e21ce0ba7d686744421826a5fe849fce1f4673d.svg]] + +#+begin_src jupyter-python + for model, data in aux.model_data_iterator(Δ_models): + inter = model.interaction_energy(data) + print(abs(inter.value).max()) +#+end_src - 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 'aa73b52f3e07edf0fa67327b5bca8bcddc7bcde8ff1727adec1303fb9a22eb78'. -#+end_example -[[file:./.ob-jupyter/46246e7d08d0a701ac2a68c2e4bf9f13274dbf05.svg]] -:END: Looks very much like a speed limit. Likely due tho lower temperature suppression.