Changeset 38367 for branches/eam_branches/ohana.20150429/src/opihi/lib.shell/VectorIO.c

Timestamp:

Jun 4, 2015, 5:44:39 PM (11 years ago)

Author:

eugene

Message:

need to provide non-swapping options to handle compression

File:

• branches/eam_branches/ohana.20150429/src/opihi/lib.shell/VectorIO.c (modified) (4 diffs)

branches/eam_branches/ohana.20150429/src/opihi/lib.shell/VectorIO.c

@@ -2 +2 @@
   
 // write a set of vectors to a FITS FTable structure (vectors names become fits column names)
-static int WriteVectorTable (FTable *ftable, char *extname, Vector **vec, int Nvec, char *format) {
+static int WriteVectorTable (FTable *ftable, char *extname, Vector **vec, int Nvec, char *format, char nativeOrder) {
   
   int j;
@@ -59 +59 @@
   for (j = 0; j < Nvec; j++) {
     if (vec[j][0].type == OPIHI_FLT) {
-      gfits_set_bintable_column_reformat (theader, ftable, vec[j][0].name, "double", vec[j][0].elements.Flt, vec[j][0].Nelements);
+      gfits_set_bintable_column_reformat (theader, ftable, vec[j][0].name, "double", vec[j][0].elements.Flt, vec[j][0].Nelements, nativeOrder);
     } else {
-      gfits_set_bintable_column_reformat (theader, ftable, vec[j][0].name, "int", vec[j][0].elements.Int, vec[j][0].Nelements);
+      gfits_set_bintable_column_reformat (theader, ftable, vec[j][0].name, "int", vec[j][0].elements.Int, vec[j][0].Nelements, nativeOrder);
     }
   }
@@ -130 +130 @@
 
   rawtable.header = &rawheader;
-  if (!WriteVectorTable (&rawtable, extname, vec, Nvec, format)) goto escape;
+  if (!WriteVectorTable (&rawtable, extname, vec, Nvec, format, (compress != NULL))) goto escape;
+  // NOTE : for compression, the table is constructed in native by order 
 
   FTable *outtable = &rawtable;
@@ -139 +140 @@
 
     dump_raw_table (&rawtable, "wd raw");
-    if (!gfits_compress_table (&rawtable, &cmptable, 10, compress)) goto escape;
+    if (!gfits_compress_table (&rawtable, &cmptable, 0, compress)) goto escape;
+
+    int i, Nfields;
+    if (!gfits_scan (&cmpheader, "TFIELDS", "%d", 1, &Nfields)) goto escape;
+    for (i = 0; i < Nfields; i++) {
+      gfits_byteswap_varlength_column (&cmptable, i+1);
+    }
     dump_cmp_table (&cmptable, "wd cmp");