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

Timestamp:

Jun 4, 2015, 8:49:25 AM (11 years ago)

Author:

eugene

Message:

rice and gzip2 both work for images; working on tables

File:

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

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

@@ +1 @@
+input tap.dvo
 
 macro test_image_all 
 
-  # foreach mode NONE GZIP_1 RICE_1
+  tapPLAN 88
+
+  foreach mode NONE GZIP_1 GZIP_2 RICE_1
   # foreach mode GZIP_1
-  foreach mode GZIP_2
+  # foreach mode GZIP_2
   # foreach mode RICE_1
     foreach bitpix 8 16 32 -32 -64
     # foreach bitpix 32
-      echo ===== bitpix = $bitpix =====
+      if (($mode == RICE_1) && ($bitpix < 0)) continue ; # RICE_1 is for int only
+      tapDIAG 2 "===== bitpix = $bitpix, mode = $mode ====="
       test_image_single $bitpix $mode
     end
   end
+
+  tapDONE
 end
 
@@ -30 +36 @@
  rd b test.cmp.fits -x 0
  set d = x - b
- stat d
+ stat -q d
+ tapOK {abs($MEAN)  < 0.001} "read our cmp $cmpmode $bitpix (MEAN  = $MEAN)"
+ tapOK {$SIGMA      < 0.001} "read our cmp $cmpmode $bitpix (SIGMA = $SIGMA)"
 
  if ($cmpmode == NONE) return
@@ -43 +51 @@
    exec fpack -r -S test.raw.fits > test.fpk.fits
  end
- echo "load fpack version"
+ # echo "load fpack version"
  rd c test.fpk.fits -x 0
  set d = x - c
- stat d
+ stat -q d
+ tapOK {abs($MEAN)  < 0.001} "read fpack $cmpmode $bitpix (MEAN  = $MEAN)"
+ tapOK {$SIGMA      < 0.001} "read fpack $cmpmode $bitpix (SIGMA = $SIGMA)"
 
  if ($cmpmode == RICE_1_0) 
@@ -53 +63 @@
  end
 
- echo "run funpack on our version"
+ # echo "run funpack on our version"
  exec funpack -S test.cmp.fits > test.fun.fits
  rd e test.fun.fits
  set d = x - e
- stat d
+ stat -q d
+ tapOK {abs($MEAN)  < 0.001} "funpack our cmp $cmpmode $bitpix (MEAN  = $MEAN)"
+ tapOK {$SIGMA      < 0.001} "funpack our cmp $cmpmode $bitpix (SIGMA = $SIGMA)"
 end