Changeset 16141 for branches/pap_branch_080117/psLib/src/fits/psFits.c

Timestamp:

Jan 18, 2008, 3:04:13 PM (19 years ago)

Author:

Paul Price

Message:

Moving all options on FITS I/O into an additional structure, psFitsOptions. The idea is that this can be carried around independently of the FITS structure, which is only created when opening a file.

File:

• branches/pap_branch_080117/psLib/src/fits/psFits.c (modified) (10 diffs)

branches/pap_branch_080117/psLib/src/fits/psFits.c

@@ -7 +7 @@
  *  @author Robert DeSonia, MHPCC
  *
- *  @version $Revision: 1.76.2.2 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-01-18 03:37:50 $
+ *  @version $Revision: 1.76.2.3 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-01-19 01:04:13 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -78 +78 @@
         fitsClose(fits);
     }
-    psFree (fits->extword);
+    psFree(fits->options);
 }
 
@@ -160 +160 @@
 
     psFits* fits = psAlloc(sizeof(psFits));
+    psMemSetDeallocator(fits, (psFreeFunc)fitsFree);
+
     fits->fd = fptr;
     fits->writable = (iomode == READWRITE);
-    fits->extword = NULL;
-    fits->conventions.compression = true;
-    fits->conventions.psBitpix = true;
-
-    fits->floatType = PS_FITS_FLOAT_NONE;
-    fits->bitpix = 0;
-    fits->scaling = PS_FITS_SCALE_NONE;
-    fits->fuzz = true;
-    fits->bscale = 1.0;
-    fits->bzero = 0.0;
-    fits->mean = NAN;
-    fits->stdev = NAN;
-    fits->stdevBits = 4;
-    fits->stdevNum = 5.0;
-
-    psMemSetDeallocator(fits,(psFreeFunc)fitsFree);
+
+    fits->options = NULL;
 
     return fits;
+}
+
+
+static void fitsOptionsFree(psFitsOptions *options)
+{
+    psFree(options->extword);
+}
+
+
+psFitsOptions *psFitsOptionsAlloc(void)
+{
+    psFitsOptions *options = psAlloc(sizeof(psFitsOptions)); // Options, to return
+    psMemSetDeallocator(options, (psFreeFunc)fitsOptionsFree);
+
+    options->extword = NULL;
+
+    options->conventions.compression = true;
+    options->conventions.psBitpix = true;
+
+    options->floatType = PS_FITS_FLOAT_NONE;
+
+    options->bitpix = 0;
+
+    options->scaling = PS_FITS_SCALE_NONE;
+    options->fuzz = true;
+    options->bscale = 1.0;
+    options->bzero = 0.0;
+    options->mean = NAN;
+    options->stdev = NAN;
+    options->stdevBits = 4;
+    options->stdevNum = 5.0;
+
+    return options;
 }
 
@@ -237 +258 @@
 }
 
-bool psFitsSetExtnameWord (psFits *fits, const char *extword)
+bool psFitsSetExtnameWord(psFits *fits, const char *extword)
 {
     PS_ASSERT_PTR_NON_NULL(fits,    false);
     PS_ASSERT_PTR_NON_NULL(extword, false);
 
-    psFree (fits->extword);
-    fits->extword = psStringCopy (extword);
+    if (!fits->options) {
+        fits->options = psFitsOptionsAlloc();
+    }
+
+    psFree(fits->options->extword);
+    fits->options->extword = psStringCopy(extword);
     return true;
 }
@@ -261 +286 @@
     int status = 0;
 
-    if (!fits->conventions.compression && !fits->extword) {
+    psFitsOptions *options = fits->options; // FITS options
+    if (options && !options->conventions.compression && !options->extword) {
         // User wants to use cfitsio.  Good luck to them!
         if (fits_movnam_hdu(fits->fd, ANY_HDU, (char*)extname, 0, &status) != 0) {
@@ -270 +296 @@
     }
 
-    bool ignoreCI = (fits->conventions.compression &&
+    bool ignoreCI = (options->conventions.compression &&
                      (strcmp(extname, "COMPRESSED_IMAGE") != 0)); // Ignore COMPRESSED_IMAGE extension name?
-    char *extword = (fits->extword ? fits->extword : "EXTNAME"); // Word to use as extension name
+    char *extword = (options->extword ? options->extword : "EXTNAME"); // Word to use as extension name
 
 #if 0
@@ -375 +401 @@
     char name[MAX_STRING_LENGTH];
 
-    char *extword = (fits->extword == NULL) ? defaultExtword : fits->extword;
+    psFitsOptions *options = fits->options; // FITS options
+    char *extword = (!options || !options->extword) ? defaultExtword : options->extword;
 
     if (fits_read_key_str(fits->fd, extword, name, NULL, &status) != 0) {
-        psError(PS_ERR_BAD_PARAMETER_NULL, true,
-                _("Header keyword %s is not found"), extword);
+        psError(PS_ERR_BAD_PARAMETER_NULL, true, _("Header keyword %s is not found"), extword);
         return NULL;
     }
@@ -392 +418 @@
     int status = 0;
 
-    char *extword = (fits->extword == NULL) ? defaultExtword : fits->extword;
+    psFitsOptions *options = fits->options; // FITS options
+    char *extword = (!options || !options->extword) ? defaultExtword : options->extword;
 
     if (fits_update_key_str(fits->fd, extword, (char*)name, NULL, &status) != 0) {
         char fitsErr[MAX_STRING_LENGTH];
         (void)fits_get_errstatus(status, fitsErr);
-        psError(PS_ERR_IO, true,
-                _("Could not write data to file. CFITSIO Error: %s"),
-                fitsErr);
+        psError(PS_ERR_IO, true, _("Could not write data to file. CFITSIO Error: %s"), fitsErr);
         return false;
     }
@@ -414 +439 @@
 
     // move to the specified HDU
-    if (! psFitsMoveExtNum(fits,extnum,relative) ) {
-        psError(PS_ERR_BAD_PARAMETER_VALUE, false,
-                "Failed to delete HDU #%d",
-                extnum);
+    if (!psFitsMoveExtNum(fits, extnum, relative) ) {
+        psError(PS_ERR_BAD_PARAMETER_VALUE, false, "Failed to delete HDU #%d", extnum);
         return false;
     }
@@ -427 +450 @@
         char fitsErr[MAX_STRING_LENGTH];
         (void)fits_get_errstatus(status, fitsErr);
-        psError(PS_ERR_IO, true,
-                _("Could not write data to file. CFITSIO Error: %s"),
-                fitsErr);
+        psError(PS_ERR_IO, true, _("Could not write data to file. CFITSIO Error: %s"), fitsErr);
         return false;
     }