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.

Reply via email to