Hello, I suspect it’s print assemble(J).M.values 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 <http://www.eonerc.rwth-aachen.de/GGE> hbuesing@eonerc.rwth-aachen.de <mailto:hbuesing@eonerc.rwth-aachen.de> ------------------------------------------------------
_______________________________________________ firedrake mailing list firedrake@imperial.ac.uk <mailto:firedrake@imperial.ac.uk> https://mailman.ic.ac.uk/mailman/listinfo/firedrake <https://mailman.ic.ac.uk/mailman/listinfo/firedrake>