getfem-users
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Getfem-users] assembling a vectorial source


From: Alex Magnes
Subject: [Getfem-users] assembling a vectorial source
Date: Sun, 30 Sep 2007 19:10:17 +0200
User-agent: Icedove 1.5.0.12 (X11/20070607)

hi everybody!

another question concerning the generic assembling procedures. i try to assemble a volumic source term "integral grad(Ni).Ts", where Ni are the basis functions of a scalar fem and Ts is a vectorial source term. so i set up a mesh with a data-fem with "FEM_NEDELEC(3)" elements and a "FEM_PK(3,1)" for setting the scalar basis functions.
then:

   assemble_source_term.push_mi(im);
   assemble_source_term.push_mf(mf_scalar);
   assemble_source_term.push_mf(mf_Ts);
assemble_source_term.push_data(Ts); <- input data, length=dof of the nedelec fem.
   assemble_source_term.push_vec(result);
assemble_source_term.set("F=data(qdim(#1),#2);" " V(#1)+=comp(Grad(#1).vBase(#2))(???????).F(???);");
   assemble_source_term.assembly();

my problem is, i do not really understand the reduction operation mentioned in the documentation, because i'm not familiar with tensor notation. the examples i found in the tests folder just showed how to handle something like "comp(Base(#1).vBase(#2))".
so how do i have to set these indices to get a result?
i hope i did a proper problem description and someone can help me.

best regards and thanks in advance, alex





reply via email to

[Prev in Thread] Current Thread [Next in Thread]