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 <[email protected]> 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
