Changeset 15801

Timestamp:

Dec 12, 2007, 6:06:59 PM (18 years ago)

Author:

eugene

Message:

added sort-by-radius

File:

• trunk/psphot/src/psphotRadialProfile.c (modified) (2 diffs)

trunk/psphot/src/psphotRadialProfile.c

@@ -1 +1 @@
 # include "psphotInternal.h"
 
+# define COMPARE(A,B) (radius->data.F32[A] < radius->data.F32[B])
+# define SWAPFUNC(TYPE,A,B) { \
+  float tmp; \
+  if (A != B) { \
+    tmp = radius->data.F32[A]; \
+    radius->data.F32[A] = radius->data.F32[B]; \
+    radius->data.F32[B] = tmp; \
+    tmp = flux->data.F32[A]; \
+    flux->data.F32[A] = flux->data.F32[B]; \
+    flux->data.F32[B] = tmp; \
+    tmp = weight->data.F32[A]; \
+    weight->data.F32[A] = weight->data.F32[B]; \
+    weight->data.F32[B] = tmp; \
+  }
+    
 bool psphotRadialProfile (pmSource *source, psMetadata *recipe, psMaskType maskVal) {
 
@@ -40 +55 @@
     flux->n = n;
 
-    // XXX need to sort here
-    // SortVectorsByRadius (radius, flux, weight);
+    // sort the vector set by the radius
+    PSSORT (radius->n, COMPARE, SWAPFUNC, NONE);
 
     return true;