Von: firedrake-bounces@imperial.ac.uk [mailto:firedrake-bounces@imperial.ac.uk] Im Auftrag von Miklos Homolya
Gesendet: 30 September 2015 10:14
An: firedrake@imperial.ac.uk
Betreff: Re: [firedrake] Inspecting values of Jacobian matrix

Hello,

I suspect it’s

print assemble(J).M.values

[Buesing, Henrik] Hmm… I get a PETSc error then (see below). Is there sth. similar for “Mat type nest”?

Traceback (most recent call last):

  File "test.py", line 151, in <module>

    print assemble(J).M.values

  File "<string>", line 2, in values

  File "/home/hb111949/Code/PyOP2/pyop2/versioning.py", line 115, in modifies

    retval = method(self, *args, **kwargs)

  File "/home/hb111949/Code/PyOP2/pyop2/petsc_base.py", line 726, in values

    return self.handle[:, :]

  File "PETSc/Mat.pyx", line 218, in petsc4py.PETSc.Mat.__getitem__ (src/petsc4py.PETSc.c:111135)

  File "PETSc/petscmat.pxi", line 962, in petsc4py.PETSc.mat_getitem (src/petsc4py.PETSc.c:30533)

  File "PETSc/petscmat.pxi", line 882, in petsc4py.PETSc.matgetvalues (src/petsc4py.PETSc.c:29285)

petsc4py.PETSc.Error: error code 56

[0] MatGetValues() line 1782 in /rwthfs/rz/cluster/home/hb111949/petsc_opt/src/mat/interface/matrix.c

[0] No support for this operation for this object type

[0] Mat type nest

what you’re looking for.

Regards,

Miklos

On 30 Sep 2015, at 10:06, Buesing, Henrik <HBuesing@eonerc.rwth-aachen.de> wrote:

Dear Firedrakers,

I would like to inspect the values of my Jacobian matrix derived from my form F. I did a

J=derivative(F,u)
print assemble(J).M

But this did not result in the actual values (see below for output). Can I somehow obtain the actual values?
Thank you!

Henrik

Output of print assemble(J).M

OP2 Mat: None_None_None_None_matrix, sparsity (OP2 Sparsity: dsets (MixedDataSet((DataSet(ExtrudedSet(Set((4, 4, 4, 4), 'None_nodes'), 5), (1,), 'None_nodes_dset'), DataSet(ExtrudedSet(Set((4, 4, 4, 4), 'None_nodes'), 5), (1,), 'None_nodes_dset'))), MixedDataSet((DataSet(ExtrudedSet(Set((4, 4, 4, 4), 'None_nodes'), 5), (1,), 'None_nodes_dset'), DataSet(ExtrudedSet(Set((4, 4, 4, 4), 'None_nodes'), 5), (1,), 'None_nodes_dset')))), rmaps (MixedMap((DecoratedMap(Map(ExtrudedSet(Set((0, 0, 0, 0), 'plexmesh_interior_facets'), 5), ExtrudedSet(Set((4, 4, 4, 4), 'None_nodes'), 5), 2, None, 'None_interior_facet_node'), frozenset(['ALL']), frozenset([]), None), DecoratedMap(Map(ExtrudedSet(Set((0, 0, 0, 0), 'plexmesh_interior_facets'), 5), ExtrudedSet(Set((4, 4, 4, 4), 'None_nodes'), 5), 2, None, 'None_interior_facet_node'), frozenset(['ALL']), frozenset([]), None))), MixedMap((DecoratedMap(Map(ExtrudedSet(Set((1, 1, 1, 1), 'plexmesh_cells'), 5), ExtrudedSet(Set((4, 4, 4, 4), 'None_nodes'), 5), 1, None, 'None_cell_node'), frozenset(['ALL', 'ON_INTERIOR_FACETS']), frozenset([]), None), DecoratedMap(Map(ExtrudedSet(Set((1, 1, 1, 1), 'plexmesh_cells'), 5), ExtrudedSet(Set((4, 4, 4, 4), 'None_nodes'), 5), 1, None, 'None_cell_node'), frozenset(['ALL', 'ON_INTERIOR_FACETS']), frozenset([]), None)))), cmaps (MixedMap((DecoratedMap(Map(ExtrudedSet(Set((0, 0, 0, 0), 'plexmesh_interior_facets'), 5), ExtrudedSet(Set((4, 4, 4, 4), 'None_nodes'), 5), 2, None, 'None_interior_facet_node'), frozenset(['ALL']), frozenset([]), None), DecoratedMap(Map(ExtrudedSet(Set((0, 0, 0, 0), 'plexmesh_interior_facets'), 5), ExtrudedSet(Set((4, 4, 4, 4), 'None_nodes'), 5), 2, None, 'None_interior_facet_node'), frozenset(['ALL']), frozenset([]), None))), MixedMap((DecoratedMap(Map(ExtrudedSet(Set((1, 1, 1, 1), 'plexmesh_cells'), 5), ExtrudedSet(Set((4, 4, 4, 4), 'None_nodes'), 5), 1, None, 'None_cell_node'), frozenset(['ALL', 'ON_INTERIOR_FACETS']), frozenset([]), None), DecoratedMap(Map(ExtrudedSet(Set((1, 1, 1, 1), 'plexmesh_cells'), 5), ExtrudedSet(Set((4, 4, 4, 4), 'None_nodes'), 5), 1, None, 'None_cell_node'), frozenset(['ALL', 'ON_INTERIOR_FACETS']), frozenset([]), None)))), name None_None_None_None_sparsity), datatype float64

--

Dipl.-Math. Henrik Büsing

Institute for Applied Geophysics and Geothermal Energy

E.ON Energy Research Center

RWTH Aachen University

------------------------------------------------------

Mathieustr. 10            |    Tel +49 (0)241 80 49907

52074 Aachen, Germany     |    Fax +49 (0)241 80 49889

------------------------------------------------------

http://www.eonerc.rwth-aachen.de/GGE

hbuesing@eonerc.rwth-aachen.de

------------------------------------------------------

_______________________________________________
firedrake mailing list
firedrake@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake