/[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 2475 by artak, Wed Jun 17 01:48:46 2009 UTC revision 2548 by jfenwick, Mon Jul 20 06:20:06 2009 UTC
# Line 1  Line 1 
1    
2  /*******************************************************  /*******************************************************
3  *  *
4  * Copyright (c) 2003-2008 by University of Queensland  * Copyright (c) 2003-2009 by University of Queensland
5  * Earth Systems Science Computational Center (ESSCC)  * Earth Systems Science Computational Center (ESSCC)
6  * http://www.uq.edu.au/esscc  * http://www.uq.edu.au/esscc
7  *  *
# 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 228  void Paso_Pattern_Aggregiation(Paso_Spar Line 228  void Paso_Pattern_Aggregiation(Paso_Spar
228    }    }
229    
230    
231      #pragma omp parallel for private(i,iptr,diag) schedule(static)      #pragma omp parallel for private(i,iptr) reduction(+:diag) schedule(static)
232        for (i=0;i<n;++i) {        for (i=0;i<n;++i) {
233          diag = 0;          diag = 0;
234          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.2475  
changed lines
  Added in v.2548

  ViewVC Help
Powered by ViewVC 1.1.26