Mesh.cpp

//
//  file    : Mesh.cpp
//  date    : 22 dec 2001
//  changed : 28 jul 2002
//  author  : jh
//


#include <qfile.h>
#include <qtextstream.h> 
#include <qptrvector.h> 
#include <qdatetime.h>


#include "Mesh.h"


extern char* __delaunay_lib_version;


Mesh::Mesh()
{
  init();
}


void
Mesh::init()
{
 _vertices.setAutoDelete(true);
 _triangles.setAutoDelete(true);
}



void
Mesh::clear()
{
  _vertices.clear();
  _triangles.clear();
}





void
Mesh::add_triangle(Triangle* t) 
{
#ifdef _DEBUG_
  fprintf(stderr, "<Mesh::addTriangle()>\n");
#endif

  _triangles.insert(t->get_ID(), t ); 

#ifdef _DEBUG_
  fprintf(stderr, "</Mesh::addTriangle()>\n");
#endif

}



void
Mesh::add_vertex(Vertex* p) 
{
  _vertices.insert(p->ID(), p ); 
}




void
Mesh::remove_triangle(Triangle* t)
{
  _triangles.remove( t->get_ID() );
}



void
Mesh::remove_vertex(Vertex* p)
{
  _vertices.remove( p->ID() );
}





void
Mesh::write_XML(const QString& filename)
{
  QFile f(filename);
  if(!f.open(IO_WriteOnly)) return;
  QTextStream* cout = new QTextStream( &f );
  write_XML(*cout);
}



void
Mesh::write_XML(QTextStream &cout)
{  
#ifdef _SOFISTIK_
  cout << "+PROG GENF M 40 \n";
  cout << "  SYST RAUM GDIV 1000 GDIR ZZ OPTI NEIN \n\n";
  cout << "  MAT 1 E 1000 \n\n";

  QDictIterator<Vertex> point_it(_vertices);

  for( ; point_it.current(); ++point_it ) {

    cout << "  KNOT " << point_it.current()->ID();
    cout << " X " << point_it.current()->x();
    cout << " Y " << point_it.current()->y();
    cout << " Z 0 \n";
    
  }

  cout << "\n\n\n";

  
  QDictIterator<Triangle> tri_it(_triangles);
  
  for(; tri_it.current(); ++tri_it) {

    cout << "  QUAD " << tri_it.current()->get_ID();
    cout << "  " << tri_it.current()->vertex_ID(0);
    cout << "  " << tri_it.current()->vertex_ID(1);
    cout << "  " << tri_it.current()->vertex_ID(2);
    cout << " MNR 1 D 1.0 \n";
  }

  cout << endl;

  cout << "ENDE\n\n\n";

  cout << "+PROG ASE M 40 \n";
  cout << "  SYST RAUM PROB LINE \n";
  cout << "  STEU SOLV 2 \n";
  cout << "  LF 1 EGY 0.0 \n";
  cout << "  #include Lasten.inp \n";
  cout << "ENDE\n";
#endif


  cout << "<?xml version=\"1.0\" encoding=\"ISO-8859-1\"?>\n";
  cout << "<!-- created by delaunay mesh generator \n";
  cout << "     " << __delaunay_lib_version << endl;
  cout << "     " << QTime::currentTime().toString() << " -->\n";
  cout << "<objects>\n";
  

  QDictIterator<Vertex> point_it(_vertices);

  for( ; point_it.current(); ++point_it )
    point_it.current()->write_XML(cout);


  QDictIterator<Triangle> tri_it(_triangles);

  for(; tri_it.current(); ++tri_it)
    tri_it.current()->write_XML(cout);

  cout << "</objects>\n";


}

---