/[escript]/trunk/paso/src/SystemMatrix.c
ViewVC logotype

Diff of /trunk/paso/src/SystemMatrix.c

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

revision 633 by dhawcroft, Thu Mar 23 05:37:00 2006 UTC revision 969 by ksteube, Tue Feb 13 23:02:23 2007 UTC
# Line 75  Paso_SystemMatrix* Paso_SystemMatrix_all Line 75  Paso_SystemMatrix* Paso_SystemMatrix_all
75          out->num_rows=out->pattern->n_index;          out->num_rows=out->pattern->n_index;
76          out->num_cols=out->pattern->n_ptr;          out->num_cols=out->pattern->n_ptr;
77          n_norm = out->num_cols * out->col_block_size;          n_norm = out->num_cols * out->col_block_size;
78         } else if (type & MATRIX_FORMAT_TRILINOS_CRS) {
79            if (type & MATRIX_FORMAT_SYM) {
80               Paso_setError(TYPE_ERROR,"Generation of matrix pattern for symmetric CSC is not implemented yet for TRILINOS.");
81               return NULL;
82            } else {
83    #ifdef TRILINOS
84               Initialize_TrilinosData(out->trilinos_data, pattern);
85           printf("ksteube paso/src/SystemMatrix.c row_block_size=%d col_block_size=%d numCPUs=%d\n", row_block_size, col_block_size, pattern->MPIInfo->size);
86               out->row_block_size=row_block_size;
87               out->col_block_size=col_block_size;
88    #else
89               Paso_setError(TYPE_ERROR,"You need to recompile with -DTRILINOS in order to use setSolverPackage(mypde.TRILINOS).");
90               return NULL;
91    #endif
92            }
93            out->pattern=Paso_SystemMatrixPattern_reference(pattern); /* Count number of references to pattern */
94            out->num_rows=pattern->numLocal;
95            out->num_cols=out->pattern->n_ptr;
96            n_norm = out->num_cols * out->col_block_size;   /* ksteube what is this used for? */
97       } else {       } else {
98          if (type & MATRIX_FORMAT_SYM) {          if (type & MATRIX_FORMAT_SYM) {
99             Paso_setError(TYPE_ERROR,"Generation of matrix pattern for symmetric CSR is not implemented yet.");             Paso_setError(TYPE_ERROR,"Generation of matrix pattern for symmetric CSR is not implemented yet.");

Legend:
Removed from v.633  
changed lines
  Added in v.969

  ViewVC Help
Powered by ViewVC 1.1.26