[petsc-users] SLEPc - control of MUMPS in Fortran

[Barry Smith]

  Please send a Fortran code that demonstrates the difficulty. In theory one can do this in Fortran as simply as C but since it is far less used there could be problems.

    Barry

> On Sep 19, 2017, at 11:02 AM, Appel, Thibaut <t.appel17 at imperial.ac.uk> wrote:
> 
> Good afternoon,
> 
> I am using SLEPc in shift-and-invert mode to solve linear generalized EVPs. I am also using MUMPS for solving the linear systems that stem from the problem. After browsing a lot of pages, I am quite surprised by the lack of documentation available regarding the MUMPS interface with SLEPc/PETSc.
> 
> - I would like, in my code, to set and print MUMPS control and output parameters (CNTL,ICNTL,INFO) to monitor the process but it seems like a tedious task in Fortran. I found this discussion https://lists.mcs.anl.gov/pipermail/petsc-users/2012-May/013437.html that suggests hard-coded statements in c++ to set MUMPS parameters. Is there a clean solution in Fortran? 
> 
> - When using MUMPS together with SLEPc, can you perform iterative refinement? I tried to set ICNTL(10) to do so but nothing happens. Can I trust SLEPc for automatically setting MUMPS parameters that will ensure best scalability? If you have further advice for using optimally SLEPc together with MUMPS I would really appreciate it and enjoy discussing it with you.
> 
> Thank you in advance,
> 
> 
> Thibaut