MATSHELL a bit broken...

Lisandro Dalcin dalcinl at gmail.com
Thu Oct 25 16:20:48 CDT 2007


It seems MATSHELL implementation is a bit broken. If I'm not wrong, if
a user set the operations MOP_ASSEMBLY_END, then the use of the
automatic MatScale/MatShift would lead to unexpected result...

Should this operation be managed the same way as
mult/multtranspose/get diagonal currently is, that is, by backupping
the user provided function pointer?

BTW, sorry if I disturb you asking for making changes that seems
trivial, but I really do not like to change any line of code after
knowing your opinion.

-- 
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




More information about the petsc-dev mailing list