Hi Wolfgang, the function achieves the desired results. Thanks a lot!
Cheers, Jose On Tuesday, March 21, 2023 at 9:10:12 PM UTC+1 Wolfgang Bangerth wrote: > > 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/95846f08-c4a0-477a-a8d0-fb38f37f6c78n%40googlegroups.com.
