Alberto,
my suspicion would be that you are computing terms for the interface matrix
that are not part of the sparsity pattern of the Trilinos matrix. (PETSc does
not require you to say up front where the nonzero entries of the matrix are
going to be, if I recall correctly.)
This can happen because for "normal" CG methods, you only get entries in the
matrix if DoFs i,j are located on the same cell. But for methods that contain
jump terms on the interfaces, you also get contributions for DoFs i,j located
on different cells (and not necessarily at the face you're integrating over).
Take a look at how step-47 builds its sparsity pattern.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/138efd45-e249-2dbf-af90-b4e2edb4e40f%40colostate.edu.