On 10/17/19 9:45 PM, Sam wrote:
> I don't really need to evaluate the function at any interior point. I
> eventually wanted to project the boundary data obtained from the
> solution to another FE_FaceQ of different degree defined on a
> different triangulation. So to pass a boundary_function map to
> vectortools::project_boundary_values, I was planning to use the
> FEFieldfunction to first define a function from solution vector
> corresponding to dof_handler_faceQ. I assume it won't work because
> FEFieldFunction uses Fevalues instead of FeFaceValues at quadrature
> points.
Yes. And quadrature points may not exactly lie on faces either.
> How else can I pass the boundary data obtained from FE_FaceQ
> to vectortools::project_boundary_values.
You can't. But you can take a look at how that function is implemented
and then copy it into your code. In all of the places where you are
evaluating the function via FEFieldFunction, you'll just have to
directly evaluate your existing function using FEFaceValues. It's really
no different than what you do in many other contexts where an existing
solution enters into the next equation you want to solve (e.g., in step-15).
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/fb36e926-a34b-5fd7-bae3-ea8d1d5cd66e%40colostate.edu.
