Simon, defining a class derived from dealii::Function sounds sensible. Whether calling VectorTools::interpolate makes sense, depends on what you want to do with the interpolated function. In may cases, it's enough to just evaluate it locally during assembly using functionality like FEValues::get_function_values.
Best, Daniel On Mon, Jul 4, 2022 at 10:26 AM Simon <simon.wieshe...@gmail.com> wrote: > Dear all: > > > I am solving a nonlinear PDE on a Triangulation T_1. > > I know the analytical representation of a scalar function of two variables > (no space co-ordinates, but two invariants of a quantity) > and my goal is to find a discretized version of the analytical function > (on a second, two-dimensional triangulation T_2). > The nodal values of T_2 should be the exact values corresponding to the > analytical solution. > As for the values between the nodes, I want to start with a bilinear > interpolation. > > The coupling to my original PDE works as follows: > For each quadrature point defined on T_1, I compute the two input > variables of the above function based on my solution vector, find the > corresponding cell of T_2 and want to retrieve the interpolated nodal > values. > > My idea is to use the function > VectorTools::interpolate , > but I do not know a suitable Function object to hand over. > > Is this an appropriate approach at all? > > > Thank you, > Simon > > -- > 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/14b1d8f6-f902-4414-a5c6-348ae3f5003en%40googlegroups.com > <https://groups.google.com/d/msgid/dealii/14b1d8f6-f902-4414-a5c6-348ae3f5003en%40googlegroups.com?utm_medium=email&utm_source=footer> > . > -- 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/CAOYDWbKNTV%3D8Bipyq_gzbMoLvxdvv-%2BmHLipn27DbOcK-eWL_A%40mail.gmail.com.
