On 2/11/25 03:17, Junxiang Wang wrote:
Therefore, Can I interpret the issue I encounter in that it is troublesome to
distribute the local_rhs when the hanging node is located on the shared edge
of two elements and also the two elements belong to different processors
respectively?
Junxiang:
I don't know what the concrete problem in your code is, but in general deal.II
can deal with hanging nodes just fine. You just need to make sure that you
distribute local contributions into the global matrix and vector objects using
the AffineConstraints::distribute_local_to_global() functions. step-40 is the
canonical place to see how that is done. I would try and follow the same
structure as shown in step-40 in your code as close as possible.
My suggestion for you is to see if you can simplify the situation (say, use a
right hand side f=1 for simplicity), output the local contributions to make
sure you get the same with one or multiple processes, and output the completed
right hand side vector (again to make sure you get the same with one and
multiple processes).
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 visit
https://groups.google.com/d/msgid/dealii/6a9ac2a8-3daa-4bbf-8073-aed5c6b2fa0d%40colostate.edu.
