Dear Nicolò,
The data has to be a vector field on
the domain. The brick is equivalent to the term (which can be
added equivalently with 'add_linear_term')
'(V.Normal).Test_u'
For 'V' the data and 'u' the vector
variable. So the dimension is 3xN with N the number of degrees of
freedom of a scalar field. If you use a vector field, this should
be just N, a priori.
Best regards,
Yves
Le 22/10/2019 à 11:48, Nicolò Cambiaso
a écrit :
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 <address@hidden> 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
-- Yves Renard (address@hidden) tel : (33) 04.72.43.87.08 INSA-Lyon 20, rue Albert Einstein 69621 Villeurbanne Cedex, FRANCE http://math.univ-lyon1.fr/~renard ---------