diff options
author | ROUVREAU Vincent <vincent.rouvreau@inria.fr> | 2020-11-10 15:48:33 +0100 |
---|---|---|
committer | ROUVREAU Vincent <vincent.rouvreau@inria.fr> | 2020-11-10 15:48:33 +0100 |
commit | 1d7c18e91a73a539e807d323a47d32bbba18dac7 (patch) | |
tree | 5455a6f508a2823a9c5b2be4e815aacdbc9a6a62 /src/Coxeter_triangulation/include/gudhi | |
parent | b118b86b7f2e7ae9ee7f080ec89d956b161aedfd (diff) |
Add this check in debug mode
Diffstat (limited to 'src/Coxeter_triangulation/include/gudhi')
-rw-r--r-- | src/Coxeter_triangulation/include/gudhi/Freudenthal_triangulation.h | 1 |
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); |