Re: [deal.II] discretization of a analytical function on a second triangulation

2022-07-06 Thread Wolfgang Bangerth
On 7/6/22 01:02, Simon Wiesheier wrote: Basically, my motivation to call VectorTools::interpolate is to compute the nodal values (scalar field --> DOFs) in the triangulation based on the formula implemented in the just derived class. During the assembly routine, I access some interpolated valu

Re: [deal.II] discretization of a analytical function on a second triangulation

2022-07-06 Thread Daniel Arndt
Simon, if you get all information you need using VectorTools::interpolate then this sounds good. Note that there also is FEFieldFunction( https://www.dealii.org/current/doxygen/deal.II/classFunctions_1_1FEFieldFunction.html ) that implements the interpolation for you. Best, Daniel On Wed, Jul

Re: [deal.II] discretization of a analytical function on a second triangulation

2022-07-06 Thread Simon Wiesheier
Dear Daniel, yes, I defined a class derived from deal::Function and it seems to work for me. Basically, my motivation to call VectorTools::interpolate is to compute the nodal values (scalar field --> DOFs) in the triangulation based on the formula implemented in the just derived class. During th

Re: [deal.II] discretization of a analytical function on a second triangulation

2022-07-05 Thread Daniel Arndt
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_f

[deal.II] discretization of a analytical function on a second triangulation

2022-07-04 Thread Simon
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-dimensi