/[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 631 by dhawcroft, Thu Mar 23 04:27:32 2006 UTC revision 969 by ksteube, Tue Feb 13 23:02:23 2007 UTC
# Line 2  Line 2 
2    
3  /*  /*
4  ********************************************************************************  ********************************************************************************
5  *               Copyright 2006 by ACcESS MNRF                                *  *               Copyright   2006 by ACcESS MNRF                                *
6  *                                                                              *  *                                                                              *
7  *                 http://www.access.edu.au                                     *  *                 http://www.access.edu.au                                     *
8  *           Primary Business: Queensland, Australia                            *  *           Primary Business: Queensland, Australia                            *
# 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.631  
changed lines
  Added in v.969

  ViewVC Help
Powered by ViewVC 1.1.26