Re: [deal.II] shape function gradient at arbitrary points in the element

2022-07-07 Thread Wolfgang Bangerth
On 7/7/22 02:34, Simon Wiesheier wrote: Just to make sure that I chose a reasonable approach: //create FEValues object whenever I need it, that is, 'number_of_cells*number_of_qps_per_cell' times FEValues<2> fe_values(dof_handler.get_fe(),    Quadrature<2>(m

Re: [deal.II] shape function gradient at arbitrary points in the element

2022-07-07 Thread Simon Wiesheier
Just to make sure that I chose a reasonable approach: //create FEValues object whenever I need it, that is, 'number_of_cells*number_of_qps_per_cell' times FEValues<2> fe_values(dof_handler.get_fe(), Quadrature<2>(my_point_in_ref_coords), update_gradients); f

Re: [deal.II] shape function gradient at arbitrary points in the element

2022-07-06 Thread Wolfgang Bangerth
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 Vecto

Re: [deal.II] shape function gradient at arbitrary points in the element

2022-07-05 Thread Simon Wiesheier
Dear Martin, thanks for pointing out to the FEPointEvaluation class. However, my intent is to compute the gradients of shape functions in real coordinates. " so in that case the recommendation would be to use FEValues, despite that being terribly inefficient." What is the approach to achieve thi

Re: [deal.II] shape function gradient at arbitrary points in the element

2022-07-05 Thread Martin Kronbichler
Dear Simon, We have the class FEPointEvalation, https://dealii.org/developer/doxygen/deal.II/classFEPointEvaluation.html , which implements the operation of evaluating a solution at arbitrary points (as handed in as an array to points in unit coordinates). To do this for positions in real coo

[deal.II] shape function gradient at arbitrary points in the element

2022-07-05 Thread Simon
Dear all: I have to compute the gradient with respect to real space co-ordinates of the shape functions belonging to a FE_Q element at arbitrary points in the element, that is, not only at quadrature points. What is the way to do this in dealii? FEValues provides shape function values, gradien