#include "petscmat.h" PetscErrorCode MatZeroRowsIS(Mat mat,IS is,PetscScalar diag)Collective on Mat
mat | - the matrix | |
is | - index set of rows to remove | |
diag | - value put in all diagonals of eliminated rows |
If the option MatSetOption(mat,MAT_KEEP_ZEROED_ROWS) the nonzero structure of the matrix is not changed (even for AIJ and BAIJ matrices) the values are merely zeroed.
The user can set a value in the diagonal entry (or for the AIJ and row formats can optionally remove the main diagonal entry from the nonzero structure as well, by passing 0.0 as the final argument).
For the parallel case, all processes that share the matrix (i.e., those in the communicator used for matrix creation) MUST call this routine, regardless of whether any rows being zeroed are owned by them.
Each processor should list the rows that IT wants zeroed
Level:intermediate
Location:src/mat/interface/matrix.c
Index of all Mat routines
Table of Contents for all manual pages
Index of all manual pages