EnsiteToXml.cpp File Reference

#include <iostream>
#include <fstream>
#include <iomanip>
#include <sstream>
#include <vector>
#include <list>
#include <tinyxml.h>

Include dependency graph for EnsiteToXml.cpp:

Go to the source code of this file.

Classes
struct	Vertex
	Represents a vertex in the mesh. More...
struct	Edge
struct	TwoDElement

Functions
void	readEnsite (const char *filename, vector< Vertex > &pVertexList, vector< TwoDElement > &pElementList)
int	GetEdge (vector< int > &vert, vector< Edge > &edges, int elm_type)
void	createEdgeList (vector< Vertex > &pVertexList, vector< Edge > &pEdgeList, vector< TwoDElement > &pElementList)
void	WriteToXMLFile (const char *outfile, const vector< Vertex > &nodes, const vector< Edge > &edges, const vector< TwoDElement > &twoDElements)
void	WriteMshFile (const char *outfile, const vector< Vertex > &nodes, const vector< Edge > &edges, const vector< TwoDElement > &twoDElements)
int	main (int argc, char *argv[])

Function Documentation

void createEdgeList	(	vector< Vertex > &	pVertexList,
		vector< Edge > &	pEdgeList,
		vector< TwoDElement > &	pElementList
	)

Definition at line 141 of file EnsiteToXml.cpp.

References GetEdge().

Referenced by main().

{
    int i,j,edgeid;

    cout << "Creating edge list..." << flush;
    for (i = 0; i < pElementList.size(); ++i)
    {
        for (j = 0; j < 3; ++j)
        {
            vector<int> verts;
            verts.push_back(pElementList[i].verts[j%3]);
            verts.push_back(pElementList[i].verts[(j+1)%3]);
            edgeid = GetEdge(verts, pEdgeList, 2);
            pElementList[i].edges.push_back(edgeid);
        }
    }
    cout << "done." << endl;
}

int GetEdge	(	vector< int > &	vert,
		vector< Edge > &	edges,
		int	elm_type
	)

Definition at line 113 of file EnsiteToXml.cpp.

Referenced by Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), createEdgeList(), GenerateMapEidsv1v2(), MapEdgeVertices(), Nektar::Utilities::ProcessSpherigon::Process(), Nektar::LocalRegions::TriExp::v_FwdTrans_BndConstrained(), and Nektar::LocalRegions::QuadExp::v_FwdTrans_BndConstrained().

{
  int i;
  int size = edges.size();
  for(i = 0; i < size; i++)
{
  if( ((vert[0] == edges[i].vert[0])&&
       (vert[1] == edges[i].vert[1]))
                ||
      ((vert[1] == edges[i].vert[0])&&
       (vert[0] == edges[i].vert[1])) )
    {
      return i;
    }
}

  // for the last edge of an element
  vector<int> edgevert;
  edgevert.push_back(vert[0]);
  edgevert.push_back(vert[1]);
  if(elm_type==9)
{
  edgevert.push_back(vert[2]);
}
  edges.push_back( Edge(size,edgevert) );
  return size;
}

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 294 of file EnsiteToXml.cpp.

References createEdgeList(), readEnsite(), WriteMshFile(), and WriteToXMLFile().

{
    if(argc != 4)
    {
        fprintf(stderr,"Usage: EnsiteToXml ensite.xml session.xml mesh.msh\n");
        exit(1);
    }

    vector<Vertex> vVertexList;
    vector<Edge> vEdgeList;
    vector<TwoDElement> vElementList;

    readEnsite(argv[1], vVertexList, vElementList);
    createEdgeList(vVertexList, vEdgeList, vElementList);
    WriteToXMLFile(argv[2], vVertexList, vEdgeList, vElementList);
    WriteMshFile(argv[3], vVertexList, vEdgeList, vElementList);
}

void readEnsite	(	const char *	filename,
		vector< Vertex > &	pVertexList,
		vector< TwoDElement > &	pElementList
	)

Definition at line 57 of file EnsiteToXml.cpp.

Referenced by main().

{
    cout << "Reading Ensite file..." << flush;
    TiXmlDocument vEnsite(filename);
    TiXmlElement * vElmt;
    stringstream vTmp;
    int n;

    if (!vEnsite.LoadFile()) {
        cout << "Unable to open file " << filename << endl;
        exit(-1);
    }

    vElmt = vEnsite.FirstChildElement("DIF");
    if (!vElmt) cout << "Unable to find DIF tag" << endl;
    vElmt = vElmt->FirstChildElement("DIFBody");
    if (!vElmt) cout << "Unable to find DIFBody tag" << endl;
    vElmt = vElmt->FirstChildElement("Volumes");
    if (!vElmt) cout << "Unable to find Volumes tag" << endl;
    vElmt = vElmt->FirstChildElement("Volume");
    if (!vElmt) cout << "Unable to find Volume tag" << endl;

    TiXmlElement * vVertices = vElmt->FirstChildElement("Vertices");
    if (!vVertices) cout << "Unable to find Vertices tag" << endl;
    int vNVertices = atoi(vVertices->Attribute("number"));
    TiXmlElement * vPolygons = vElmt->FirstChildElement("Polygons");
    if (!vPolygons) cout << "Unable to find Polygons tag" << endl;
    int vNPolygons = atoi(vPolygons->Attribute("number"));

    // Read vertices
    vTmp.write(vVertices->GetText(), strlen(vVertices->GetText()));
    double x, y, z;
    for (n = 0; n < vNVertices; ++n)
    {
        vTmp >> x >> y >> z;
        pVertexList.push_back(Vertex(n,x,y,z));
    }

    // Read polygons
    vTmp.clear();
    vTmp.write(vPolygons->GetText(), strlen(vPolygons->GetText()));
    int a, b, c;
    for (n = 0; n < vNPolygons; ++n)
    {
        vTmp >> a >> b >> c;
        vector<int> vert;
        vert.push_back(a-1);
        vert.push_back(b-1);
        vert.push_back(c-1);
        pElementList.push_back(TwoDElement(n,2,vert));
    }
    cout << "done." << endl;
}

