Add this check in debug mode

1 files changed, 1 insertions, 0 deletions

diff --git a/src/Coxeter_triangulation/include/gudhi/Freudenthal_triangulation.h b/src/Coxeter_triangulation/include/gudhi/Freudenthal_triangulation.h
index dd059184..845e4883 100644
--- a/src/Coxeter_triangulation/include/gudhi/Freudenthal_triangulation.h
+++ b/src/Coxeter_triangulation/include/gudhi/Freudenthal_triangulation.h
@@ -149,6 +149,7 @@ class Freudenthal_triangulation {
     } else {
       Eigen::VectorXd p_vect(d);
       for (std::size_t i = 0; i < d; i++) p_vect(i) = point[i];
+      assert(p_vect.size() == offset_.size());
       Eigen::VectorXd x_vect = colpivhouseholderqr_.solve(p_vect - offset_);
       for (std::size_t i = 0; i < d; i++) {
         double x_i = scale * x_vect(i);