Changeset 34490 for branches/eam_branches/ipp-20120905/Ohana/src/opihi/lib.data/PlotVectors.c

Timestamp:

Sep 30, 2012, 8:03:14 AM (14 years ago)

Author:

eugene

Message:

add "where" option to plot, zplot, zcplot

File:

• branches/eam_branches/ipp-20120905/Ohana/src/opihi/lib.data/PlotVectors.c (modified) (4 diffs)

branches/eam_branches/ipp-20120905/Ohana/src/opihi/lib.data/PlotVectors.c

@@ -1 +1 @@
 # include "display.h"
 
-int PlotVectorSingle (int kapa, Vector *vec, char *mode) {
-
-  int i, Npts;
+int PlotVectorSingle (int kapa, Vector *vec, char *mask, char *mode) {
+
+  int i, Npts, Nout;
   float *temp;
 
@@ -9 +9 @@
   ALLOCATE (temp, float, Npts);
 
+  Nout = 0;
   if (vec->type == OPIHI_FLT) {
     opihi_flt *value = vec->elements.Flt;
     for (i = 0; i < Npts; i++) {
-      temp[i] = value[i];
+      if (mask && mask[i]) continue;
+      temp[Nout] = value[i];
+      Nout ++;
     }
   } else {
     opihi_int *value = vec->elements.Int;
     for (i = 0; i < Npts; i++) {
-      temp[i] = value[i];
-    }
-  }
-  KapaPlotVector (kapa, Npts, temp, mode);
-
-  free (temp);
-
-  return (TRUE);
-}
-
-int PlotVectorPair (int kapa, Vector *xVec, Vector *yVec, Graphdata *graphmode) {
-
-  int i, Npts;
+      if (mask && mask[i]) continue;
+      temp[Nout] = value[i];
+      Nout ++;
+    }
+  }
+  KapaPlotVector (kapa, Nout, temp, mode);
+
+  free (temp);
+
+  return (TRUE);
+}
+
+int PlotVectorPair (int kapa, Vector *xVec, Vector *yVec, char *mask, Graphdata *graphmode) {
+
+  int i, Npts, Nout;
   float *temp;
 
@@ -37 +42 @@
   ALLOCATE (temp, float, Npts);
 
-  KapaPrepPlot (kapa, Npts, graphmode);
-
+  Nout = 0;
   if (xVec->type == OPIHI_FLT) {
     opihi_flt *value = xVec->elements.Flt;
     for (i = 0; i < Npts; i++) {
-      temp[i] = value[i];
+      if (mask && mask[i]) continue;
+      temp[Nout] = value[i];
+      Nout ++;
     }
   } else {
     opihi_int *value = xVec->elements.Int;
     for (i = 0; i < Npts; i++) {
-      temp[i] = value[i];
-    }
-  }
-  KapaPlotVector (kapa, Npts, temp, "x");
-
+      if (mask && mask[i]) continue;
+      temp[Nout] = value[i];
+      Nout ++;
+    }
+  }
+  KapaPrepPlot (kapa, Nout, graphmode);
+  KapaPlotVector (kapa, Nout, temp, "x");
+
+  Nout = 0;
   if (yVec->type == OPIHI_FLT) {
     opihi_flt *value = yVec->elements.Flt;
     for (i = 0; i < Npts; i++) {
-      temp[i] = value[i];
+      if (mask && mask[i]) continue;
+      temp[Nout] = value[i];
+      Nout ++;
     }
   } else {
     opihi_int *value = yVec->elements.Int;
     for (i = 0; i < Npts; i++) {
-      temp[i] = value[i];
-    }
-  }
-  KapaPlotVector (kapa, Npts, temp, "y");
-
-  free (temp);
-
-  return (TRUE);
-}
-
-int PlotVectorTriplet (int kapa, Vector *xVec, Vector *yVec, Vector *zVec, Graphdata *graphmode) {
-
-  int i, Npts;
+      if (mask && mask[i]) continue;
+      temp[Nout] = value[i];
+      Nout ++;
+    }
+  }
+  KapaPlotVector (kapa, Nout, temp, "y");
+
+  free (temp);
+
+  return (TRUE);
+}
+
+int PlotVectorTriplet (int kapa, Vector *xVec, Vector *yVec, Vector *zVec, char *mask, Graphdata *graphmode) {
+
+  int i, Npts, Nout;
   float *temp;
 
@@ -81 +95 @@
   ALLOCATE (temp, float, Npts);
 
-  KapaPrepPlot (kapa, Npts, graphmode);
-
+  Nout = 0;
   if (xVec->type == OPIHI_FLT) {
     opihi_flt *value = xVec->elements.Flt;
     for (i = 0; i < Npts; i++) {
-      temp[i] = value[i];
+      if (mask && mask[i]) continue;
+      temp[Nout] = value[i];
+      Nout ++;
     }
   } else {
     opihi_int *value = xVec->elements.Int;
     for (i = 0; i < Npts; i++) {
-      temp[i] = value[i];
-    }
-  }
-  KapaPlotVector (kapa, Npts, temp, "x");
-
+      if (mask && mask[i]) continue;
+      temp[Nout] = value[i];
+      Nout ++;
+    }
+  }
+  KapaPrepPlot (kapa, Nout, graphmode);
+  KapaPlotVector (kapa, Nout, temp, "x");
+
+  Nout = 0;
   if (yVec->type == OPIHI_FLT) {
     opihi_flt *value = yVec->elements.Flt;
     for (i = 0; i < Npts; i++) {
-      temp[i] = value[i];
+      if (mask && mask[i]) continue;
+      temp[Nout] = value[i];
+      Nout ++;
     }
   } else {
     opihi_int *value = yVec->elements.Int;
     for (i = 0; i < Npts; i++) {
-      temp[i] = value[i];
-    }
-  }
-  KapaPlotVector (kapa, Npts, temp, "y");
-
+      if (mask && mask[i]) continue;
+      temp[Nout] = value[i];
+      Nout ++;
+    }
+  }
+  KapaPlotVector (kapa, Nout, temp, "y");
+
+  Nout = 0;
   if (zVec->type == OPIHI_FLT) {
     opihi_flt *value = zVec->elements.Flt;
     for (i = 0; i < Npts; i++) {
-      temp[i] = value[i];
+      if (mask && mask[i]) continue;
+      temp[Nout] = value[i];
+      Nout ++;
     }
   } else {
     opihi_int *value = zVec->elements.Int;
     for (i = 0; i < Npts; i++) {
-      temp[i] = value[i];
-    }
-  }
-  KapaPlotVector (kapa, Npts, temp, "z");
+      if (mask && mask[i]) continue;
+      temp[Nout] = value[i];
+      Nout ++;
+    }
+  }
+  KapaPlotVector (kapa, Nout, temp, "z");
 
   free (temp);