Changeset 42967 for trunk/Ohana/src/opihi/lib.data/spline.c

Timestamp:

Mar 9, 2026, 11:26:15 AM (4 months ago)

Author:

eugene

Message:

add spline options to extrapolate or saturate at end points; add vweave function to interpolate a vector with varying bin sizes

File:

• trunk/Ohana/src/opihi/lib.data/spline.c (modified) (4 diffs)

trunk/Ohana/src/opihi/lib.data/spline.c

@@ -1 +1 @@
 # include "data.h"
+
+static SplineMode MODE = SPLINE_SATURATE; // SATURATE or EXTRAPOLATE?
+
+int spline_contruct_mode (SplineMode mode) {
+  
+  switch (mode) {
+  case SPLINE_SATURATE:
+  case SPLINE_EXTRAPOLATE:
+    MODE = mode;
+    break;
+  default:
+    return FALSE;
+  }
+  return TRUE;
+}
 
 /* construct the natural spline for x, y in y2 */
@@ -36 +51 @@
   
   // saturate correction at high and low ends
-  if (X < x[0]) return y[0];
-  if (X > x[N-1]) return y[N-1];
+  if (MODE == SPLINE_SATURATE) {
+    if (X < x[0]) return y[0];
+    if (X > x[N-1]) return y[N-1];
+  }
 
   /* find correct element in array (x must be sorted) */
@@ -113 +130 @@
   opihi_flt dx, a, b, value;
   
-  // linear extrapolation past endpoints
-  if (X < x[0]) {
-    hi = 1;
-    lo = 0;
-    goto evaluate;
-    // alternative: saturate correction at high and low ends
-    // return y[0]; 
-  }
-  if (X > x[N-1]) {
-    hi = N - 1;
-    lo = N - 2;
-    goto evaluate;
-    // alternative: saturate correction at high and low ends
-    // return y[N-1];
+  // saturate correction at high and low ends
+  if (MODE == SPLINE_SATURATE) {
+    if (X < x[0]) return y[0];
+    if (X > x[N-1]) return y[N-1];
   }
 
@@ -140 +147 @@
     }
   }
-
-evaluate:
 
   /* error condition: duplicate abssisca */