Obtain the Schur complement from eliminating part of the matrix in another part.


#include "petscksp.h" 
PetscErrorCode MatGetSchurComplement(Mat A, IS isrow0, IS iscol0, IS isrow1, IS iscol1, MatReuse mreuse, Mat *S, MatSchurComplementAinvType ainvtype, MatReuse preuse, Mat *Sp)


Input Parameters#

Output Parameters#

  • S - exact Schur complement, often of type MATSCHURCOMPLEMENT which is difficult to use for preconditioning

  • Sp - approximate Schur complement from which a preconditioner can be built A11 - A10 inv(DIAGFORM(A00)) A01


Since the real Schur complement is usually dense, providing a good approximation to Sp usually requires application-specific information.

Sometimes users would like to provide problem-specific data in the Schur complement, usually only for special row and column index sets. In that case, the user should call PetscObjectComposeFunction() on the *S matrix and pass mreuse of MAT_REUSE_MATRIX to set “MatGetSchurComplement_C” to their function. If their function needs to fall back to the default implementation, it should call MatGetSchurComplement_Basic().

MatCreateSchurComplement() takes as arguments the four submatrices and returns the virtual Schur complement (what this function returns in S).

MatSchurComplementGetPmat() takes the virtual Schur complement and returns an explicit approximate Schur complement (what this returns in Sp).

In other words calling MatCreateSchurComplement() followed by MatSchurComplementGetPmat() produces the same output as this function but with slightly different inputs. The actually submatrices of the original block matrix instead of index sets to the submatrices.

Developer Note#

The API that includes MatGetSchurComplement(), MatCreateSchurComplement(), MatSchurComplementGetPmat() should be refactored to remove redundancy and be clearer and simpler.

See Also#

KSP: Linear System Solvers, MatCreateSubMatrix(), PCFIELDSPLIT, MatCreateSchurComplement(), MatSchurComplementAinvType





Edit on GitLab

Index of all KSP routines
Table of Contents for all manual pages
Index of all manual pages