Still a variation on the same question. Following the documentation,
interpolation(expr,
*args)Synopsis: V = Model.interpolation(self, string expr, {MeshFem mf | MeshImd mimd | vec pts, Mesh m}[, int region[, int extrapolation[, int rg_source]]])
I tried interpolation on MeshImd:
#-------------------------------------------------------------------------------------
degp = 2; NX = 30
tmesh = gf.Mesh('regular_simplices', np.arange(0,1+1./NX,1./NX),
np.arange(0,1+1./NX,1./NX))
tmeshfem = gf.MeshFem(tmesh, 1)
tmeshfem.set_fem(gf.Fem('FEM_PK(2,' + str(degp) + ')'))
tmeshim = gf.MeshIm(tmesh, gf.Integ('IM_TRIANGLE(6)'))
U = np.random.rand(tmeshfem.nb_basic_dof())
md = gf.Model('real')
md.add_fem_variable('u', tmeshfem)
md.to_variables(U)
GU = md.interpolation('Grad_u', tmeshim)
#-------------------------------------------------------------------------------------
and get
Traceback (most recent call last):
File "testinterpgetfem.py", line 21, in <module>
GU = md.interpolation('Grad_u', tmeshim) #pts.T, tmeshi)
File ".../Conda/deps/usr/lib/python2.7/site-packages/getfem/getfem.py", line 2800, in interpolation
return self.get("interpolation", expr, *args)
File ".../Conda/deps/usr/lib/python2.7/site-packages/getfem/getfem.py", line 2708, in get
return getfem('model_get',self.id, *args)
RuntimeError: (Getfem::InterfaceError) -- Argument 4 should be a DOUBLE REAL data array
It is not the right way to call it perhaps..?
Sorry for the noise,
best regards, Edouard.
--
Edouard Oudet : http://www-ljk.imag.fr/membres/Edouard.Oudet/
IMAG - Bureau 164
700 avenue Centrale
38400 Saint Martin d'Hères
+33 (0)4 57 42 17 71 (office LJK)
+33 (0)4 79 68 82 06 (home)