Changeset 37006 for branches/eam_branches/ipp-20140610/Ohana/src/relastro/src/UpdateObjects.c

Timestamp:

Jul 13, 2014, 11:06:04 AM (12 years ago)

Author:

eugene

Message:

refColor changed to refColorBlue,refColorRed

File:

• branches/eam_branches/ipp-20140610/Ohana/src/relastro/src/UpdateObjects.c (modified) (7 diffs)

branches/eam_branches/ipp-20140610/Ohana/src/relastro/src/UpdateObjects.c

@@ -10 +10 @@
 static double *T;
 static double *dT;
-static double *C;
+static double *C_blue;
+static double *C_red;
 
 void initObjectData (Catalog *catalog, int Ncatalog) {
@@ -38 +39 @@
   ALLOCATE (pY, double, MAX (1, Nmax));
 
-  ALLOCATE (C, double, MAX (1, Nmax));
+  ALLOCATE (C_blue, double, MAX (1, Nmax));
+  ALLOCATE (C_red,  double, MAX (1, Nmax));
 }  
 
@@ -58 +60 @@
   free (pY);
 
-  free (C);
+  free (C_blue);
+  free (C_red);
 }  
 
@@ -66 +69 @@
 
   off_t j, k, m;
-  int i, N, Nc, Nsecfilt, mode, result, status, XVERB;
+  int i, N, NcBlue, NcRed, Nsecfilt, mode, result, status, XVERB;
   StatType statsR, statsD;
   Coords coords;
@@ -130 +133 @@
       }
 
-      Nc = 0;
+      NcBlue = 0;
+      NcRed = 0;
       N = 0;
       m = catalog[i].average[j].measureOffset;
@@ -225 +229 @@
 
         if (setRefColor) {
-          float color = getColor (m+k, i);
-          if (!isnan(color)) {
-            C[Nc] = color;
-            Nc++;
+          float colorBlue = getColorBlue (m+k, i);
+          if (!isnan(colorBlue)) {
+            C_blue[NcBlue] = colorBlue;
+            NcBlue++;
+          }
+          float colorRed = getColorRed (m+k, i);
+          if (!isnan(colorRed)) {
+            C_red[NcRed] = colorRed;
+            NcRed++;
           }
         }
@@ -360 +369 @@
 
       if (setRefColor) {
-        dsort (C, Nc);
-        float colorMedian = (Nc > 0) ? C[(int)(0.5*Nc)] : NAN;
-        catalog[i].average[j].refColor = colorMedian;
+        float colorMedian;
+        dsort (C_blue, NcBlue);
+        colorMedian = (NcBlue > 0) ? C_blue[(int)(0.5*NcBlue)] : NAN;
+        catalog[i].average[j].refColorBlue = colorMedian;
+        dsort (C_red, NcRed);
+        colorMedian = (NcRed > 0) ? C_red[(int)(0.5*NcRed)] : NAN;
+        catalog[i].average[j].refColorRed = colorMedian;
       }