Changeset 20859 for branches/eam_branch_20081124/Ohana/src/opihi/cmd.data/integrate.c

Timestamp:

Nov 30, 2008, 11:13:15 AM (18 years ago)

Author:

eugene

Message:

allow functions to use FLT or INT vectors as appropriate; document valid types

File:

• branches/eam_branch_20081124/Ohana/src/opihi/cmd.data/integrate.c (modified) (2 diffs)

branches/eam_branch_20081124/Ohana/src/opihi/cmd.data/integrate.c

@@ -3 +3 @@
 int integrate (int argc, char **argv) {
   
-  int i, N, VERBOSE;
+  int i, N, VERBOSE, isFloat;
   opihi_flt *X, *Y;
-  double start, end, value, range;
+  double start, end, value, range, dvalue;
   Vector *vecx, *vecy;
 
@@ -25 +25 @@
   end   = atof (argv[4]);
 
-  X = vecx[0].elements.Flt;
-  Y = vecy[0].elements.Flt;
+  isFloat = vecy[0].type == OPIHI_FLT;
+  Yf = vecy[0].elements.Flt;
+  Yi = vecy[0].elements.Int;
 
   value = 0;
   range = 0;
-  for (i = 0; i < vecx[0].Nelements-1; i++, X++, Y++) {
-    if ((*X >= start) && (*X <= end)) {
-      value += *Y * (X[1] - X[0]);
-      range += (X[1] - X[0]);
-    }
-  }      
+
+  if (vecx[0].type == OPIHI_FLT) {
+    opihi_flt *X = vecx[0].elements.Flt;
+    for (i = 0; i < vecx[0].Nelements-1; i++, X++, Yi++, Yf++) {
+      if ((*X >= start) && (*X <= end)) {
+        dvalue = (isFloat) ? *Yf : *Yi;
+        value += *Y * (X[1] - X[0]);
+        range += (X[1] - X[0]);
+      }
+    }      
+  } else {
+    opihi_int *X = vecx[0].elements.Int;
+    for (i = 0; i < vecx[0].Nelements-1; i++, X++, Yi++, Yf++) {
+      if ((*X >= start) && (*X <= end)) {
+        dvalue = (isFloat) ? *Yf : *Yi;
+        value += *Y * (X[1] - X[0]);
+        range += (X[1] - X[0]);
+      }
+    }      
+  }
 
   set_variable ("sum", value);