make coupling stronger

python/otto_motor/subprojects/cycle_length_coupling_strength/project.org

@@ -303,96 +303,26 @@ coupling-change/cycle time.
 #+RESULTS:
 
 #+begin_src jupyter-python :results none :tangle tangle/weak_coupling_model.py
-  weak_coupling_model = make_model(50, .01)
-  weak_coupling_model.num_cycles = 15
-  weak_coupling_model.k_max = 2
+  weak_coupling_model = make_model(50, .2)
+  weak_coupling_model.num_cycles = 5
 #+end_src
 
 
 ** Integrate
 #+begin_src jupyter-python  :tangle tangle/weak_coupling_model.py
-  ot.integrate_online(weak_coupling_model, 100, analyze_kwargs=dict(every=10_000))
+  ot.integrate_online(weak_coupling_model, 50_000, analyze_kwargs=dict(every=10_000))
 #+end_src
 
 #+RESULTS:
-:RESULTS:
-#+begin_example
-  [INFO    hops.core.integration     290962] Choosing the nonlinear integrator.
-  [INFO    root                      290962] Starting analysis process.
-  [INFO    root                      290962] Started analysis process with pid 293338.
-  [INFO    hops.core.hierarchy_data  290962] Creating the streaming fifo at: /home/hiro/Documents/Projects/UNI/master/eflow_paper/python/otto_motor/subprojects/cycle_length_coupling_strength/results_ef8e3679c35345f6c603306d9c116e7f72af929efd9037b1f3858e1f3d306a87.fifo
-  [INFO    hops.core.integration     290962] Using 16 integrators.
-  [INFO    hops.core.integration     290962] Some 100 trajectories have to be integrated.
-  [INFO    hops.core.integration     290962] Using 66 hierarchy states.
-    0% 0/100 [00:00<?, ?it/s][INFO    hops.core.signal_delay    290962] caught sig 'SIGINT'
-
-    0% 0/100 [00:10<?, ?it/s]
-  2023-07-10 14:21:28,916       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.
-  2023-07-10 14:21:28,917       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.
-  2023-07-10 14:21:28,919       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.
-  2023-07-10 14:21:28,920       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.
-  [INFO    hops.core.signal_delay    290962] caught 1 signal(s)
-  2023-07-10 14:21:28,921       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.
-  [INFO    hops.core.signal_delay    290962] emit signal 'SIGINT'
-  [INFO    hops.core.signal_delay    290962] caught sig 'SIGINT'
-  2023-07-10 14:21:28,924       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.
-  2023-07-10 14:21:28,927       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.
-  2023-07-10 14:21:28,928       WARNING worker.py:1404 -- A worker died or was killed while executing a task by an unexpected system error. To troubleshoot the problem, check the logs for the dead worker. RayTask ID: 42eb0de981b10f9ded44f707c2238576ac033f1501000000 Worker ID: ab3bee80797ba4d55a20a02edba9fc252c2f48eed15636ec6543f53b Node ID: 6f18af98e8f1029efc0c6aab222a211f324ccc8890393f00385bece8 Worker IP address: 10.121.151.166 Worker port: 37335 Worker PID: 291261
-  2023-07-10 14:21:28,961       WARNING worker.py:1404 -- A worker died or was killed while executing a task by an unexpected system error. To troubleshoot the problem, check the logs for the dead worker. RayTask ID: c73273ba2de9e743dd3864a96c16b15a37f8a0bb01000000 Worker ID: 47c18d4d67e434163287de7d8bd0dad4d633903880d6e64886a1b419 Node ID: 6f18af98e8f1029efc0c6aab222a211f324ccc8890393f00385bece8 Worker IP address: 10.121.151.166 Worker port: 46419 Worker PID: 291262
-  2023-07-10 14:21:28,962       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.
-  2023-07-10 14:21:29,022       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.
-  2023-07-10 14:21:29,045       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.
-  2023-07-10 14:21:29,058       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.
-  2023-07-10 14:21:29,091       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.
-  2023-07-10 14:21:29,129       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.
-  2023-07-10 14:21:29,160       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.
-  2023-07-10 14:21:29,206       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.
-  2023-07-10 14:21:29,247       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.
-  [INFO    hops.core.signal_delay    290962] caught 1 signal(s)
-  [INFO    hops.core.signal_delay    290962] emit signal 'SIGINT'
-#+end_example
-# [goto error]
-#+begin_example
-  ---------------------------------------------------------------------------
-  KeyboardInterrupt                         Traceback (most recent call last)
-  Cell In[86], line 1
-  ----> 1 ot.integrate_online(weak_coupling_model, 100, analyze_kwargs=dict(every=10_000))
-
-  File ~/Documents/Projects/UNI/master/eflow_paper/python/otto_motor/subprojects/cycle_length_coupling_strength/otto_utilities.py:363, in integrate_online(model, n, stream_folder, **kwargs)
-      362 def integrate_online(model, n, stream_folder=None, **kwargs):
-  --> 363     aux.integrate(
-      364         model,
-      365         n,
-      366         stream_file=("" if stream_folder is None else stream_folder)
-      367         + f"results_{model.hexhash}.fifo",
-      368         analyze=True,
-      369         **kwargs,
-      370     )
-
-  File ~/src/two_qubit_model/hiro_models/model_auxiliary.py:201, in integrate(model, n, data_path, clear_pd, single_process, stream_file, analyze, results_path, analyze_kwargs)
-      199         supervisor.integrate_single_process(clear_pd)
-      200     else:
-  --> 201         supervisor.integrate(clear_pd)
-      203 cleanup(0)
-
-  File /nix/store/fai1b55231rnk4jyj0kjicdnqcgdf9ph-python3-3.9.15-env/lib/python3.9/site-packages/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 /nix/store/fai1b55231rnk4jyj0kjicdnqcgdf9ph-python3-3.9.15-env/lib/python3.9/site-packages/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 /nix/store/fai1b55231rnk4jyj0kjicdnqcgdf9ph-python3-3.9.15-env/lib/python3.9/site-packages/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:
+: [INFO    hops.core.integration     290962] Choosing the nonlinear integrator.
+: [INFO    root                      290962] Starting analysis process.
+: [INFO    root                      290962] Started analysis process with pid 304855.
+: [INFO    hops.core.hierarchy_data  290962] Creating the streaming fifo at: /home/hiro/Documents/Projects/UNI/master/eflow_paper/python/otto_motor/subprojects/cycle_length_coupling_strength/results_f9d1bda1cf5d0f7419d0cc1b1b5db082c13604538255d3dfbc05af9f0d354f49.fifo
+: [INFO    hops.core.integration     290962] Using 16 integrators.
+: [INFO    hops.core.integration     290962] Some 16 trajectories have to be integrated.
+: [INFO    hops.core.integration     290962] Using 1001 hierarchy states.
+: 100% 16/16 [01:28<00:00,  5.52s/it]
+:
 
 ** Analysis
 #+begin_src jupyter-python
