Re: [Getfem-users] periodic boundary conditions example and source term question

[Kajetan Sikorski]

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