I have a vector-valued solution, call it H, on a codimension-one geometry
(described by a MappingQEulerian object) and I am very interested in
obtaining H*nu, the dot product of H with the normal vector nu.
I have already obtained H*H using a class derived from
DataPostProcessorScalar, and that makes sense to me - the function
signature of compute_derived_quantities_vector() looks like this:
template <int spacedim>
void
VectorValuedSolutionSquared<spacedim>::compute_derived_quantities_vector
(const std::vector<Vector<double> > &uh,
const
std::vector<std::vector<Tensor<1, spacedim> > > & /*duh*/,
const
std::vector<std::vector<Tensor<2, spacedim> > > & /*dduh*/,
const
std::vector<Point<spacedim> > & /*normals*/,
const
std::vector<Point<spacedim> > & /*evaluation_points*/,
std::vector<Vector<double> >
&computed_quantities) const
since H*H doesn't require anything other than uh, but now I would like to
take advantage of that vector of normals (which, for dealii 8.44.1, should
be of type const std::vector<Tensor<1,spacedim> >, see here - but I digress)
At any rate, I don't know how to access this normals vector, and if I were
to compute it ahead of time using fe_values.get_all_normals() then I am not
sure how to pass it to compute_derived_quantities(), since my code never
actually calls that function.
At the end of the day I think this comes down to me fighting with a
solution H of type Vector<double> of size dof_handler.n_dofs(), where that
dof_handler was built with vector-valued finite elements, and now I want to
look at essentially a scalar function on my geometry.
Any help would be great,
Tom
--
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.
