On 10/21/19 2:55 AM, Mauro Murer wrote:
>
> I tried to take constraints into account in the sparsity pattern, but I don't
> know how to do it for a rectangular matrix such in my case.
> I have rows equal to the number of velocity dofs and columns equal to
> pressure
> dofs.
> I didn't found any function for sparsity pattern in the documentation that
> take in account different constraints (in my case velocity and pressure
> constraints).
Then just build the sparsity pattern as you do and use
AffineConstraints::compress() after the fact -- I believe that there is such a
function that can work on rectangular matrices.
In general, you are discussing why it is simpler to just keep all components
of the solution in the same DoFHandler object. In that case, the rectangular
matrix is just one block in a square matrix, and all of the constraint
handling things work as expected.
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/146c2446-0d06-1cec-a908-0b474e4747dd%40colostate.edu.
