I've tried to use the interpolator function without
success; here what I've done:
getfem::mesh_fem mf_p0_3;
getfem::pfem pf_p0;
pf_p0 = getfem::fem_descriptor("FEM_PK(3,0)");
mf_p0_3.init_with_mesh(mesh);
mf_p0_3.set_qdim(N);
mf_p0_3.set_finite_element(mesh.convex_index(), pf_p0);
std::vector<scalar_type> U_v_int;
U_v_int.resize(mf_p0_3.nb_basic_dof());
getfem::interpolation(mf_v,mf_p0_3,U_v,U_v_int);
mf_v is the RT0 fem and the computed velocity in RT0 is
stored in U_v. I have the following error:
terminate called after throwing an instance of
'gmm::gmm_error'
what(): Error in
/usr/local/include/getfem/getfem_fem.h, line 728 :
Wrong size for coeff vector
Thank you for the availability,
Calogero B. Rizzo