'Regression test for https://github.com/spacetelescope/PyFITS/issues/71'

def test_compression_with_gzip_column(self):

    arr = np.zeros((2, 7000), dtype='float32')
    arr[0] = np.linspace(0, 1, 7000)
    arr[1] = np.random.normal(size=7000)
    hdu = fits.CompImageHDU(data=arr)
    hdu.writeto(self.temp('test.fits'))
    with fits.open(self.temp('test.fits')) as hdul:
        comp_hdu = hdul[1]
        assert np.all((comp_hdu.data[0] == arr[0]))
