// test_nav1.java    

import java.io.*;

public class test_nav1
{
  nav1 N = new nav1();

  public test_nav1() // constructor
  {
    System.out.println("test_nav1 running");
    N.test();
    write_bound();
    System.out.println("test_nav1 finished");
  }

  void write_bound()
  {

    try
    {
      FileWriter fout = new FileWriter("nav1_boundary.dat");
      BufferedWriter fileout = new BufferedWriter(fout);

      System.out.println("writing nav1_boundary.dat");

    double xmin = 0.0;
    double xmax = 1.0;
    double ymin = 0.0;
    double ymax = 1.0;
    double zmin = 0.0;
    double zmax = 1.0;
    double tmin = 0.0;
    double tmax = 1.0;
    double x, y, z, t;
    double hx, hy, hz, ht;
    int nx = 2; // just boundary
    int ny = 2;
    int nz = 2;
    int nt = 2;

    hx = (xmax-xmin)/(nx-1);
    hy = (ymax-ymin)/(ny-1);
    hz = (zmax-zmin)/(nz-1);
    ht = (tmax-tmin)/(nt-1);
    for(int i=0; i<nx; i++)
    {
      x = xmin+hx*i;
      for(int j=0; j<nx; j++)
      {
        y = ymin+hy*j;
        for(int k=0; k<nx; k++)
        {
          z = zmin+hz*k;
          for(int m=0; m<nx; m++)
          {
            t = tmin+ht*m;
            fileout.write(" "+N.u(x,y,z,t)+" "+x+" "+y+" "+z+" "+t+" \n");
	  } // end t
	} // end z
      } // end y
    } // end x
    fileout.close();
    System.out.println("finished writing");
    }
    catch(FileNotFoundException exception)
    {
    System.out.println("nav1_boundary.dat  not found");
    }
    catch(Exception exception) // catches all exceptions
    {
      System.out.println("test_nav1.java output-file error");
    } // end try
  } // end write_bound

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