1 // { dg-do compile }
2 // { dg-additional-options "-Wno-return-type" }
3
4 namespace std {
5 template<class,class>struct pair{};
6 template<typename _Tp> struct _Vector_base {
7 struct _Vector_impl
8 {
9 _Tp* _M_start;
10 _Tp* _M_finish;
11 _Tp* _M_end_of_storage;
12 };
13 _Vector_impl _M_impl;
14 };
15 template<typename _Tp >
16 struct vector : _Vector_base<_Tp>
17 {
18 vector(const vector& __x);
19 };
20 }
21 namespace boost {
22 struct G{};
23 template <class T, class U >
24 struct modable2
25 : G
26 { };
27 }
28 namespace CGAL {
29 struct Rep { };
30 struct Handle
31 {
32 Handle() ;
33 Handle(const Handle& x) ;
34 Rep* PTR;
35 };
36 template <typename ET_> class Lazy_exact_nt
37 : Handle
38 , boost::modable2< Lazy_exact_nt<ET_>, int >
39 , boost::modable2< Lazy_exact_nt<ET_>, double >
40 { };
41 struct CC_iterator { };
42 struct Triangulation_data_structure_3 {
43 typedef CC_iterator Vertex_handle;
44 typedef CC_iterator Cell_handle;
45 typedef std::pair<Cell_handle, int> Facet;
46 };
47 template < class GT, class Tds_ > struct Triangulation_3 {
48 typedef Tds_ Tds;
49 typedef typename GT::Point_3 Point;
50 typedef typename Tds::Facet Facet;
51 typedef typename Tds::Vertex_handle Vertex_handle;
52 enum Locate_type { VERTEX=0, EDGE, FACET, CELL, OUTSIDE_CONVEX_HULL, OUTSIDE_AFFINE_HULL };
53 Tds _tds;
54 bool is_infinite(const Facet & f) const ;
55 };
56 template < class Gt, class Tds_ > struct Delaunay_triangulation_3 : public Triangulation_3<Gt, Tds_> { };
57 namespace Surface_mesher { enum Verbose_flag { VERBOSE, NOT_VERBOSE }; }
58 enum Mesher_level_conflict_status { NO_CONFLICT = 0, CONFLICT_BUT_ELEMENT_CAN_BE_RECONSIDERED, CONFLICT_AND_ELEMENT_SHOULD_BE_DROPPED };
59 struct Null_mesher_level {
60 template <typename P, typename Z> Mesher_level_conflict_status test_point_conflict_from_superior(P, Z) ;
61 };
62 template < class Tr, class Derived, class Element, class Previous, class Triangulation_traits > struct Mesher_level {
63 typedef Tr Triangulation;
64 typedef typename Triangulation::Point Point;
65 typedef typename Triangulation::Vertex_handle Vertex_handle;
66 typedef typename Triangulation_traits::Zone Zone;
67 typedef Previous Previous_level;
derivedMesher_level68 Derived& derived() { return static_cast<Derived&>(*this); }
69 Previous& previous_level;
Mesher_levelMesher_level70 Mesher_level(Previous_level& previous)
71 : previous_level(previous)
72 { }
73 Vertex_handle insert(Point p, Zone& z) ; // { dg-warning "used but never defined" }
74 Zone conflicts_zone(const Point& p, Element e) ; // { dg-warning "used but never defined" }
75 Element get_next_element() ; // { dg-warning "used but never defined" }
before_insertionMesher_level76 template <class Mesh_visitor> void before_insertion(Element& e, const Point& p, Zone& zone, Mesh_visitor visitor) {
77 visitor.before_insertion(e, p, zone);
78 }
after_insertionMesher_level79 template <class Mesh_visitor> void after_insertion(Vertex_handle vh, Mesh_visitor visitor) {
80 derived().after_insertion_impl(vh);
81 }
after_no_insertionMesher_level82 template <class Mesh_visitor> void after_no_insertion(const Element& e, const Point& p, Zone& zone, Mesh_visitor visitor) {
83 visitor.after_no_insertion(e, p, zone);
84 }
refineMesher_level85 template <class Mesh_visitor> void refine(Mesh_visitor visitor)
86 {
87 Element e = get_next_element();
88 const Mesher_level_conflict_status result = try_to_refine_element(e, visitor);
89 }
try_to_refine_elementMesher_level90 template <class Mesh_visitor> Mesher_level_conflict_status try_to_refine_element(Element e, Mesh_visitor visitor)
91 {
92 Point p ;
93 Zone zone = conflicts_zone(p, e);
94 const Mesher_level_conflict_status result = test_point_conflict(p, zone);
95 before_insertion(e, p, zone, visitor);
96 Vertex_handle v = insert(p, zone);
97 after_insertion(v, visitor);
98 after_no_insertion(e, p, zone, visitor);
99 }
test_point_conflictMesher_level100 Mesher_level_conflict_status test_point_conflict(const Point& p, Zone& zone)
101 {
102 return previous_level.test_point_conflict_from_superior(p, zone);
103 }
104 };
105 struct Null_mesh_visitor {
before_insertionNull_mesh_visitor106 template <typename E, typename P, typename Z> void before_insertion(E, P, Z) const {}
after_no_insertionNull_mesh_visitor107 template <typename E, typename P, typename Z> void after_no_insertion(E, P, Z) const {}
108 };
109 template <class Tr> struct Triangulation_ref_impl {
110 Triangulation_ref_impl(Tr& t);
111 };
112 template <typename Tr> struct Triangulation_mesher_level_traits_3
113 : public Triangulation_ref_impl<Tr>
114 {
115 typedef typename Tr::Facet Facet;
Triangulation_mesher_level_traits_3Triangulation_mesher_level_traits_3116 Triangulation_mesher_level_traits_3(Tr& t)
117 : Triangulation_ref_impl<Tr>(t)
118 { }
119 struct Zone {
120 typedef std::vector<int*> Cells;
121 typedef std::vector<Facet> Facets;
122 typedef typename Tr::Locate_type Locate_type;
123 Locate_type locate_type;
124 Cells cells;
125 Facets boundary_facets;
126 Facets internal_facets;
127 };
128 };
129 namespace Surface_mesher {
130 namespace details {
131 template <typename Base> struct Triangulation_generator {
132 typedef typename Base::Complex_2_in_triangulation_3 C2T3;
133 typedef typename C2T3::Triangulation Triangulation;
134 typedef Triangulation Type;
135 typedef Type type;
136 };
137 template <typename Base> struct Facet_generator {
138 typedef typename Triangulation_generator<Base>::type Tr;
139 typedef typename Tr::Facet Type;
140 typedef Type type;
141 };
142 template <typename Base, typename Self, typename Element, typename PreviousLevel = Null_mesher_level> struct Mesher_level_generator {
143 typedef typename Base::Complex_2_in_triangulation_3 C2T3;
144 typedef typename C2T3::Triangulation Triangulation;
145 typedef Triangulation_mesher_level_traits_3<Triangulation> Tr_m_l_traits_3;
146 typedef Mesher_level <Triangulation, Self, Element, PreviousLevel, Tr_m_l_traits_3> Type;
147 };
148 }
149 template < class C2T3, class Surface_, class SurfaceMeshTraits, class Criteria_ > struct Surface_mesher_base
150 : public Triangulation_mesher_level_traits_3<typename C2T3::Triangulation>
151 {
152 typedef C2T3 Complex_2_in_triangulation_3;
153 typedef Surface_ Surface;
154 typedef SurfaceMeshTraits Surface_mesh_traits;
155 typedef Criteria_ Criteria;
156 typedef typename C2T3::Triangulation Tr;
157 typedef typename Tr::Vertex_handle Vertex_handle;
158 typedef typename Tr::Facet Facet;
Surface_mesher_baseSurface_mesher_base159 Surface_mesher_base (C2T3& co, const Surface& s, const Surface_mesh_traits& mesh_traits, const Criteria& c)
160 : Triangulation_mesher_level_traits_3<Tr>(co.triangulation()), c2t3(co), tr(co.triangulation()), surf(s), meshtraits(mesh_traits), criteria(c)
161 { }
162 C2T3& c2t3;
163 Tr& tr;
164 const Surface& surf;
165 const Surface_mesh_traits& meshtraits;
166 const Criteria& criteria;
after_insertion_implSurface_mesher_base167 void after_insertion_impl(const Vertex_handle& v) {
168 after_insertion_handle_opposite_facet (Facet ());
169 after_insertion_handle_incident_facet (Facet ());
170 }
after_insertion_handle_incident_facetSurface_mesher_base171 void after_insertion_handle_incident_facet (const Facet& f) {
172 tr.is_infinite(f) ;
173 new_facet<false>(f);
174 }
175 template <bool remove_from_complex_if_not_in_restricted_Delaunay> void new_facet (const Facet& f) ; // { dg-warning "used but never defined" }
after_insertion_handle_opposite_facetSurface_mesher_base176 void after_insertion_handle_opposite_facet (const Facet& f) {
177 after_insertion_handle_incident_facet (f);
178 }
179 };
180 template < typename Base, typename Element = typename details::Facet_generator<Base>::type, typename PreviousLevel = Null_mesher_level, Verbose_flag verbose = NOT_VERBOSE > struct Surface_mesher
181 : public Base , public details::Mesher_level_generator< Base, Surface_mesher<Base, Element, PreviousLevel, verbose>, Element, PreviousLevel >::Type
182 {
183 typedef typename Base::Complex_2_in_triangulation_3 C2T3;
184 typedef typename Base::Surface Surface;
185 typedef typename Base::Criteria Criteria;
186 typedef typename Base::Surface_mesh_traits Surface_mesh_traits;
187 typedef typename details::Mesher_level_generator< Base, Surface_mesher<Base, Element, PreviousLevel, verbose>, Element, PreviousLevel >::Type Mesher_lvl;
188 using Mesher_lvl::refine;
189 Null_mesher_level null_mesher_level;
190 Null_mesh_visitor null_visitor;
191 bool initialized;
Surface_mesherSurface_mesher192 Surface_mesher(C2T3& c2t3, const Surface& surface, const Surface_mesh_traits& mesh_traits, const Criteria& criteria)
193 : Base(c2t3, surface, mesh_traits, criteria), Mesher_lvl(null_mesher_level), initialized(false)
194 { }
refine_meshSurface_mesher195 void refine_mesh () {
196 refine(null_visitor);
197 }
198 };
199 }
200 template <typename Surface> struct Surface_mesh_traits_generator_3 {
201 typedef typename Surface::Surface_mesher_traits_3 Type;
202 typedef Type type;
203 };
204 template < class Tr, typename Edge_info_ = void > struct Complex_2_in_triangulation_3 {
205 typedef Tr Triangulation;
206 Triangulation& triangulation();
207 };
208 template <class Tr> struct Surface_mesh_complex_2_in_triangulation_3
209 : public Complex_2_in_triangulation_3<Tr>
210 { };
211 struct Non_manifold_tag {};
212 template < typename C2T3, typename SurfaceMeshTraits_3, typename Criteria, typename Tag > struct Make_surface_mesh_helper {
213 typedef Surface_mesher::Surface_mesher_base< C2T3, typename SurfaceMeshTraits_3::Surface_3, SurfaceMeshTraits_3, Criteria> Mesher_base;
214 };
215 template <typename C2T3, typename SurfaceMeshTraits_3, typename Criteria, typename Tag, Surface_mesher::Verbose_flag verbosity = Surface_mesher::NOT_VERBOSE > struct Surface_mesher_generator {
216 typedef typename Make_surface_mesh_helper< C2T3, SurfaceMeshTraits_3, Criteria, Tag>::Mesher_base Mesher_base;
217 typedef Surface_mesher::Surface_mesher< Mesher_base, typename Surface_mesher::details::Facet_generator<Mesher_base>::type, Null_mesher_level, verbosity> Mesher;
218 typedef Mesher type;
219 };
make_surface_mesh(C2T3 & c2t3,const typename SurfaceMeshTraits_3::Surface_3 & surface,const SurfaceMeshTraits_3 & surface_mesh_traits,const Criteria & criteria)220 template <typename C2T3, typename SurfaceMeshTraits_3, typename Criteria> void make_surface_mesh(C2T3& c2t3, const typename SurfaceMeshTraits_3::Surface_3& surface, const SurfaceMeshTraits_3& surface_mesh_traits, const Criteria& criteria) {
221 typedef typename Surface_mesher_generator< C2T3, SurfaceMeshTraits_3, Criteria, Non_manifold_tag, Surface_mesher::NOT_VERBOSE >::type Mesher;
222 Mesher mesher(c2t3, surface, surface_mesh_traits, criteria);
223 mesher.refine_mesh();
224 }
225 template <class Kernel> struct Surface_mesh_triangulation_generator_3 {
226 typedef CGAL::Triangulation_data_structure_3 Tds;
227 typedef CGAL::Delaunay_triangulation_3<Kernel, Tds> Type;
228 };
229 namespace Surface_mesher {
230 namespace { struct Return_min { }; }
231 template < class GT, class Surface, class Unused = Return_min > struct Implicit_surface_oracle_3 {
232 typedef Surface Surface_3;
233 };
234 }
235 template< typename GT> struct Implicit_surface_3 {
236 typedef GT Geom_traits;
237 typedef Implicit_surface_3<Geom_traits > Self;
238 typedef Surface_mesher::Implicit_surface_oracle_3< Geom_traits, Self> Surface_mesher_traits_3;
239 };
240 }
241 struct K {
242 struct Point_3 {
243 CGAL::Lazy_exact_nt<double> a[3];
244 };
245 };
246 typedef CGAL::Surface_mesh_triangulation_generator_3<K>::Type Tr;
247 typedef CGAL::Surface_mesh_complex_2_in_triangulation_3<Tr> C2T3;
248 typedef CGAL::Implicit_surface_3<K > Surface;
249 typedef CGAL::Surface_mesh_traits_generator_3<Surface>::type Traits;
f()250 void f() {
251 C2T3 c2t3 ;
252 CGAL::make_surface_mesh(c2t3, Surface(), Traits(), 3);
253 }
254