// time_simeq_thread.java  not a great benchmark

import java.text.*;

class time_simeq_thread
{
  int N; // can use command line
  int NP;
  DecimalFormat f1 = new DecimalFormat("0.000E00");

  time_simeq_thread(int N, int NP)
  {
    this.N = N;
    this.NP = NP;
    double A[] = new double[N*N];
    double X[] = new double[N];
    double Y[] = new double[N];
    double t, total;
    double start, now, next;
  
    System.out.println("time_simeq_thread.java running, N="+N+", NP="+NP);
    start = System.currentTimeMillis();

    init_matrix(N, A, Y);
    now = System.currentTimeMillis();
    System.out.println("matrix initialized in "+((now-start)/1000)+" sec");
    simeq_thread ST = new simeq_thread(NP);
    ST.psimeq(NP, N, A, Y, X); // new
    // new simeq(N, A, Y, X); // old
    next = System.currentTimeMillis();
    total = 0.0;
    for(int i=0; i<N; i++)
    {
      t = 0.0;
      for(int j=0; j<N; j++) t = t + A[i*N+j]*X[j];
        total = total + Math.abs(t-Y[i]);
    }
    System.out.println("N="+N+", NP="+NP+", error="+total+
                       ", "+((next-now)/1000)+" seconds");
  } // end time_simeq_thread

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

  public static void main (String[] args)
  {
    int n = 200; // default, can use command line
    int np = 4;
    if(args.length>1) np = Integer.parseInt(args[1]);
    if(args.length>0) n = Integer.parseInt(args[0]);
    new time_simeq_thread(n, np);
  } // end main

} // end class time_simeq_thread.java