/[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 1902 by artak, Wed Oct 22 03:54:14 2008 UTC revision 1907 by phornby, Wed Oct 22 11:41:18 2008 UTC
# Line 60  void Paso_Pattern_coup(Paso_SparseMatrix Line 60  void Paso_Pattern_coup(Paso_SparseMatrix
60       /* is there any vertex available ?*/       /* is there any vertex available ?*/
61       while (Paso_Util_isAny(n,mis_marker,IS_AVAILABLE)) {       while (Paso_Util_isAny(n,mis_marker,IS_AVAILABLE)) {
62    
63             #pragma omp parallel for private(naib,i,iptr,flag) schedule(static)             #pragma omp parallel for private(i,iptr,flag) schedule(static)
64             for (i=0;i<n;++i) {             for (i=0;i<n;++i) {
65                if (mis_marker[i]==IS_AVAILABLE) {                if (mis_marker[i]==IS_AVAILABLE) {
66                   flag=IS_IN_MIS;                   flag=IS_IN_MIS;
# Line 87  void Paso_Pattern_coup(Paso_SparseMatrix Line 87  void Paso_Pattern_coup(Paso_SparseMatrix
87                  }                  }
88              }              }
89                        
90                #pragma omp parallel for private(naib,i,iptr) schedule(static)                #pragma omp parallel for private(i,iptr) schedule(static)
91                for (i=0;i<n;i++) {                for (i=0;i<n;i++) {
92                 if (mis_marker[i]==IS_AVAILABLE) {                 if (mis_marker[i]==IS_AVAILABLE) {
93                   diagptr=A->pattern->ptr[i];                   diagptr=A->pattern->ptr[i];
# Line 136  void Paso_Pattern_RS(Paso_SparseMatrix* Line 136  void Paso_Pattern_RS(Paso_SparseMatrix*
136  {  {
137    index_t index_offset=(A->pattern->type & PATTERN_FORMAT_OFFSET1 ? 1:0);    index_t index_offset=(A->pattern->type & PATTERN_FORMAT_OFFSET1 ? 1:0);
138    dim_t i,j;    dim_t i,j;
139    index_t naib,iptr;    index_t iptr;
140    double threshold,min_offdiagonal;    double threshold,min_offdiagonal;
141    bool_t flag;    bool_t flag;
142    dim_t n=A->pattern->numOutput;    dim_t n=A->pattern->numOutput;

Legend:
Removed from v.1902  
changed lines
  Added in v.1907

  ViewVC Help
Powered by ViewVC 1.1.26