/[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 2652 by artak, Mon Sep 7 05:04:45 2009 UTC revision 2662 by artak, Tue Sep 15 03:05:23 2009 UTC
# Line 245  void Paso_Pattern_Aggregiation(Paso_Spar Line 245  void Paso_Pattern_Aggregiation(Paso_Spar
245       for (i=0;i<n;++i) {       for (i=0;i<n;++i) {
246         if (mis_marker[i]==IS_AVAILABLE) {         if (mis_marker[i]==IS_AVAILABLE) {
247          eps_Aii = theta*theta*diags[i];          eps_Aii = theta*theta*diags[i];
         val=0.;  
248          for (iptr=A->pattern->ptr[i];iptr<A->pattern->ptr[i+1]; ++iptr) {          for (iptr=A->pattern->ptr[i];iptr<A->pattern->ptr[i+1]; ++iptr) {
249              j=A->pattern->index[iptr];              j=A->pattern->index[iptr];
250              val=A->val[iptr];              val=A->val[iptr];
251              if(j!= i) {              if(j!= i) {
252                if(val*val>=eps_Aii * diags[j]) {                if((val*val)>=(eps_Aii*diags[j])) {
253                 Paso_IndexList_insertIndex(&(index_list[i]),j);                 Paso_IndexList_insertIndex(&(index_list[i]),j);
254                }                }
255              }              }

Legend:
Removed from v.2652  
changed lines
  Added in v.2662

  ViewVC Help
Powered by ViewVC 1.1.26