On 7/5/22 09:28, Simon Wiesheier wrote:
What is the approach to achieve this?
My only idea is to create a new FEValues object for each quadrature point and
pass an appropriate Quadrature object. But seems to be very inefficient.
This is basically what the VectorTools::point_value() and
VectorTools::point_gradient() functions do. You might want to read through
their implementations to see how this is done.
The only other approach is with classes like FEPointEvaluation. If you read
through their implementation, you might be able to understand how they can
also be used to evaluate just a single shape function. In the end, the values
and gradients of a single shape function correspond to the values and
gradients of a finite element field that corresponds to a solution vector with
one 1.0 and the rest of the elements set to zero.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
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 dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/6527dccb-aebc-f424-5217-1c19a9b8b3c2%40colostate.edu.