@@ -402,28 +332,36 @@ coupling-change/cycle time.
 
 #+RESULTS:
 :RESULTS:
-| array | ((-0.03997036 -0.03997036 -0.03996733 ... -0.21910526 -0.21910526 -0.21910526)) |
-[[file:./.ob-jupyter/c1241b127ae3595367627adfdc9a1f6114cbc6af.svg]]
+| array | ((-0.18507824 -0.18507818 -0.18506804 ... -0.87570649 -0.87570633 -0.87570633)) |
+[[file:./.ob-jupyter/ed60941fac1fe754ef3b56d5f56b7b83a3d91291.svg]]
 :END:
+
+
 #+begin_src jupyter-python
-  weak_coupling_model.total_power().value[weak_coupling_model.strobe[1][7]:-1].mean()
 
   #weak_coupling_model.strobe[1]
 #+end_src
 
 #+RESULTS:
-: 0.0
+:RESULTS:
+# [goto error]
+: ---------------------------------------------------------------------------
+: IndexError                                Traceback (most recent call last)
+: Cell In[92], line 1
+: ----> 1 weak_coupling_model.total_power().value[weak_coupling_model.strobe[1][7]:-1].mean()
+:       3 #weak_coupling_model.strobe[1]
+:
+: IndexError: index 7 is out of bounds for axis 0 with size 6
+:END:
 
 
 
 #+begin_src jupyter-python
-weak_coupling_model.efficiency(steady_idx=7)
-weak_coupling_model.steady_bath_energy_change(steady_idx=7, bath=1)
-weak_coupling_model.steady_total_energy_change(steady_idx=7)
+weak_coupling_model.efficiency(steady_idx=-3)
 #+end_src
 
 #+RESULTS:
-: EnsembleValue([(100, 8.326672684688674e-17, 0.010617712466080408)])
+: EnsembleValue([(16, 0.3649132926357113, 0.21802876349189884)])
 
 
 * Things to Look At

python/otto_motor/subprojects/cycle_length_coupling_strength/tangle/coupling_strength.py

@@ -202,10 +202,6 @@ fs.export_fig("coupling_speed_scan_power_efficiency_uncertainty")
 ot.plot_energy(weak_coupling_model)
 weak_coupling_model.total_energy_from_power().slice(slice(weak_coupling_model.strobe[1][1], None, 1)).values
 
-weak_coupling_model.total_power().value[weak_coupling_model.strobe[1][7]:-1].mean()
-
 #weak_coupling_model.strobe[1]
 
-weak_coupling_model.efficiency(steady_idx=7)
-weak_coupling_model.steady_bath_energy_change(steady_idx=7, bath=1)
-weak_coupling_model.steady_total_energy_change(steady_idx=7)
+weak_coupling_model.efficiency(steady_idx=-3)

python/otto_motor/subprojects/cycle_length_coupling_strength/tangle/weak_coupling_model.py

@@ -44,8 +44,7 @@ def make_model(Θ, δ):
           L_shift=(0, 0),
       )
 
-weak_coupling_model = make_model(50, .01)
-weak_coupling_model.num_cycles = 15
-weak_coupling_model.k_max = 2
+weak_coupling_model = make_model(50, .2)
+weak_coupling_model.num_cycles = 5
 
-ot.integrate_online(weak_coupling_model, 100, analyze_kwargs=dict(every=10_000))
+ot.integrate_online(weak_coupling_model, 50_000, analyze_kwargs=dict(every=10_000))