mirror of
https://github.com/vale981/master-thesis
synced 2025-03-05 10:01:43 -05:00
try with 9 HO levels
This commit is contained in:
parent
0c165279f3
commit
710326b8c4
2 changed files with 176 additions and 15 deletions
|
|
@ -114,18 +114,25 @@
|
|||
return params
|
||||
|
||||
|
||||
model_keys = [dict(bcf_scale=0.2, wc=ω_c, γ=.5, T=.6, k_max=4, bcf_terms=5, max_HO_level=5) for ω_c in [2]]
|
||||
model_keys = [dict(bcf_scale=0.2, wc=ω_c, γ=.5, T=.6, k_max=3, bcf_terms=5, max_HO_level=9) for ω_c in [2]]
|
||||
multi_params = [ho_duo(**keys) for keys in model_keys]
|
||||
multi_params = multi_params
|
||||
#+end_src
|
||||
|
||||
#+RESULTS:
|
||||
: stocproc.stocproc - INFO - Loaded instance from cache.
|
||||
: INFO:stocproc.stocproc:Loaded instance from cache.
|
||||
: stocproc.stocproc - INFO - Loaded instance from cache.
|
||||
: INFO:stocproc.stocproc:Loaded instance from cache.
|
||||
: stocproc.stocproc - INFO - Loaded instance from cache.
|
||||
: INFO:stocproc.stocproc:Loaded instance from cache.
|
||||
#+begin_example
|
||||
INFO:stocproc.stocproc:Loaded instance from cache.
|
||||
stocproc.stocproc - INFO - Loaded instance from cache.
|
||||
INFO:stocproc.stocproc:Loaded instance from cache.
|
||||
stocproc.stocproc - INFO - Loaded instance from cache.
|
||||
INFO:stocproc.stocproc:Loaded instance from cache.
|
||||
stocproc.stocproc - INFO - Loaded instance from cache.
|
||||
INFO:stocproc.stocproc:Loaded instance from cache.
|
||||
stocproc.stocproc - INFO - Loaded instance from cache.
|
||||
INFO:stocproc.stocproc:Loaded instance from cache.
|
||||
stocproc.stocproc - INFO - Loaded instance from cache.
|
||||
INFO:stocproc.stocproc:Loaded instance from cache.
|
||||
#+end_example
|
||||
|
||||
|
||||
#+begin_src jupyter-python
|
||||
|
|
@ -135,7 +142,10 @@
|
|||
#+end_src
|
||||
|
||||
#+RESULTS:
|
||||
:RESULTS:
|
||||
: 2022-08-04 09:20:43,308 INFO worker.py:956 -- Connecting to existing Ray cluster at address: 141.30.17.225:48665
|
||||
: RayContext(dashboard_url='', python_version='3.9.13', ray_version='1.13.0', ray_commit='e4ce38d001dbbe09cd21c497fedd03d692b2be3e', address_info={'node_ip_address': '141.30.17.225', 'raylet_ip_address': '141.30.17.225', 'redis_address': None, 'object_store_address': '/tmp/ray/session_2022-08-03_16-51-50_330688_106699/sockets/plasma_store', 'raylet_socket_name': '/tmp/ray/session_2022-08-03_16-51-50_330688_106699/sockets/raylet', 'webui_url': '', 'session_dir': '/tmp/ray/session_2022-08-03_16-51-50_330688_106699', 'metrics_export_port': 40539, 'gcs_address': '141.30.17.225:48665', 'address': '141.30.17.225:48665', 'node_id': '2b485e9daf283f489d86e8256bde8e8a015de5563d8e1923adda5db4'})
|
||||
:END:
|
||||
|
||||
|
||||
#+begin_src jupyter-python
|
||||
|
|
@ -143,7 +153,7 @@
|
|||
for params in multi_params:
|
||||
supervisor = HOPSSupervisor(
|
||||
params,
|
||||
10_000,
|
||||
2_000,
|
||||
data_path="ho_data",
|
||||
data_name="two_baths",
|
||||
)
|
||||
|
|
@ -154,8 +164,111 @@
|
|||
#+end_src
|
||||
|
||||
#+RESULTS:
|
||||
: ho_data/two_baths/_f/two_baths_fbeb4100da4ff052f0c2153129535231_1.h5
|
||||
: 0it [00:00, ?it/s]
|
||||
:RESULTS:
|
||||
: 0% 1/10000 [02:21<394:12:56, 141.93s/it]
|
||||
: 2022-08-04 09:48:05,978 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information.
|
||||
: 2022-08-04 09:48:05,981 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information.
|
||||
: 2022-08-04 09:48:05,984 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information.
|
||||
: 2022-08-04 09:48:05,986 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information.
|
||||
# [goto error]
|
||||
#+begin_example
|
||||
[0;31m---------------------------------------------------------------------------[0m
|
||||
[0;31mKeyboardInterrupt[0m Traceback (most recent call last)
|
||||
Input [0;32mIn [47][0m, in [0;36m<cell line: 2>[0;34m()[0m
|
||||
[1;32m 3[0m supervisor [38;5;241m=[39m HOPSSupervisor(
|
||||
[1;32m 4[0m params,
|
||||
[1;32m 5[0m [38;5;241m10_000[39m,
|
||||
[1;32m 6[0m data_path[38;5;241m=[39m[38;5;124m"[39m[38;5;124mho_data[39m[38;5;124m"[39m,
|
||||
[1;32m 7[0m data_name[38;5;241m=[39m[38;5;124m"[39m[38;5;124mtwo_baths[39m[38;5;124m"[39m,
|
||||
[1;32m 8[0m )
|
||||
[1;32m 9[0m [38;5;28mprint[39m(supervisor[38;5;241m.[39mget_data([38;5;28;01mTrue[39;00m)[38;5;241m.[39mhdf5_name)
|
||||
[0;32m---> 10[0m [43msupervisor[49m[38;5;241;43m.[39;49m[43mintegrate[49m[43m([49m[43m)[49m
|
||||
[1;32m 11[0m supervisors[38;5;241m.[39mappend(supervisor)
|
||||
|
||||
File [0;32m~/src/hops/hops/core/integration.py:1288[0m, in [0;36mHOPSSupervisor.integrate[0;34m(self, clear_pd)[0m
|
||||
[1;32m 1285[0m [38;5;28;01mbreak[39;00m
|
||||
[1;32m 1287[0m integration[38;5;241m.[39mupdate()
|
||||
[0;32m-> 1288[0m data[38;5;241m.[39mnew_samples(
|
||||
[1;32m 1289[0m idx[38;5;241m=[39mindex,
|
||||
[1;32m 1290[0m incomplete[38;5;241m=[39mincomplete,
|
||||
[1;32m 1291[0m psi0[38;5;241m=[39mpsi0,
|
||||
[1;32m 1292[0m aux_states[38;5;241m=[39maux_states,
|
||||
[1;32m 1293[0m stoc_proc[38;5;241m=[39mstoc_proc,
|
||||
[1;32m 1294[0m result_type[38;5;241m=[39m[38;5;28mself[39m[38;5;241m.[39mparams[38;5;241m.[39mHiP[38;5;241m.[39mresult_type,
|
||||
[1;32m 1295[0m normed[38;5;241m=[39m[38;5;28mself[39m[38;5;241m.[39m_normed_average,
|
||||
[1;32m 1296[0m rng_seed[38;5;241m=[39mseed,
|
||||
[1;32m 1297[0m )
|
||||
|
||||
File [0;32m~/src/hops/hops/core/signal_delay.py:87[0m, in [0;36msig_delay.__exit__[0;34m(self, exc_type, exc_val, exc_tb)[0m
|
||||
[1;32m 84[0m [38;5;28;01mif[39;00m [38;5;28mlen[39m([38;5;28mself[39m[38;5;241m.[39msigh[38;5;241m.[39msigs_caught) [38;5;241m>[39m [38;5;241m0[39m [38;5;129;01mand[39;00m [38;5;28mself[39m[38;5;241m.[39mhandler [38;5;129;01mis[39;00m [38;5;129;01mnot[39;00m [38;5;28;01mNone[39;00m:
|
||||
[1;32m 85[0m [38;5;28mself[39m[38;5;241m.[39mhandler([38;5;28mself[39m[38;5;241m.[39msigh[38;5;241m.[39msigs_caught)
|
||||
[0;32m---> 87[0m [38;5;28;43mself[39;49m[38;5;241;43m.[39;49m[43m_restore[49m[43m([49m[43m)[49m
|
||||
|
||||
File [0;32m~/src/hops/hops/core/signal_delay.py:68[0m, in [0;36msig_delay._restore[0;34m(self)[0m
|
||||
[1;32m 66[0m [38;5;28;01mfor[39;00m i, s [38;5;129;01min[39;00m [38;5;28menumerate[39m([38;5;28mself[39m[38;5;241m.[39msigs):
|
||||
[1;32m 67[0m signal[38;5;241m.[39msignal(s, [38;5;28mself[39m[38;5;241m.[39mold_handlers[i])
|
||||
[0;32m---> 68[0m [38;5;28;43mself[39;49m[38;5;241;43m.[39;49m[43msigh[49m[38;5;241;43m.[39;49m[43memit[49m[43m([49m[43m)[49m
|
||||
|
||||
File [0;32m~/src/hops/hops/core/signal_delay.py:42[0m, in [0;36mSigHandler.emit[0;34m(self)[0m
|
||||
[1;32m 40[0m [38;5;28;01mfor[39;00m s [38;5;129;01min[39;00m [38;5;28mself[39m[38;5;241m.[39msigs_caught:
|
||||
[1;32m 41[0m log[38;5;241m.[39minfo([38;5;124m"[39m[38;5;124memit signal [39m[38;5;124m'[39m[38;5;132;01m{}[39;00m[38;5;124m'[39m[38;5;124m"[39m[38;5;241m.[39mformat(SIG_MAP[s]))
|
||||
[0;32m---> 42[0m [43mos[49m[38;5;241;43m.[39;49m[43mkill[49m[43m([49m[43mos[49m[38;5;241;43m.[39;49m[43mgetpid[49m[43m([49m[43m)[49m[43m,[49m[43m [49m[43ms[49m[43m)[49m
|
||||
|
||||
[0;31mKeyboardInterrupt[0m:
|
||||
ho_data/two_baths/_4/two_baths_449cdb9aeeb1a354762f4a6b2f69bcee_1.h5
|
||||
2% 44/2000 [09:32<7:04:09, 13.01s/it]
|
||||
2022-08-04 09:57:41,769 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information.
|
||||
2022-08-04 09:57:41,771 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information.
|
||||
2022-08-04 09:57:41,780 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information.
|
||||
2022-08-04 09:57:41,782 ERROR worker.py:94 -- Unhandled error (suppress with 'RAY_IGNORE_UNHANDLED_ERRORS=1'): The worker died unexpectedly while executing this task. Check python-core-worker-*.log files for more information.
|
||||
#+end_example
|
||||
# [goto error]
|
||||
#+begin_example
|
||||
[0;31m---------------------------------------------------------------------------[0m
|
||||
[0;31mKeyboardInterrupt[0m Traceback (most recent call last)
|
||||
Input [0;32mIn [50][0m, in [0;36m<cell line: 2>[0;34m()[0m
|
||||
[1;32m 3[0m supervisor [38;5;241m=[39m HOPSSupervisor(
|
||||
[1;32m 4[0m params,
|
||||
[1;32m 5[0m [38;5;241m2_000[39m,
|
||||
[1;32m 6[0m data_path[38;5;241m=[39m[38;5;124m"[39m[38;5;124mho_data[39m[38;5;124m"[39m,
|
||||
[1;32m 7[0m data_name[38;5;241m=[39m[38;5;124m"[39m[38;5;124mtwo_baths[39m[38;5;124m"[39m,
|
||||
[1;32m 8[0m )
|
||||
[1;32m 9[0m [38;5;28mprint[39m(supervisor[38;5;241m.[39mget_data([38;5;28;01mTrue[39;00m)[38;5;241m.[39mhdf5_name)
|
||||
[0;32m---> 10[0m [43msupervisor[49m[38;5;241;43m.[39;49m[43mintegrate[49m[43m([49m[43m)[49m
|
||||
[1;32m 11[0m supervisors[38;5;241m.[39mappend(supervisor)
|
||||
|
||||
File [0;32m~/src/hops/hops/core/integration.py:1288[0m, in [0;36mHOPSSupervisor.integrate[0;34m(self, clear_pd)[0m
|
||||
[1;32m 1285[0m [38;5;28;01mbreak[39;00m
|
||||
[1;32m 1287[0m integration[38;5;241m.[39mupdate()
|
||||
[0;32m-> 1288[0m data[38;5;241m.[39mnew_samples(
|
||||
[1;32m 1289[0m idx[38;5;241m=[39mindex,
|
||||
[1;32m 1290[0m incomplete[38;5;241m=[39mincomplete,
|
||||
[1;32m 1291[0m psi0[38;5;241m=[39mpsi0,
|
||||
[1;32m 1292[0m aux_states[38;5;241m=[39maux_states,
|
||||
[1;32m 1293[0m stoc_proc[38;5;241m=[39mstoc_proc,
|
||||
[1;32m 1294[0m result_type[38;5;241m=[39m[38;5;28mself[39m[38;5;241m.[39mparams[38;5;241m.[39mHiP[38;5;241m.[39mresult_type,
|
||||
[1;32m 1295[0m normed[38;5;241m=[39m[38;5;28mself[39m[38;5;241m.[39m_normed_average,
|
||||
[1;32m 1296[0m rng_seed[38;5;241m=[39mseed,
|
||||
[1;32m 1297[0m )
|
||||
|
||||
File [0;32m~/src/hops/hops/core/signal_delay.py:87[0m, in [0;36msig_delay.__exit__[0;34m(self, exc_type, exc_val, exc_tb)[0m
|
||||
[1;32m 84[0m [38;5;28;01mif[39;00m [38;5;28mlen[39m([38;5;28mself[39m[38;5;241m.[39msigh[38;5;241m.[39msigs_caught) [38;5;241m>[39m [38;5;241m0[39m [38;5;129;01mand[39;00m [38;5;28mself[39m[38;5;241m.[39mhandler [38;5;129;01mis[39;00m [38;5;129;01mnot[39;00m [38;5;28;01mNone[39;00m:
|
||||
[1;32m 85[0m [38;5;28mself[39m[38;5;241m.[39mhandler([38;5;28mself[39m[38;5;241m.[39msigh[38;5;241m.[39msigs_caught)
|
||||
[0;32m---> 87[0m [38;5;28;43mself[39;49m[38;5;241;43m.[39;49m[43m_restore[49m[43m([49m[43m)[49m
|
||||
|
||||
File [0;32m~/src/hops/hops/core/signal_delay.py:68[0m, in [0;36msig_delay._restore[0;34m(self)[0m
|
||||
[1;32m 66[0m [38;5;28;01mfor[39;00m i, s [38;5;129;01min[39;00m [38;5;28menumerate[39m([38;5;28mself[39m[38;5;241m.[39msigs):
|
||||
[1;32m 67[0m signal[38;5;241m.[39msignal(s, [38;5;28mself[39m[38;5;241m.[39mold_handlers[i])
|
||||
[0;32m---> 68[0m [38;5;28;43mself[39;49m[38;5;241;43m.[39;49m[43msigh[49m[38;5;241;43m.[39;49m[43memit[49m[43m([49m[43m)[49m
|
||||
|
||||
File [0;32m~/src/hops/hops/core/signal_delay.py:42[0m, in [0;36mSigHandler.emit[0;34m(self)[0m
|
||||
[1;32m 40[0m [38;5;28;01mfor[39;00m s [38;5;129;01min[39;00m [38;5;28mself[39m[38;5;241m.[39msigs_caught:
|
||||
[1;32m 41[0m log[38;5;241m.[39minfo([38;5;124m"[39m[38;5;124memit signal [39m[38;5;124m'[39m[38;5;132;01m{}[39;00m[38;5;124m'[39m[38;5;124m"[39m[38;5;241m.[39mformat(SIG_MAP[s]))
|
||||
[0;32m---> 42[0m [43mos[49m[38;5;241;43m.[39;49m[43mkill[49m[43m([49m[43mos[49m[38;5;241;43m.[39;49m[43mgetpid[49m[43m([49m[43m)[49m[43m,[49m[43m [49m[43ms[49m[43m)[49m
|
||||
|
||||
[0;31mKeyboardInterrupt[0m:
|
||||
#+end_example
|
||||
:END:
|
||||
|
||||
* Flow
|
||||
#+begin_src jupyter-python :results none
|
||||
|
|
@ -202,8 +315,8 @@
|
|||
#+RESULTS:
|
||||
:RESULTS:
|
||||
: WARNING:matplotlib.legend:No artists with labels found to put in legend. Note that artists whose label start with an underscore are ignored when legend() is called with no argument.
|
||||
: <matplotlib.legend.Legend at 0x7ff3f665a1f0>
|
||||
[[file:./.ob-jupyter/eae302d599fab252f5be9982c699f16f4bf0d274.svg]]
|
||||
: <matplotlib.legend.Legend at 0x7ff3f665a130>
|
||||
[[file:./.ob-jupyter/b0eb1c317bdf1184cddc45ce8f2d1b5a7e1d661c.svg]]
|
||||
:END:
|
||||
|
||||
* Analytic
|
||||
|
|
@ -250,6 +363,32 @@
|
|||
exact_flows.append(exact_flow)
|
||||
#+end_src
|
||||
|
||||
#+begin_src jupyter-python
|
||||
params = multi_params[0]
|
||||
tt = params.IntP.t
|
||||
α_0 = util.BCF(
|
||||
params.IntP.t_max,
|
||||
hops.util.bcf.OhmicBCF_zeroTemp(s=1, eta=1, w_c=keys["wc"]),
|
||||
num_terms=7,
|
||||
resolution=0.01,
|
||||
)
|
||||
α_alt = gf.BCF(
|
||||
params.IntP.t_max,
|
||||
factors=params.SysP.g[0] / params.SysP.bcf_scale[0],
|
||||
exponents=params.SysP.w[0],
|
||||
)
|
||||
|
||||
pu.plot_complex(
|
||||
tt,
|
||||
hops.util.bcf.OhmicBCF_zeroTemp(s=1, eta=1, w_c=keys["wc"])(tt) - α_alt.approx(tt),
|
||||
)
|
||||
#+end_src
|
||||
|
||||
#+RESULTS:
|
||||
:RESULTS:
|
||||
| <Figure | size | 640x440 | with | 1 | Axes> | <AxesSubplot:> |
|
||||
[[file:./.ob-jupyter/8359947a3fe7de0ee3390e3803949d14be4af47f.svg]]
|
||||
:END:
|
||||
|
||||
#+begin_src jupyter-python
|
||||
fig, ax = plt.subplots()
|
||||
|
|
@ -283,7 +422,10 @@
|
|||
#+end_src
|
||||
|
||||
#+RESULTS:
|
||||
[[file:./.ob-jupyter/6f4f530dfd0f24c23fbf75506a97fb37981c8475.svg]]
|
||||
:RESULTS:
|
||||
: WARNING:matplotlib.legend:No artists with labels found to put in legend. Note that artists whose label start with an underscore are ignored when legend() is called with no argument.
|
||||
[[file:./.ob-jupyter/7ad76b73dba0620873d9a9ceb124cffd520ecafd.svg]]
|
||||
:END:
|
||||
|
||||
#+begin_src jupyter-python
|
||||
pu.plot_convergence(
|
||||
|
|
|
|||
|
|
@ -109,7 +109,7 @@ def ho_duo(
|
|||
return params
|
||||
|
||||
|
||||
model_keys = [dict(bcf_scale=0.2, wc=ω_c, γ=.5, T=.6, k_max=4, bcf_terms=5, max_HO_level=5) for ω_c in [2]]
|
||||
model_keys = [dict(bcf_scale=0.2, wc=ω_c, γ=.5, T=.6, k_max=3, bcf_terms=5, max_HO_level=9) for ω_c in [2]]
|
||||
multi_params = [ho_duo(**keys) for keys in model_keys]
|
||||
multi_params = multi_params
|
||||
|
||||
|
|
@ -121,7 +121,7 @@ supervisors = []
|
|||
for params in multi_params:
|
||||
supervisor = HOPSSupervisor(
|
||||
params,
|
||||
10_000,
|
||||
2_000,
|
||||
data_path="ho_data",
|
||||
data_name="two_baths",
|
||||
)
|
||||
|
|
@ -200,6 +200,25 @@ for params, keys in zip(multi_params, model_keys):
|
|||
exact_flow = [C.flow(params.IntP.t, 0, False), C.flow(params.IntP.t, 1, False)]
|
||||
exact_flows.append(exact_flow)
|
||||
|
||||
params = multi_params[0]
|
||||
tt = params.IntP.t
|
||||
α_0 = util.BCF(
|
||||
params.IntP.t_max,
|
||||
hops.util.bcf.OhmicBCF_zeroTemp(s=1, eta=1, w_c=keys["wc"]),
|
||||
num_terms=7,
|
||||
resolution=0.01,
|
||||
)
|
||||
α_alt = gf.BCF(
|
||||
params.IntP.t_max,
|
||||
factors=params.SysP.g[0] / params.SysP.bcf_scale[0],
|
||||
exponents=params.SysP.w[0],
|
||||
)
|
||||
|
||||
pu.plot_complex(
|
||||
tt,
|
||||
hops.util.bcf.OhmicBCF_zeroTemp(s=1, eta=1, w_c=keys["wc"])(tt) - α_alt.approx(tt),
|
||||
)
|
||||
|
||||
fig, ax = plt.subplots()
|
||||
for params, flow, ex_flow, keys in zip(multi_params, flow_hops, exact_flows, model_keys):
|
||||
consistency = (-1 * flow.for_bath(0)).consistency(ex_flow[0])
|
||||
|
|
|
|||
Loading…
Add table
Reference in a new issue