/* nav1_deriv.c  compute coefficients for PDE solution        */
/*   construct for each order, for each dimension, xc, shown  */
/*   xc[point][solution_subscript]                            */

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

void nav1_deriv(int order, int nx, double xmin, double xmax, double xc[])
{
  int debug1 = 1; /* basic derivative terms */
  double hx;
  int i, ii;

  printf("nav1_deriv running \n");
  hx = (xmax-xmin)/(nx-1);
  printf("nx=%d, xmin=%6.3f, xmax=%6.3f, hx=%8.5f \n",nx,xmin,xmax,hx);

  for(i=0; i<nx; i++) /* set up descrete derivative in x direction */
  {
    rderiv(order,nx,i,hx,&xc[i*nx]);
    if(debug1)
      for(ii=0; ii<nx; ii++)
        printf("order=%d, xc[%d][%d]=%15.9f\n",order,i,ii,xc[i*nx+ii]);
  }
  if(debug1) printf("\n");
} /* end nav1_deriv */
/* end  nav1_deriv.c */