[Getfem-users] Problem assembling convection contribution

[Cédric Venet]

Hi,

I am trying to assemble the term "v.Grad u" where v is a field of vector 
data on my mesh. u is a scalar field, and the unknow of my equation. I 
am trying to do this:

  getfem::generic_assembly assem;
  assem.push_mi(fm.GetIntegrationMesh());
  assem.push_mf(fm.GetShapeMesh());
  std::vector<T> tmpve(1,m_coef);
  assem.push_data(tmpve);
  assem.push_data(m_speed);
  assem.push_mat(m_A);
  
assem.set("F=data(1);M$1(#1,#1)+=comp(Base(#1).Grad(#1).Base(#1))(i,:,k,:).data$2(k,i).F(p)");
  assem.assembly();

however, it seems we cannot use two indexes on a data (here is the the 
node/DOF number and k the dimension number).
can I do this? should I push three data vector? should I define a new 
fem mesh for the data?

regards,
Cédric Venet
Laboratoire MAS