Thank you, Simon and David. Both answers were helpful. I will try implementing it both ways for practice.
Best, Corbin On Thursday, April 29, 2021 at 8:40:48 AM UTC-4 d.arnd...@gmail.com wrote: > 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 < > corbin....@gmail.com>: > >> 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 dealii+un...@googlegroups.com. >> 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 dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/602b1076-c4ec-4451-bd86-9c391b7644d6n%40googlegroups.com.
