first comments after reviewing and trying the new matrix factorization code
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' b) F->preallocated sould be PETSC_TRUE c) F->assembled sould be PETSC_TRUE Is this OK? 2) What should be the output of MatView(F, viewer) ? -- 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
participants (1)
-
Lisandro Dalcin