DelaunayMeshGenerator.h

#ifndef _DELAUNAYMESHGENERATOR_H_
#define _DELAUNAYMESHGENERATOR_H_

/*
    This library is free software; you can redistribute it and/or
    modify it under the terms of the GNU Library General Public
    License as published by the Free Software Foundation; either
    version 2 of the License, or (at your option) any later version.

    This library is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    Library General Public License for more details.

    You should have received a copy of the GNU Library General Public License
    along with this library; see the file COPYING.LIB.  If not, write to
    the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
    Boston, MA 02111-1307, USA.

    file    : DelaunayMeshGenerator.h

*/


#include <qptrlist.h>
#include <qxml.h>

#include "Vertex.h"
#include "Edge.h"
#include "Triangle.h"

#include "Mesh.h"





#define CEps 1e-8




class DelaunayMeshGenerator : public QXmlDefaultHandler
{

public:


  DelaunayMeshGenerator(const QString& filename, bool verbose = false);


  
  DelaunayMeshGenerator(QIODevice& iodevice, bool verbose = false);



  ~DelaunayMeshGenerator();



  void generate();



  Mesh* mesh() const { return _mesh; }



  bool startElement(const QString&, const QString&, 
                    const QString& qName, const QXmlAttributes& );
   


  void write_XML(QTextStream& cout);



protected:


  void add_vertex(Vertex* v);


  void add_edge(Edge* e) { _edges.append(e); }

  
  Vertex* vertex(const QString& id) 
  {
    for(Vertex* v = _vertices.first(); v; v = _vertices.next()) {
      if(v->ID() == id) return v;
    }
    return 0;
  }



private:

  QPtrList<Vertex> _vertices;
  QPtrList<Vertex> _active_vertices;

  QPtrList<Edge>  _edges;
  QPtrList<Edge>  _active_edges;

  QPtrList<Triangle> _triangles;
 
  bool _verbose;
  
  double minX, minY, maxX, maxY;
  double box_minX, box_minY, box_maxX, box_maxY;
  
  double u_alt;

  int _count_ID;


  QString createID() {
    QString str;
    str.sprintf("id_%d", _count_ID++);
    return str;
  }


  Mesh* _mesh;


  void init();

  void reset();


  void open(QIODevice& iodevice) 
  {
    QXmlInputSource  source( &iodevice );
    QXmlSimpleReader reader;
    reader.setContentHandler(this);
    reader.parse( source ); 
  }


  
  bool test_intersection(Vertex* v1, Vertex* v2);


  bool is_voronoi_triangle(Edge* edge, Vertex* v);


  void create_triangle(Edge* basisedge, Vertex* vertex);

};


#endif

---