Changeset 38340 for branches/eam_branches/ohana.20150429/src/libfits/test/compress.sh

Timestamp:

May 30, 2015, 7:59:48 PM (11 years ago)

Author:

eugene

Message:

more work on compression : cfitsio assumes gzip, not zlib (I have made uncompress general )

File:

• branches/eam_branches/ohana.20150429/src/libfits/test/compress.sh (modified) (2 diffs)

branches/eam_branches/ohana.20150429/src/libfits/test/compress.sh

@@ +1 @@
+
+macro testall
+
+  foreach bitpix 8 16 32 -32 -64
+    echo ===== bitpix = $bitpix =====
+    test1 $bitpix
+  end
+end
 
 macro test1
- mcreate z 500 500
+ if ($0 != 2)
+   echo "USAGE: test_image (bitpix)"
+   break
+ end
+
+ mcreate z 8 8
  set x = xramp(z)
- wd x test.raw.fits -bitpix -32 -bzero 0 -bscale 1
- wd x test.cmp.fits -bitpix -32 -bzero 0 -bscale 1 -compress
+ wd x test.raw.fits -bitpix $1 -bzero 0 -bscale 1
+ wd x test.cmp.fits -bitpix $1 -bzero 0 -bscale 1 -compress
  rd b test.cmp.fits -x 0
  set d = x - b
  stat d
- exec fpack -g -S test.raw.fits > test.fpk.fits
- exec funpack -S test.cmp.fits > test.fun.fits
+
+ if ($bitpix == -64)
+   echo "*** bitpix $bitpix fails fpack / funpack ***"
+ else
+   exec fpack -g -S test.raw.fits > test.fpk.fits
+   echo "load fpack version"
+   rd c test.fpk.fits -x 0
+   set d = x - c
+   stat d
+ end
+
+ if ($bitpix == -64)
+   echo "*** bitpix $bitpix fails fpack / funpack ***"
+ else
+   exec funpack -S test.cmp.fits > test.fun.fits
+   rd e test.fun.fits
+   set d = x - e
+   stat d
+ end
 end
 
@@ -25 +55 @@
 macro test3
 
+ $NPT = 10
+ create x 0 $NPT
+ set y = x^2
+ create ID 1 {$NPT + 1} -int
+
+ write -fits test test.raw.tbl x y ID
+ write -fits test test.cmp.tbl x y ID -compress-mode NONE
+
+ foreach f x y ID
+   set $f\_raw = $f
+   delete $f
+ end
+
+ echo "===== raw ====="
+ data test.raw.tbl
+ read -fits test x y ID
+
+ foreach f x y ID
+   set dv = $f - $f\_raw
+   vstat dv
+ end
+ delete -q x y ID
+
+ echo "===== cmp ====="
+ data test.cmp.tbl
+ read -fits test x y ID
+
+ foreach f x y ID
+   set dv = $f - $f\_raw
+   vstat dv
+ end
+end
+
+macro test4
+
  create x 0 100
  set y = x^2
+ create ID 1 101 -int
 
- create ID 1 101 -int
  write -fits test test.raw.tbl x y ID
+ write -fits test test.cmp.tbl x y ID -compress
 
- write -fits test test.cmp.tbl x y ID -compress
+ foreach f x y ID
+   set $f\_raw = $f
+   delete $f
+ end
+
+ echo "===== raw ====="
+ data test.raw.tbl
+ read -fits test x y ID
+
+ foreach f x y ID
+   set dv = $f - $f\_raw
+   vstat dv
+ end
+ delete -q x y ID
+
+ echo "===== cmp ====="
+ data test.cmp.tbl
+ read -fits test x y ID
+
+ foreach f x y ID
+   set dv = $f - $f\_raw
+   vstat dv
+ end
 end