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

example implicit sphere

This commit is contained in:
Pierre Alliez 2023-12-23 19:25:43 +01:00
parent f782daba86
commit 5239cac781
1 changed files with 7 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
Isosurfacing_3/examples/Isosurfacing_3/marching_cubes_implicit_sphere.cpp
View File

@ -24,25 +24,24 @@ int main(int, char**)
{ {
// box domain and spacing vector // box domain and spacing vector
const CGAL::Bbox_3 bbox{ -1.0, -1.0, -1.0, 1.0, 1.0, 1.0 }; const CGAL::Bbox_3 bbox{ -1.0, -1.0, -1.0, 1.0, 1.0, 1.0 };
const FT spacing = 0.04; const FT spacing = 0.05;
const Vector vec_spacing(spacing, spacing, spacing); const Vector vec_spacing(spacing, spacing, spacing);
// create domain with sphere function // create domain with sphere function
auto domain = CGAL::Isosurfacing::create_implicit_Cartesian_grid_domain<Kernel>(bbox, vec_spacing, sphere_function); auto domain = CGAL::Isosurfacing::create_implicit_Cartesian_grid_domain<Kernel>
(bbox, vec_spacing, sphere_function);
// points and triangles for the output indexed mesh // points and triangles for the output indexed mesh
Point_range points; Point_range points;
Triangle_range triangles; Triangle_range triangles;
const tbb::tick_count start = tbb::tick_count::now(); // run marching cubes with given isovalue
std::cout << "marching cubes...";
// execute marching cubes with a given isovalue
const FT isovalue = 0.8; const FT isovalue = 0.8;
const tbb::tick_count start = tbb::tick_count::now();
CGAL::Isosurfacing::marching_cubes(domain, isovalue, points, triangles); CGAL::Isosurfacing::marching_cubes(domain, isovalue, points, triangles);
const tbb::tick_count end = tbb::tick_count::now(); const tbb::tick_count end = tbb::tick_count::now();
std::cout << "done (" << (end - start).seconds() << "s)" << std::endl;
std::cout << (end - start).seconds() << std::endl;
// save ouput indexed mesh to a file, in the OFF format // save ouput indexed mesh to a file, in the OFF format
CGAL::IO::write_OFF("output.off", points, triangles); CGAL::IO::write_OFF("output.off", points, triangles);