def test_pqw_returns_dimensionless(): r_0 = ([1, 0, 0] * u.au).to(u.km) v_0 = ([0, 6, 0] * u.au / u.year).to(u.km / u.day) ss = RVState(Sun, r_0, v_0) p, q, w = ss.pqw() assert p.unit == u.one assert q.unit == u.one assert w.unit == u.one