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
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
