tols

Author: kaelyndunnell

src/hisp/model.py

@@ -217,7 +217,7 @@ def B_stepsize(self, t: float) -> float:
             # the stepsize is 1/10 of the duration of the pulse
             if pulse.pulse_type == "FP":
                 if relative_time_within_sub_pulse < pulse.duration_no_waiting:
-                    value = 0.01  # s # usually 0.01
+                    value = 0.0001  # s # usually 0.01
                 else:
                     value = pulse.duration_no_waiting / 10
             elif pulse.pulse_type == "BAKE":
@@ -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 / 10000  # usually /100
         return periodic_step_function(
             relative_time,
             period_on=pulse.duration_no_waiting,
@@ -331,17 +331,17 @@ 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-13  # 1e-12 for just glow scenario
+            rtol = 1e-11  # 1e-12 for just glow scenario
         elif pulse.pulse_type == "BAKE":
-            rtol = 1e-12  # 1e-15 for just glow scenario
+            rtol = 1e-11  # 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:
                 rtol = 1e-12  # 1e-12
             elif pulse.pulse_type == "FP_D":
                 rtol = 1e-9
             else:
-                rtol = 1e-14  # 1e-6 or 1e-8 most recently
+                rtol = 1e-12  # 1e-6 or 1e-8 most recently
         elif pulse.pulse_type == "RISP":
             rtol = 1e-7
         else: