Jose,
I am working on a problem with several subdomains. At the interface
between them a boundary integral is to be evaluated. I am identifying
the interface by comparing the material_id of neighboring cells (or
their active_fe_index as I am using a different FESystem per subdomain).
In order to speed up the search during assembly, a Vector<float> is
previously filled with 1.0 at the cells where the
material_id/active_fe_index differ. This approach works in serial but in
parallel the material_id() call of a neighbor cell outside the locally
owned subdomain always returns 0 (An assertion is missing here). As
such, not only the interface between subdomains is marked but also the
interface between locally owned subdomains, as shown in the attached picture
If I understand you right, then you want to have the correct material_id
set also on ghost cells? If so, take a look at this function:
https://dealii.org/developer/doxygen/deal.II/namespaceGridTools.html#a9565dbf2f8e45fee28e40806870e2c98
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].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/ae02e863-ffb4-5b4b-6a13-61fbb6b31157%40colostate.edu.
