conformalmesher.h

Go to the documentation of this file.

#ifndef CONFORMALMESHER_H
#define CONFORMALMESHER_H

#include "util.h"
#include "cell.h"
#include <vector>
#include "hash_classes.h"
#include "hash_map.h"
#include "hash_set.h"

/* Predeclaration of Classes */
class BezierVertex;
class BezierEdge;
class BezierMesh;
class BoundaryMesh;

class ConformalMesher {
    private:
        /* the mesh we're creating, in two parts */
        BezierMesh *bezier;
        BoundaryMesh *bdry;

    public:
        class Subspline; // just to allow having operator<< be static
    private:
        /* data structures needed during the course of the algorithm */
        std::list<Subspline*> subsplines; /* all splines, for the dtor */
        std::vector<Subspline*> spline_stack; /* of segments to split */
        hashers::hash_set<Subspline*, CastHash<Subspline*> > is_encroached;
        hashers::hash_set<Subspline*, CastHash<Subspline*> > in_stack;
        hashers::hash_map<BezierVertex*, std::vector<Subspline*> >
                          incident_edges;

        std::vector<Point2D> remove_points;
        std::vector<std::pair<BoundaryFace*,Point2D> > face_points;

        /* the functions invoked by meshing */
        void create_bbox();
        void insert_d0_vertices();
        void create_subsplines();
        void conform();
        void reticulate_splines();
        void set_face_pointers();
        void excavate_holes();

        /* subfunctions of those */
        void link_edge(Subspline*);
        void unlink_edge(Subspline *s);
        bool segment_appears(Subspline*);
        bool angle_is_positive(Subspline*,BezierVertex*,BezierEdge*);
        BezierVertex *split_destroy_segment(Subspline*,
                                            std::vector<Subspline*>&);
        void collect_faces(Point2D, hashers::hash_set<BezierTriangle*>&);
        void set_face_pointer(BoundaryFace *face);
        void remove_hole(Point2D);

    public:
        ConformalMesher(BezierMesh*, BoundaryMesh*);
        ~ConformalMesher();

        /* Hacks to deal with how to link geometry and topology. */
        void add_face_point(BoundaryFace *f, Point2D p) {
          face_points.push_back(std::pair<BoundaryFace*,Point2D>(f,p));
        }
        void add_hole(Point2D p) {
          remove_points.push_back(p);
        }

        void mesh();
};

#endif /* MESHCONSTRUCTOR_H */

---