Changeset 42482 for trunk/Ohana/src/opihi/cmd.data

Timestamp:

Aug 2, 2023, 3:05:46 PM (3 years ago)

Author:

eugene

Message:

uniq: handle (ignore) nan values in input vector

File:

• trunk/Ohana/src/opihi/cmd.data/uniq.c (modified) (2 diffs)

trunk/Ohana/src/opihi/cmd.data/uniq.c

@@ -45 +45 @@
     opihi_flt *indata = NULL;
     ALLOCATE (indata, opihi_flt, ivec[0].Nelements);
-    memcpy (indata, ivec->elements.Flt, ivec[0].Nelements*sizeof(opihi_flt));
 
-    dsort (indata, ivec->Nelements);
+    // manually copy and exclude any nan values
+    int Nindata = 0;
+    for (i = 0; i < ivec->Nelements; i++) {
+      if (!isfinite(ivec[0].elements.Flt[i])) continue;
+      indata[Nindata] = ivec[0].elements.Flt[i];
+      Nindata ++;
+    }
+    // XXX do not do this: memcpy (indata, ivec->elements.Flt, ivec[0].Nelements*sizeof(opihi_flt));
+
+    dsort (indata, Nindata);
 
     Nnew = 0;
@@ -54 +62 @@
     opihi_flt *vsrc = indata;
 
-    int onePercent = ivec->Nelements / 100;
+    int onePercent = Nindata / 100;
 
     struct sigaction *old_sigaction = SetInterrupt();
-    for (i = 0; (i < ivec->Nelements) && !interrupt; Nnew++) {
+    for (i = 0; (i < Nindata) && !interrupt; Nnew++) {
       vtgt[Nnew] = *vsrc;
       int Ndup = 0;