tols

Author: kaelyndunnell

src/hisp/model.py

@@ -228,7 +228,7 @@ def B_stepsize(self, t: float) -> float:
                 else:
                     value = pulse.duration_no_waiting / 10
             else:
-                value = pulse.duration_no_waiting / 1000  # usually /100
+                value = pulse.duration_no_waiting / 100  # usually /100
         return periodic_step_function(
             relative_time,
             period_on=pulse.duration_no_waiting,
@@ -331,9 +331,9 @@ def make_custom_rtol(self, t: float) -> float:
         pulse = self.scenario.get_pulse(t)
         relative_time = t - self.scenario.get_time_start_current_pulse(t)
         if pulse.pulse_type == "GDC" or pulse.pulse_type == "ICWC":
-            rtol = 1e-11  # 1e-12 for just glow scenario
+            rtol = 1e-12  # 1e-12 for just glow scenario
         elif pulse.pulse_type == "BAKE":
-            rtol = 1e-11  # 1e-15 for just glow scenario
+            rtol = 1e-12  # 1e-15 for just glow scenario
         elif pulse.pulse_type == "FP" or pulse.pulse_type == "FP_D":
             # rtol = 1e-10
             if relative_time % pulse.total_duration > pulse.duration_no_waiting:
@@ -343,7 +343,7 @@ def make_custom_rtol(self, t: float) -> float:
             else:
                 rtol = 1e-8  # 1e-6 or 1e-8 most recently
         elif pulse.pulse_type == "RISP":
-            rtol = 1e-6 # usually 1e-7
+            rtol = 1e-5 # usually 1e-7
         else:
             rtol = 1e-8 # usually 1e-10
         return rtol