an small and perhaps unimportant gotcha
Lisandro Dalcin
dalcinl at gmail.com
Tue Dec 2 09:12:27 CST 2008
Other (really nice!) change I've noticed is that now we can call
{Vec|Mat}Create() and next {Vec|Mat}SetType() and it works!
However, see this.
In [1]: from petsc4py import PETSc
In [2]: x = PETSc.Vec().create()
In [3]: x.setType('seq')
In [4]: print x.getType()
None
In [5]: A = PETSc.Mat().create()
In [6]: A.setType('seqaij')
In [7]: print A.getType()
None
Internally, the "type_name" field is never set, it still is NULL (then
petsc4y returns None).
The problem is that forcing the "type_name" to be set will be really
dangerous (rmember PetscValidType macro).
Any comments?
--
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