Changeset 100

Timestamp:

Jan 12, 2010, 1:54:07 PM (14 years ago)

Author:

(none)

Message:

PPPP - jacobi_vector for parallel

File:

• proiecte/PPPP/eigenface/new/jacobi_vector.c (modified) (2 diffs)

proiecte/PPPP/eigenface/new/jacobi_vector.c

@@ -19 +19 @@
         z = (float *) malloc(n * sizeof(float));
         
+        /* Fork a team of threads giving them their own copies of variables */
+        #pragma omp parallel private(nthreads, tid, i, j, ip, iq) shared(tresh, theta, tau, sm, s, h, g, c, b, z, a, n, d, v, nrot)
+        {
+
+        /* Obtain thread number */
+        tid = omp_get_thread_num();
+        printf("Hello World from thread = %d\n", tid);
+
+        /* Only master thread does this */
+        if (tid == 0) 
+    {
+        nthreads = omp_get_num_threads();
+        printf("Number of threads = %d\n", nthreads);
+    }
+
+        
+        #pragma omp parallel for
         for (ip=0;ip<n;ip++) {                                                                                                                  // Initialize to the identity matrix.
                 for (iq=0;iq<n;iq++) v[ip*n + iq]=0.0;
                         v[ip*n + ip]=1.0;
         }
+        #pragma omp parallel for
         for (ip=0;ip<n;ip++) {                                                                                                                  // Initialize b and d to the diagonal of a.
                 b[ip]=d[ip]=a[ip*n + ip];                                                                                                               
                 z[ip]=0.0;                                                                                                                              // This vector will accumulate terms of the form t*a[pq] as in equation (11.1.14).
         }
-        
+        }       
         *nrot=0;
         for (i=1;i<=50;i++) {
                 sm=0.0;
+        #pragma omp parallel for reduction(+:sm)
                 for (ip=0;ip<n-1;ip++) {                                                                                                                // Sum off-diagonal elements.
                         for (iq=ip+1;iq<n;iq++)
@@ -70 +89 @@
                                         a[ip*n + iq]=0.0;
                                         
+                                        #pragma omp parallel for
                                         for (j=0;j<=ip-1;j++) {                                                                                         // Case of rotations 1 <= j < p.
                                                 ROTATE(a,j,ip,j,iq)
                                         }
+                                        #pragma omp parallel for
                                         for (j=ip+1;j<=iq-1;j++) {                                                                              // Case of rotations p < j < q.
                                                 ROTATE(a,ip,j,j,iq)
                                         }
+                                        #pragma omp parallel for
                                         for (j=iq+1;j<n;j++) {                                                                                  // Case of rotations q < j <= n.
                                                 ROTATE(a,ip,j,iq,j)
                                         }
+                                        #pragma omp parallel for
                                         for (j=0;j<n;j++) {
                                                 ROTATE(v,j,ip,j,iq)