CGAL users discussion list

["Amir Vaxman" <>]

I enter a constraint triangulation by vertices, and then operate the mesh generator, but I see no change. A little debugging indicated that there were no bad faces. I tried it on several outputs which apparently had bad faces, so it must be something in my definition. My definitions are:

 

 

 

typedef CGAL::Exact_predicates_inexact_constructions_kernel Kernel;

typedef CGAL::Triangulation_vertex_base_2<Kernel> FaceVertexBase;

typedef CGAL::Triangulation_vertex_base_with_info_2<MeshVertex*, Kernel, FaceVertexBase> IFaceVertexBase;  //MeshVertex is some object defined earlier

typedef CGAL::Delaunay_mesh_face_base_2<Kernel> FaceFaceBase;

 

template <class Gt,

          class Fb >

class Enriched_face_base_2 : public Fb {

public:

  typedef Gt Geom_traits;

  typedef typename Fb::Vertex_handle Vertex_handle;

  typedef typename Fb::Face_handle Face_handle;

 

  template < typename TDS2 >

  struct Rebind_TDS {

    typedef typename Fb::template Rebind_TDS<TDS2>::Other Fb2;

    typedef Enriched_face_base_2<Gt,Fb2> Other;

  };

 

protected:

 

public:

      bool Inside;

      bool Inquired;

  Enriched_face_base_2(): Fb(), Inside(false), Inquired(false){};

 

  Enriched_face_base_2(Vertex_handle v0,

                   Vertex_handle v1,

                   Vertex_handle v2)

    : Fb(v0,v1,v2), Inside(false), Inquired(false) {};

 

  Enriched_face_base_2(Vertex_handle v0,

                   Vertex_handle v1,

                   Vertex_handle v2,

                   Face_handle n0,

                   Face_handle n1,

                   Face_handle n2)

    : Fb(v0,v1,v2,n0,n1,n2), Inside(false), Inquired(false) {};

 

};

typedef Enriched_face_base_2<Kernel, FaceFaceBase> IFaceBase;

typedef CGAL::Triangulation_data_structure_2<IFaceVertexBase, IFaceBase> Trids;

typedef CGAL::Constrained_Delaunay_triangulation_2<Kernel, Trids> MeshCDT;

typedef CGAL::Delaunay_mesh_size_criteria_2<MeshCDT> Criteria;

typedef CGAL::Delaunay_mesher_2<MeshCDT, Criteria> TriMesher;

 

…entering input as vertices in a MeshCDT object named “mCDT”

 

CGAL::refine_Delaunay_mesh_2(mCDT, Criteria(0.125, 0.5));

 

 

 

Thanks,

Amir.