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. 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 -- 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.
