Upon reading the documentation closer, I think I've found the way to access
u^k (+/-) across each interface
> - I can't seem to be able to access u^k(+) and u^k(-) using a
> FEInterfaceValues fiv, because it's vector-valued
> - this same procedure works fine for scalars. On any face, I can use
> something like fiv.get_face_values(1).get_function_values(scalar_field,
> qpt_vector) to get scalar(-) at the face quadrature points (and use
> get_face_values(0) for the scalar(+) value)
>
> Am I missing something? There must be an easy way to access jumps and
> averages of vector-valued DG data on the right-hand side. Step-47 almost
> does this, but all the vector-valued jumps appear in the bilinear form
> rather than the right-hand side.
>
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"?
Thank you,
Corbin
--
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/789d59fa-2bc9-4dca-a577-504713b29e87n%40googlegroups.com.
