def test_compute_phi(): rx, x = compute_phi(traj, opt=True) ry, y = compute_phi(traj, opt=False) eq(rx, ry) eq(x, y)