corrects interpretation of sort order

nasa · Jul 27, 2022 · 820cc1a · 820cc1a
1 parent fdea519
commit 820cc1a
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/ref_inflate.c b/src/ref_inflate.c
@@ -433,19 +433,19 @@ REF_FCN static REF_STATUS ref_inflate_sort_rail(REF_INT n, REF_DBL *x,
   ref_malloc(tmp, n, REF_DBL);
   RSS(ref_sort_heap_dbl(n, x, order), "heap");
   for (node = 0; node < n; node++) {
-    tmp[order[node]] = x[node];
+    tmp[node] = x[order[node]];
   }
   for (node = 0; node < n; node++) {
     x[node] = tmp[node];
   }
   for (node = 0; node < n; node++) {
-    tmp[order[node]] = yz[0 + 2 * node];
+    tmp[node] = yz[0 + 2 * order[node]];
   }
   for (node = 0; node < n; node++) {
     yz[0 + 2 * node] = tmp[node];
   }
   for (node = 0; node < n; node++) {
-    tmp[order[node]] = yz[1 + 2 * node];
+    tmp[node] = yz[1 + 2 * order[node]];
   }
   for (node = 0; node < n; node++) {
     yz[1 + 2 * node] = tmp[node];