Question #240117 on Yade changed:
https://answers.launchpad.net/yade/+question/240117
Status: Answered => Open
burak er is still having a problem:
Hi Anton,
I created a class[1] which I use to import mesh file than create
elements using gmsh by its releated structure GModel. When I use this
class gmshModel variable does not load mesh file(without giving error).
If I compile the same lines [2] without yade the mesh loading goes very
well.
[1]class GMshMesh:public Serializable{
public:
shared_ptr<GModel> gmshModel;
void initfunc(void){
gmshModel=shared_ptr<GModel>(new GModel);
}
void LoadMesh(std::string filename)
{
boost::filesystem::path mypath(filename);
if(mypath.extension()!=boost::filesystem::path(".msh")){
std::string errormessage ="A .msh file must be
specified for the input";
throw std::invalid_argument(errormessage);
}
else
{
std::cout<<"model is created "<<std::endl;
gmshModel->readMSH("/home/burak/boxcly.msh");
std::cout<<"model is read "<<std::endl;
gmshModel->getMeshStatus();
// Create domain
elems.clear();// clear our vector
std::cout<<"elements are cleared "<<std::endl;
interfacingpairs.clear();
std::cout<<"interfacing pairs are cleared "<<std::endl;
size_t numberofelements=gmshModel->getNumMeshElements();//get
number of mesh elements
std::cout<<"getting number of elements
as"<<numberofelements<<std::endl;
elems.resize(numberofelements); // resize it as the
new number of elements
std::cout<<"resizing elements "<<std::endl;
for(size_t iter=0;iter<numberofelements;iter++){
//We can check element type here if required but for
now we assume it to be tetrahedral element
//if(Ng_GetElementType(iter)==NG_ELEMENT_TYPE.NG_TET)
{std::cout<<"getting element "<< iter<<std::endl;
int vertindex[4];
MElement*
element=gmshModel->getMeshElementByTag(iter);
std::vector<MVertex*> vertices;
element->getVertices(vertices);
//Create new elements
Vector3r vertex0;
vertex0<<vertices[0]->x(),vertices[0]->y(),vertices[0]->z();
Vector3r vertex1;
vertex1<<vertices[1]->x(),vertices[1]->y(),vertices[1]->z();
Vector3r vertex2;
vertex2<<vertices[2]->x(),vertices[2]->y(),vertices[2]->z();
Vector3r vertex3;
vertex3<<vertices[3]->x(),vertices[3]->y(),vertices[3]->z();
//share_ptr handles deletion of the newly created
element
elems[iter]=shared_ptr<MeshElement>(new
Linear4NodeTetMeshElem(vertex0,vertex1,vertex2,vertex3));
}
}
}
}
public:
YADE_CLASS_BASE_DOC_ATTRS_INIT_CTOR_PY(GMshMesh,Serializable,"An
interfacing pair",
,
initfunc();
,
/* ctor */,
/* py */
.def("doMeshing",&GMshMesh::LoadMesh,"Does meshing")
.def("getElements",&GMshMesh::getElements,"Return a
reference to element object")
.def("getPairs",&GMshMesh::getPairs,"Return a reference
to element object")
);
private:
MeshElementContainer elems;//Mesh Elements
ElementPairContainer interfacingpairs;//Interacing pairs
};
REGISTER_SERIALIZABLE(GMshMesh);
[2]#include <iostream>
#include<set>
#include<boost/shared_ptr.hpp>
int main(void)
{
boost::shared_ptr<GModel> gmshModel=boost::shared_ptr<GModel>(new GModel());
std::cout<<"model is created "<<std::endl;
gmshModel->readMSH("/home/burak/boxcly.msh");
std::cout<<"model is read "<<std::endl;
gmshModel->getMeshStatus();
// Create domain
size_t numberofelements=gmshModel->getMaxElementNumber();//get number
of mesh elements
std::cout<<numberofelements;
return 0;
}
--
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
