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
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
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
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
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