Hello everybody, is there a way to get the Jacobian of a mapping (or its determinant) for an arbitrary point in a cell? I know that usually this quantity is can be obtained via the FEValues Object for the quadrature points. But I need the Jacobian at different points in each cell.
I see two ways that might work, but both are a bit annoying: - Wrap the points that are required into a (Pseudo)Quadrature. As these points are different for each cell, this requires a new FEValues object for each cell. A variation might be to use the hp variant of FEValues and use a collection of quadratures, one for each cell. This would allow to at least recycle the FEValues object. - Calculate the Jacobian on my own from the information I have available (probably Mapping::support_point_value is everything I need). But I do not want to reinvent the wheel, Mapping can already calculate that. Is there another, better way that I missed? If not, which way do you think is better? Thank you very much in advance Johannes Reinhardt _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
