/[escript]/trunk/esys2/finley/src/finleyC/System_add.c
ViewVC logotype

Diff of /trunk/esys2/finley/src/finleyC/System_add.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 102 by jgs, Wed Dec 15 07:08:39 2004 UTC revision 117 by jgs, Fri Apr 1 05:48:57 2005 UTC
# Line 72  void  Finley_SystemMatrix_add(Finley_Sys Line 72  void  Finley_SystemMatrix_add(Finley_Sys
72                  for (k=in->pattern->ptr[iptr]-PTR_OFFSET;k<in->pattern->ptr[iptr+1]-PTR_OFFSET;++k) {                  for (k=in->pattern->ptr[iptr]-PTR_OFFSET;k<in->pattern->ptr[iptr+1]-PTR_OFFSET;++k) {
73                      if (in->pattern->index[k]==index) {                      if (in->pattern->index[k]==index) {
74                            for (ic=0;ic<col_block_size;++ic) {                            for (ic=0;ic<col_block_size;++ic) {
75                                  i_Sol=ic+num_subblocks_Sol*l_col;                                  i_Sol=ic+col_block_size*l_col;
76                                  for (ir=0;ir<row_block_size;++ir) {                                  for (ir=0;ir<row_block_size;++ir) {
77                                     i_Equa=ir+num_subblocks_Equa*l_row;                                     i_Equa=ir+row_block_size*l_row;
78                             in->val[k*block_size+ir+row_block_size*ic]+=                             in->val[k*block_size+ir+row_block_size*ic]+=
79                                             array[INDEX4(i_Equa,i_Sol,k_Equa,k_Sol,num_Equa,num_Sol,NN_Equa)];                                             array[INDEX4(i_Equa,i_Sol,k_Equa,k_Sol,num_Equa,num_Sol,NN_Equa)];
80                                  }                                  }
# Line 88  void  Finley_SystemMatrix_add(Finley_Sys Line 88  void  Finley_SystemMatrix_add(Finley_Sys
88           }           }
89     }     }
90  }  }
 /*  
  * $Log$  
  * Revision 1.4  2004/12/15 07:08:33  jgs  
  * *** empty log message ***  
  *  
  *  
  *  
  */  

Legend:
Removed from v.102  
changed lines
  Added in v.117

  ViewVC Help
Powered by ViewVC 1.1.26