mirror of https://github.com/CGAL/cgal
Working on the hierarchy
This commit is contained in:
Nico Kruithof
parent
7b1ab8eb5e
commit
62ee5054c2
|
|
@ -0,0 +1,31 @@
|
|||
#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
|
||||
#include <CGAL/Delaunay_triangulation_2.h>
|
||||
#include <CGAL/Triangulation_hierarchy_2.h>
|
||||
#include <CGAL/point_generators_2.h>
|
||||
#include <CGAL/algorithm.h>
|
||||
#include <cassert>
|
||||
|
||||
typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
|
||||
|
||||
typedef CGAL::Triangulation_vertex_base_2<K> Vbb;
|
||||
typedef CGAL::Triangulation_hierarchy_vertex_base_2<Vbb> Vb;
|
||||
typedef CGAL::Triangulation_face_base_2<K> Fb;
|
||||
typedef CGAL::Triangulation_data_structure_2<Vb,Fb> Tds;
|
||||
typedef CGAL::Delaunay_triangulation_2<K,Tds> Dt;
|
||||
typedef CGAL::Triangulation_hierarchy_2<Dt> Triangulation;
|
||||
typedef Triangulation::Point Point;
|
||||
typedef CGAL::Creator_uniform_2<double,Point> Creator;
|
||||
|
||||
int main( )
|
||||
{
|
||||
std::cout << "insertion of 1000 random points" << std::endl;
|
||||
Triangulation t;
|
||||
CGAL::Random_points_in_square_2<Point,Creator> g(1.);
|
||||
CGAL::cpp11::copy_n( g, 1000, std::back_inserter(t));
|
||||
|
||||
//verbose mode of is_valid ; shows the number of vertices at each level
|
||||
std::cout << "The number of vertices at successive levels" << std::endl;
|
||||
assert(t.is_valid(true));
|
||||
|
||||
return 0;
|
||||
}
|
||||
|
|
@ -136,6 +136,13 @@ public:
|
|||
hints[level] = v->face();
|
||||
|
||||
v->set_down (prev);
|
||||
// NGHK: Added
|
||||
if (hierarchy[level]->number_of_sheets()[0] != 1) {
|
||||
std::vector<Vertex_handle> vtc
|
||||
= hierarchy[level]->periodic_copies(v);
|
||||
for (unsigned int i=0 ; i<vtc.size() ; i++) vtc[i]->set_down(prev);
|
||||
}
|
||||
|
||||
prev->set_up (v);
|
||||
prev = v;
|
||||
}
|
||||
|
|
@ -202,7 +209,7 @@ Periodic_2_triangulation_hierarchy_2(const Periodic_2_triangulation_hierarchy_2<
|
|||
// create an empty triangulation to be able to delete it !
|
||||
hierarchy[0] = this;
|
||||
for(int i=1;i<maxlevel;++i)
|
||||
hierarchy[i] = new PTr_Base(tr.geom_traits());
|
||||
hierarchy[i] = new PTr_Base(tr.domain(), tr.geom_traits());
|
||||
copy_triangulation(tr);
|
||||
}
|
||||
|
||||
|
|
@ -242,8 +249,6 @@ copy_triangulation(const Periodic_2_triangulation_hierarchy_2<PTr> &tr)
|
|||
}
|
||||
}
|
||||
|
||||
add_hidden_vertices_into_map(Weighted_tag(), V);
|
||||
|
||||
{
|
||||
for(int i=1;i<maxlevel;++i) {
|
||||
for( Finite_vertices_iterator it=hierarchy[i]->finite_vertices_begin();
|
||||
|
|
|
|||
|
|
@ -0,0 +1,70 @@
|
|||
// Copyright (c) 1998 INRIA Sophia-Antipolis (France).
|
||||
// All rights reserved.
|
||||
//
|
||||
// This file is part of CGAL (www.cgal.org).
|
||||
// You can redistribute it and/or modify it under the terms of the GNU
|
||||
// General Public License as published by the Free Software Foundation,
|
||||
// either version 3 of the License, or (at your option) any later version.
|
||||
//
|
||||
// Licensees holding a valid commercial license may use this file in
|
||||
// accordance with the commercial license agreement provided with the software.
|
||||
//
|
||||
// This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
|
||||
// WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
|
||||
//
|
||||
// $URL$
|
||||
// $Id$
|
||||
//
|
||||
//
|
||||
// Author(s) : Olivier Devillers <Olivivier.Devillers@sophia.inria.fr>
|
||||
// Mariette Yvinec <Mariette.Yvinec@sophia.inria.fr>
|
||||
|
||||
#ifndef CGAL_PERIODIC_2_TRIANGULATION_HIERARCHY_VERTEX_BASE_2_H
|
||||
#define CGAL_PERIODIC_2_TRIANGULATION_HIERARCHY_VERTEX_BASE_2_H
|
||||
|
||||
#include <CGAL/basic.h>
|
||||
|
||||
namespace CGAL {
|
||||
|
||||
template < class Gt, class Vb = CGAL::Periodic_2_triangulation_vertex_base_2<Gt> >
|
||||
class Periodic_2_triangulation_hierarchy_vertex_base_2
|
||||
: public Vb
|
||||
{
|
||||
typedef typename Vb::Triangulation_data_structure Tds;
|
||||
typedef Gt Geom_traits;
|
||||
typedef typename Gt::Point_2 Point;
|
||||
typedef Tds Triangulation_data_structure;
|
||||
typedef typename Tds::Face_handle Face_handle;
|
||||
typedef typename Tds::Vertex_handle Vertex_handle;
|
||||
|
||||
template < typename Tds2 >
|
||||
struct Rebind_Tds {
|
||||
typedef typename Vb::template Rebind_Tds<Tds2>::Other Vb2;
|
||||
typedef Triangulation_vertex_base_2<Gt, Vb2> Other;
|
||||
};
|
||||
typedef typename Vb::Offset Offset;
|
||||
|
||||
Periodic_2_triangulation_hierarchy_vertex_base_2()
|
||||
: Base(), _up(0), _down(0)
|
||||
{}
|
||||
Periodic_2_triangulation_hierarchy_vertex_base_2(const Point & p, Face_handle f)
|
||||
: Base(p,f), _up(0), _down(0)
|
||||
{}
|
||||
Periodic_2_triangulation_hierarchy_vertex_base_2(const Point & p)
|
||||
: Base(p), _up(0), _down(0)
|
||||
{}
|
||||
|
||||
Vertex_handle up() {return _up;}
|
||||
Vertex_handle down() {return _down;}
|
||||
void set_up(Vertex_handle u) {_up=u;}
|
||||
void set_down(Vertex_handle d) { _down=d;}
|
||||
|
||||
|
||||
private:
|
||||
Vertex_handle _up; // same vertex one level above
|
||||
Vertex_handle _down; // same vertex one level below
|
||||
};
|
||||
|
||||
} //namespace CGAL
|
||||
|
||||
#endif // CGAL_PERIODIC_2_TRIANGULATION_HIERARCHY_VERTEX_BASE_2_H
|
||||
|
|
@ -0,0 +1,84 @@
|
|||
// Copyright (c) 1998 INRIA Sophia-Antipolis (France).
|
||||
// All rights reserved.
|
||||
//
|
||||
// This file is part of CGAL (www.cgal.org).
|
||||
// You can redistribute it and/or modify it under the terms of the GNU
|
||||
// General Public License as published by the Free Software Foundation,
|
||||
// either version 3 of the License, or (at your option) any later version.
|
||||
//
|
||||
// Licensees holding a valid commercial license may use this file in
|
||||
// accordance with the commercial license agreement provided with the software.
|
||||
//
|
||||
// This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
|
||||
// WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
|
||||
//
|
||||
// $URL: svn+ssh://mcaroli@scm.gforge.inria.fr/svn/cgal/trunk/Periodic_2_triangulation_2/test/Periodic_2_triangulation_2/test_periodic_2_delaunay_2.cpp $
|
||||
// $Id: test_periodic_2_delaunay_2.cpp 48874 2009-04-22 12:54:28Z mcaroli $
|
||||
//
|
||||
//
|
||||
// Author(s) : Nico Kruithof
|
||||
// Manuel Caroli
|
||||
|
||||
#include <iostream>
|
||||
#include <fstream>
|
||||
|
||||
#include <CGAL/Periodic_2_Delaunay_triangulation_2.h>
|
||||
#include <CGAL/Periodic_2_triangulation_hierarchy_2.h>
|
||||
#include <CGAL/Periodic_2_triangulation_traits_2.h>
|
||||
|
||||
#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
|
||||
typedef CGAL::Exact_predicates_inexact_constructions_kernel K1;
|
||||
typedef CGAL::Periodic_2_triangulation_traits_2<K1> PTT1;
|
||||
typedef CGAL::Periodic_2_triangulation_vertex_base_2<PTT1> PVB1;
|
||||
typedef CGAL::Periodic_2_triangulation_face_base_2<PTT1> PFB1;
|
||||
typedef CGAL::Triangulation_hierarchy_vertex_base_2<PVB1> PHVB1;
|
||||
typedef CGAL::Triangulation_data_structure_2<PHVB1,PFB1> Tds1;
|
||||
typedef CGAL::Periodic_2_Delaunay_triangulation_2<PTT1,Tds1> PDT1;
|
||||
// Explicit instantiation of the whole class :
|
||||
template class CGAL::Periodic_2_triangulation_hierarchy_2<PDT1>;
|
||||
|
||||
// #include <CGAL/Exact_predicates_exact_constructions_kernel.h>
|
||||
// typedef CGAL::Exact_predicates_exact_constructions_kernel K2;
|
||||
// typedef CGAL::Periodic_2_triangulation_traits_2<K2> PTT2;
|
||||
// typedef CGAL::Periodic_2_triangulation_vertex_base_2<PTT2> DSVB2;
|
||||
// typedef CGAL::Periodic_2_triangulation_face_base_2<PTT2> DSFB2;
|
||||
// typedef CGAL::Triangulation_vertex_base_2<PTT2,DSVB2> VBB2;
|
||||
// typedef CGAL::Triangulation_hierarchy_vertex_base_2<VBB2> VB2;
|
||||
// typedef CGAL::Triangulation_face_base_2<PTT2,DSFB2> FB2;
|
||||
// typedef CGAL::Triangulation_data_structure_2<VB2,FB2> Tds2;
|
||||
// typedef CGAL::Periodic_2_Delaunay_triangulation_2<PTT2,Tds2> PDT2;
|
||||
// // Explicit instantiation of the whole class :
|
||||
// template class CGAL::Periodic_2_triangulation_hierarchy_2<PDT2>;
|
||||
|
||||
// #include <CGAL/MP_Float.h>
|
||||
// #include <CGAL/Simple_homogeneous.h>
|
||||
// typedef CGAL::Simple_homogeneous<CGAL::MP_Float> K3;
|
||||
// typedef CGAL::Periodic_2_triangulation_traits_2<K3> PTT3;
|
||||
// typedef CGAL::Periodic_2_triangulation_vertex_base_2<PTT3> DSVB3;
|
||||
// typedef CGAL::Periodic_2_triangulation_face_base_2<PTT3> DSFB3;
|
||||
// typedef CGAL::Triangulation_vertex_base_2<PTT3,DSVB3> VBB3;
|
||||
// typedef CGAL::Triangulation_hierarchy_vertex_base_2<VBB3> VB3;
|
||||
// typedef CGAL::Triangulation_face_base_2<PTT3,DSFB3> FB3;
|
||||
// typedef CGAL::Triangulation_data_structure_2<VB3,FB3> Tds3;
|
||||
// typedef CGAL::Periodic_2_Delaunay_triangulation_2<PTT3,Tds3> PDT3;
|
||||
// // Explicit instantiation of the whole class :
|
||||
// template class CGAL::Periodic_2_triangulation_hierarchy_2<PDT3>;
|
||||
|
||||
// #include <CGAL/_test_cls_periodic_3_delaunay_3.h>
|
||||
int main()
|
||||
{
|
||||
// TODO(NGHK): Re-enable
|
||||
|
||||
typedef CGAL::Periodic_2_triangulation_hierarchy_2< PDT1 > P3T3_1;
|
||||
//_test_cls_periodic_2_delaunay_2( P3T3_1() );
|
||||
|
||||
// typedef CGAL::Periodic_2_triangulation_hierarchy_2< PDT2 > P3T3_2;
|
||||
// // this takes too much time for the test suite.
|
||||
// //_test_cls_periodic_2_delaunay_2( P3T3_2() );
|
||||
|
||||
// typedef CGAL::Periodic_2_triangulation_hierarchy_2< PDT3 > P3T3_3;
|
||||
// // this takes too much time for the test suite.
|
||||
// //_test_cls_periodic_2_delaunay_2( P3T3_3(), true );
|
||||
|
||||
return 0;
|
||||
}
|
||||
|
|
@ -1,496 +1,7 @@
|
|||
#include "./types.h"
|
||||
#include "types.h"
|
||||
#include "interface_test.h"
|
||||
#include <CGAL/Periodic_2_triangulation_2.h>
|
||||
|
||||
template <class T>
|
||||
void test_constructor() {
|
||||
typedef typename T::Iso_rectangle Iso_rectangle;
|
||||
typedef typename T::Geom_traits Geom_traits;
|
||||
|
||||
T t;
|
||||
T t2(t);
|
||||
CGAL_assertion(t == t2);
|
||||
T t3 = t2;
|
||||
CGAL_assertion(t == t3);
|
||||
T t4(Iso_rectangle(0,0,2,2));
|
||||
T t5(Iso_rectangle(0,0,2,2), Geom_traits());
|
||||
t5.clear();
|
||||
|
||||
t.insert(Point(0.5, 0.5));
|
||||
CGAL_assertion(t != t2);
|
||||
CGAL_assertion(t != t3);
|
||||
|
||||
T t6(t);
|
||||
CGAL_assertion(t == t6);
|
||||
T t7 = t6;
|
||||
CGAL_assertion(t == t7);
|
||||
|
||||
t.clear();
|
||||
CGAL_assertion(t != t6);
|
||||
CGAL_assertion(t != t7);
|
||||
|
||||
std::vector<Point> pts;
|
||||
pts.push_back(Point(0.5, 0.5));
|
||||
pts.push_back(Point(0.25, 0.5));
|
||||
pts.push_back(Point(0.5, 0.25));
|
||||
pts.push_back(Point(0.25, 0.25));
|
||||
T t8(pts.begin(), pts.end());
|
||||
|
||||
t7.swap(t8);
|
||||
}
|
||||
|
||||
template <class T>
|
||||
void test_global_access() {
|
||||
T t;
|
||||
const T &t_const = t;
|
||||
|
||||
typename T::Geom_traits gt;
|
||||
gt = t.geom_traits();
|
||||
|
||||
t_const.tds();
|
||||
t.tds();
|
||||
t_const.domain();
|
||||
t_const.number_of_sheets();
|
||||
t_const.dimension();
|
||||
|
||||
size_t number_of_vertices = t_const.number_of_vertices();
|
||||
CGAL_assertion(number_of_vertices == t.number_of_vertices());
|
||||
size_t number_of_faces = t_const.number_of_faces();
|
||||
CGAL_assertion(number_of_faces == t.number_of_faces());
|
||||
size_t number_of_stored_vertices = t_const.number_of_stored_vertices();
|
||||
CGAL_assertion(number_of_stored_vertices == t.number_of_stored_vertices());
|
||||
size_t number_of_stored_faces = t_const.number_of_stored_faces();
|
||||
CGAL_assertion(number_of_stored_faces == t.number_of_stored_faces());
|
||||
|
||||
size_t number_of_edges = t_const.number_of_edges();
|
||||
CGAL_assertion(number_of_edges == t.number_of_edges());
|
||||
size_t number_of_stored_edges = t_const.number_of_stored_edges();
|
||||
CGAL_assertion(number_of_stored_edges == t.number_of_stored_edges());
|
||||
|
||||
bool ext1 = t_const.is_extensible_triangulation_in_1_sheet_h1();
|
||||
CGAL_assertion(ext1 == t.is_extensible_triangulation_in_1_sheet_h1());
|
||||
bool ext2 = t_const.is_extensible_triangulation_in_1_sheet_h2();
|
||||
CGAL_assertion(ext2 == t.is_extensible_triangulation_in_1_sheet_h2());
|
||||
bool is_triang1 = t_const.is_triangulation_in_1_sheet();
|
||||
CGAL_assertion(is_triang1 == t.is_triangulation_in_1_sheet());
|
||||
t.convert_to_1_sheeted_covering();
|
||||
t.convert_to_9_sheeted_covering();
|
||||
}
|
||||
|
||||
template <class T>
|
||||
void test_geometric_access() {
|
||||
typedef typename T::Point Point;
|
||||
typedef typename T::Segment Segment;
|
||||
typedef typename T::Triangle Triangle;
|
||||
typedef typename T::Periodic_point Periodic_point;
|
||||
typedef typename T::Periodic_segment Periodic_segment;
|
||||
typedef typename T::Periodic_triangle Periodic_triangle;
|
||||
typedef typename T::Vertex_iterator Vertex_iterator;
|
||||
typedef typename T::Face_iterator Face_iterator;
|
||||
|
||||
T t;
|
||||
const T &t_const = t;
|
||||
|
||||
t.insert(Point(0.5, 0.5), Face_handle());
|
||||
t.insert(Point(0.7, 0.5));
|
||||
t.insert(Point(0.7, 0.7));
|
||||
|
||||
Vertex_iterator vit = t.vertices_begin();
|
||||
Face_iterator fit = t.faces_begin();
|
||||
|
||||
Periodic_point pp0 = t_const.periodic_point(vit);
|
||||
Periodic_point pp1 = t_const.periodic_point(fit, 0);
|
||||
CGAL_USE(pp1);
|
||||
Periodic_segment ps0 = t_const.periodic_segment(fit, 0);
|
||||
CGAL_USE(ps0);
|
||||
Periodic_segment ps1 = t_const.periodic_segment(*t.edges_begin());
|
||||
CGAL_USE(ps1);
|
||||
Periodic_triangle pt0 = t_const.periodic_triangle(fit);
|
||||
CGAL_USE(pt0);
|
||||
|
||||
Point p0 = t_const.point(pp0);
|
||||
CGAL_USE(p0);
|
||||
Segment s0 = t_const.segment(ps0);
|
||||
CGAL_USE(s0);
|
||||
Triangle t0 = t_const.triangle(pt0);
|
||||
CGAL_USE(t0);
|
||||
}
|
||||
|
||||
template <class T>
|
||||
void test_predicates() {
|
||||
typedef typename T::Vertex_handle Vertex_handle;
|
||||
typedef typename T::Face_handle Face_handle;
|
||||
|
||||
T t;
|
||||
|
||||
Vertex_handle vh0 = t.insert(Point(0.5, 0.5));
|
||||
Vertex_handle vh1 = t.insert(Point(0.7, 0.5));
|
||||
Vertex_handle vh2 = t.insert(Point(0.7, 0.7));
|
||||
|
||||
t.is_edge(vh0, vh1);
|
||||
Face_handle fh; int i;
|
||||
t.is_edge(vh0, vh1, fh, i);
|
||||
|
||||
t.is_face(vh0, vh1, vh2);
|
||||
t.is_face(vh0, vh1, vh2, fh);
|
||||
}
|
||||
|
||||
template <class T>
|
||||
void test_queries() {
|
||||
typedef typename T::Vertex_handle Vertex_handle;
|
||||
typedef typename T::Face_handle Face_handle;
|
||||
|
||||
T t;
|
||||
const T &t_const = t;
|
||||
CGAL_USE(t_const);
|
||||
|
||||
Point p0(0.5, 0.5);
|
||||
Vertex_handle vh0 = t.insert(p0);
|
||||
CGAL_USE(vh0);
|
||||
Vertex_handle vh1 = t.insert(Point(0.7, 0.5));
|
||||
CGAL_USE(vh1);
|
||||
Vertex_handle vh2 = t.insert(Point(0.7, 0.7));
|
||||
CGAL_USE(vh2);
|
||||
|
||||
Face_handle fh = t_const.locate(p0);
|
||||
fh = t_const.locate(Point(0.5, 0.5), fh);
|
||||
|
||||
typename T::Locate_type lt; int li;
|
||||
fh = t_const.locate(p0, lt, li);
|
||||
fh = t_const.locate(p0, lt, li, fh);
|
||||
|
||||
t_const.oriented_side(fh, p0);
|
||||
}
|
||||
|
||||
template <class T>
|
||||
void test_iterators() {
|
||||
typedef typename T::Vertex_handle Vertex_handle;
|
||||
typedef typename T::Face_handle Face_handle;
|
||||
|
||||
T t;
|
||||
const T &t_const = t;
|
||||
|
||||
Point p0(0.5, 0.5);
|
||||
Vertex_handle vh0 = t.insert(p0);
|
||||
CGAL_USE(vh0);
|
||||
Vertex_handle vh1 = t.insert(Point(0.7, 0.5));
|
||||
CGAL_USE(vh1);
|
||||
Vertex_handle vh2 = t.insert(Point(0.7, 0.7));
|
||||
CGAL_USE(vh2);
|
||||
|
||||
for (typename T::Vertex_iterator vit = t_const.vertices_begin();
|
||||
vit != t_const.vertices_end(); ++vit) {
|
||||
}
|
||||
for (typename T::Edge_iterator eit = t_const.edges_begin();
|
||||
eit != t_const.edges_end(); ++eit) {
|
||||
}
|
||||
for (typename T::Face_iterator fit = t_const.faces_begin();
|
||||
fit != t_const.faces_end(); ++fit) {
|
||||
}
|
||||
|
||||
for (typename T::Periodic_point_iterator ppit = t_const.periodic_points_begin();
|
||||
ppit != t_const.periodic_points_end(); ++ppit) {
|
||||
}
|
||||
for (typename T::Periodic_point_iterator ppit =
|
||||
t_const.periodic_points_begin(T::STORED);
|
||||
ppit != t_const.periodic_points_end(T::STORED); ++ppit) {
|
||||
}
|
||||
for (typename T::Periodic_point_iterator ppit =
|
||||
t_const.periodic_points_begin(T::UNIQUE);
|
||||
ppit != t_const.periodic_points_end(T::UNIQUE); ++ppit) {
|
||||
}
|
||||
for (typename T::Periodic_point_iterator ppit =
|
||||
t_const.periodic_points_begin(T::STORED_COVER_DOMAIN);
|
||||
ppit != t_const.periodic_points_end(T::STORED_COVER_DOMAIN); ++ppit) {
|
||||
}
|
||||
for (typename T::Periodic_point_iterator ppit =
|
||||
t_const.periodic_points_begin(T::UNIQUE_COVER_DOMAIN);
|
||||
ppit != t_const.periodic_points_end(T::UNIQUE_COVER_DOMAIN); ++ppit) {
|
||||
}
|
||||
|
||||
for (typename T::Periodic_segment_iterator psit = t_const.periodic_segments_begin();
|
||||
psit != t_const.periodic_segments_end(); ++psit) {
|
||||
}
|
||||
for (typename T::Periodic_segment_iterator psit =
|
||||
t_const.periodic_segments_begin(T::STORED);
|
||||
psit != t_const.periodic_segments_end(T::STORED); ++psit) {
|
||||
}
|
||||
for (typename T::Periodic_segment_iterator psit =
|
||||
t_const.periodic_segments_begin(T::UNIQUE);
|
||||
psit != t_const.periodic_segments_end(T::UNIQUE); ++psit) {
|
||||
}
|
||||
for (typename T::Periodic_segment_iterator psit =
|
||||
t_const.periodic_segments_begin(T::STORED_COVER_DOMAIN);
|
||||
psit != t_const.periodic_segments_end(T::STORED_COVER_DOMAIN); ++psit) {
|
||||
}
|
||||
for (typename T::Periodic_segment_iterator psit =
|
||||
t_const.periodic_segments_begin(T::UNIQUE_COVER_DOMAIN);
|
||||
psit != t_const.periodic_segments_end(T::UNIQUE_COVER_DOMAIN); ++psit) {
|
||||
}
|
||||
|
||||
for (typename T::Periodic_triangle_iterator ptit = t_const.periodic_triangles_begin();
|
||||
ptit != t_const.periodic_triangles_end(); ++ptit) {
|
||||
}
|
||||
for (typename T::Periodic_triangle_iterator ptit =
|
||||
t_const.periodic_triangles_begin(T::STORED);
|
||||
ptit != t_const.periodic_triangles_end(T::STORED); ++ptit) {
|
||||
}
|
||||
for (typename T::Periodic_triangle_iterator ptit =
|
||||
t_const.periodic_triangles_begin(T::UNIQUE);
|
||||
ptit != t_const.periodic_triangles_end(T::UNIQUE); ++ptit) {
|
||||
}
|
||||
for (typename T::Periodic_triangle_iterator ptit =
|
||||
t_const.periodic_triangles_begin(T::STORED_COVER_DOMAIN);
|
||||
ptit != t_const.periodic_triangles_end(T::STORED_COVER_DOMAIN); ++ptit) {
|
||||
}
|
||||
for (typename T::Periodic_triangle_iterator ptit =
|
||||
t_const.periodic_triangles_begin(T::UNIQUE_COVER_DOMAIN);
|
||||
ptit != t_const.periodic_triangles_end(T::UNIQUE_COVER_DOMAIN); ++ptit) {
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
template <class T>
|
||||
void test_circulators() {
|
||||
typedef typename T::Vertex_handle Vertex_handle;
|
||||
typedef typename T::Face_handle Face_handle;
|
||||
|
||||
T t;
|
||||
const T &t_const = t;
|
||||
|
||||
Point p0(0.5, 0.5);
|
||||
Vertex_handle vh0 = t.insert(p0);
|
||||
CGAL_USE(vh0);
|
||||
Vertex_handle vh1 = t.insert(Point(0.7, 0.5));
|
||||
CGAL_USE(vh1);
|
||||
Vertex_handle vh2 = t.insert(Point(0.7, 0.7));
|
||||
CGAL_USE(vh2);
|
||||
|
||||
typename T::Face_circulator fcir = t_const.incident_faces(vh0);
|
||||
fcir = t_const.incident_faces(vh0, fcir);
|
||||
|
||||
typename T::Edge_circulator ecir = t_const.incident_edges(vh0);
|
||||
ecir = t_const.incident_edges(vh0, fcir);
|
||||
|
||||
typename T::Vertex_circulator vcir = t_const.adjacent_vertices(vh0);
|
||||
vcir = t_const.adjacent_vertices(vh0, fcir);
|
||||
|
||||
Vertex_handle v_mirror = t_const.mirror_vertex(fcir, 0);
|
||||
CGAL_USE(v_mirror);
|
||||
t_const.mirror_index(fcir, 0);
|
||||
}
|
||||
|
||||
template <class T>
|
||||
void test_modifiers() {
|
||||
typedef typename T::Vertex_handle Vertex_handle;
|
||||
typedef typename T::Face_handle Face_handle;
|
||||
|
||||
T t;
|
||||
const T &t_const = t;
|
||||
|
||||
Point p0(0.5, 0.5);
|
||||
Point p1(0.8, 0.6);
|
||||
Point p2(0.7, 0.7);
|
||||
|
||||
Vertex_handle vh0 = t.insert(p0);
|
||||
Face_handle fh = t.faces_begin();
|
||||
Vertex_handle vh1 = t.insert(Point(0.7, 0.5), fh);
|
||||
Vertex_handle vh2 = t.insert(Point(0.7, 0.7));
|
||||
|
||||
|
||||
t.flip(fh, 0);
|
||||
|
||||
t.clear();
|
||||
vh0 = t.insert(p0);
|
||||
vh1 = t.insert(p1);
|
||||
vh2 = t.insert(p2);
|
||||
|
||||
typename T::Locate_type lt; int li;
|
||||
fh = t_const.locate(p0, lt, li);
|
||||
t.insert(p0, lt, fh, li);
|
||||
t.push_back(p0);
|
||||
|
||||
t.clear();
|
||||
std::vector<Point> pts;
|
||||
pts.push_back(p0);
|
||||
pts.push_back(p1);
|
||||
pts.push_back(p0);
|
||||
pts.push_back(p2);
|
||||
t.insert(pts.begin(), pts.end());
|
||||
|
||||
t.clear();
|
||||
vh0 = t.insert(p0);
|
||||
t.remove_first(vh0);
|
||||
|
||||
t.clear();
|
||||
vh0 = t.insert_first(p0);
|
||||
vh1 = t.insert(p1);
|
||||
if (t.degree(vh1) == 3) {
|
||||
// The vertex has degree 6 in case we are testing the Delaunay triangulation
|
||||
t.remove_degree_3(vh1);
|
||||
}
|
||||
|
||||
t.clear();
|
||||
vh0 = t.insert_first(p0);
|
||||
|
||||
fh = t_const.locate(p1, lt, li);
|
||||
CGAL_assertion(lt == T::FACE);
|
||||
t.insert_in_face(p1, fh);
|
||||
fh = t_const.locate(p2, lt, li);
|
||||
CGAL_assertion(lt == T::EDGE);
|
||||
t.insert_in_edge(p2, fh, li);
|
||||
|
||||
t.clear();
|
||||
t.insert(pts.begin(), pts.end());
|
||||
for (typename T::Vertex_iterator vit = t_const.vertices_begin();
|
||||
vit != t_const.vertices_end(); ++vit) {
|
||||
if (t_const.degree(vit) == 3) {
|
||||
t.remove_degree_3(vit);
|
||||
vit = t_const.vertices_begin();
|
||||
}
|
||||
}
|
||||
|
||||
t.clear();
|
||||
vh0 = t.insert(p0);
|
||||
t.remove_first(vh0);
|
||||
|
||||
// star_hole is not tested
|
||||
}
|
||||
|
||||
template <class T>
|
||||
void test_miscellaneous() {
|
||||
typedef typename T::Vertex_handle Vertex_handle;
|
||||
typedef typename T::Face_handle Face_handle;
|
||||
|
||||
T t;
|
||||
|
||||
Point p0(0.5, 0.5);
|
||||
Point p1(0.8, 0.6);
|
||||
Point p2(0.7, 0.7);
|
||||
Vertex_handle vh0, vh1, vh2;
|
||||
vh0 = t.insert(p0);
|
||||
vh1 = t.insert(p1);
|
||||
vh2 = t.insert(p2);
|
||||
|
||||
t.set_domain(typename T::Iso_rectangle(0,0,2,2));
|
||||
int i = t.ccw(0);
|
||||
int j = t.cw(0);
|
||||
CGAL_assertion(i+j == 3);
|
||||
|
||||
t = T();
|
||||
vh0 = t.insert(p0);
|
||||
vh1 = t.insert(p1);
|
||||
vh2 = t.insert(p2);
|
||||
Face_handle fh = t.faces_begin();
|
||||
t.flippable(fh, 0);
|
||||
t.degree(vh0);
|
||||
|
||||
t.is_valid();
|
||||
t.is_valid(true);
|
||||
t.is_valid(false);
|
||||
t.is_valid(true, 0);
|
||||
t.is_valid(false, 0);
|
||||
}
|
||||
|
||||
template <class T>
|
||||
void test_io(T &pt1) {
|
||||
bool ex = false; // Exact predicates
|
||||
|
||||
std::cout << "I/O" << std::endl;
|
||||
std::cout << " ascii" << std::endl;
|
||||
|
||||
std::stringstream ss1;
|
||||
ss1 << pt1;
|
||||
|
||||
T pt1r;
|
||||
ss1 >> pt1r;
|
||||
|
||||
assert(CGAL::is_ascii(ss1));
|
||||
if (!ex) assert(pt1 == pt1r);
|
||||
|
||||
std::cout << " binary" << std::endl;
|
||||
pt1r.clear();
|
||||
// There are problems with the IO of exact number types in binary mode.
|
||||
if (!ex) {
|
||||
std::stringstream ss1b;
|
||||
CGAL::set_binary_mode(ss1b);
|
||||
ss1b << pt1;
|
||||
|
||||
ss1b >> pt1r;
|
||||
assert(CGAL::is_binary(ss1b));
|
||||
|
||||
assert(pt1 == pt1r);
|
||||
}
|
||||
|
||||
std::cout << " pretty" << std::endl;
|
||||
|
||||
pt1r.clear();
|
||||
std::stringstream ss1p;
|
||||
CGAL::set_pretty_mode(ss1p);
|
||||
ss1p << pt1;
|
||||
|
||||
assert(CGAL::is_pretty(ss1p));
|
||||
}
|
||||
|
||||
template <class T>
|
||||
void test_io() {
|
||||
T t;
|
||||
test_io(t);
|
||||
|
||||
t.insert(Point(0.5, 0.5));
|
||||
test_io(t);
|
||||
|
||||
t.insert(Point(0.6, 0.8));
|
||||
t.insert(Point(0.2, 0.6));
|
||||
test_io(t);
|
||||
|
||||
// One cover for the Delaunay triangulation
|
||||
t.clear();
|
||||
for (int x=0; x<5; ++x)
|
||||
for (int y=0; y<5; ++y)
|
||||
t.insert(Point(x/5.0, y/5.0));
|
||||
test_io(t);
|
||||
}
|
||||
|
||||
template <class T>
|
||||
void test() {
|
||||
test_constructor<T>();
|
||||
test_global_access<T>();
|
||||
test_geometric_access<T>();
|
||||
test_predicates<T>();
|
||||
test_queries<T>();
|
||||
test_iterators<T>();
|
||||
test_circulators<T>();
|
||||
test_modifiers<T>();
|
||||
test_miscellaneous<T>();
|
||||
test_io<T>();
|
||||
}
|
||||
|
||||
|
||||
template <class T>
|
||||
void test_nearest() {
|
||||
typedef typename T::Vertex_handle Vertex_handle;
|
||||
|
||||
Point p0(0.5, 0.5);
|
||||
Point p1(0.8, 0.6);
|
||||
Point p2(0.7, 0.7);
|
||||
Vertex_handle vh0, vh1, vh2;
|
||||
|
||||
T t;
|
||||
CGAL_assertion(t.nearest_vertex(p0) == Vertex_handle());
|
||||
|
||||
vh0 = t.insert(p0);
|
||||
CGAL_assertion(t.get_original_vertex(t.nearest_vertex(p0)) == vh0);
|
||||
CGAL_assertion(t.get_original_vertex(t.nearest_vertex(p1)) == vh0);
|
||||
CGAL_assertion(t.get_original_vertex(t.nearest_vertex(p2)) == vh0);
|
||||
|
||||
vh1 = t.insert(p1);
|
||||
vh2 = t.insert(p2);
|
||||
CGAL_assertion(t.get_original_vertex(t.nearest_vertex(p0)) == vh0);
|
||||
CGAL_assertion(t.get_original_vertex(t.nearest_vertex(p1)) == vh1);
|
||||
CGAL_assertion(t.get_original_vertex(t.nearest_vertex(p2)) == vh2);
|
||||
|
||||
}
|
||||
|
||||
int main() {
|
||||
typedef Periodic_2_triangulation_2<Gt> P2T2;
|
||||
typedef Periodic_2_Delaunay_triangulation_2<Gt> DP2T2;
|
||||
|
|
|
|||
Loading…
Reference in New Issue