// time_psimeq.java  multiple runs

import java.text.*;

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

  time_psimeq(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_psimeq.java running, N="+N+", NP="+NP);
    start = System.currentTimeMillis();

    init_matrix(N, A, Y);
    now = System.currentTimeMillis();
    System.out.println("matrix initialized at "+((now-start)/1000));
    new psimeq(NP, N, A, Y, X);
    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_psimeq

  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();
    }
  } // 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_psimeq(n, np);
  } // end main

} // end class time_psimeq.java