Changeset 40751

Timestamp:

May 25, 2019, 8:56:56 AM (7 years ago)

Author:

eugene

Message:

add dsort_indexonly

Location:

branches/eam_branches/ohana.20190329/src/libohana

Files:

• include/ohana_sort.h (modified) (1 diff)
• src/sorts.c (modified) (1 diff)

branches/eam_branches/ohana.20190329/src/libohana/include/ohana_sort.h

@@ -83 +83 @@
 void isortfour (int *X, int *Y, int *Z, int *W, int N);
 
+void dsort_indexonly (double *X, off_t *S, off_t N);
+
 #endif

branches/eam_branches/ohana.20190329/src/libohana/src/sorts.c

@@ -196 +196 @@
 }
 
+/* sort the index of a vector (vector stays unsorted) */
+void dsort_indexonly (double *X, off_t *S, off_t N) {
+  
+# define SWAPFUNC(A,B){ off_t itmp; \
+  itmp = S[A]; S[A] = S[B]; S[B] = itmp; \
+}
+# define COMPARE(A,B)((!isfinite(X[S[A]]) && isfinite(X[S[B]])) || (X[S[A]] < X[S[B]]))
+
+  OHANA_SORT (N, COMPARE, SWAPFUNC);
+
+# undef SWAPFUNC
+# undef COMPARE
+
+}
+
+