petsc-3.10.5 2019-03-28
Report Typos and Errors


Solves A X = B, given a factored matrix.


#include "petscmat.h" 
PetscErrorCode MatMatSolve(Mat A,Mat B,Mat X)
Neighbor-wise Collective on Mat

Input Parameters

A - the factored matrix
B - the right-hand-side matrix (dense matrix)

Output Parameter

X -the result matrix (dense matrix)


The matrices b and x cannot be the same. I.e., one cannot call MatMatSolve(A,x,x).


Most users should usually employ the simplified KSP interface for linear solvers instead of working directly with matrix algebra routines such as this. See, e.g., KSPCreate(). However KSP can only solve for one vector (column of X) at a time.

When using SuperLU_Dist as a parallel solver PETSc will use the SuperLU_Dist functionality to solve multiple right hand sides simultaneously. For MUMPS it calls a separate solve for each right hand side since MUMPS does not yet support distributed right hand sides.

Since the resulting matrix X must always be dense we do not support sparse representation of the matrix B.

See Also

MatMatSolveTranspose(), MatLUFactor(), MatCholeskyFactor()






MatMatSolve_MKL_CPARDISO in src/mat/impls/aij/mpi/mkl_cpardiso/mkl_cpardiso.c
MatMatSolve_MUMPS in src/mat/impls/aij/mpi/mumps/mumps.c
MatMatSolve_STRUMPACK in src/mat/impls/aij/mpi/strumpack/strumpack.c
MatMatSolve_SuperLU_DIST in src/mat/impls/aij/mpi/superlu_dist/superlu_dist.c
MatMatSolve_SeqAIJ_inplace in src/mat/impls/aij/seq/aijfact.c
MatMatSolve_SeqAIJ in src/mat/impls/aij/seq/aijfact.c
MatMatSolve_MKL_PARDISO in src/mat/impls/aij/seq/mkl_pardiso/mkl_pardiso.c
MatMatSolve_SuperLU in src/mat/impls/aij/seq/superlu/superlu.c
MatMatSolve_SeqDense in src/mat/impls/dense/seq/dense.c
MatMatSolve_Elemental in src/mat/impls/elemental/matelem.cxx

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