def test_coma_equivalent_radius(self): sigma = 1 * u.arcsec aper = GaussianAperture(sigma) assert aper.coma_equivalent_radius() == np.sqrt(np.pi / 2) * u.arcsec