def test_logprior(self):
        for theta in self.prior.sample(5, 0):
            g1 = spop.approx_fprime(theta, self.prior.get_logprior, 1e-8)
            _, g2 = self.prior.get_logprior(theta, True)
            nt.assert_allclose(g1, g2, rtol=1e-6)