On 10/17/19 6:41 PM, sam the dude wrote:
> Thanks for the response. I think FE_FaceQ is exactly what I am looking for,
> at least in the case of representing the normal trace of RT_0.
> But I have trouble using the FEFieldFunction with the FE_FaceQ element. I
> tried the following:
>
> /Functions::FEFieldFunction<dim,DoFHandler<dim>,Vector<double>>
> fe_interface_data(dof_handler_face,solution_vector);
> /
> / MappingQGeneric<dim> mapping_generic(1);/
> /
> DoFTools::map_dofs_to_support_points(mapping_generic,dof_handler_face,support_points);
> fe_interface_data.value_list(support_points,values);/
>
> And got all values to be NaN. Am I doing something wrong?
No -- a finite element that only lives on the faces simply can't be evaluated
in the interior of cells (i.e., at arbitrary points). If that's what you need
to do, then you ought to rethink how exactly your shape functions are defined.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
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 dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/3137b188-ad69-a57a-ac87-0893d7cdca02%40colostate.edu.
