[petsc-users] Using iterative refinement of MUMPS from PETSC

[Barry Smith]

  Evan,

   Just use PETSc to do the iterative refinement for you. For example -ksp_type gmres -pc_type lu -ksp_max_its 2 

   Barry

On Sep 11, 2014, at 3:24 PM, Evan Um <evanum at gmail.com> wrote:

> Dear PETSC and MUMPS users,
> 
> I try to use an iterative refinement option (ICNTL(10)=max # of iterative refinement) of MUMPS in my PETSC application. MUMPS manual says that if the solution is kept distributed (ICNTL(21)=1), the iterative refinement option is disabled. When a problem is solved using KSPSolve() with multiple cores, the solution is automatically kept distributed over the processors. Does this mean that iterative refinement of MUMPS is available from PETSC only when we run our application on a single core? I tried to solve my problem on multiple cores with ICNTL(10)=10 and INCTL(21)=0, but the program  crashed. Does anyone know how to use iterative refinement option of MUMPS from PETSC on a parallel computer? In advance, thanks for your comments.
> 
> Regards,
> Evan
> 
> 
>