/[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 2881 by jfenwick, Thu Jan 28 02:03:15 2010 UTC revision 3259 by jfenwick, Mon Oct 11 01:48:14 2010 UTC
# Line 38  void  Paso_SystemMatrix_MatrixVector(con Line 38  void  Paso_SystemMatrix_MatrixVector(con
38                                       double* out) {                                       double* out) {
39    
40    /*double *snd_buffer=NULL, *rcv_buffer=NULL;*/    /*double *snd_buffer=NULL, *rcv_buffer=NULL;*/
41    Paso_MPIInfo *mpi_info=A->mpi_info;    Esys_MPIInfo *mpi_info=A->mpi_info;
42    
43    if (A->type & MATRIX_FORMAT_CSC) {    if (A->type & MATRIX_FORMAT_CSC) {
44       if ( mpi_info->size>1) {       if ( mpi_info->size>1) {
45             Paso_setError(SYSTEM_ERROR,"Paso_SystemMatrix_MatrixVector: CSC is not supported by MPI.");             Esys_setError(SYSTEM_ERROR,"Paso_SystemMatrix_MatrixVector: CSC is not supported by MPI.");
46             return;             return;
47       } else {       } else {
48         if (A->type & MATRIX_FORMAT_OFFSET1) {         if (A->type & MATRIX_FORMAT_OFFSET1) {
# Line 52  void  Paso_SystemMatrix_MatrixVector(con Line 52  void  Paso_SystemMatrix_MatrixVector(con
52         }         }
53       }       }
54    } else if (A->type & MATRIX_FORMAT_TRILINOS_CRS) {    } else if (A->type & MATRIX_FORMAT_TRILINOS_CRS) {
55             Paso_setError(SYSTEM_ERROR,"Paso_SystemMatrix_MatrixVector: TRILINOS is not supported with MPI.");             Esys_setError(SYSTEM_ERROR,"Paso_SystemMatrix_MatrixVector: TRILINOS is not supported with MPI.");
56             return;             return;
57    } else {    } else {
58       if (A->type & MATRIX_FORMAT_OFFSET1) {       if (A->type & MATRIX_FORMAT_OFFSET1) {
59             if ( mpi_info->size>1) {             if ( mpi_info->size>1) {
60                Paso_setError(SYSTEM_ERROR,"Paso_SystemMatrix_MatrixVector: CSR with index 1 is not supported by MPI.");                Esys_setError(SYSTEM_ERROR,"Paso_SystemMatrix_MatrixVector: CSR with index 1 is not supported by MPI.");
61                return;                return;
62             } else {             } else {
63                Paso_SparseMatrix_MatrixVector_CSR_OFFSET1(alpha,A->mainBlock,in,beta,out);                Paso_SparseMatrix_MatrixVector_CSR_OFFSET1(alpha,A->mainBlock,in,beta,out);
64             }             }
65       } else {       } else {
66           if (Paso_noError()) {           if (Esys_noError()) {
67              Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(alpha,A,in,beta,out);              Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(alpha,A,in,beta,out);
68           }           }
69       }       }

Legend:
Removed from v.2881  
changed lines
  Added in v.3259

  ViewVC Help
Powered by ViewVC 1.1.26