/[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 1882 by artak, Tue Oct 14 04:39:02 2008 UTC revision 1883 by artak, Tue Oct 14 04:46:27 2008 UTC
# Line 64  void Paso_Pattern_coup(Paso_SparseMatrix Line 64  void Paso_Pattern_coup(Paso_SparseMatrix
64                   flag=IS_AVAILABLE;                   flag=IS_AVAILABLE;
65                   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) {
66                       naib=A->pattern->index[iptr]-index_offset;                       naib=A->pattern->index[iptr]-index_offset;
67                       if (naib!=i && A->val[naib]<threshold*A->val[i]) {                       if (naib!=i && A->val[iptr]>=threshold*A->val[i]) {
68                          flag=IS_IN_MIS;                          flag=IS_IN_MIS;
69                          break;                          break;
70                       }                       }
# Line 78  void Paso_Pattern_coup(Paso_SparseMatrix Line 78  void Paso_Pattern_coup(Paso_SparseMatrix
78                if (mis_marker[i]==IS_AVAILABLE) {                if (mis_marker[i]==IS_AVAILABLE) {
79                   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) {
80                       naib=A->pattern->index[iptr]-index_offset;                       naib=A->pattern->index[iptr]-index_offset;
81                        if (naib!=i && mis_marker[naib]==IS_IN_MIS && A->val[iptr]/A->val[i]>=-threshold)                       if (naib!=i && mis_marker[naib]==IS_IN_MIS && A->val[iptr]/A->val[i]>=-threshold){
82                           mis_marker[i]=IS_IN_MIS;                           mis_marker[i]=IS_IN_MIS;
83                       else                       }
84                         else {
85                           mis_marker[i]=IS_CONNECTED_TO_MIS;                           mis_marker[i]=IS_CONNECTED_TO_MIS;
86                         }
87                   }                   }
88                }                }
89             }             }
Legend:
Removed from v.1882  
changed lines
  Added in v.1883
  ViewVC Help
Powered by ViewVC 1.1.26