def t_multigamma():
xvals = list(map(np.atleast_1d, [0, .1, 2, 100]))
x = t.dvector('x')
x.tag.test_value = xvals[0]
p = t.iscalar('p')
p.tag.test_value = 1
multigammaln = function([x, p], ps.multigammaln(x, p))
def ssmultigammaln(a, b):
return ss.multigammaln(a[0], b)
for p in [0, 1, 2, 3, 4, 100]:
for x in xvals:
yield check_vals, multigammaln, ssmultigammaln, x, p
def t_multigamma():
xvals = map(np.atleast_1d, [0, .1, 2, 100])
x = t.dvector('x')
x.tag.test_value = xvals[0]
p = t.iscalar('p')
p.tag.test_value = 1
multigammaln = function([x, p], ps.multigammaln(x, p))
def ssmultigammaln(a, b):
return ss.multigammaln(a[0], b)
for p in [0, 1, 2, 3, 4, 100]:
for x in xvals:
yield check_vals, multigammaln, ssmultigammaln, x, p
def test_functions():
xvals = list(map(np.atleast_1d, [.01, .1, 2, 100, 10000]))
x = t.dvector('x')
x.tag.test_value = xvals[0]
p = t.iscalar('p')
p.tag.test_value = 1
gammaln = function([x], ps.gammaln(x))
psi = function([x], ps.psi(x))
multigammaln = function([x, p], ps.multigammaln(x, p))
for x in xvals:
yield check_vals, gammaln, ss.gammaln, x
for x in xvals[1:]:
yield check_vals, psi, ss.psi, x
def test_functions():
xvals = map(np.atleast_1d, [.01, .1, 2, 100, 10000])
x = t.dvector('x')
x.tag.test_value = xvals[0]
p = t.iscalar('p')
p.tag.test_value = 1
gammaln = function([x], ps.gammaln(x))
psi = function([x], ps.psi(x))
multigammaln = function([x, p], ps.multigammaln(x, p))
for x in xvals:
yield check_vals, gammaln, ss.gammaln, x
for x in xvals[1:]:
yield check_vals, psi, ss.psi, x
