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
>