Zeros all entries (except possibly the main diagonal) of a set of rows and columns of a matrix.
mat - the matrix
numRows - the number of rows/columns to zero
rows - the global row indices
diag - value put in the diagonal of the eliminated rows
x - optional vector of the solution for zeroed rows (other entries in vector are not used), these must be set before this call
b - optional vector of the right hand side, that will be adjusted by provided solution entries
This routine, along with
MatZeroRows(), is typically used to eliminate known Dirichlet boundary conditions from a linear system.
For each zeroed row, the value of the corresponding
b is set to diag times the value of the corresponding
The other entries of
b will be adjusted by the known values of
x times the corresponding matrix entries in the columns that are being eliminated
If the resulting linear system is to be solved with
KSP then one can (but does not have to) call
KSPSetInitialGuessNonzero() to allow the
Krylov method to take advantage of the known solution on the zeroed rows.
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.
MatZeroRows() this does not change the nonzero structure of the matrix, it merely zeros those entries in the matrix.
Each processor can indicate any rows in the entire matrix to be zeroed (i.e. each process does NOT have to list only rows local to itself).
MAT_NO_OFF_PROC_ZERO_ROWS does not apply to this routine.
MatZeroRowsColumns_MPIAIJ in src/mat/impls/aij/mpi/mpiaij.c
MatZeroRowsColumns_SeqAIJ in src/mat/impls/aij/seq/aij.c
MatZeroRowsColumns_MPIBAIJ in src/mat/impls/baij/mpi/mpibaij.c
MatZeroRowsColumns_SeqBAIJ in src/mat/impls/baij/seq/baij.c
MatZeroRowsColumns_SeqDense in src/mat/impls/dense/seq/dense.c
MatZeroRowsColumns_HYPRE in src/mat/impls/hypre/mhypre.c
MatZeroRowsColumns_IS in src/mat/impls/is/matis.c
MatZeroRowsColumns_SeqSBAIJ in src/mat/impls/sbaij/seq/sbaij.c
MatZeroRowsColumns_Shell in src/mat/impls/shell/shell.c