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

Diff of /trunk/esys2/finley/src/finleyC/System_saveMM.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 108 by jgs, Thu Jan 27 06:21:59 2005 UTC
# Line 43  void Finley_SystemMatrix_saveMM(Finley_S Line 43  void Finley_SystemMatrix_saveMM(Finley_S
43           for (ir = 0; ir< A_p->row_block_size; ir++) {           for (ir = 0; ir< A_p->row_block_size; ir++) {
44          for (iPtr = A_p->pattern->ptr[iRow] - PTR_OFFSET;iPtr < A_p->pattern->ptr[iRow+1] - PTR_OFFSET; iPtr++) {          for (iPtr = A_p->pattern->ptr[iRow] - PTR_OFFSET;iPtr < A_p->pattern->ptr[iRow+1] - PTR_OFFSET; iPtr++) {
45                 for (ic = 0; ic< A_p->col_block_size; ic++) {                 for (ic = 0; ic< A_p->col_block_size; ic++) {
46                fprintf(fileHandle_p, "%12d %12d %22.15g\n",                fprintf(fileHandle_p, "%12d %12d %e\n",
47                    iRow*A_p->row_block_size+ir+ 1,                    iRow*A_p->row_block_size+ir+ 1,
48                    (A_p->pattern->index[iPtr]-INDEX_OFFSET)*A_p->col_block_size+ic+1,                    (A_p->pattern->index[iPtr]-INDEX_OFFSET)*A_p->col_block_size+ic+1,
49                    A_p->val[iPtr*A_p->block_size+ir+A_p->row_block_size*ic]);                    A_p->val[iPtr*A_p->block_size+ir+A_p->row_block_size*ic]);
# Line 57  void Finley_SystemMatrix_saveMM(Finley_S Line 57  void Finley_SystemMatrix_saveMM(Finley_S
57           for (ic = 0; ic< A_p->col_block_size; ic++) {           for (ic = 0; ic< A_p->col_block_size; ic++) {
58          for (iPtr = A_p->pattern->ptr[iCol] - PTR_OFFSET;iPtr < A_p->pattern->ptr[iCol+1] - PTR_OFFSET; iPtr++) {          for (iPtr = A_p->pattern->ptr[iCol] - PTR_OFFSET;iPtr < A_p->pattern->ptr[iCol+1] - PTR_OFFSET; iPtr++) {
59                 for (ir = 0; ir< A_p->row_block_size; ir++) {                 for (ir = 0; ir< A_p->row_block_size; ir++) {
60                fprintf(fileHandle_p, "%12d %12d %22.15g\n",                fprintf(fileHandle_p, "%12d %12d %e\n",
61                    (A_p->pattern->index[iPtr]-INDEX_OFFSET)*A_p->row_block_size+ir+1,                    (A_p->pattern->index[iPtr]-INDEX_OFFSET)*A_p->row_block_size+ir+1,
62                    iCol*A_p->col_block_size+ic+ 1,                    iCol*A_p->col_block_size+ic+ 1,
63                    A_p->val[iPtr*A_p->block_size+ir+A_p->row_block_size*ic]);                    A_p->val[iPtr*A_p->block_size+ir+A_p->row_block_size*ic]);
# Line 72  void Finley_SystemMatrix_saveMM(Finley_S Line 72  void Finley_SystemMatrix_saveMM(Finley_S
72        
73    return;    return;
74  }  }
 /*  
  * $Log$  
  * Revision 1.4  2004/12/15 07:08:33  jgs  
  * *** empty log message ***  
  *  
  *  
  *  
  */  

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

  ViewVC Help
Powered by ViewVC 1.1.26