@@ -0,0 +1,556 @@
+from __future__ import print_function, division%0A%0Aimport tempfile%0Aimport os%0A%0Aimport numpy as np%0Aimport vip_hci as vip%0A%0A%0Adef test_saveable_dataset():%0A    cube = np.zeros((5, 10, 10))%0A    angles = np.linspace(1, 2, 5)%0A    fwhm = 4  # test non-numpy type saving/loading%0A%0A    ds = vip.HCIDataset(cube=cube, angles=angles, fwhm=fwhm)%0A%0A    fd, fn = tempfile.mkstemp(prefix=%22vip_%22)%0A%0A    ds.save(fn)%0A%0A    ds2 = vip.HCIDataset.load(fn)%0A%0A    assert np.allclose(ds2.cube, cube)%0A    assert np.allclose(ds2.angles, angles)%0A    assert ds2.fwhm == fwhm%0A%0A    os.remove(fn)%0A

add test