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

[Roman Putanowicz]

> Dear All, 
> 
> Apparently a simple question but I could not find a simple answer for it:
> I would like to evaluate solution to my problem a given point in space.

Answering my own post: the class that solves the above problem is
getfem::interpolator_on_mesh_fem.

I have overlooked it in the reference documentation. 

Its usage is roughly equivalent to:

 getfem::base_node A; bgeot::sc(A)=0.01,0.01;
 plain_vector T;
 gmm::resize(T, mf_T.nb_dof());
 gmm::copy(model.real_variable("T"), T);
 getfem::interpolator_on_mesh_fem itpr(mf_T, T);
 getfem::base_vector solA;
 getfem::base_matrix gradA;
 itpr.eval(A, solA, gradA); 

> I tired to use some slicer object but cannot figure out the pattern of its 
> usage. I would appreciate any hint.

The question still remains about slicers. Is it possible to use them for 
this purpose (especially if we have a cloud of points instead of a
single point)?

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