Merge pull request #76 from edger-asiimwe/kgut-parameter-refactor

Refactored kgut parameters alpha and beta
jxx123 · Mar 16, 2024 · 61919d7 · 61919d7
2 parents 80abe07 + b07a5d7
commit 61919d7
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/simglucose/patient/t1dpatient.py b/simglucose/patient/t1dpatient.py
@@ -133,8 +133,8 @@ def model(t, x, action, params, last_Qsto, last_foodtaken):
         dxdt[0] = -params.kmax * x[0] + d
 
         if Dbar > 0:
-            aa = 5 / 2 / (1 - params.b) / Dbar
-            cc = 5 / 2 / params.d / Dbar
+            aa = 5 / (2 * Dbar * (1 - params.b))
+            cc = 5 / (2 * Dbar * params.d)
             kgut = params.kmin + (params.kmax - params.kmin) / 2 * (
                 np.tanh(aa * (qsto - params.b * Dbar)) -
                 np.tanh(cc * (qsto - params.d * Dbar)) + 2)