On 09/27/2016 11:18 AM, Timo Heister wrote:
That raises the question if there is ever a valid reason to call
set_boundary_indicator on an interior face
No.
or if we should put an
Assert() into place.
That's difficult because a face itself doesn't know whether it's at the
boundary or not. Only cells know whether they're at the boundary. In other
words, once you write
cell->face(f)->set_boundary(42)
you've lost the information at the first ->. We would need to add a flag to
all faces that says whether they're at the boundary. We could then also write a
face->at_boundary()
function, and do
Assert (!at_boundary(), ...)
in the implementation of set_boundary. But we will have to have that flag first.
Feel free to copy this into a github issue if you think it would be a
worthwhile project.
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.
