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

Diff of /temp_trunk_copy/paso/src/SystemMatrix_nullifyRowsAndCols.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 53  void Paso_SystemMatrix_nullifyRowsAndCol Line 53  void Paso_SystemMatrix_nullifyRowsAndCol
53        }        }
54      }      }
55        }        }
56        } else if (A->type & MATRIX_FORMAT_TRILINOS_CRS) {
57          /*
58            Plan: Collect non-zeros on each CPU,
59        use MPI to concatenate,
60        sort with qsort,
61        and then use binary search to check if a row is non-zero
62          */
63          fprintf(stderr, "SystemMatrix_nullifyRowsAndCols: need to implement MATRIX_FORMAT_TRILINOS_CRS (and for block sizes != 1)\n");
64          exit(1);
65      } else {      } else {
66        #pragma omp parallel for private(irow, iptr,icol) schedule(static)        #pragma omp parallel for private(irow, iptr,icol) schedule(static)
67        for (irow=0;irow< A->pattern->n_ptr;irow++) {        for (irow=0;irow< A->pattern->n_ptr;irow++) {
# Line 90  void Paso_SystemMatrix_nullifyRowsAndCol Line 99  void Paso_SystemMatrix_nullifyRowsAndCol
99        }        }
100      }      }
101        }        }
102        } else if (A->type & MATRIX_FORMAT_TRILINOS_CRS) {
103          fprintf(stderr, "SystemMatrix_nullifyRowsAndCols 2: need to implement MATRIX_FORMAT_TRILINOS_CRS\n");
104          exit(1);
105      } else {      } else {
106        #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)
107        for (ir=0;ir< A->pattern->n_ptr;ir++) {        for (ir=0;ir< A->pattern->n_ptr;ir++) {

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

  ViewVC Help
Powered by ViewVC 1.1.26