petsc-3.13.6 2020-09-29
Report Typos and Errors

MatMatSolve

Solves A X = B, given a factored matrix.

Synopsis

#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 MATDENSE (or sparse -- when using MUMPS)

Output Parameter

X -the result matrix (dense matrix)

Notes

If B is a MATDENSE matrix then one can call MatMatSolve(A,B,B); otherwise, B and X cannot be the same.

Notes

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.

See Also

MatMatSolveTranspose(), MatLUFactor(), MatCholeskyFactor()

Level

developer

Location

src/mat/interface/matrix.c

Examples

src/ksp/ksp/tutorials/ex74.c.html

Implementations

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_SeqDenseCUDA in src/mat/impls/dense/seq/cuda/densecuda.cu
MatMatSolve_SeqDense in src/mat/impls/dense/seq/dense.c
MatMatSolve_Elemental in src/mat/impls/elemental/matelem.cxx
MatMatSolve_CHOLMOD in src/mat/impls/sbaij/seq/cholmod/sbaijcholmod.c
MatMatSolve_SeqSBAIJ_1_NaturalOrdering in src/mat/impls/sbaij/seq/sbaijfact2.c

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