songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

extract polylines from example and add corresponding data file

This commit is contained in:
Jane Tournois 2017-01-09 12:03:37 +01:00
parent 89b0f864d4
commit f6bb4822fb
2 changed files with 25 additions and 41 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
Mesh_3/examples/Mesh_3/data/hybrid_example.polylines.txt Executable file
View File

@ -0,0 +1,15 @@
2 -1 -1 -1 1 -1 -1
2 -1 -1 -1 -1 1 -1
2 -1 -1 -1 -1 -1 1
2 1 -1 -1 1 1 -1
2 1 -1 -1 1 -1 1
2 -1 1 -1 1 1 -1
2 -1 1 -1 -1 1 1
2 -1 -1 1 1 -1 1
2 -1 -1 1 -1 1 1
2 1 1 -1 1 1 0
2 -1 1 1 0 1 1
2 1 -1 1 1 0 1
21 0 1 1 0.05 0.68775 1 0.1 0.56411 1 0.15 0.473217 1 0.2 0.4 1 0.25 0.338562 1 0.3 0.285857 1 0.35 0.240066 1 0.4 0.2 1 0.45 0.164835 1 0.5 0.133975 1 0.55 0.106971 1 0.6 0.0834849 1 0.65 0.0632503 1 0.7 0.0460608 1 0.75 0.0317542 1 0.8 0.0202041 1 0.85 0.011314 1 0.9 0.00501256 1 0.95 0.00125078 1 1 0 1
21 1 1 0 1 0.68775 0.05 1 0.56411 0.1 1 0.473217 0.15 1 0.4 0.2 1 0.338562 0.25 1 0.285857 0.3 1 0.240066 0.35 1 0.2 0.4 1 0.164835 0.45 1 0.133975 0.5 1 0.106971 0.55 1 0.0834849 0.6 1 0.0632503 0.65 1 0.0460608 0.7 1 0.0317542 0.75 1 0.0202041 0.8 1 0.011314 0.85 1 0.00501256 0.9 1 0.00125078 0.95 1 0 1
21 0 1 1 0.05 1 0.68775 0.1 1 0.56411 0.15 1 0.473217 0.2 1 0.4 0.25 1 0.338562 0.3 1 0.285857 0.35 1 0.240066 0.4 1 0.2 0.45 1 0.164835 0.5 1 0.133975 0.55 1 0.106971 0.6 1 0.0834849 0.65 1 0.0632503 0.7 1 0.0460608 0.75 1 0.0317542 0.8 1 0.0202041 0.85 1 0.011314 0.9 1 0.00501256 0.95 1 0.00125078 1 1 0

51
Mesh_3/examples/Mesh_3/mesh_hybrid_mesh_domain.cpp
View File

@ -17,6 +17,9 @@
// Ouput
#include <CGAL/Mesh_3/Dump_c3t3.h>
// Read 1D features from input file
#include "read_polylines.h"
// Sphere Domain
typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef K::FT FT;
@ -215,50 +218,16 @@ int main()
Domain domain(sphere_domain, polyhedron_domain);
// Read polyline features
const char* lines_fname = "data/hybrid_example.polylines.txt";
std::vector<std::vector<Point> > featured_curves;
featured_curves.resize(15);
featured_curves[ 0].push_back(Point(-1, -1, -1));
featured_curves[ 0].push_back(Point( 1, -1, -1));
featured_curves[ 1].push_back(Point(-1, -1, -1));
featured_curves[ 1].push_back(Point(-1, 1, -1));
featured_curves[ 2].push_back(Point(-1, -1, -1));
featured_curves[ 2].push_back(Point(-1, -1, 1));
featured_curves[ 3].push_back(Point( 1, -1, -1));
featured_curves[ 3].push_back(Point( 1, 1, -1));
featured_curves[ 4].push_back(Point( 1, -1, -1));
featured_curves[ 4].push_back(Point( 1, -1, 1));
featured_curves[ 5].push_back(Point(-1, 1, -1));
featured_curves[ 5].push_back(Point( 1, 1, -1));
featured_curves[ 6].push_back(Point(-1, 1, -1));
featured_curves[ 6].push_back(Point(-1, 1, 1));
featured_curves[ 7].push_back(Point(-1, -1, 1));
featured_curves[ 7].push_back(Point( 1, -1, 1));
featured_curves[ 8].push_back(Point(-1, -1, 1));
featured_curves[ 8].push_back(Point(-1, 1, 1));
featured_curves[ 9].push_back(Point( 1, 1, -1));
featured_curves[ 9].push_back(Point( 1, 1, 0));
featured_curves[10].push_back(Point(-1, 1, 1));
featured_curves[10].push_back(Point( 0, 1, 1));
featured_curves[11].push_back(Point( 1, -1, 1));
featured_curves[11].push_back(Point( 1, 0, 1));
featured_curves[12].push_back(Point(0, 1, 1));
featured_curves[13].push_back(Point(1, 1, 0));
featured_curves[14].push_back(Point(0, 1, 1));
for(double x = 0.05; x < 0.98; x+=0.05) {
const double y = 1 - CGAL::sqrt(1 - CGAL::square(x - 1));
featured_curves[12].push_back(Point(x, y, 1));
featured_curves[13].push_back(Point(1, y, x));
featured_curves[14].push_back(Point(x, 1, y));
if (!read_polylines(lines_fname, featured_curves))
{ // see file "read_polylines.h"
std::cerr << "Error: Cannot read file " << lines_fname << std::endl;
return false;
}
featured_curves[12].push_back(Point(1, 0, 1));
featured_curves[13].push_back(Point(1, 0, 1));
featured_curves[14].push_back(Point(1, 1, 0));
// Add features for protection
domain.add_features(featured_curves.begin(), featured_curves.end());
// Criteria