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

Example DC on octree (still uses the octree wrapper)

This commit is contained in:
Pierre Alliez 2024-01-05 14:10:41 +01:00
parent d2c94538ae
commit 26848f3d67
1 changed files with 5 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
Isosurfacing_3/examples/Isosurfacing_3/dual_contouring_octree.cpp
View File

@ -47,7 +47,6 @@ struct Refine_one_eighth
bool operator()(const Octree_wrapper::Octree::Node& n) const bool operator()(const Octree_wrapper::Octree::Node& n) const
{ {
// n.depth()
if(n.depth() < min_depth_) if(n.depth() < min_depth_)
return true; return true;
@ -62,30 +61,27 @@ struct Refine_one_eighth
if(leaf_coords[1] >= octree_dim_ / 2) if(leaf_coords[1] >= octree_dim_ / 2)
return false; return false;
// if(leaf_coords[2] >= octree_dim_ / 2)
// return false;
return true; return true;
} }
}; };
int main(int, char**) int main(int, char**)
{ {
const CGAL::Bbox_3 bbox{-1., -1., -1., 1., 1., 1.}; const CGAL::Bbox_3 bbox{-1., -1., -1., 1., 1., 1.};
Octree_wrapper octree_wrap { bbox }; Octree_wrapper octree_wrap { bbox };
Refine_one_eighth split_predicate(3, 4); Refine_one_eighth split_predicate(3, 5);
octree_wrap.refine(split_predicate); octree_wrap.refine(split_predicate);
auto sphere_function = [&](const Point& p) -> FT auto sphere_function = [&](const Point& p) -> FT
{ {
return sqrt(p.x() * p.x() + p.y() * p.y() + p.z() * p.z()); return std::sqrt(p.x() * p.x() + p.y() * p.y() + p.z() * p.z());
}; };
auto sphere_gradient = [&](const Point& p) -> Vector auto sphere_gradient = [&](const Point& p) -> Vector
{ {
const Vector g = p - CGAL::ORIGIN; const Vector g = p - CGAL::ORIGIN;
return g / sqrt(g.squared_length()); return g / std::sqrt(g.squared_length());
}; };
auto domain = CGAL::Isosurfacing::create_implicit_octree_domain(octree_wrap, auto domain = CGAL::Isosurfacing::create_implicit_octree_domain(octree_wrap,
@ -97,7 +93,7 @@ int main(int, char**)
CGAL::Isosurfacing::dual_contouring(domain, 0.8, points, polygons); CGAL::Isosurfacing::dual_contouring(domain, 0.8, points, polygons);
CGAL::IO::write_OFF("result.off", points, polygons); CGAL::IO::write_OFF("output.off", points, polygons);
return EXIT_SUCCESS; return EXIT_SUCCESS;
} }