On 6/10/21 2:07 AM, Corbin Foucart wrote:
Indeed,
std::vector<Vector<double>> velocity_qpt_vals(n_face_qpts,
Vector<double>(dim));
fiv.get_fe_face_values(1).get_function_values(velocity_field,
velocity_qpt_vals);
accesses the values; the issue was I hadn't explicitly instantiated the
Vector<double>s with length dim, which was causing a segfault. However, I'm
curious -- is there a better way to compute the jumps and averages of DG data
other than requesting them via fe_face_values and computing them "by hand"?
Not right now. The FEInterfaceValues class should eventually acquire functions
get_function_jumps
get_function_averages
get_function_gradient_jumps
get_function_gradient_averages
etc, similar to the FEValues::get_function_values function and friends. But
nobody has written these functions yet. It would of course be fantastic if you
could give this a try -- it would solve your problem, and would be very useful
to others as well! We'd love walking you through the steps to make that happen!
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
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 [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/11001043-6a0f-de69-2b46-48fa0ef62817%40colostate.edu.
