def METHOD_NAME():
    x = pt.vector("x", shape=(1,))
    p = pt.scalar("p")
    xvals = [np.array([v], dtype=config.floatX) for v in [0.1, 2, 5, 10, 50, 100]]
    multigammaln_ = function([x, p], multigammaln(x, p), mode="FAST_COMPILE")
    def ref_multigammaln(a, b):
        return np.array(scipy.special.multigammaln(a[0], b), config.floatX)
    for p in [0, 1, 2, 3, 4, 100]:
        for x in xvals:
            if np.all(x > 0.5 * (p - 1)):
                check_vals(multigammaln_, ref_multigammaln, x, p)