def test_probit(): p = np.array([0.01, 0.25, 0.5, 0.75, 0.99]) np.testing.assert_allclose(invprobit(probit(p)).eval(), p, atol=1e-5)