Changeset 42682

Timestamp:

Jun 1, 2024, 3:17:30 PM (2 years ago)

Author:

eugene

Message:

linear extrapolation at endpoints

File:

• branches/eam_branches/ipp-20230313/Ohana/src/opihi/lib.data/spline.c (modified) (2 diffs)

branches/eam_branches/ipp-20230313/Ohana/src/opihi/lib.data/spline.c

@@ -113 +113 @@
   opihi_flt dx, a, b, value;
   
-  // saturate correction at high and low ends
-  if (X < x[0]) return y[0];
-  if (X > x[N-1]) return y[N-1];
+  // 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];
+  }
 
   /* find correct element in array (x must be sorted) */
@@ -128 +140 @@
     }
   }
+
+evaluate:
 
   /* error condition: duplicate abssisca */