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

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

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

revision 2551 by gross, Thu Jul 23 09:19:15 2009 UTC revision 2552 by artak, Fri Jul 24 00:42:10 2009 UTC
# Line 178  void Paso_Pattern_RS(Paso_SparseMatrix* Line 178  void Paso_Pattern_RS(Paso_SparseMatrix*
178              if(A->val[iptr]<=threshold) {              if(A->val[iptr]<=threshold) {
179                 if(j!=i) {                 if(j!=i) {
180                  Paso_IndexList_insertIndex(&(index_list[i]),j);                  Paso_IndexList_insertIndex(&(index_list[i]),j);
181                  Paso_IndexList_insertIndex(&(index_list[j]),i);                  /*Paso_IndexList_insertIndex(&(index_list[j]),i);*/
182                  }                  }
183              }              }
184          }          }
# Line 229  void Paso_Pattern_Aggregiation(Paso_Spar Line 229  void Paso_Pattern_Aggregiation(Paso_Spar
229    }    }
230    
231    
232      #pragma omp parallel for private(i,iptr) reduction(+:diag) schedule(static)      #pragma omp parallel for private(i,iptr,diag) schedule(static)
233        for (i=0;i<n;++i) {        for (i=0;i<n;++i) {
234          diag = 0;          diag = 0;
235          for (iptr=A->pattern->ptr[i];iptr<A->pattern->ptr[i+1]; ++iptr) {          for (iptr=A->pattern->ptr[i];iptr<A->pattern->ptr[i+1]; ++iptr) {

Legend:
Removed from v.2551  
changed lines
  Added in v.2552

  ViewVC Help
Powered by ViewVC 1.1.26