/* test_mpf_nuderivd.c  just simple test, needs full */
/*                possible 2D, 3D and 4D for scattered points */

#include "mpf_nuderivd.h" 
#include <stdlib.h>
#include <stdio.h>

#define npoints 20
#define digits 100
int main(int argc, char *argv[])
{
  int i, norder, point;
  double xg[npoints];
  double cx[npoints];

  printf("test_mpf_nuderivd.c running\n");

  for(i=0; i<npoints; i++)
    xg[i] = (double)i/(double)(npoints-1);

  for(i=0; i<npoints; i++)
  {
    printf("xg[%2d]=%g \n", i, xg[i]);
  }

  for(norder=10; norder<npoints-1; norder++)
  {
    for(point=0; point<npoints; point++)
    {
      printf("norder=%2d, npoints=%2d, point=%2d \n",
                 norder, npoints, point);
      mpf_nuderivd(digits, norder, npoints, point, xg, cx);
      for(i=0; i<npoints; i++)
      {
        printf("cx[%2d]=%g \n", i, cx[i]);
      }
    }
  }

  gmp_printf("test_mpf_nuderivd finished \n");
  return 0;
} /* end test_mpf_nuderivd.c */