using MUMPS with mg

[Lisandro Dalcin]

On 6/23/08, Thomas Geenen <geenen at gmail.com> 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:

MatConvert(A, MATSAME, MAT_REUSE_MATRIX,&A)

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