def test_vectorization(self):
a = flib.gev([4,5,6], xi=.2,mu=4,sigma=1)
b = flib.gev([4,5,6], xi=[.2,.2,.2],mu=4,sigma=1)
c = flib.gev([4,5,6], xi=.2,mu=[4,4,4],sigma=1)
d = flib.gev([4,5,6], xi=.2,mu=4,sigma=[1,1,1])
assert_equal(a,b)
assert_equal(b,c)
assert_equal(c,d)
def test_vectorization(self):
a = flib.gev([4, 5, 6], xi=.2, mu=4, sigma=1)
b = flib.gev([4, 5, 6], xi=[.2, .2, .2], mu=4, sigma=1)
c = flib.gev([4, 5, 6], xi=.2, mu=[4, 4, 4], sigma=1)
d = flib.gev([4, 5, 6], xi=.2, mu=4, sigma=[1, 1, 1])
assert_equal(a, b)
assert_equal(b, c)
assert_equal(c, d)
def test_with_scipy(self):
if not SP:
raise nose.SkipTest("SciPy not installed.")
x = [1,2,3,4]
scipy_y = log(genextreme.pdf(x, -.3, 4, 2))
flib_y = []
for i in x:
flib_y.append(flib.gev(i, .3, 4, 2))
assert_array_almost_equal(scipy_y,flib_y,5)
def test_with_scipy(self):
if not SP:
raise nose.SkipTest("SciPy not installed.")
x = [1, 2, 3, 4]
scipy_y = log(genextreme.pdf(x, -.3, 4, 2))
flib_y = []
for i in x:
flib_y.append(flib.gev(i, .3, 4, 2))
assert_array_almost_equal(scipy_y, flib_y, 5)
def test_limit(self):
x = [1,2,3]
a = flib.gev(x, 0.00001, 0, 1)
b = flib.gev(x, 0, 0, 1)
assert_almost_equal(a,b,4)
def test_limit(self):
x = [1, 2, 3]
a = flib.gev(x, 0.00001, 0, 1)
b = flib.gev(x, 0, 0, 1)
assert_almost_equal(a, b, 4)
