yade-users team mailing list archive

[burak er <question240117@xxxxxxxxxxxxxxxxxxxxx>]

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.