/* nag_1d_spline_deriv (e02bcc) Example Program.
 *
 * NAGPRODCODE Version.
 *
 * Copyright 2016 Numerical Algorithms Group.
 *
 * Mark 26, 2016.
 *
 */

#include <nag.h>
#include <stdio.h>
#include <nag_stdlib.h>
#include <nage02.h>

int main(void)
{
  Integer exit_status = 0, i, j, l, m, ncap, ncap7;
  NagError fail;
  Nag_DerivType derivs;
  Nag_Spline spline;
  double s[4], x;

  INIT_FAIL(fail);

  /* Initialize spline */
  spline.lamda = 0;
  spline.c = 0;

  printf("nag_1d_spline_deriv (e02bcc) Example Program Results\n");
  scanf("%*[^\n]"); /* Skip heading in data file */
  while (scanf("%" NAG_IFMT "%" NAG_IFMT "", &ncap, &m) != EOF)
  {
    if (m <= 0) {
      printf("Invalid m.\n");
      exit_status = 1;
      return exit_status;
    }
    if (ncap > 0) {
      ncap7 = ncap + 7;
      spline.n = ncap7;
      if (!(spline.c = NAG_ALLOC(ncap7, double)) ||
          !(spline.lamda = NAG_ALLOC(ncap7, double)))
      {
        printf("Allocation failure\n");
        exit_status = -1;
        goto END;
      }
    }
    else {
      printf("Invalid ncap.\n");
      exit_status = 1;
      return exit_status;
    }
    for (j = 0; j < ncap7; j++)
      scanf("%lf", &(spline.lamda[j]));
    for (j = 0; j < ncap + 3; j++)
      scanf("%lf", &(spline.c[j]));
    printf("        x             Spline    1st deriv  "
           "2nd deriv   3rd deriv");
    for (i = 1; i <= m; i++) {
      scanf("%lf", &x);
      derivs = Nag_LeftDerivs;
      for (j = 1; j <= 2; j++) {
        /* nag_1d_spline_deriv (e02bcc).
         * Evaluation of fitted cubic spline, function and
         * derivatives
         */
        nag_1d_spline_deriv(derivs, x, s, &spline, &fail);
        if (fail.code != NE_NOERROR) {
          printf("Error from nag_1d_spline_deriv (e02bcc).\n%s\n",
                 fail.message);
          exit_status = 1;
          goto END;
        }

        if (derivs == Nag_LeftDerivs) {
          printf("\n\n%11.4f   Left", x);
          for (l = 0; l < 4; l++)
            printf("%11.4f", s[l]);
        }
        else {
          printf("\n%11.4f  Right", x);
          for (l = 0; l < 4; l++)
            printf("%11.4f", s[l]);
        }
        derivs = Nag_RightDerivs;
      }
    }
    printf("\n");
  END:
    NAG_FREE(spline.c);
    NAG_FREE(spline.lamda);
  }
  return exit_status;
}