Changeset 1096 in Sophya

Timestamp:

Jul 26, 2000, 5:02:55 PM (25 years ago)

Author:

ansari

Message:

Extension test OpenMP (tompCXX.cc) - Reza 26/7/2000

File:

• trunk/Eval/Speed/tompCXX.cc (modified) (4 diffs)

trunk/Eval/Speed/tompCXX.cc

@@ -9 +9 @@
 void Mult(int n, double *v1, double *v2, double *v3);
 void MultOMP(int n, double *v1, double *v2, double *v3);
+void MultOMP2(int N, int n, double *v1, double *v2, double *v3);
 
 int main (int narg, char *arg[])
@@ -18 +19 @@
 if (narg < 2) { 
   printf("\n Usage tompCXX P/x [N M] \n \n");
-  printf(" P --> Calling MultOMP, x ---> Calling Mult \n");
+  printf(" P -> Calling MultOMP2 p -> Calling MultOMP, x -> Calling Mult \n");
   printf(" N (=100): External loop number   M : Array Size (500000) \n");
   exit(0);
@@ -36 +37 @@
 PrtTim("End of Init ");
 
-if (*arg[1] == 'P') { 
+if (*arg[1] == 'p') { 
+  printf("Calling MultOMP2(N=%d, Size= %d) \n", N, M);
+  for(i=0; i<N; i++) MultOMP2(N, M, v1, v2, v3);
+}
+else 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);
@@ -68 +73 @@
   v3[k] = v1[k] * v2[k];
 }
+void MultOMP2(int N, int n, double *v1, double *v2, double *v3)
+{
+int k;
+int i;
+#pragma omp parallel default(shared)
+#pragma omp parallel for
+for(i=0; i<N; i++)
+#pragma omp parallel for
+  for(k=0; k<n; k++)
+     v3[k] = v1[k] * v2[k];
+}