How to deal with the negative value appearing in the np.log() and np.sqrt() function when calculating density about pymle HOT 1 CLOSED

Thanks for your great question! I have rewritten the Ozaki density function to make clear that each of these terms is indeed well-defined (this form is also more efficient :). You can see that (exp(x*t) -1)/x is always positive for positive t ):

    sig = self._model.diffusion(x0, t)
    mu = self._model.drift(x0, t)
    mu_x = self._model.drift_x(x0, t)
    temp = mu * (np.exp(mu_x * t) - 1) / mu_x

    Mt = x0 + temp
    Kt = (2 / t) * np.log(1 + temp / x0)
    Vt = sig * np.sqrt((np.exp(Kt * t) - 1) / Kt)

    return np.exp(-0.5 * ((xt - Mt) / Vt) ** 2) / (np.sqrt(2 * np.pi) * Vt)

Also, I removed the if z <=0 expression from above, as you mention this is not going to be optimal, and its better that if something like this happens we should return a nan value, allowing the optimizer to adapt the parameters accordingly

from pymle.