From a62aefa7ff1143a1a35e5fd4ce2e1c8ff5235892 Mon Sep 17 00:00:00 2001
From: Thomas Krijnen <t.krijnen@gmail.com>
Date: Mon, 11 Dec 2023 21:00:04 +0100
Subject: [PATCH] Disable automatic density substitution in cgal kernel

---
 src/ifcgeom/kernels/cgal/CgalKernel.cpp | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/ifcgeom/kernels/cgal/CgalKernel.cpp b/src/ifcgeom/kernels/cgal/CgalKernel.cpp
index 4fa2f456db..cfb6cabce0 100644
--- a/src/ifcgeom/kernels/cgal/CgalKernel.cpp
+++ b/src/ifcgeom/kernels/cgal/CgalKernel.cpp
@@ -153,7 +153,7 @@ CGAL::Nef_polyhedron_3<Kernel_> ifcopenshell::geometry::utils::create_nef_polyhe
 #endif
 
 bool CgalKernel::convert(const taxonomy::shell::ptr l, cgal_shape_t& shape) {
-	if (l->children.size() > 100) {
+	if (false && l->children.size() > 100) {
 		static double inf = 1.e9; //  std::numeric_limits<double>::infinity();
 		std::pair<Eigen::Vector3d, Eigen::Vector3d> minmax(
 			Eigen::Vector3d(+inf, +inf, +inf),
@@ -174,6 +174,7 @@ bool CgalKernel::convert(const taxonomy::shell::ptr l, cgal_shape_t& shape) {
 		});
 		auto diag = minmax.second - minmax.first;
 		double volume = diag(0) * diag(1) * diag(2);
+		// @todo volume van be zero also..
 		double density = num_points / volume;
 		Logger::Notice("Density " + boost::lexical_cast<std::string>(density), l->instance);
 		if (density > 5000) {