On 11/08/2016 05:42 PM, Hamed Babaei wrote:
Since I use Newton method and I want to apply boundary displacement as
constraint, I should apply constraint only in zero Newton iteration. So I can
not do so in setup_system but after assembly. I followed step-44 for that.
Then keep two ConstraintMatrix objects around: one that only contains hanging
node constraints, and the other that has both hanging node constraints and
boundary values. You would use the first for creating the sparsity pattern,
and the second to distribute local contributions into the global objects in
assembly.
That said, you *do* need to apply boundary values in subsequent Newton
iterations. They just have to be zero displacement. Take a look at step-15.
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.
