Changeset 33944 for branches/eam_branches/ipp-20120405/Ohana/src/opihi/cmd.data/threshold.c

Timestamp:

May 29, 2012, 4:18:58 PM (14 years ago)

Author:

eugene

Message:

split out time from csystem inits in dbExtractAverageInit to dbExtractAverageInitAve; allow option to skip result output for parallel ops; fix remote -reload; flush, fsync, fsetlockfile on remote command to ensure all machines see it (and report errors caught); handle equal value bins in threshold function; set global variable to hold the scrip name

File:

• branches/eam_branches/ipp-20120405/Ohana/src/opihi/cmd.data/threshold.c (modified) (2 diffs)

branches/eam_branches/ipp-20120405/Ohana/src/opihi/cmd.data/threshold.c

@@ -69 +69 @@
     x0 = isFltX ? vecx[0].elements.Flt[Nhi] : vecy[0].elements.Int[Nhi];
     x1 = isFltX ? vecx[0].elements.Flt[Nhi+1]   : vecy[0].elements.Int[Nhi+1];
-    Xvalue = (value - y0) * (x1 - x0) / (y1 - y0) + x0;
+    if (y0 == y1) {
+      Xvalue = 0.5*(x0 + x1);
+    } else {
+      Xvalue = (value - y0) * (x1 - x0) / (y1 - y0) + x0;
+    }    
   } else {
     // interpolate to value:
@@ -76 +80 @@
     x0 = isFltX ? vecx[0].elements.Flt[Nhi-1] : vecy[0].elements.Int[Nhi-1];
     x1 = isFltX ? vecx[0].elements.Flt[Nhi]   : vecy[0].elements.Int[Nhi];
-    Xvalue = (value - y0) * (x1 - x0) / (y1 - y0) + x0;
+    if (y0 == y1) {
+      Xvalue = 0.5*(x0 + x1);
+    } else {
+      Xvalue = (value - y0) * (x1 - x0) / (y1 - y0) + x0;
+    }
   }