def nll(self, p):
    ll = np.sum((np.power((self.model(*p) - self.epsf), 2) * self.epsf))
    return (ll if np.isfinite(ll) else 1e+25)
def nll(self, p):
    ll = np.sum((np.power((self.model(*p) - self.epsf), 2) * self.epsf))
    return (ll if np.isfinite(ll) else 1e+25)