#include <iostream>
#include<time.h>
#include<omp.h>
#include<sys/time.h>
using namespace std;


void lu_factor(double **A, int n){

  int i,j,k;
  i=j=k=0;
  for(k=0;k<n-1;k++){
    for(i=k+1;i<n;i++)
	A[i][k] = A[i][k]/A[k][k];
    for(i=k+1;i<n;i++)
       for(j=k+1;j<n;j++)
	   A[i][j] = A[i][j]-A[i][k]*A[k][j];
  }

}

void lu(double **A,int n){

  int  i, j, count, count_rate, count_max,T1,T2;

  char buffer[30];
  struct timeval tv,hh;
  time_t curtime;
  long double StartTime, EndTime;
  suseconds_t st,end;
  gettimeofday(&tv, NULL); 
  	lu_factor(A,n);
  gettimeofday(&hh, NULL); 
  
  cout<< "\n THE TOTAL TIME IN LU: "<<hh.tv_sec-tv.tv_sec<<"sec : "<<hh.tv_usec-tv.tv_usec<<"usec.\n";

}


int main(int argc, char *argv[]) {
  int n,i,j,p;
  double **A;

  double processors;

	n = atoi( argv[1] );
	

# ifdef _OPENMP
	if(argv[2]!=NULL)
		processors = atoi( argv[2] );	
	else 
		{cout<<"Error: Number of Processors Unspecified. Setting Default = 1!!"; processors = 1;}
	omp_set_num_threads(processors);
	cout<<"\n THE NUMBER OF PROCESSORS: "<<processors;
# endif

  cout<<"\n ARRAY SIZE: "<<n;

  A = new double *[n];
  for(i=0;i<n;i++)
      *(A+i) = new double[n];

  for(i=0;i<n;i++)
     for(j=0;j<n;j++)
	 A[i][j] = rand();
     
	lu(A,n);

}