def test_gamma_visual(self): x = np.linspace(0.1, 10) y = gamma(x) z = gamma_(x) plt.figure(figsize=(10, 6)) sb.lineplot(x, y, label="Local") sb.lineplot(x, z, label="Scipy") plt.show()
def gamma(x): #if(x <= 0): return float("nan") return gamma_(x)
def test_running_gamma(self): x = np.linspace(0.1, 10) y = gamma(x) z = gamma_(x) self.assertAlmostEqual(np.square(y - z).mean(), 0.0, delta=1e-2)
def factorial(value): """Uses scipy.special.gamma.""" from scipy.special import gamma as gamma_ return gamma_(value + 1)