On 12/07/2016 10:50 PM, RAJAT ARORA wrote:
I am using deal.ii to solve a 3D solid mechanics problem which uses
Parallel::Distributed::Triangulation.
I have a question related to setting the boundary indicators on the
triangulation when no automatic partitioning is allowed.
To set the boundary ids,
1) Should I only loop over active_cells (triangulation.begin_active() to
triangulation.end()) or all cells (triangulation.begin() to
triangulation.end()) ?
2) Should the boundary ids be set only for locally_owned_cells or for all
cells (locally_owned + ghost ) ?
You need to set boundary_ids for at least locally owned and ghost cells. But
there is no harm setting them for all cells, so I might just go with 1/ above.
It's easiest to set them for *all* cells (not just the active ones) because
then boundary ids can be inherited upon mesh refinement from parent face to
child face.
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.
For more options, visit https://groups.google.com/d/optout.
