Corbin, if you only ever need the DoF mapping periodic DoFTools::make_periodicity_constrainst ( https://www.dealii.org/current/doxygen/deal.II/namespaceDoFTools.html#a03324e6e060c6a55191b2c60ad871eab) sound like the right tool for you to use. Just create a separate AffineConstraints object to be filled and hand it to that function. Then, you can loop through the filled object to either use it as a map directly or you can create a separate data structure. Note that the constraints might not be in the order (bottom -> top) you need them to be. So you might need to adjust that yourself.
Best, Daniel Am Do., 29. Apr. 2021 um 04:48 Uhr schrieb Corbin Foucart < [email protected]>: > Hi all, > > I'm looking to map DoF from the bottom of an extruded volumetric mesh to > the corresponding DOF on the surface of the mesh. Right now, I'm saving the > cell "number"/ID of each cell in the bottom layer of the volume mesh, and > want to match it with a surface cell by comparing their centers > > - I could loop over all the volume cells until I find the one matching > the surface cell center, but as I have the cell numbers, I would rather > access the cell center from the cell id (which looks to me like an unsigned > int) > - something like get_cell(cell_id)->center --- is this possible? > - Alternatively, I was thinking that a periodic boundary condition > linking the bottom of the volume mesh with the top might be able to achieve > the same thing, but I wasn't sure how to go about setting it up > > Any pointers would be appreciated! > Corbin > > -- > 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/122e1b62-e672-446b-8c25-ce6012c1fe3fn%40googlegroups.com > <https://groups.google.com/d/msgid/dealii/122e1b62-e672-446b-8c25-ce6012c1fe3fn%40googlegroups.com?utm_medium=email&utm_source=footer> > . > -- 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/CAOYDWbJvkw_0yW2KtsMntKMt7eSEYBUq_jJOoP2EvSjWT5O%2B4w%40mail.gmail.com.
