/[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 2446 by artak, Fri May 29 00:39:44 2009 UTC revision 2447 by artak, Fri May 29 04:32:53 2009 UTC
# Line 222  void Paso_Pattern_Aggregiation(Paso_Spar Line 222  void Paso_Pattern_Aggregiation(Paso_Spar
222      }      }
223            
224    if (A->pattern->type & PATTERN_FORMAT_SYM) {    if (A->pattern->type & PATTERN_FORMAT_SYM) {
225      Paso_setError(TYPE_ERROR,"Paso_Pattern_RS: symmetric matrix pattern is not supported yet");      Paso_setError(TYPE_ERROR,"Paso_Pattern_Aggregiation: symmetric matrix pattern is not supported yet");
226      return;      return;
227    }    }
228    
# Line 231  void Paso_Pattern_Aggregiation(Paso_Spar Line 231  void Paso_Pattern_Aggregiation(Paso_Spar
231        for (i=0;i<n;++i) {        for (i=0;i<n;++i) {
232          diag = 0;          diag = 0;
233          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) {
234              if(A->pattern->index[iptr] != i){              if(A->pattern->index[iptr] == i){
235                  diag+=A->val[iptr];                  diag+=A->val[iptr];
236              }              }
237          }          }

Legend:
Removed from v.2446  
changed lines
  Added in v.2447

  ViewVC Help
Powered by ViewVC 1.1.26