/[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 2686 by artak, Tue Sep 29 03:39:36 2009 UTC revision 2699 by artak, Wed Sep 30 05:43:20 2009 UTC
# Line 162  void Paso_Pattern_RS(Paso_SparseMatrix* Line 162  void Paso_Pattern_RS(Paso_SparseMatrix*
162      Paso_setError(TYPE_ERROR,"Paso_Pattern_RS: symmetric matrix pattern is not supported yet");      Paso_setError(TYPE_ERROR,"Paso_Pattern_RS: symmetric matrix pattern is not supported yet");
163      return;      return;
164    }    }
165      #pragma omp parallel for private(i,iptr,max_offdiagonal,threshold,j) schedule(static)      /*#pragma omp parallel for private(i,iptr,max_offdiagonal,threshold,j) schedule(static)*/
166      for (i=0;i<n;++i) {      for (i=0;i<n;++i) {
167        if(mis_marker[i]==IS_AVAILABLE) {        if(mis_marker[i]==IS_AVAILABLE) {
168          max_offdiagonal = DBL_MIN;          max_offdiagonal = DBL_MIN;
# Line 287  void Paso_Pattern_greedy(Paso_Pattern* p Line 287  void Paso_Pattern_greedy(Paso_Pattern* p
287      return;      return;
288    }    }
289        
290     /* We do not need this loop if we set IS_IN_MIS=IS_AVAILABLE. */     /* We do not need this loop if we set IS_REMOVED=IS_AVAILABLE. */
291     #pragma omp parallel for private(i) schedule(static)     #pragma omp parallel for private(i) schedule(static)
292     for (i=0;i<n;++i)     for (i=0;i<n;++i)
293          if(mis_marker[i]==IS_AVAILABLE)          if(mis_marker[i]==IS_AVAILABLE)
# Line 318  void Paso_Pattern_greedy(Paso_Pattern* p Line 318  void Paso_Pattern_greedy(Paso_Pattern* p
318                      break;                      break;
319                  }                  }
320                }                }
            }  
321             if (passed) mis_marker[i]=IS_REMOVED;             if (passed) mis_marker[i]=IS_REMOVED;
322               }
323          }          }
324    
325       /* swap to TRUE/FALSE in mis_marker */       /* swap to TRUE/FALSE in mis_marker */

Legend:
Removed from v.2686  
changed lines
  Added in v.2699

  ViewVC Help
Powered by ViewVC 1.1.26