void WriteMshFile	(	const char *	outfile,
		const vector< Vertex > &	nodes,
		const vector< Edge > &	edges,
		const vector< TwoDElement > &	twoDElements
	)

Definition at line 264 of file EnsiteToXml.cpp.

Referenced by main().

{
    cout << "Write GMSH file..." << flush;
    int n;
    ofstream out(outfile);

    out << "$MeshFormat" << endl;
    out << "2.1 0 8" << endl;
    out << "$EndMeshFormat" << endl;
    out << "$Nodes" << endl;
    out << nodes.size() << endl;
    for (n = 0; n < nodes.size(); ++n)
    {
        out << n+1 << " " << nodes[n].x << " " << nodes[n].y << " "
                << nodes[n].z << endl;
    }
    out << "$EndNodes" << endl;
    out << "$Elements" << endl;
    out << twoDElements.size() << endl;
    for (n = 0; n < twoDElements.size(); ++n)
    {
        out << n+1 << " 2 3 1 0 0 " << twoDElements[n].verts[0] + 1 << " "
                << twoDElements[n].verts[1] + 1
                << " " << twoDElements[n].verts[2] + 1 << endl;
    }
    out << "$EndElements" << endl;
    cout << "done." << endl;
}

void WriteToXMLFile	(	const char *	outfile,
		const vector< Vertex > &	nodes,
		const vector< Edge > &	edges,
		const vector< TwoDElement > &	twoDElements
	)

Definition at line 162 of file EnsiteToXml.cpp.

References Nektar::iterator.

Referenced by main().

{
    cout << "Writing XML file..." << flush;

    TiXmlDocument doc;
    TiXmlDeclaration* decl = new TiXmlDeclaration( "1.0", "utf-8", "");
    doc.LinkEndChild( decl );

    TiXmlElement * root = new TiXmlElement( "NEKTAR" );
    doc.LinkEndChild( root );

    //---------------------------------------------
    // Write DIM and SPACE
    TiXmlElement * geomTag = new TiXmlElement( "GEOMETRY" );
    geomTag->SetAttribute("DIM", "2");
    geomTag->SetAttribute("SPACE", "3");
    root->LinkEndChild( geomTag );
    //---------------------------------------------

    //---------------------------------------------
    // Write VERTEX
    TiXmlElement* verTag = new TiXmlElement( "VERTEX" );

    for( int i = 0; i < nodes.size(); ++i ) {
        stringstream s;
        s << scientific << setprecision(3) <<  nodes[i].x << " "
          << nodes[i].y << " " << nodes[i].z;
        TiXmlElement * v = new TiXmlElement( "V" );
        v->SetAttribute("ID",nodes[i].id);
        v->LinkEndChild( new TiXmlText(s.str()) );
        verTag->LinkEndChild(v);
    }
    geomTag->LinkEndChild( verTag );
    //---------------------------------------------


    //---------------------------------------------
    // Write EDGE
    int edgecnt = 0;
    verTag = new TiXmlElement( "EDGE" );

    for( int i = 0; i < edges.size(); ++i ) {
        stringstream s;

        s << setw(5) << edges[i].vert[0] << "  " << edges[i].vert[1] << "   ";
        TiXmlElement * e = new TiXmlElement( "E" );
        e->SetAttribute("ID",edgecnt++);
        e->LinkEndChild( new TiXmlText(s.str()) );
        verTag->LinkEndChild(e);
    }
    geomTag->LinkEndChild( verTag );
    //---------------------------------------------


    //--------------------------------------------
    // Write ELEMENT
    verTag = new TiXmlElement( "ELEMENT" );
    int elmcnt = 0;

    for(int i=0; i<twoDElements.size(); ++i){
        stringstream st;
        for(int j=0; j<twoDElements[i].edges.size(); ++j){
            st << setw(5) << twoDElements[i].edges[j] << " ";
        }

        TiXmlElement *elm_tag;
        elm_tag = new TiXmlElement("T");
        elm_tag->SetAttribute("ID", elmcnt++);
        elm_tag->LinkEndChild( new TiXmlText(st.str()) );
        verTag->LinkEndChild(elm_tag);
    }
    geomTag->LinkEndChild( verTag );
    //--------------------------------------------------


    //--------------------------------------------------
    // Write COMPOSITE
    verTag = new TiXmlElement("COMPOSITE");

    list<int>::iterator it;
    list<int> compList;

    stringstream st;
    st << " T[0-" << twoDElements.size()-1 << "]";

    TiXmlElement *comp_tag = new TiXmlElement("C"); // Composite
    comp_tag->SetAttribute("ID", "0");
    comp_tag->LinkEndChild( new TiXmlText(st.str()) );
    verTag->LinkEndChild(comp_tag);
    geomTag->LinkEndChild( verTag );
    //--------------------------------------------------

    TiXmlElement * domain = new TiXmlElement ("DOMAIN" );
    domain->LinkEndChild( new TiXmlText( " C[0] " ));
    geomTag->LinkEndChild( domain );

    doc.SaveFile(outfile );

    cout << "done." << endl;
} // end of function WriteToXMLFile