Jose David Bermeol <[email protected]> writes: > So for which solvers MatMatSolve does not work??
There is a default implementation (MatMatSolve_Basic) that just solves single vector problems in sequence. So you can use the interface, but the execution might vary between solvers. The MUMPS interface is actually a liability here because they insist on having the entire vector on rank 0. This is a terrible interface and not memory-scalable. If you naively put many RHS vectors on rank 0, you would likely run out of memory.
pgpc417bx9hOT.pgp
Description: PGP signature
