On Thu, Sep 11, 2014 at 3:24 PM, Evan Um <[email protected]> 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. >
Just use Newton to solve your system, which for a linear operator is iterative refinement. Matt > Regards, > Evan > > > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener
