def test_coma_equivalent_radius(self): shape = (0.8, 2) * u.arcsec aper = RectangularAperture(shape) r = aper.coma_equivalent_radius() assert np.isclose(r.value, 0.66776816346357259)