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