#include <stdlib.h>
#include <stdio.h>

void InitTim();
void PrtTim(char *);

void Mult(int n, double *v1, double *v2, double *v3);
void MultOMP(int n, double *v1, double *v2, double *v3);

int main (int narg, char *arg[])
{
int i,j, N, M;

double *v1, *v2, *v3;

if (narg < 2) { 
  printf("\n Usage tompC P/x [N M] \n \n");
  printf(" P --> Calling MultOMP, x ---> Calling Mult \n");
  printf(" N (=100): External loop number   M : Array Size (500000) \n");
  exit(0);
}
N = 100;
M = 500000;
if (narg > 2) N = atoi(arg[2]);
if (narg > 3) M = atoi(arg[3]);

InitTim();

v1 = (double *)malloc(M*sizeof(double));
v2 = (double *)malloc(M*sizeof(double));
v3 = (double *)malloc(M*sizeof(double));
 for(i=0; i<M; i++) { v1[i] = random()%1000;  v2[i] = random()%5000; }
PrtTim("End of Init ");

if (*arg[1] == 'P') { 
  printf("Calling N=%d times MultOMP(Size= %d) \n", N, M);
  for(i=0; i<N; i++) MultOMP(M, v1, v2, v3);
}
else { 
  printf("Calling N=%d times Mult(Size= %d) \n", N, M);
  for(i=0; i<N; i++) Mult(M, v1, v2, v3);
}

PrtTim("End of Multiplication ");

free(v1);
free(v2);
free(v3);

exit(0);
}


void Mult(int n, double *v1, double *v2, double *v3)
{
int k;
for(k=0; k<n; k++)
  v3[k] = v1[k] * v2[k];
}

void MultOMP(int n, double *v1, double *v2, double *v3)
{
int k;
#pragma omp parallel for
for(k=0; k<n; k++)
  v3[k] = v1[k] * v2[k];
}