-----Ursprüngliche Nachricht----- Von: firedrake-bounces@imperial.ac.uk [mailto:firedrake- bounces@imperial.ac.uk] Im Auftrag von Lawrence Mitchell Gesendet: 30 September 2015 10:31 An: firedrake@imperial.ac.uk Betreff: Re: [firedrake] Inspecting values of Jacobian matrix
-----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. [Buesing, Henrik] Great! That worked. Unfortunately the values are wrong... I will post another question concerning this. Thank you!
Henrik
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-----
_______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake