On 09/29/2017 06:43 AM, [email protected] wrote:
The system matrix needs to have 3x3 blocks since the block for (n_u, n_u) no
longer corresponds to a vector-valued FE space (RT) but to Qn x Qn. Hence, I
have to manually generate blocks of the form (n_u_x,n_u_x), (n_u_x, n_u_y),
etc. This would mean that I have to rewrite many parts of the Schur complement
solver, which hard codes the 2x2 block structure of the matrix.
The issue is how you compute n_u. As Timo already says, how you block things
is your choice. The problem is just that in step-20, n_u needs to be computed
differently from step-22 because in the former case, the finite element used
for the velocity is not primitive whereas in step-22 it is.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
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 [email protected].
For more options, visit https://groups.google.com/d/optout.
