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