23 #ifndef WITNESS_COMPLEX_ALL_FACES_IN_H_ 24 #define WITNESS_COMPLEX_ALL_FACES_IN_H_ 30 template <
typename SimplicialComplexForWitness,
32 bool all_faces_in(Simplex& simplex,
33 double* filtration_value,
34 SimplicialComplexForWitness& sc) {
35 typedef typename SimplicialComplexForWitness::Simplex_handle Simplex_handle;
37 if (simplex.size() == 1)
41 for (
typename Simplex::iterator not_it = simplex.begin(); not_it != simplex.end(); ++not_it) {
43 for (
typename Simplex::iterator it = simplex.begin(); it != simplex.end(); ++it)
46 Simplex_handle facet_sh = sc.find(facet);
47 if (facet_sh == sc.null_simplex())
49 else if (sc.filtration(facet_sh) > *filtration_value)
50 *filtration_value = sc.filtration(facet_sh);
55 #endif // WITNESS_COMPLEX_ALL_FACES_IN_H_