/[escript]/branches/doubleplusgood/paso/src/SystemMatrix_MatrixVector.cpp
ViewVC logotype

Diff of /branches/doubleplusgood/paso/src/SystemMatrix_MatrixVector.cpp

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

revision 3369 by gross, Fri Nov 19 06:26:11 2010 UTC revision 3445 by gross, Wed Jan 19 06:02:15 2011 UTC
# Line 84  void  Paso_SystemMatrix_MatrixVector_CSR Line 84  void  Paso_SystemMatrix_MatrixVector_CSR
84    /* start exchange */    /* start exchange */
85    Paso_SystemMatrix_startCollect(A,in);    Paso_SystemMatrix_startCollect(A,in);
86    /* process main block */    /* process main block */
87    Paso_SparseMatrix_MatrixVector_CSR_OFFSET0(alpha,A->mainBlock,in,beta,out);    if (A->type & MATRIX_FORMAT_DIAGONAL_BLOCK) {
88         Paso_SparseMatrix_MatrixVector_CSR_OFFSET0_DIAG(alpha,A->mainBlock,in,beta,out);
89      } else {
90         Paso_SparseMatrix_MatrixVector_CSR_OFFSET0(alpha,A->mainBlock,in,beta,out);
91      }
92    /* finish exchange */    /* finish exchange */
93    remote_values=Paso_SystemMatrix_finishCollect(A);    remote_values=Paso_SystemMatrix_finishCollect(A);
94    /* process couple block */    /* process couple block */
95    if (A->col_coupleBlock->pattern->ptr!=NULL) {    if (A->col_coupleBlock->pattern->ptr!=NULL) {
96        Paso_SparseMatrix_MatrixVector_CSR_OFFSET0(alpha,A->col_coupleBlock,remote_values,1.,out);        if (A->type & MATRIX_FORMAT_DIAGONAL_BLOCK) {
97             Paso_SparseMatrix_MatrixVector_CSR_OFFSET0_DIAG(alpha,A->col_coupleBlock,remote_values,1.,out);
98          } else {
99             Paso_SparseMatrix_MatrixVector_CSR_OFFSET0(alpha,A->col_coupleBlock,remote_values,1.,out);
100          }
101    }    }
102        
103  }  }

Legend:
Removed from v.3369  
changed lines
  Added in v.3445

  ViewVC Help
Powered by ViewVC 1.1.26