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

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

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

revision 969 by ksteube, Tue Feb 13 23:02:23 2007 UTC revision 971 by ksteube, Wed Feb 14 04:40:49 2007 UTC
# Line 53  void Paso_SystemMatrix_nullifyRowsAndCol Line 53  void Paso_SystemMatrix_nullifyRowsAndCol
53        }        }
54      }      }
55        }        }
     } else if (A->type & MATRIX_FORMAT_TRILINOS_CRS) {  
       /*  
         Plan: Collect non-zeros on each CPU,  
     use MPI to concatenate,  
     sort with qsort,  
     and then use binary search to check if a row is non-zero  
       */  
       fprintf(stderr, "SystemMatrix_nullifyRowsAndCols: need to implement MATRIX_FORMAT_TRILINOS_CRS (and for block sizes != 1)\n");  
       exit(1);  
56      } else {      } else {
57        #pragma omp parallel for private(irow, iptr,icol) schedule(static)        #pragma omp parallel for private(irow, iptr,icol) schedule(static)
58        for (irow=0;irow< A->pattern->n_ptr;irow++) {        for (irow=0;irow< A->pattern->n_ptr;irow++) {
# Line 99  void Paso_SystemMatrix_nullifyRowsAndCol Line 90  void Paso_SystemMatrix_nullifyRowsAndCol
90        }        }
91      }      }
92        }        }
     } else if (A->type & MATRIX_FORMAT_TRILINOS_CRS) {  
       fprintf(stderr, "SystemMatrix_nullifyRowsAndCols 2: need to implement MATRIX_FORMAT_TRILINOS_CRS\n");  
       exit(1);  
93      } else {      } else {
94        #pragma omp parallel for private(l,irow, iptr,icol,ir,irb,icb) schedule(static)        #pragma omp parallel for private(l,irow, iptr,icol,ir,irb,icb) schedule(static)
95        for (ir=0;ir< A->pattern->n_ptr;ir++) {        for (ir=0;ir< A->pattern->n_ptr;ir++) {

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

  ViewVC Help
Powered by ViewVC 1.1.26