Re: [Getfem-users] How to get a solution at given point

[Roman Putanowicz]

> Hello,
> 
> thanks for posting the code for obtaining the value of solution in a
> given point.
> I am going to give it a try.
> Nevertheless, I would like to ask, if there is some way how to get a 
> derivative
> of the solution in a given point. Just for a verification purposes, I need to
> model single tetrahedron (quadratic) and to get the solution and its
> derivative in
> arbitrary point inside the tetra.
In the code I have posted the method eval of getfem::interpolator_on_mesh_fem
evaluates both the value of the field and its gradiet (the third argument).
> >  getfem::interpolator_on_mesh_fem itpr(mf_T, T);
> >  getfem::base_vector solA;
> >  getfem::base_matrix gradA;
> >  itpr.eval(A, solA, gradA);

Regards,

Roman
-- 
Roman Putanowicz, PhD  < address@hidden  >
Institute for Computational Civil Engng (L-5)
Dept. of Civil Engng, Cracow Univ. of Technology
www.l5.pk.edu.pl, tel. +48 12 628 2569, fax 2034