Felix,

also, if the point you want is not a quadrature point, you can exploit the 
following trick:

Point<dim> my_point_in_reference_cell;

std::vector<Point<dim> > quad(1, my_point_in_reference_cell);

and then you can initialize an 

FEValues<dim> fev(fe, quad, update_inverse_jacobians);

This will give you access to what you need by


fev.reinit(cell);
auto &jac_invs = fev.get_inverse_jacobians();

Best,
Luca.

> On 14 Jun 2017, at 2:20, Weixiong Zheng <[email protected]> wrote:
> 
> Felix,
> 
> Would you take a look at if [1] is what you want to have?
> 
> [1]: 
> https://www.dealii.org/8.4.0/doxygen/deal.II/classFEValuesBase.html#a665ab681011424cb1b0ba11c4155a121
> 
> Best,
> Weixiong
> 
> 在 2017年6月13日星期二 UTC-7上午9:35:09,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:
> 
> 
> 
> 
> 
> 
> 
> How do I get the inverse Jacobian at this point?
> 
> 
> 
> 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 
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an 
> email to [email protected].
> For more options, visit https://groups.google.com/d/optout.

-- 
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 [email protected].
For more options, visit https://groups.google.com/d/optout.

Reply via email to