Gabriel Peters
Endenicher Str. 310
53121 Bonn
00491525/5478185
gabriel.pet...@koeln.de
Am 05.04.19 um 15:55 schrieb Daniel Arndt
> Gabriel,
>
> to extend on my last post. Finally, I think that your approach with the two
> DoFHandler objects will give you some more problems.
> What is likely happening here is that the AffineConstraints object here is
> used both for rows and columns although row dofs and column dofs are
> unrelated.
> What you wanted here is to shift the pressure dofs such that they don't
> overlap with the velocoty dofs. However, you can't just do this using a
> renumbering since it is currently not allowed.
That would habe Been my question to your First post, how to apply 2 constraints
in the „constraints.add... function, but your described Problem seems Even
Harder to solve.
> Even if you get this to work here without shifting, you will have the same
> problem when assembling your matrix and applying constraints.
> I agree with David that the better approach here is to use a FESystem that
> accounts both for velocity and pressure and follow step-32.
>
> Best,
> Daniel
>
> --
This Sounds like but Problem.
But thanks a lot again. I will try next week to implement the code using
Blockmatrices as prescribed. I‘ll let you know if it works or if further
questions come up.
Best,
Gabriel
--
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.
For more options, visit https://groups.google.com/d/optout.
