/[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

branches/more_shared_ptrs_from_1812/paso/src/Pattern_coupling.c revision 1851 by jfenwick, Mon Oct 6 03:16:43 2008 UTC trunk/paso/src/Pattern_coupling.c revision 1890 by artak, Fri Oct 17 00:14:22 2008 UTC
# Line 42  Line 42 
42  #define IS_IN_MIS -3  #define IS_IN_MIS -3
43  #define IS_CONNECTED_TO_MIS -4  #define IS_CONNECTED_TO_MIS -4
44    
45  void Paso_Pattern_cop(Paso_SparseMatrix* A, index_t* mis_marker) {  
46    
47    void Paso_Pattern_coup(Paso_SparseMatrix* A, index_t* mis_marker) {
48    
49    index_t index_offset=(A->pattern->type & PATTERN_FORMAT_OFFSET1 ? 1:0);    index_t index_offset=(A->pattern->type & PATTERN_FORMAT_OFFSET1 ? 1:0);
50    dim_t i;    dim_t i,j;
51    double threshold=0.05;    double threshold=0.05;
52    index_t naib,iptr;    index_t naib,iptr;
53    bool_t flag;    bool_t flag;
# Line 61  void Paso_Pattern_cop(Paso_SparseMatrix* Line 63  void Paso_Pattern_cop(Paso_SparseMatrix*
63             #pragma omp parallel for private(naib,i,iptr,flag) schedule(static)             #pragma omp parallel for private(naib,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_NOW;                   flag=IS_AVAILABLE;
67                   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) {
68                       naib=A->pattern->index[iptr]-index_offset;                       naib=A->pattern->index[iptr]-index_offset;
69                       if (naib!=i && A->val[naib]<threshold*A->val[i]) {                       if (naib!=i && A->val[iptr]>=threshold*A->val[i]) {
70                          flag=IS_AVAILABLE;                          flag=IS_IN_MIS;
71                          break;                          break;
72                       }                       }
73                   }                   }
# Line 78  void Paso_Pattern_cop(Paso_SparseMatrix* Line 80  void Paso_Pattern_cop(Paso_SparseMatrix*
80                if (mis_marker[i]==IS_AVAILABLE) {                if (mis_marker[i]==IS_AVAILABLE) {
81                   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) {
82                       naib=A->pattern->index[iptr]-index_offset;                       naib=A->pattern->index[iptr]-index_offset;
83                       if (naib!=i && mis_marker[i]==IS_IN_MIS_NOW && A->val[naib]/A->val[i]>=-threshold)                       if (naib!=i && mis_marker[naib]==IS_IN_MIS && A->val[iptr]/A->val[i]>=-threshold){
84                           mis_marker[naib]=IS_IN_MIS;                           mis_marker[i]=IS_IN_MIS;
85                         }
86                         else {
87                             mis_marker[i]=IS_CONNECTED_TO_MIS;
88                         }
89                   }                   }
                  mis_marker[i]=IS_CONNECTED_TO_MIS;  
90                }                }
91             }             }
92       }       }
93       /* swap to TRUE/FALSE in mis_marker */       /* swap to TRUE/FALSE in mis_marker */
94       #pragma omp parallel for private(i) schedule(static)       #pragma omp parallel for private(i) schedule(static)
95       for (i=0;i<n;i++) mis_marker[i]=(mis_marker[i]==IS_IN_MIS);       for (i=0;i<n;i++) mis_marker[i]=(mis_marker[i]==IS_IN_MIS);
96    }
97    
98    /*
99     *
100     * Return a strength of connection mask using the classical
101     * strength of connection measure by Ruge and Stuben.
102     *
103     * Specifically, an off-diagonal entry A[i.j] is a strong
104     * connection if:
105     *  
106     *      -A[i,j] >= theta * max( -A[i,k] )   where k != i
107     *
108     * Otherwise, the connection is weak.
109     *  
110     */  
111    void Paso_Pattern_RS(Paso_SparseMatrix* A, index_t* mis_marker, double theta)
112    {
113      index_t index_offset=(A->pattern->type & PATTERN_FORMAT_OFFSET1 ? 1:0);
114      dim_t i,j;
115      index_t naib,iptr;
116      double threshold,min_offdiagonal;
117      bool_t flag;
118      dim_t n=A->pattern->numOutput;
119      if (A->pattern->type & PATTERN_FORMAT_SYM) {
120        Paso_setError(TYPE_ERROR,"Paso_Pattern_mis: symmetric matrix pattern is not supported yet");
121        return;
122      }
123    
124    /* is there any vertex available ?*/
125         while (Paso_Util_isAny(n,mis_marker,IS_AVAILABLE)) {
126    
127          #pragma omp parallel for private(i,iptr,min_offdiagonal,threshold) schedule(static)
128          for (i=0;i<n;++i) {
129            min_offdiagonal = A->val[A->pattern->ptr[i]];
130            for (iptr=A->pattern->ptr[i]-index_offset;iptr<A->pattern->ptr[i+1]-index_offset; ++iptr) {
131                if(A->pattern->index[iptr] != i){
132                    min_offdiagonal = MIN(min_offdiagonal,A->val[iptr]);
133                }
134            }
135    
136            threshold = theta*min_offdiagonal;
137            #pragma omp parallel for private(iptr) schedule(static)
138            for (iptr=A->pattern->ptr[i]-index_offset;iptr<A->pattern->ptr[i+1]-index_offset; ++iptr) {
139                mis_marker[i]=IS_CONNECTED_TO_MIS;
140                if(A->val[iptr] <= threshold){
141                    if(A->pattern->index[iptr] != i){
142                        mis_marker[i]=IS_IN_MIS;
143                    }
144                }
145            }
146          }
147        }
148         /* swap to TRUE/FALSE in mis_marker */
149         #pragma omp parallel for private(i) schedule(static)
150         for (i=0;i<n;i++) mis_marker[i]=(mis_marker[i]==IS_IN_MIS);
151  }  }
152  #undef IS_AVAILABLE  #undef IS_AVAILABLE
153  #undef IS_IN_MIS_NOW  #undef IS_IN_MIS_NOW
Legend:
Removed from v.1851  
changed lines
  Added in v.1890
  ViewVC Help
Powered by ViewVC 1.1.26