Hi Simon,
This could be helpful:
https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#how-to-get-the-mapped-position-of-support-points-of-my-element
Here's a code snippet for Option-1, which I usually use.
const auto & feSystem = dofHandler.get_fe();
dealii::MappingQ<dim, dim> cellMapping(feSystem.degree);
const auto & referenceCellSupportPoints =
feSystem.get_unit_support_points();
std::vector<dealii::Point<dim>>
realCellSupportPoints(referenceCellSupportPoints.size());
unsigned int supportPointCtr = 0;
for (const auto &supportPoint : referenceCellSupportPoints)
{
realCellSupportPoints[supportPointCtr] =
cellMapping.transform_unit_to_real_cell(cellIter, supportPoint);
++supportPointCtr;
}
Best regards,
Paras
On Mon, Apr 19, 2021 at 12:14 PM Simon <simon.wieshe...@gmail.com> wrote:
> Dear All,
>
> I am implementing a local postprocessing SPR-Approach for a mechanical
> problem, in which I sloppy speaking need the global coordinates of all
> nodes.
>
> I use Elements of type FE_Q<dim>(p).
> For p=1 there is no problem: I can simply loop over all cells and vertices
> and use cell->vertex(...) in order to get the global coordinates of
> vertex(...).
> Being p=2 for instance, this approach does not work anymore. In this case
> is there an easy way to figure out the coordinates of midnodes as well?
>
> Thanks for help!
>
> Best
> Simon
>
> --
> 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/3d87b2c4-61a4-4125-982d-f7341f44492an%40googlegroups.com
> <https://groups.google.com/d/msgid/dealii/3d87b2c4-61a4-4125-982d-f7341f44492an%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/CAEU6zmSQ-9squp34HNynpjsc02VcdpNxqiNPAq%2Bq4LQTmgwckQ%40mail.gmail.com.
