This is strange. I just tried KSP ex2 and after MatSetType() the type_name is<br>set. What does *.setType() call?<br><br>  Matt<br><br><div class="gmail_quote">On Tue, Dec 2, 2008 at 9:12 AM, Lisandro Dalcin <span dir="ltr"><<a href="mailto:dalcinl@gmail.com">dalcinl@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Other (really nice!) change I've noticed is that now we can call<br>
{Vec|Mat}Create() and next {Vec|Mat}SetType() and it works!<br>
<br>
However, see this.<br>
<br>
In [1]: from petsc4py import PETSc<br>
<br>
In [2]: x = PETSc.Vec().create()<br>
<br>
In [3]: x.setType('seq')<br>
<br>
In [4]: print x.getType()<br>
<br>
None<br>
<br>
In [5]: A = PETSc.Mat().create()<br>
<br>
In [6]: A.setType('seqaij')<br>
<br>
In [7]: print A.getType()<br>
<br>
None<br>
<br>
<br>
Internally, the "type_name" field is never set, it still is NULL (then<br>
petsc4y returns None).<br>
The problem is that forcing the "type_name" to be set will be really<br>
dangerous (rmember PetscValidType macro).<br>
Any comments?<br>
<font color="#888888"><br>
<br>
<br>
--<br>
Lisandro Dalcín<br>
---------------<br>
Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC)<br>
Instituto de Desarrollo Tecnológico para la Industria Química (INTEC)<br>
Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)<br>
PTLC - Güemes 3450, (3000) Santa Fe, Argentina<br>
Tel/Fax: +54-(0)342-451.1594<br>
<br>
</font></blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>