-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 On 30/09/15 09:19, Buesing, Henrik wrote:
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"?
You can look at the values of the individual blocks, corresponding to the subspaces in your mixed space, like so: M = assemble(J).M print M[i, j].values where i and j are the block row and column indices. Lawrence -----BEGIN PGP SIGNATURE----- Version: GnuPG v1 iQEcBAEBAgAGBQJWC53IAAoJECOc1kQ8PEYvrUQIAL2cbn+iDOcC1KItC07ZuI93 UuzrbHBDoGqXZLhhcdddJkk3XVqDr47xuu4UFOSBLyEROimi/X5lyLiE+h+d48yg /OtSp0r2+iwUUQcxKhqKSumWaLr62qc71NV61nTcwHW/iJAM4FhxUlmCfe+lpht6 01UkufcgKwt7qBomXMXgQjFkV7qXp0tIaS+YlDLaRelD4xV3/paWmjNhzNrZxtqx KUWC82bUhakv3M2sCdsHiBXtiHvOrPaC3LgFBQHhqnKk83Nxia/ols4JE1+LzBTH 3WUuLzU+Cw8ERB/vQkeMcBb0XYTKANWqmiCQZYcUW1W2KHLhBL+LIZP1yoMI+io= =8MAd -----END PGP SIGNATURE-----