Returns a matrix suitable to calls to MatXXFactorSymbolic()


#include "petscmat.h" 
PetscErrorCode MatGetFactor(Mat mat, MatSolverType type, MatFactorType ftype, Mat *f)


Input Parameters#

Output Parameter#

  • f - the factor matrix used with MatXXFactorSymbolic() calls. Can be NULL in some cases, see notes below.

Options Database Key#

  • -mat_factor_bind_factorization <host, device> - Where to do matrix factorization? Default is device (might consume more device memory. One can choose host to save device memory). Currently only supported with MATSEQAIJCUSPARSE matrices.


The return matrix can be NULL if the requested factorization is not available, since some combinations of matrix types and factorization types registered with MatSolverTypeRegister() cannot be fully tested if not at runtime.

Users usually access the factorization solvers via KSP

Some PETSc matrix formats have alternative solvers available that are contained in alternative packages such as pastix, superlu, mumps etc.

PETSc must have been ./configure to use the external solver, using the option –download-package

Some of the packages have options for controlling the factorization, these are in the form -prefix_mat_packagename_packageoption where prefix is normally obtained from the calling KSP/PC. If MatGetFactor() is called directly one can set call MatSetOptionsPrefixFactor() on the originating matrix or MatSetOptionsPrefix() on the resulting factor matrix.

Developer Notes#

This should actually be called MatCreateFactor() since it creates a new factor object

See Also#

Matrices, Mat, Matrix Factorization, KSP, MatSolverType, MatFactorType, MatCopy(), MatDuplicate(), MatGetFactorAvailable(), MatFactorGetCanUseOrdering(), MatSolverTypeRegister(), MAT_FACTOR_LU, MAT_FACTOR_CHOLESKY, MAT_FACTOR_ICC, MAT_FACTOR_ILU, MAT_FACTOR_QR







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