Re: [deal.II] Compute jacobian at point (no quadrature point)

2017-06-14 Thread Felix Lorenz
Thank you very much. This is exactly what I wanted. Best regards, Felix -- 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

Re: [deal.II] Compute jacobian at point (no quadrature point)

2017-06-14 Thread luca.heltai
Felix, also, if the point you want is not a quadrature point, you can exploit the following trick: Point my_point_in_reference_cell; std::vector quad(1, my_point_in_reference_cell); and then you can initialize an FEValues fev(fe, quad, update_inverse_jacobians); This will give you

[deal.II] Compute jacobian at point (no quadrature point)

2017-06-13 Thread Felix Lorenz
Hi everyone, I want to compute the transformation of a gradient at a point (no quadrature point) on the reference cell to the real space. Mathematically formulated I want to compute: