<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">

</head>

<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">

Good afternoon,

<div class=""><br class="">

</div>

<div class="">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.</div>

<div class=""><br class="">

</div>

<div class="">- 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

<a href="https://lists.mcs.anl.gov/pipermail/petsc-users/2012-May/013437.html" class="">

https://lists.mcs.anl.gov/pipermail/petsc-users/2012-May/013437.html</a> that suggests hard-coded statements in c++ to set MUMPS parameters. Is there a clean solution in Fortran? </div>

<div class=""><br class="">

</div>

<div class="">- 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.</div>

<div class=""><br class="">

</div>

<div class="">Thank you in advance,</div>

<div class=""><br class="">

</div>

<div class=""><br class="">

</div>

<div class="">Thibaut </div>

</body>

</html>