Closed FReina closed 3 years ago
I found out I made a terrible typo at some point while writing the pseudocode, and never checked again. The equations for the Confined:
def __call__(self, t, D_micro, delta, tau): return D_micro * (tau/t) * (1 - np.exp(-tau/t)) + \ delta ** 2 / (2 * t * (1 - 2 * self.R * self.dt / t))
and Hop model
def __call__(self, t, D_macro, D_micro, delta, tau): return D_macro + \ D_micro * (tau/t) * (1 - np.exp(-tau/t)) + \ delta ** 2 / (2 * t * (1 - 2 * self.R * self.dt / t))
should instead have (1 - np.exp(-t/tau)).
(1 - np.exp(-t/tau))
I found out while discussing the paper I am finishing with Chris. This should solve many problems of convergence. I will open another issue for the bounds.
I found out I made a terrible typo at some point while writing the pseudocode, and never checked again. The equations for the Confined:
def __call__(self, t, D_micro, delta, tau): return D_micro * (tau/t) * (1 - np.exp(-tau/t)) + \ delta ** 2 / (2 * t * (1 - 2 * self.R * self.dt / t))
and Hop model
def __call__(self, t, D_macro, D_micro, delta, tau): return D_macro + \ D_micro * (tau/t) * (1 - np.exp(-tau/t)) + \ delta ** 2 / (2 * t * (1 - 2 * self.R * self.dt / t))
should instead have
(1 - np.exp(-t/tau))
.I found out while discussing the paper I am finishing with Chris. This should solve many problems of convergence. I will open another issue for the bounds.