Add the sandbox example for cgal dependent mesh boolean

3 files changed, 91 insertions, 2 deletions

diff --git a/sandboxes/CMakeLists.txt b/sandboxes/CMakeLists.txt
index 3372698c3..91d6ca225 100644
--- a/sandboxes/CMakeLists.txt
+++ b/sandboxes/CMakeLists.txt
@@ -1,2 +1,3 @@
-add_subdirectory(slasupporttree)
-add_subdirectory(openvdb)
+#add_subdirectory(slasupporttree)
+#add_subdirectory(openvdb)
+add_subdirectory(meshboolean)
diff --git a/sandboxes/meshboolean/CMakeLists.txt b/sandboxes/meshboolean/CMakeLists.txt
new file mode 100644
index 000000000..3306d7ed0
--- /dev/null
+++ b/sandboxes/meshboolean/CMakeLists.txt
@@ -0,0 +1,5 @@
+
+find_package(CGAL REQUIRED)
+add_executable(meshboolean MeshBoolean.cpp)
+
+target_link_libraries(meshboolean libslic3r CGAL::CGAL)
\ No newline at end of file
diff --git a/sandboxes/meshboolean/MeshBoolean.cpp b/sandboxes/meshboolean/MeshBoolean.cpp
new file mode 100644
index 000000000..f10c54d7f
--- /dev/null
+++ b/sandboxes/meshboolean/MeshBoolean.cpp
@@ -0,0 +1,83 @@
+#include <libslic3r/TriangleMesh.hpp>
+#undef PI
+#include <igl/readOFF.h>
+//#undef IGL_STATIC_LIBRARY
+#include <igl/copyleft/cgal/mesh_boolean.h>
+
+#include <Eigen/Core>
+#include <iostream>
+
+#include <admesh/stl.h>
+
+#include <boost/nowide/cstdio.hpp>
+#include <boost/log/trivial.hpp>
+
+namespace Slic3r {
+
+bool its_write_obj(const Eigen::MatrixXd &V, Eigen::MatrixXi &F, const char *file)
+{
+
+  	FILE *fp = boost::nowide::fopen(file, "w");
+  	if (fp == nullptr) {
+		BOOST_LOG_TRIVIAL(error) << "stl_write_obj: Couldn't open " << file << " for writing";
+    	return false;
+  	}
+
+	for (size_t i = 0; i < V.rows(); ++ i)
+    	fprintf(fp, "v %lf %lf %lf\n", V(i, 0), V(i, 1), V(i, 2));
+  	for (size_t i = 0; i < F.rows(); ++ i)
+    	fprintf(fp, "f %d %d %d\n", F(i, 0) + 1, F(i, 1) + 1, F(i, 2) + 1);
+  	fclose(fp);
+  	return true;
+}
+
+void mesh_boolean_test(const std::string &fname)
+{
+  using namespace Eigen;
+  using namespace std;
+//  igl::readOFF(TUTORIAL_SHARED_PATH "/cheburashka.off",VA,FA);
+//  igl::readOFF(TUTORIAL_SHARED_PATH "/decimated-knight.off",VB,FB);
+  // Plot the mesh with pseudocolors
+//  igl::opengl::glfw::Viewer viewer;
+
+  // Initialize
+//  update(viewer);
+
+  //igl::copyleft::cgal::mesh_boolean(VA,FA,VB,FB,boolean_type,VC,FC,J);
+
+
+	TriangleMesh mesh;
+	mesh.ReadSTLFile(fname.c_str());
+	mesh.repair(true);
+    its_write_obj(mesh.its, (fname + "-imported0.obj").c_str());
+
+	Eigen::MatrixXd VA,VB,VC;
+	Eigen::VectorXi J,I;
+	Eigen::MatrixXi FA,FB,FC;
+	igl::MeshBooleanType boolean_type(igl::MESH_BOOLEAN_TYPE_UNION);
+
+
+  	typedef Eigen::Map<const Eigen::Matrix<float, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor | Eigen::DontAlign>> MapMatrixXfUnaligned;
+    typedef Eigen::Map<const Eigen::Matrix<int,   Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor | Eigen::DontAlign>> MapMatrixXiUnaligned;
+
+	Eigen::MatrixXd V = MapMatrixXfUnaligned(mesh.its.vertices.front().data(), mesh.its.vertices.size(), 3).cast<double>();
+    Eigen::MatrixXi F = MapMatrixXiUnaligned(mesh.its.indices.front().data(), mesh.its.indices.size(), 3);
+
+    its_write_obj(V, F, (fname + "-imported.obj").c_str());
+
+    // Self-union to clean up 
+    igl::copyleft::cgal::mesh_boolean(V, F, Eigen::MatrixXd(), Eigen::MatrixXi(), boolean_type, VC, FC);
+
+    its_write_obj(VC, FC, (fname + "-fixed.obj").c_str());
+}
+
+} // namespace Slic3r
+
+int main(const int argc, const char * argv[])
+{
+    if (argc < 1) return -1;
+    
+    Slic3r::mesh_boolean_test(argv[1]);
+    
+    return 0;
+}