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)
def factorial(value):
"""Uses scipy.special.gamma."""
from scipy.special import gamma as gamma_
return gamma_(value + 1)
