Files
OrcaSlicer/deps_src/libigl/igl/copyleft/cgal/points_inside_component.cpp
Donovan Baarda dc5897d7b5 Update eigen to v5.0.1 and libigl to v2.6.0. (#11311)
* Update eigen from v3.3.7 to v5.0.1.

This updates eigen from v3.3.7 released on  December 11, 2018-12-11 to v5.0.1
released on 2025-11-11. There have be a large number of bug-fixes,
optimizations, and improvements between these releases. See the details at;

https://gitlab.com/libeigen/eigen/-/releases

It retains the previous custom minimal `CMakeLists.txt`, and adds a
README-OrcaSlicer.md that explains what version and parts of the upstream
eigen release have been included, and where the full release can be found.

* Update libigl from v2.0.0 (or older) to v2.6.0.

This updates libigl from what was probably v2.0.0 released on 2018-10-16 to
v2.6.0 released on 2025-05-15. It's possible the old version was even older
than that but there is no version indicators in the code and I ran out of
patience identifying missing changes and only went back as far as v2.0.0.

There have been a large number of bug-fixes, optimizations, and improvements
between these versions. See the following for details;

https://github.com/libigl/libigl/releases

I retained the minimal custom `CMakeLists.txt`, added `README.md` from the
libigl distribution which identifies the version, and added a
README-OrcaSlicer.md that details the version and parts that have been
included.

* Update libslic3r for libigl v2.6.0 changes.

This updates libslic3r for all changes moving to eigen v5.0.1 and libigl
v2.6.0. Despite the large number of updates to both dependencies, no changes
were required for the eigen update, and only one change was required for the
libigl update.

For libigl, `igl::Hit` was changed to a template taking the Scalar type to
use. Previously it was hard-coded to `float`, so to minimize possible impact
I've updated all places it is used from `igl::Hit` to `igl::Hit<float>`.

* Add compiler option `-DNOMINMAX` for libigl with MSVC.

MSVC by default defines `min(()` and `max()` macros that break
`std::numeric_limits<>::max()`. The upstream cmake that we don't include
adds `-DNOMINMAX` for the libigl module when compiling with MSVC, so we need
to add the same thing here.

* Fix src/libslic3r/TriangleMeshDeal.cpp for the unmodified upstream libigl.

This fixes `TriangleMeshDeal.cpp` to work with the unmodified upstream
libigl v2.6.0. loop.{h,cpp} implementation.

This file and feature was added in PR "BBS Port: Mesh Subdivision" (#12150)
which included changes to `loop.{h,cpp}` in the old version of libigl. This PR
avoids modifying the included dependencies, and uses the updated upstream
versions of those files without any modifications, which requires fixing
TriangleMeshDeal.cpp to work with them.

In particular, the modifications made to `loop.{h,cpp}` included changing the
return type from void to bool, adding additional validation checking of the
input meshes, and returning false if they failed validation. These added
checks looked unnecessary and would only have caught problems if the input
mesh was very corrupt.

To make `TriangleMeshDeal.cpp` work without this built-in checking
functionality, I removed checking/handling of any `false` return value.

There was also a hell of a lot of redundant copying and casting back and forth
between float and double, so I cleaned that up. The input and output meshs use
floats for the vertexes, and there would be no accuracy benefits from casting
to and from doubles for the simple weighted average operations done by
igl::loop(). So this just uses `Eigen:Map` to use the original input mesh
vertex data directly without requiring any copy or casting.

* Move eigen from included `deps_src` to externaly fetched `deps`.

This copys what PrusaSlicer did and moved it from an included dependency under
`deps_src` to an externaly fetched dependency under `deps`. This requires
updating some `CMakeList.txt` configs and removing the old and obsolete
`cmake/modules/FindEigen3.cmake`. The details of when this was done in
PrusaSlicer and the followup fixes are at;

* 21116995d7
* https://github.com/prusa3d/PrusaSlicer/issues/13608
* https://github.com/prusa3d/PrusaSlicer/pull/13609
* e3c277b9ee

For some reason I don't fully understand this also required fixing
`src/slic3r/GUI/GUI_App.cpp` by adding `#include <boost/nowide/cstdio.hpp>` to
fix an `error: ‘remove’ is not a member of ‘boost::nowide'`. The main thing I
don't understand is how it worked before. Note that this include is in the
PrusaSlicer version of this file, but it also significantly deviates from what
is currently in OrcaSlicer in many other ways.

* Whups... I missed adding the deps/Eigen/Eigen.cmake file...

* Tidy some whitespace indenting in CMakeLists.txt.

* Ugh... tabs indenting needing fixes.

* Change the include order of deps/Eigen.

It turns out that although Boost includes some references to Eigen, Eigen also
includes some references to Boost for supporting some of it's additional
numeric types.

I don't think it matters much since we are not using these features, but I
think technically its more correct to say Eigen depends on Boost than the
other way around, so I've re-ordered them.

* Add source for Eigen 5.0.1 download to flatpak yml config.

* Add explicit `DEPENDS dep_Boost to deps/Eigen.

I missed this before. This ensures we don't rely on include orders to make
sure Boost is installed before we configure Eigen.

* Add `DEPENDS dep_Boost dep_GMP dep_MPFR` to deps/Eigen.

It turns out Eigen can also use GMP and MPFR for multi-precision and
multi-precision-rounded numeric types if they are available.

Again, I don't think we are using these so it doesn't really matter, but it is
technically correct and ensures they are there if we ever do need them.

* Fix deps DEPENDENCY ordering for GMP, MPFR, Eigen, and CGAL.

I think this is finally correct. Apparently CGAL also optionally depends on
Eigen, so the correct dependency order from lowest to highest is GMP, MPFR, Eigen, and CGAL.

---------

Co-authored-by: Donovan Baarda <dbaarda@google.com>
Co-authored-by: Noisyfox <timemanager.rick@gmail.com>
2026-05-12 15:09:13 +08:00

357 lines
20 KiB
C++

// This file is part of libigl, a simple c++ geometry processing library.
//
// Copyright (C) 2015 Qingnan Zhou <qnzhou@gmail.com>
//
// This Source Code Form is subject to the terms of the Mozilla Public License
// v. 2.0. If a copy of the MPL was not distributed with this file, You can
// obtain one at http://mozilla.org/MPL/2.0/.
#include "points_inside_component.h"
#include "../../LinSpaced.h"
#include "../../parallel_for.h"
#include "../../PlainMatrix.h"
#include "order_facets_around_edge.h"
#include "assign_scalar.h"
#include <CGAL/AABB_tree.h>
#include <CGAL/AABB_traits.h>
#include <CGAL/AABB_triangle_primitive.h>
#include <CGAL/Exact_predicates_exact_constructions_kernel.h>
#include <cassert>
#include <list>
#include <limits>
#include <vector>
namespace igl {
namespace copyleft
{
namespace cgal {
namespace points_inside_component_helper {
typedef CGAL::Exact_predicates_exact_constructions_kernel Kernel;
typedef Kernel::Ray_3 Ray_3;
typedef Kernel::Point_3 Point_3;
typedef Kernel::Vector_3 Vector_3;
typedef Kernel::Triangle_3 Triangle;
typedef Kernel::Plane_3 Plane_3;
typedef std::vector<Triangle>::iterator Iterator;
typedef CGAL::AABB_triangle_primitive<Kernel, Iterator> Primitive;
typedef CGAL::AABB_traits<Kernel, Primitive> AABB_triangle_traits;
typedef CGAL::AABB_tree<AABB_triangle_traits> Tree;
template<typename DerivedF, typename DerivedI>
void extract_adj_faces(
const Eigen::MatrixBase<DerivedF>& F,
const Eigen::MatrixBase<DerivedI>& I,
const size_t s, const size_t d,
std::vector<int>& adj_faces) {
const size_t num_faces = I.rows();
for (size_t i=0; i<num_faces; i++) {
Eigen::Matrix<typename DerivedF::Scalar,3,1> f = F.row(I(i, 0));
if (((size_t)f[0] == s && (size_t)f[1] == d) ||
((size_t)f[1] == s && (size_t)f[2] == d) ||
((size_t)f[2] == s && (size_t)f[0] == d)) {
adj_faces.push_back((I(i, 0)+1) * -1);
continue;
}
if (((size_t)f[0] == d && (size_t)f[1] == s) ||
((size_t)f[1] == d && (size_t)f[2] == s) ||
((size_t)f[2] == d && (size_t)f[0] == s)) {
adj_faces.push_back(I(i, 0)+1);
continue;
}
}
}
template<typename DerivedF, typename DerivedI>
void extract_adj_vertices(
const Eigen::MatrixBase<DerivedF>& F,
const Eigen::MatrixBase<DerivedI>& I,
const size_t v, std::vector<int>& adj_vertices) {
std::set<size_t> unique_adj_vertices;
const size_t num_faces = I.rows();
for (size_t i=0; i<num_faces; i++) {
Eigen::Matrix<typename DerivedF::Scalar,3,1> f = F.row(I(i, 0));
if ((size_t)f[0] == v) {
unique_adj_vertices.insert(f[1]);
unique_adj_vertices.insert(f[2]);
} else if ((size_t)f[1] == v) {
unique_adj_vertices.insert(f[0]);
unique_adj_vertices.insert(f[2]);
} else if ((size_t)f[2] == v) {
unique_adj_vertices.insert(f[0]);
unique_adj_vertices.insert(f[1]);
}
}
adj_vertices.resize(unique_adj_vertices.size());
std::copy(unique_adj_vertices.begin(),
unique_adj_vertices.end(),
adj_vertices.begin());
}
template<typename DerivedV, typename DerivedF, typename DerivedI>
bool determine_point_edge_orientation(
const Eigen::MatrixBase<DerivedV>& V,
const Eigen::MatrixBase<DerivedF>& F,
const Eigen::MatrixBase<DerivedI>& I,
const Point_3& query, size_t s, size_t d) {
// Algorithm:
//
// Order the adj faces around the edge (s,d) clockwise using
// query point as pivot. (i.e. The first face of the ordering
// is directly after the pivot point, and the last face is
// directly before the pivot.)
//
// The point is outside if the first and last faces of the
// ordering forms a convex angle. This check can be done
// without any construction by looking at the orientation of the
// faces. The angle is convex iff the first face contains (s,d)
// as an edge and the last face contains (d,s) as an edge.
//
// The point is inside if the first and last faces of the
// ordering forms a concave angle. That is the first face
// contains (d,s) as an edge and the last face contains (s,d) as
// an edge.
//
// In the special case of duplicated faces. I.e. multiple faces
// sharing the same 3 corners, but not necessarily the same
// orientation. The ordering will always rank faces containing
// edge (s,d) before faces containing edge (d,s).
//
// Therefore, if there are any duplicates of the first faces,
// the ordering will always choose the one with edge (s,d) if
// possible. The same for the last face.
//
// In the very degenerated case where the first and last face
// are duplicates, but with different orientations, it is
// equally valid to think the angle formed by them is either 0
// or 360 degrees. By default, 0 degree is used, and thus the
// query point is outside.
std::vector<int> adj_faces;
extract_adj_faces(F, I, s, d, adj_faces);
const size_t num_adj_faces = adj_faces.size();
assert(num_adj_faces > 0);
PlainMatrix<DerivedV,1> pivot_point(1, 3);
igl::copyleft::cgal::assign_scalar(query.x(), pivot_point(0, 0));
igl::copyleft::cgal::assign_scalar(query.y(), pivot_point(0, 1));
igl::copyleft::cgal::assign_scalar(query.z(), pivot_point(0, 2));
using VectorXI = Eigen::Matrix<typename DerivedF::Scalar, Eigen::Dynamic, 1>;
VectorXI order;
order_facets_around_edge(V, F, s, d,
adj_faces, pivot_point, order);
assert((size_t)order.size() == num_adj_faces);
if (adj_faces[order[0]] > 0 &&
adj_faces[order[num_adj_faces-1] < 0]) {
return true;
} else if (adj_faces[order[0]] < 0 &&
adj_faces[order[num_adj_faces-1] > 0]) {
return false;
} else {
throw "The input mesh does not represent a valid volume";
}
throw "The input mesh does not represent a valid volume";
return false;
}
template<typename DerivedV, typename DerivedF, typename DerivedI>
bool determine_point_vertex_orientation(
const Eigen::MatrixBase<DerivedV>& V,
const Eigen::MatrixBase<DerivedF>& F,
const Eigen::MatrixBase<DerivedI>& I,
const Point_3& query, size_t s) {
std::vector<int> adj_vertices;
extract_adj_vertices(F, I, s, adj_vertices);
const size_t num_adj_vertices = adj_vertices.size();
std::vector<Point_3> adj_points;
for (size_t i=0; i<num_adj_vertices; i++) {
const size_t vi = adj_vertices[i];
adj_points.emplace_back(V(vi,0), V(vi,1), V(vi,2));
}
// A plane is on the exterior if all adj_points lies on or to
// one side of the plane.
auto is_on_exterior = [&](const Plane_3& separator) -> bool{
size_t positive=0;
size_t negative=0;
for (const auto& point : adj_points) {
switch(separator.oriented_side(point)) {
case CGAL::ON_POSITIVE_SIDE:
positive++;
break;
case CGAL::ON_NEGATIVE_SIDE:
negative++;
break;
case CGAL::ON_ORIENTED_BOUNDARY:
break;
default:
throw "Unknown plane-point orientation";
}
}
auto query_orientation = separator.oriented_side(query);
bool r =
(positive == 0 && query_orientation == CGAL::POSITIVE)
||
(negative == 0 && query_orientation == CGAL::NEGATIVE);
return r;
};
size_t d = std::numeric_limits<size_t>::max();
Point_3 p(V(s,0), V(s,1), V(s,2));
for (size_t i=0; i<num_adj_vertices; i++) {
const size_t vi = adj_vertices[i];
for (size_t j=i+1; j<num_adj_vertices; j++) {
Plane_3 separator(p, adj_points[i], adj_points[j]);
if (separator.is_degenerate()) {
throw "Input mesh contains degenerated faces";
}
if (is_on_exterior(separator)) {
d = vi;
assert(!CGAL::collinear(p, adj_points[i], query));
break;
}
}
if (d < (size_t)V.rows()) break;
}
if (d > (size_t)V.rows()) {
// All adj faces are coplanar, use the first edge.
d = adj_vertices[0];
}
return determine_point_edge_orientation(V, F, I, query, s, d);
}
template<typename DerivedV, typename DerivedF, typename DerivedI>
bool determine_point_face_orientation(
const Eigen::MatrixBase<DerivedV>& V,
const Eigen::MatrixBase<DerivedF>& F,
const Eigen::MatrixBase<DerivedI>& I,
const Point_3& query, size_t fid) {
// Algorithm: A point is on the inside of a face if the
// tetrahedron formed by them is negatively oriented.
Eigen::Matrix<typename DerivedF::Scalar,3,1> f = F.row(I(fid, 0));
const Point_3 v0(V(f[0], 0), V(f[0], 1), V(f[0], 2));
const Point_3 v1(V(f[1], 0), V(f[1], 1), V(f[1], 2));
const Point_3 v2(V(f[2], 0), V(f[2], 1), V(f[2], 2));
auto result = CGAL::orientation(v0, v1, v2, query);
if (result == CGAL::COPLANAR) {
throw "Cannot determine inside/outside because query point lies exactly on the input surface.";
}
return result == CGAL::NEGATIVE;
}
}
}
}
}
template<typename DerivedV, typename DerivedF, typename DerivedI,
typename DerivedP, typename DerivedB>
IGL_INLINE void igl::copyleft::cgal::points_inside_component(
const Eigen::MatrixBase<DerivedV>& V,
const Eigen::MatrixBase<DerivedF>& F,
const Eigen::MatrixBase<DerivedI>& I,
const Eigen::MatrixBase<DerivedP>& P,
Eigen::PlainObjectBase<DerivedB>& inside) {
using namespace igl::copyleft::cgal::points_inside_component_helper;
if (F.rows() <= 0 || I.rows() <= 0) {
throw "Inside check cannot be done on empty facet component.";
}
const size_t num_faces = I.rows();
std::vector<Triangle> triangles;
for (size_t i=0; i<num_faces; i++) {
const Eigen::Matrix<typename DerivedF::Scalar,3,1> f = F.row(I(i, 0));
triangles.emplace_back(
Point_3(V(f[0], 0), V(f[0], 1), V(f[0], 2)),
Point_3(V(f[1], 0), V(f[1], 1), V(f[1], 2)),
Point_3(V(f[2], 0), V(f[2], 1), V(f[2], 2)));
if (triangles.back().is_degenerate()) {
throw "Input facet components contains degenerated triangles";
}
}
Tree tree(triangles.begin(), triangles.end());
tree.accelerate_distance_queries();
enum ElementType { VERTEX, EDGE, FACE };
auto determine_element_type = [&](
size_t fid, const Point_3& p, size_t& element_index) -> ElementType{
const Eigen::Matrix<typename DerivedF::Scalar,3,1> f = F.row(I(fid, 0));
const Point_3 p0(V(f[0], 0), V(f[0], 1), V(f[0], 2));
const Point_3 p1(V(f[1], 0), V(f[1], 1), V(f[1], 2));
const Point_3 p2(V(f[2], 0), V(f[2], 1), V(f[2], 2));
if (p == p0) { element_index = 0; return VERTEX; }
if (p == p1) { element_index = 1; return VERTEX; }
if (p == p2) { element_index = 2; return VERTEX; }
if (CGAL::collinear(p0, p1, p)) { element_index = 2; return EDGE; }
if (CGAL::collinear(p1, p2, p)) { element_index = 0; return EDGE; }
if (CGAL::collinear(p2, p0, p)) { element_index = 1; return EDGE; }
element_index = 0;
return FACE;
};
const size_t num_queries = P.rows();
inside.resize(num_queries, 1);
//for (size_t i=0; i<num_queries; i++) {
igl::parallel_for(num_queries, [&](const int i) {
const Point_3 query(P(i,0), P(i,1), P(i,2));
auto projection = tree.closest_point_and_primitive(query);
auto closest_point = projection.first;
size_t fid = projection.second - triangles.begin();
size_t element_index;
switch (determine_element_type(fid, closest_point, element_index)) {
case VERTEX:
{
const size_t s = F(I(fid, 0), element_index);
inside(i) = determine_point_vertex_orientation(
V, F, I, query, s);
}
break;
case EDGE:
{
const size_t s = F(I(fid, 0), (element_index+1)%3);
const size_t d = F(I(fid, 0), (element_index+2)%3);
inside(i) = determine_point_edge_orientation(
V, F, I, query, s, d);
}
break;
case FACE:
inside(i) = determine_point_face_orientation(V, F, I, query, fid);
break;
default:
throw "Unknown closest element type!";
}
}, 1000);
}
template<typename DerivedV, typename DerivedF, typename DerivedP,
typename DerivedB>
IGL_INLINE void igl::copyleft::cgal::points_inside_component(
const Eigen::MatrixBase<DerivedV>& V,
const Eigen::MatrixBase<DerivedF>& F,
const Eigen::MatrixBase<DerivedP>& P,
Eigen::PlainObjectBase<DerivedB>& inside) {
using VectorXI = Eigen::Matrix<typename DerivedF::Scalar, Eigen::Dynamic, 1>;
VectorXI I = igl::LinSpaced<VectorXI>(F.rows(), 0, F.rows()-1);
igl::copyleft::cgal::points_inside_component(V, F, I, P, inside);
}
#ifdef IGL_STATIC_LIBRARY
// Explicit template instantiation
// generated by autoexplicit.sh
template void igl::copyleft::cgal::points_inside_component<Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, -1, 0, -1, -1>, Eigen::Matrix<double, -1, 3, 0, -1, 3>, Eigen::Matrix<int, -1, 1, 0, -1, 1>>(Eigen::MatrixBase<Eigen::Matrix<double, -1, -1, 0, -1, -1>> const&, Eigen::MatrixBase<Eigen::Matrix<int, -1, -1, 0, -1, -1>> const&, Eigen::MatrixBase<Eigen::Matrix<double, -1, 3, 0, -1, 3>> const&, Eigen::PlainObjectBase<Eigen::Matrix<int, -1, 1, 0, -1, 1>>&);
// generated by autoexplicit.sh
template void igl::copyleft::cgal::points_inside_component<Eigen::Matrix<CGAL::Epeck::FT, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, -1, 0, -1, -1>, Eigen::Matrix<CGAL::Epeck::FT, -1, 3, 0, -1, 3>, Eigen::Matrix<int, -1, 1, 0, -1, 1>>(Eigen::MatrixBase<Eigen::Matrix<CGAL::Epeck::FT, -1, -1, 0, -1, -1>> const&, Eigen::MatrixBase<Eigen::Matrix<int, -1, -1, 0, -1, -1>> const&, Eigen::MatrixBase<Eigen::Matrix<CGAL::Epeck::FT, -1, 3, 0, -1, 3>> const&, Eigen::PlainObjectBase<Eigen::Matrix<int, -1, 1, 0, -1, 1>>&);
template void igl::copyleft::cgal::points_inside_component< Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix< int, -1, -1, 0, -1, -1>, Eigen::Matrix< int, -1, -1, 0, -1, -1>, Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix< int, -1, -1, 0, -1, -1> > ( Eigen::MatrixBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix< int, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix< int, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> > const&, Eigen::PlainObjectBase<Eigen::Matrix< int, -1, -1, 0, -1, -1> >&);
template void igl::copyleft::cgal::points_inside_component< Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix< int, -1, -1, 0, -1, -1>, Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix< int, -1, -1, 0, -1, -1> > ( Eigen::MatrixBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix< int, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> > const&, Eigen::PlainObjectBase<Eigen::Matrix< int, -1, -1, 0, -1, -1> >&);
template void igl::copyleft::cgal::points_inside_component<Eigen::Matrix<CGAL::Epeck::FT, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, -1, 0, -1, -1>, Eigen::Matrix<CGAL::Epeck::FT, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, 1, 0, -1, 1> >(Eigen::MatrixBase<Eigen::Matrix<CGAL::Epeck::FT, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix<int, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix<CGAL::Epeck::FT, -1, -1, 0, -1, -1> > const&, Eigen::PlainObjectBase<Eigen::Matrix<int, -1, 1, 0, -1, 1> >&);
template void igl::copyleft::cgal::points_inside_component<Eigen::Matrix<CGAL::Epeck::FT, -1, 3, 0, -1, 3>, Eigen::Matrix<int, -1, 3, 0, -1, 3>, Eigen::Matrix<CGAL::Epeck::FT, -1, 3, 0, -1, 3>, Eigen::Matrix<int, -1, 1, 0, -1, 1> >(Eigen::MatrixBase<Eigen::Matrix<CGAL::Epeck::FT, -1, 3, 0, -1, 3> > const&, Eigen::MatrixBase<Eigen::Matrix<int, -1, 3, 0, -1, 3> > const&, Eigen::MatrixBase<Eigen::Matrix<CGAL::Epeck::FT, -1, 3, 0, -1, 3> > const&, Eigen::PlainObjectBase<Eigen::Matrix<int, -1, 1, 0, -1, 1> >&);
template void igl::copyleft::cgal::points_inside_component<Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, -1, 0, -1, -1>, Eigen::Matrix<CGAL::Epeck::FT, -1, 3, 0, -1, 3>, Eigen::Array<bool, -1, 1, 0, -1, 1> >(Eigen::MatrixBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix<int, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix<CGAL::Epeck::FT, -1, 3, 0, -1, 3> > const&, Eigen::PlainObjectBase<Eigen::Array<bool, -1, 1, 0, -1, 1> >&);
template void igl::copyleft::cgal::points_inside_component<Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, -1, 0, -1, -1>, Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, 1, 0, -1, 1> >(Eigen::MatrixBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix<int, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> > const&, Eigen::PlainObjectBase<Eigen::Matrix<int, -1, 1, 0, -1, 1> >&);
template void igl::copyleft::cgal::points_inside_component<Eigen::Matrix<double, -1, 3, 0, -1, 3>, Eigen::Matrix<int, -1, 3, 0, -1, 3>, Eigen::Matrix<double, -1, 3, 0, -1, 3>, Eigen::Matrix<int, -1, 1, 0, -1, 1> >(Eigen::MatrixBase<Eigen::Matrix<double, -1, 3, 0, -1, 3> > const&, Eigen::MatrixBase<Eigen::Matrix<int, -1, 3, 0, -1, 3> > const&, Eigen::MatrixBase<Eigen::Matrix<double, -1, 3, 0, -1, 3> > const&, Eigen::PlainObjectBase<Eigen::Matrix<int, -1, 1, 0, -1, 1> >&);
#endif