On Jun 23, 2008, at 1:16 PM, Lisandro Dalcin wrote:
On 6/23/08, Thomas Geenen <[email protected]> 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