'Test that np.func(jd, **kw) is the same as t1.func(**kw) where
jd is a similarly shaped array and the same integral values.'

@pytest.mark.parametrize('kw,    func', itertools.product(kwargs, functions))
def test_funcs(self, kw, func):

    t1v = getattr(self.t1, func)(**kw)
    jdv = getattr(np, func)(self.jd, **kw)
    assert np.all((t1v.value == jdv))
    assert (t1v.shape == jdv.shape)
