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.