/[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 1935 by artak, Mon Oct 27 01:31:28 2008 UTC revision 1936 by artak, Mon Oct 27 23:30:57 2008 UTC
# Line 47  void Paso_Pattern_coup(Paso_SparseMatrix Line 47  void Paso_Pattern_coup(Paso_SparseMatrix
47    dim_t i,j;    dim_t i,j;
48    /*double threshold=0.05;*/    /*double threshold=0.05;*/
49    index_t iptr,*index,*where_p,diagptr;    index_t iptr,*index,*where_p,diagptr;
50    bool_t flag,fail;    bool_t fail;
51    dim_t n=A->pattern->numOutput;    dim_t n=A->pattern->numOutput;
52    if (A->pattern->type & PATTERN_FORMAT_SYM) {    if (A->pattern->type & PATTERN_FORMAT_SYM) {
53      Paso_setError(TYPE_ERROR,"Paso_Pattern_mis: symmetric matrix pattern is not supported yet");      Paso_setError(TYPE_ERROR,"Paso_Pattern_mis: symmetric matrix pattern is not supported yet");
# Line 162  void Paso_Pattern_RS(Paso_SparseMatrix* Line 162  void Paso_Pattern_RS(Paso_SparseMatrix*
162          for (iptr=A->pattern->ptr[i]-index_offset;iptr<A->pattern->ptr[i+1]-index_offset; ++iptr) {          for (iptr=A->pattern->ptr[i]-index_offset;iptr<A->pattern->ptr[i+1]-index_offset; ++iptr) {
163              if(A->val[iptr-index_offset] < threshold){              if(A->val[iptr-index_offset] < threshold){
164                      mis_marker[i]=IS_CONNECTED_TO_MIS;                      mis_marker[i]=IS_CONNECTED_TO_MIS;
165                      break;                      #ifndef _OPENMP    
166                         break;
167                        #endif
168              }              }
169          }          }
170        }        }
Legend:
Removed from v.1935  
changed lines
  Added in v.1936
  ViewVC Help
Powered by ViewVC 1.1.26