Changeset 42946 for branches/eam_branches/ipp-ppmerge-20241229/Ohana/src/opihi/lib.data/spline.c

Timestamp:

Nov 23, 2025, 4:33:16 PM (8 months ago)

Author:

eugene

Message:

add spline option to extrapolate (vs saturate); add vweave function (like drizzle for a vector)

File:

• branches/eam_branches/ipp-ppmerge-20241229/Ohana/src/opihi/lib.data/spline.c (modified) (3 diffs)

branches/eam_branches/ipp-ppmerge-20241229/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) */
@@ -114 +131 @@
   
   // 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) */