Fix bug in effTime time integration scheme

nicodirkes · Feb 26, 2024 · 2814668 · 2814668
1 parent 5668d7c
commit 2814668
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/hemtracer/hemolysis_model/hemolysis_model.py b/hemtracer/hemolysis_model/hemolysis_model.py
@@ -200,9 +200,10 @@ def _compute_HI5(self, t: NDArray, G: NDArray) -> NDArray:
         dt = np.diff(t)
         C_tau_beta = self._C * tau**self._beta
 
-        for i in range(1, len(t)):
+        for i in range(1, len(t)-1):
             t_eff = (IH[i-1] / C_tau_beta[i-1])**(1/self._alpha)
             IH[i] = C_tau_beta[i] * (t_eff + dt[i-1])**self._alpha
+        IH[-1] = IH[-2]
 
         return IH