/[escript]/trunk/paso/src/AMG.c
ViewVC logotype

Diff of /trunk/paso/src/AMG.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3484 by gross, Wed Mar 23 04:06:52 2011 UTC revision 3485 by gross, Thu Mar 24 22:44:40 2011 UTC
# Line 717  void Paso_Preconditioner_AMG_transposeSt Line 717  void Paso_Preconditioner_AMG_transposeSt
717     index_t i, j;     index_t i, j;
718     dim_t p;     dim_t p;
719     dim_t len=0;     dim_t len=0;
720     #pragma omp parallel private(i) schedule(static)     #pragma omp parallel for private(i) schedule(static)
721     for (i=0; i<nT ;++i) {     for (i=0; i<nT ;++i) {
722        degree_ST[i]=0;        degree_ST[i]=0;
723     }     }
# Line 857  printf("%d is in D\n",i); Line 857  printf("%d is in D\n",i);
857       */       */
858       /* w is updated  for local rows only */       /* w is updated  for local rows only */
859       {       {
860          #pragma omp parallel for private(j, jptr)          #pragma omp parallel for private(i, jptr)
861          for (i=0; i< my_n; ++i) {          for (i=0; i< my_n; ++i) {
862    
863             for (jptr=0; jptr<degree_ST[i]; ++jptr) {             for (jptr=0; jptr<degree_ST[i]; ++jptr) {
# Line 870  printf("%d reduced by %d\n",i,j); Line 870  printf("%d reduced by %d\n",i,j);
870             }             }
871                        
872          }          }
873          #pragma omp parallel for private(j)          #pragma omp parallel for private(i, jptr)
874          for (i=my_n; i< n; ++i) {          for (i=my_n; i< n; ++i) {
875             for (jptr=0; jptr<degree_ST[i]; ++jptr) {             for (jptr=0; jptr<degree_ST[i]; ++jptr) {
876            const index_t j = ST[offset_ST[i]+jptr];            const index_t j = ST[offset_ST[i]+jptr];

Legend:
Removed from v.3484  
changed lines
  Added in v.3485

  ViewVC Help
Powered by ViewVC 1.1.26