def test_J_with_lam_and_zero_theta(): l = Logysterical(tr) theta = np.zeros((l.tr["X"].shape[1], 1)) assert l.j(theta) - 0.69314718056 < epsilon