[petsc-users] Using factored complex matrices from MUMPS as a preconditioner in PETSC

[Barry Smith]

> On Sep 23, 2017, at 8:38 PM, Evan Um <evanum at gmail.com> wrote:
> 
> Dear PETSC Users,
> 
> My system matrix comes from finite element modeling and is complex and unstructured. Its typical size is a few millions-by a few millions. I wondering if I can use MUMPS parallel direct solver as a preconditioner in PETSC. For example, I want to pass factored matrices to Krylov iterative solvers such as QMR. Is there any PETSC+MUMPS example code for the purpose?

  You don't pass factored matrices you just pass the original matrix and use -pc_type lu -pc_factor_mat_solver_package mumps

> Can PETSC call the latest MUMPS that supports block low rank approximation? 

  No, send us info on it and we'll see if we can add an interface


> 
> In advance, thank you very much for your comments.
> 
> Best,
> Evan
> 
> 
> 
> 
>