using MUMPS with mg

Lisandro Dalcin dalcinl at
Mon Jun 23 13:16:48 CDT 2008

On 6/23/08, Thomas Geenen <geenen at> wrote:

> i will check if just adding
>  ierr = MatConvert(A, MATAIJMUMPS, MAT_REUSE_MATRIX,&A);CHKERRQ(ierr);
>  will do the trick

Indeed. That would do the trick. But perhaps better is to add:


Then, you can pass '-matconvert_type aijmumps' (check for the actual
option name, to busy right now to look at the source) to actually use
MATAIJMUMPS. If you do not pass the option, the the MatConvert() call
is just a non-op .

Perhaps this is a candidate for petsc-dev? Or MatSetSolverType() will
handle this in the near future?

Lisandro Dalcín
Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC)
Instituto de Desarrollo Tecnológico para la Industria Química (INTEC)
Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)
PTLC - Güemes 3450, (3000) Santa Fe, Argentina
Tel/Fax: +54-(0)342-451.1594

More information about the petsc-dev mailing list