htswanalysis/MACS/lib/gsl/gsl-1.11/doc/examples/interp.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 i;
  11   double xi, yi, x[10], y[10];
  12
  13   printf ("#m=0,S=2\n");
  14
  15   for (i = 0; i < 10; i++)
  16     {
  17       x[i] = i + 0.5 * sin (i);
  18       y[i] = i + cos (i * i);
  19       printf ("%g %g\n", x[i], y[i]);
  20     }
  21
  22   printf ("#m=1,S=0\n");
  23
  24   {
  25     gsl_interp_accel *acc
  26       = gsl_interp_accel_alloc ();
  27     gsl_spline *spline
  28       = gsl_spline_alloc (gsl_interp_cspline, 10);
  29
  30     gsl_spline_init (spline, x, y, 10);
  31
  32     for (xi = x[0]; xi < x[9]; xi += 0.01)
  33       {
  34         yi = gsl_spline_eval (spline, xi, acc);
  35         printf ("%g %g\n", xi, yi);
  36       }
  37     gsl_spline_free (spline);
  38     gsl_interp_accel_free (acc);
  39   }
  40   return 0;
  41 }