From 2a4f60b822b15c34058220beffd311f46f11d3b3 Mon Sep 17 00:00:00 2001 From: Hind-M Date: Tue, 2 Aug 2022 12:04:08 +0200 Subject: Get filtration from face when it's the same as simplex --- src/Cech_complex/include/gudhi/Cech_complex_blocker.h | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) (limited to 'src') diff --git a/src/Cech_complex/include/gudhi/Cech_complex_blocker.h b/src/Cech_complex/include/gudhi/Cech_complex_blocker.h index 22b99c5c..6f082494 100644 --- a/src/Cech_complex/include/gudhi/Cech_complex_blocker.h +++ b/src/Cech_complex/include/gudhi/Cech_complex_blocker.h @@ -106,7 +106,10 @@ class Cech_blocker { #if CGAL_VERSION_NR >= 1050000000 if(exact_) CGAL::exact(sph.second); #endif - radius = std::sqrt(cast_to_fv(sph.second)); + if(k != sc_ptr_->null_key()) + radius = sc_ptr_->filtration(face_opposite_vertex.first); + else + radius = std::sqrt(cast_to_fv(sph.second)); #ifdef DEBUG_TRACES std::clog << "center: " << sph.first << ", radius: " << radius << std::endl; #endif // DEBUG_TRACES -- cgit v1.2.3