Dear Yves,
Sorry for my disturb. Right now, I can accomplish what I want by the code I
in the last of the email.
Second, you have memtioned that do not use the copy constructor for
getfem::mesh. Then, what will you suggest for the situation in my code II?
I tried the code III which however has some problems.
sincerely
Tsai
code I:
dal::bit_vector MeshRegion_points_index(const getfem::mesh_fem&
meshFem,const getfem::mesh_im& meshIm,unsigned int rg_id)
{
unsigned int nb_dof=meshFem.nb_dof();
dal::bit_vector temp;
vector<REAL> Btemp(nb_dof);
gmm::resize(Btemp,nb_dof); gmm::clear(Btemp);
getfem::generic_assembly assem;
assem.push_mi(meshIm);
assem.push_mf(meshFem);
assem.push_vec(Btemp);
assem.set("V(#1)+=comp(Base(#1))(:)");
assem.assembly(rg_id);
for(unsigned int i=0;i<nb_dof;i++)
if(Btemp[i])temp[i]=true;
return temp;
}
code II:
class FemModel
{
mesh msh;
mesh_fem mf;
mesh_im mim;
FemModel(mesh MSH,mesh_fem MF, mesh_im MIM):msh(MSH),mf(MF),min(MIM){}
}
code III:
class FemModel
{
mesh msh;
mesh_fem mf;
mesh_im mim;
FemModel(mesh MSH,mesh_fem MF, mesh_im MIM) {
msh.copy_from(MSH);
mim.init_with_mesh(msh);
mf_u.init_with_mesh(msh);
}
}
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
