Sorry, I forgot to define you the region, that is a side of the cube,
gf_mesh_set( Mesh_M , 'region', RIGHT_BOUND, fb1);
fb1 = gf_mesh_get( Mesh_M , 'outer faces with direction', [ 1 0 0], 0.01);
Thank you very much,
Have a nice day
Nicolò Cambiaso
Il giorno mar 22 ott 2019 alle ore 11:43 Nicolò Cambiaso <
cambiaso.nic...@gmail.com> ha scritto:
> Dear GetFem users,
>
> I'm trying to apply a normal boundary condition, using the command
> gf_model_set(model M, 'add normal source term brick', mesh_im mim,string
> varname, string dataname, int region),
> using matlab interface.
>
> and I previously defined the components as follows:
> mim = gf_mesh_im( mf, IT);
> IT = gf_integ('IM_TETRAHEDRON(1)')
> mf = gf_mesh_fem(Mesh_M, 3);
> Mesh_M is a cube, [0,1]x [0,1]x [0,1]
> string varname is 'ut', and it was defined as: gf_model_set(md, 'add fem
> variable', 'ut', mf); (md is the model).
>
> but the problem is the "string dataname". The dataname was initialized by
> gf_model_set(md, 'add initialized data', 'dataname', Data);
> and previously Data was defined as a 3xN matrix, where 3 is the dimension
> of the vectors (3D). Hoping right now that everything is correct, what is N
> ??
> Because it's not the number of degrees of freedom of mf, as i thought at
> the beginning.
>
> Thank you very much,
> Have a nice day
>
> Nicolò Cambiaso
>
