Once again, apologies for asking questions which are more oriented towards vector calculus than the deal.II library/code itself.
In Step-20, under the introduction section on solvers, it says: "If we sort our degrees of freedom so that all velocity come before all pressure variables, then we can subdivide the linear system Ax=b into the following blocks:" (MBTB0)(UP)=(FG), I am a bit confused about the origins of this block-matrix, particular the "B" term in the top-right block (and to a lesser extent to its transpose term in the bottom-left block). I shall explain more below. *The confusing aspects* Recalling that the above system matrix is the discretized version of the original (continuous) PDEs, this corresponds to the assembly of the bilinear form: (v,Kâ1u) â(div v,p) â(q,div u) - Now, the top-left block M of the matrix operator is clear; We are multiplying two vector shape functions which leads to the so-called "mass matrix". - But the entry in the top-right block, B is not clear. This is just a product of a scalar shape function and the negative divergence of a vector shape function. *How is this operator, the gradient?* - Yes, the bottom-left block is actually the same as #2 above, right? ie. the *product* of a scalar shape function and the negative divergence of a vector shape function? *Why is this entry indicated as just the -div operator?* Apologies for asking all these vector-calculus questions in this forum. Regards, Krishna -- 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/32c21917-6e37-43a2-b47b-63b139c6a920%40googlegroups.com.