htswanalysis/MACS/lib/gsl/gsl-1.11/doc/examples/interpp.c

   1 #include <stdlib.h>
   2 #include <stdio.h>
   3 #include <math.h>
   4 #include <gsl/gsl_errno.h>
   5 #include <gsl/gsl_spline.h>
   6
   7 int
   8 main (void)
   9 {
  10   int N = 4;
  11   double x[4] = {0.00, 0.10,  0.27,  0.30};
  12   double y[4] = {0.15, 0.70, -0.10,  0.15}; /* Note: first = last
  13                                                for periodic data */
  14
  15   gsl_interp_accel *acc = gsl_interp_accel_alloc ();
  16   const gsl_interp_type *t = gsl_interp_cspline_periodic;
  17   gsl_spline *spline = gsl_spline_alloc (t, N);
  18
  19   int i; double xi, yi;
  20
  21   printf ("#m=0,S=5\n");
  22   for (i = 0; i < N; i++)
  23     {
  24       printf ("%g %g\n", x[i], y[i]);
  25     }
  26
  27   printf ("#m=1,S=0\n");
  28   gsl_spline_init (spline, x, y, N);
  29
  30   for (i = 0; i <= 100; i++)
  31     {
  32       xi = (1 - i / 100.0) * x[0] + (i / 100.0) * x[N-1];
  33       yi = gsl_spline_eval (spline, xi, acc);
  34       printf ("%g %g\n", xi, yi);
  35     }
  36
  37   gsl_spline_free (spline);
  38   gsl_interp_accel_free (acc);
  39   return 0;
  40 }