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
