Setting MUMPS option

Question

Hi

When I used Dolfin linear algebra backend Petsc, I used MUMPS as linear solver. However, there seems no control switch to set up INCTL manually through Parameters interface. Is it possible to do so ?

The other question is that is it possible to set Petsc command line options through Parameters ?

BR
Bui

Answer 1

PETSc options can be set on the command line (you'll need to add a few lines to parse them), or you can use, for example:

PETScOptions.set("mat_mumps_icntl_14", 40.0)

The above is for Python. From the C++ interface:

PETScOptions::set("mat_mumps_icntl_14", 40.0);

Comments

It somehow works but not perfectly, for example when I set
PETScOptions.set("ksp_max_it", 600)
It does not stop the iterations after 600 steps. Is there a reason for it ?