// nuderiv_test.java  non uniformly spaced derivative coefficients test 

import java.text.*;

class nuderiv_test
{
  nuderiv_test()
  {
    int npoint;
    double x[] = {1.0, 2.0, 4.0, 7.0, 8.5};
    double y[] = new double[5];
    double yp[] = new double[5];
    double ypp[] = new double[5];
    double yppp[] = new double[5];
    double ypppp[] = new double[5];
    double c[] = new double[5];
    double nuyp; // any order 
    DecimalFormat f6= new DecimalFormat("0.0000");
    DecimalFormat fe= new DecimalFormat("0.00E00");

    System.out.println("nuderiv_test.java running ");
    npoint = 4;
    System.out.println(" testing  npoint="+npoint);
    for(int i=0; i<npoint; i++)
    {
      y[i] = yf4(x[i]);
      yp[i] = ypf4(x[i]);
      ypp[i] = yppf4(x[i]);
      yppp[i] = ypppf4(x[i]);
      System.out.println("x["+i+"]="+f6.format(x[i])+
                         ", y="+f6.format(y[i])+
                         ", y'="+f6.format(yp[i])+
                         ", y''="+f6.format(ypp[i])+
                         ", y'''="+f6.format(yppp[i]));
    }
    for(int order=1; order<npoint; order++)
    {
      for(int point=0; point<npoint; point++)
      {
        System.out.println("testing order="+order+", npoint="+npoint+
                           ", point="+point);
        new nuderiv(order, npoint, point, x, c);
        for(int i=0; i<npoint; i++)
	    System.out.println("nu c["+i+"]="+f6.format(c[i]));
        nuyp = 0.0;
        for(int i=0; i<npoint; i++)
        {
          nuyp = nuyp + c[i]*y[i];
        }
        if(order==1)
          System.out.println("order="+order+", npoint="+npoint+
                             ", point="+point+", yp="+f6.format(yp[point])+
                             ", nuyp="+f6.format(nuyp)+
	                     ", err="+fe.format(yp[point]-nuyp));
        else if(order==2)
          System.out.println("order="+order+", npoint="+npoint+
                             ", point="+point+", ypp="+f6.format(ypp[point])+
                             ", nuyp="+f6.format(nuyp)+
	                     ", err="+fe.format(ypp[point]-nuyp));
        else if(order==3)
          System.out.println("order="+order+", npoint="+npoint+
                             ", point="+point+", yppp="+f6.format(yppp[point])+
                             ", nuyp="+f6.format(nuyp)+
	                     ", err="+fe.format(yppp[point]-nuyp));
      } // end point loop
    } // end order loop

    npoint = 5;
    System.out.println(" ");
    System.out.println(" testing  npoint="+npoint);
    for(int i=0; i<npoint; i++)
    {
      y[i] = yf5(x[i]);
      yp[i] = ypf5(x[i]);
      ypp[i] = yppf5(x[i]);
      yppp[i] = ypppf5(x[i]);
      ypppp[i] = yppppf5(x[i]);
      System.out.println("x["+i+"]="+f6.format(x[i])+
                         ", y="+f6.format(y[i])+
                         ", y'="+f6.format(yp[i])+
                         ", y''="+f6.format(ypp[i])+
                         ", y'''="+f6.format(yppp[i])+
                         ", y''''="+f6.format(ypppp[i]));
    }
    for(int order=1; order<npoint; order++)
    {
      for(int point=0; point<npoint; point++)
      {
        System.out.println("testing order="+order+", npoint="+npoint+
                           ", point="+point);
        new nuderiv(order, npoint, point, x, c);
        for(int i=0; i<npoint; i++) 
	  System.out.println("nu c["+i+"]="+f6.format(c[i]));
        nuyp = 0.0;
        for(int i=0; i<npoint; i++)
        {
          nuyp = nuyp + c[i]*y[i];
        }
        if(order==1)
          System.out.println("order="+order+", npoint="+npoint+
                             ", point="+point+", yp="+f6.format(yp[point])+
                             ", nuyp="+f6.format(nuyp)+
	                     ", err="+fe.format(yp[point]-nuyp));
        else if(order==2)
          System.out.println("order="+order+", npoint="+npoint+
                             ", point="+point+", ypp="+f6.format(ypp[point])+
                             ", nuyp="+f6.format(nuyp)+
	                     ", err="+fe.format(ypp[point]-nuyp));
        else if(order==3)
          System.out.println("order="+order+", npoint="+npoint+
                             ", point="+point+", yppp="+f6.format(yppp[point])+
                             ", nuyp="+f6.format(nuyp)+
	                     ", err="+fe.format(yppp[point]-nuyp));
        else if(order==4)
          System.out.println("order="+order+", npoint="+npoint+
                             ", point="+point+", ypppp="+f6.format(ypppp[point])+
                             ", nuyp="+f6.format(nuyp)+
	                     ", err="+fe.format(ypppp[point]-nuyp));
      } // end point loop
    } // end order loop
  } // end nuderiv_test

  double yf4(double x) // test function 4 
  {
    return 3.0 -0.5* x + 0.2*x*x - 0.1*x*x*x;
  }

  double ypf4(double x) // derivative of test function 4
  {
    return -0.5 + 0.4*x - 0.3*x*x; 
  }

  double yppf4(double x) // second derivative of test function 4 
  {
    return 0.4 - 0.6*x; 
  }

  double ypppf4(double x) // third derivative of test function 4 
  {
    return -0.6; 
  }

  double yf5(double x) // test function 5 
  {
    return 3.0 -0.5* x + 0.1*x*x - 0.05*x*x*x + 0.002*x*x*x*x; 
  }

  double ypf5(double x) // derivative of test function 5
  {
    return -0.5 + 0.2*x - 0.15*x*x + 0.008*x*x*x; 
  }

  double yppf5(double x) // second derivative of test function 5 
  {
    return 0.2 - 0.3*x + 0.024*x*x; 
  }

  double ypppf5(double x) // third derivative of test function 5 
  {
    return -0.3 + 0.048*x; 
  }

  double yppppf5(double x) // fourth derivative of test function 5 
  {
    return 0.048; 
  }

  public static void main (String[] args)
  {
    new nuderiv_test();
  }
} // end class nuderiv_test