Changeset 38406 for branches/eam_branches/ohana.20150429/src/opihi/cmd.data

Timestamp:

Jun 6, 2015, 10:03:56 AM (11 years ago)

Author:

eugene

Message:

add options to cast to force single-precision floats or long long int (stuffed into double); fits table I/O puts int64 into opihi double vector

Location:

branches/eam_branches/ohana.20150429/src/opihi/cmd.data

Files:

• cast.c (modified) (4 diffs)
• read_vectors.c (modified) (2 diffs)
• vstats.c (modified) (1 diff)

branches/eam_branches/ohana.20150429/src/opihi/cmd.data/cast.c

@@ -17 +17 @@
   // out type 
   char outType = OPIHI_NOTYPE;
-  if (!strcasecmp(argv[5], "int")) outType = OPIHI_INT;
-  if (!strcasecmp(argv[5], "float")) outType = OPIHI_FLT;
+  if (!strcasecmp(argv[5], "int"))    outType = OPIHI_INT;
+  if (!strcasecmp(argv[5], "float"))  outType = OPIHI_FLT;
+  if (!strcasecmp(argv[5], "single")) outType = OPIHI_FLT;
+  if (!strcasecmp(argv[5], "int64"))  outType = OPIHI_FLT;
   if (outType == OPIHI_NOTYPE) goto usage;
+
+  int ForceSingle = FALSE;
+  if (!strcasecmp(argv[5], "single")) ForceSingle = TRUE;
+
+  int ForceI64 = FALSE;
+  if (!strcasecmp(argv[5], "int64")) ForceI64 = TRUE;
 
   if ((ovec = SelectVector (argv[1], ANYVECTOR, TRUE)) == NULL) goto error;
@@ -32 +40 @@
     opihi_flt *vo = ovec[0].elements.Flt;
     for (i = 0; i < ovec[0].Nelements; i++, vi++, vo++) {
-      *vo = *vi;
+      if (ForceSingle) {
+        float tmp = *vi;
+        *vo = tmp;
+      } else if (ForceI64) {
+        int64_t tmp = *vi;
+        *vo = tmp;
+      } else {
+        *vo = *vi;
+      }
     }
     return (TRUE);
@@ -48 +64 @@
     opihi_flt *vo = ovec[0].elements.Flt;
     for (i = 0; i < ovec[0].Nelements; i++, vi++, vo++) {
-      *vo = *vi;
+      if (ForceSingle) {
+        float tmp = *vi;
+        *vo = tmp;
+      } else {
+        *vo = *vi;
+      }
     }
     return (TRUE);
@@ -66 +87 @@
 
 usage:
-  gprint (GP_ERR, "SYNTAX: cast vec = vec as (int/float)\n");
+  gprint (GP_ERR, "SYNTAX: cast vec = vec as (int/float/single)\n");
+  gprint (GP_ERR, "  note: (single) forces vector to have single precision\n");
   return (FALSE);
 }

branches/eam_branches/ohana.20150429/src/opihi/cmd.data/read_vectors.c

@@ -487 +487 @@
   }
   if (start < 0) ESCAPE ("invalid range: start < 0\n");
-  if (start >= header.Naxis[1]) ESCAPE ("invalid range: start >= Ny (%d)\n", header.Naxis[1]);
+  if (start >= header.Naxis[1]) ESCAPE ("invalid range: start >= Ny ("OFF_T_FMT")\n", header.Naxis[1]);
   if (Nrows < 0) ESCAPE ("invalid range: Nrows < 0\n");
 
   // just a warning:
   if (start + Nrows > header.Naxis[1]) {
-    if (VERBOSE) gprint (GP_ERR, "NOTE: reading last block will return only %d rows\n", header.Naxis[1] - start);
+    if (VERBOSE) gprint (GP_ERR, "NOTE: reading last block will return only "OFF_T_FMT" rows\n", header.Naxis[1] - start);
   }
 
@@ -555 +555 @@
     
     vecType = OPIHI_INT;
-    if (!strcmp (type, "double") || !strcmp (type, "float")) {
+    if (!strcmp (type, "double") || !strcmp (type, "float") || !strcmp (type, "int64_t")) {
       vecType = OPIHI_FLT;
     }

branches/eam_branches/ohana.20150429/src/opihi/cmd.data/vstats.c

@@ -121 +121 @@
       goto skip;
     }
+    // we can hit inf if pmax and pmin are close to the max range
+    if (isinf(dx)) dx = DBL_MAX / Nbin;
 
     ALLOCATE (Nval, int, Nbin + 2);