Yves,
Thanks a lot. I've beet trying to implement what you suggested but I
just can't seem to figure it out. If I could get some help on getting
this code snippet working I think I could go from there. Right now
it's throwing an error that about vGrad having the wrong number of
arguments. Thanks in advance!
gf_workspace('clear all');
m = gf_mesh('cartesian',linspace(0,1,10),linspace(0,1,10));
femu = gf_mesh_fem(m,2);
femd = gf_mesh_fem(m,2,2);
gf_mesh_fem_set(femu,'fem',gf_fem('FEM_QK(2,2)'));
gf_mesh_fem_set(femd,'fem',gf_fem('FEM_QK(2,2)'));
mim = gf_mesh_im(m, gf_integ('IM_GAUSS_PARALLELEPIPED(2,10)'));
S = gf_mesh_fem_get(femd, 'eval', {'9.*x.^3+cos(y)' '8.*y.^2'
'exp(1).^(x.^2+y.^2)' 'x.*y'});
V = gf_asm('volumic' , 'S=data(qdim(#1),#1);V(#1)+=comp(vGrad(#1))
(i,j,i,j).S(i,j)', mim, femu, S);
Kai
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
