using MUMPS with mg

[Barry Smith]

On Jun 23, 2008, at 1:16 PM, Lisandro Dalcin wrote:

> 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?

    Even better than this setsolvertype nonsense.

    Barry

>
>
>
>
> -- 
> 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
>