Changeset 42332 for trunk/Ohana/src/opihi/cmd.data

Timestamp:

Jan 30, 2023, 9:45:29 AM (3 years ago)

Author:

eugene

Message:

adding option to spline1d to specify boundary 1st derivatives (default is 2nd derive = 0.0)

Location:

trunk/Ohana/src/opihi/cmd.data

Files:

• imspline_apply.c (modified) (1 diff)
• imspline_construct.c (modified) (1 diff)
• spline_commands.c (modified) (2 diffs)

trunk/Ohana/src/opihi/cmd.data/imspline_apply.c

@@ -77 +77 @@
       Tx1[i] = spline_apply_flt (Tyc, Ty1, Ty2, Ny, y);
     }
-    spline_construct_flt (Txc, Tx1, Nx, Tx2);
+    spline_construct_flt (Txc, Tx1, Nx, Tx2, NAN, NAN);
 
     /* apply x-dir spline to new image */

trunk/Ohana/src/opihi/cmd.data/imspline_construct.c

@@ -55 +55 @@
     }
   
-    spline_construct_flt (Tx, Ty, Ny, Ty2);
+    spline_construct_flt (Tx, Ty, Ny, Ty2, NAN, NAN);
   
     /* copy derivatives to output buffer */

trunk/Ohana/src/opihi/cmd.data/spline_commands.c

@@ -15 +15 @@
 int spline_create (int argc, char **argv) {
 
-  int i;
+  int i, N;
   Vector *xvec, *yvec;
+
+  double dyLower = NAN;
+  double dyUpper = NAN;
+  if ((N = get_argument (argc, argv, "-dyLower"))) {
+    remove_argument (N, &argc, argv);
+    dyLower = atof(argv[N]);
+    remove_argument (N, &argc, argv);
+  }
+  if ((N = get_argument (argc, argv, "-dyUpper"))) {
+    remove_argument (N, &argc, argv);
+    dyUpper = atof(argv[N]);
+    remove_argument (N, &argc, argv);
+  }
 
   if (argc != 4) {
@@ -43 +56 @@
   }    
 
-  spline_construct_dbl (myspline->xk, myspline->yk, myspline->Nknots, myspline->y2);
+  spline_construct_dbl (myspline->xk, myspline->yk, myspline->Nknots, myspline->y2, dyLower, dyUpper);
   return TRUE;
 }