petscmaster 20200805
MatReorderForNonzeroDiagonal
Changes matrix ordering to remove zeros from diagonal. This may help in the LU factorization to prevent a zero pivot.
Synopsis
#include "petscmat.h"
PetscErrorCode MatReorderForNonzeroDiagonal(Mat mat,PetscReal abstol,IS ris,IS cis)
Collective on Mat
Input Parameters
 mat   matrix to reorder

 rmap,cmap   row and column permutations. Usually obtained from
MatGetOrdering().

Notes
This is not intended as a replacement for pivoting for matrices that
have ``bad'' structure. It is only a stopgap measure. Should be called
after a call to MatGetOrdering(), this routine changes the column
ordering defined in cis.
Only works for SeqAIJ matrices
Options Database Keys (When using KSP)
pc_factor_nonzeros_along_diagonal 
Algorithm Notes:
Column pivoting is used.
1) Choice of column is made by looking at the
nonzero elements in the troublesome row for columns that are not yet
included (moving from left to right).
2) If (1) fails we check all the columns to the left of the current row
and see if one of them has could be swapped. It can be swapped if
its corresponding row has a nonzero in the column it is being
swapped with; to make sure the previous nonzero diagonal remains
nonzero
Level
intermediate
Location
src/mat/utils/zerodiag.c
Index of all Mat routines
Table of Contents for all manual pages
Index of all manual pages