On Aug 7, 2008, at 2:24 PM, Lisandro Dalcin wrote:
I've started to wrap the functionality of the new code in petsc4py. However, I'm detecting some inconsistencies and possible improvements/simplifications and then I need some directions in order to push good fixes in petsc-dev.
Then, here is my tow first questions.
1) Suppose I do
MatGetFactor(A, solver_package, factor_type, &F)
Then, if the call is successful, the following should happen
a) F->factor should be equal to 'factor_type'
Yes.
b) F->preallocated sould be PETSC_TRUE
I guess this is ok.
c) F->assembled sould be PETSC_TRUE
Hmm, not sure, should we wait until after the numerical factorization to mark this as true?
Is this OK?
2) What should be the output of MatView(F, viewer) ?
ASCII info should print matrix size, number nonzeros and solver type. Otherwise I am fine with not printing anything and generating an unsupported viewer type or format.
-- 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