// time_simeq_lup.java 

public class time_simeq_lup
{

  int N = 4096;
  double A[] = new double[N*N];
  double X[] = new double[N];
  double X_soln[] = new double[N];
  double Y[] = new double[N];

  void init_matrix(int n)
  {
    double sum;

    for(int i=0; i<n; i++)
    {
      X_soln[i] = Math.random();
      for(int j=0; j<n; j++)
        A[i*n+j] = Math.random();
    }
    for(int i=0; i<n; i++)
    {
      sum = 0.0;
      for(int j=0; j<n; j++)
        sum = sum + A[i*n+j]*X_soln[j];
      Y[i] = sum;
    }
    
    if(n<5)
    {
      for(int i=0; i<n; i++)
      {
        for(int j=0; j<n; j++)
	  System.out.println("A["+i+"]["+j+"]="+A[i*n+j]);
        System.out.println("X_soln["+i+"]="+X_soln[i]+", Y["+i+"]="+Y[i]);
        System.out.println(" ");
      }
    }
  } // end init_matrix

  time_simeq_lup()
  {
    double sum, sumerr, maxerr;
    double start, now, next;
  
    System.out.println("time_simeq_lup.java running");
    start  = System.currentTimeMillis()/1000.0;
    for(int n=4; n<N; n=n*2)
    {
      init_matrix(n);
      now  = System.currentTimeMillis()/1000.0;
      System.out.println("matrix initialized in "+(now-start)+" seconds");
      new simeq_lup(n, A, Y, X);
      next  = System.currentTimeMillis()/1000.0;
      sumerr = 0.0;
      maxerr = 0.0;
      for(int i=0; i<n; i++)
      {
        sum = 0.0;
        for(int j=0; j<n; j++) sum = sum + A[i*n+j]*X[j];
        sumerr = sumerr + Math.abs(sum-Y[i]);
        maxerr = Math.max(maxerr, Math.abs(X_soln[i]-X[i]));
        if(n<5) System.out.println("X["+i+"]="+X[i]+", err="+
                                   (X_soln[i]-X[i]));
      }
      System.out.println("n="+n+", sumerr="+sumerr+", maxerr="+maxerr);
      System.out.println("simeq_lup took "+(next-now)+" seconds");
    } // n
  } // end time_simeq_lup

  static public void main(String[] args)
  {   
    new time_simeq_lup();
  } // end main
} // end main time_simeq_lup.java