def METHOD_NAME(self):
    """
    Regression test for https://github.com/astropy/astropy/issues/4587
    """
    n = np.ones((10, 10), dtype="float32")
    comp_hdu = fits.CompImageHDU(n)
    comp_hdu.writeto(self.temp("tmp.fits"), checksum=True)
    with fits.open(self.temp("tmp.fits"), checksum=True) as hdul:
        assert np.all(hdul[1].data == comp_hdu.data)