54 using Simplex_cell_map = std::map<Simplex_handle, Hasse_cell*, Simplex_comparator<Simplex_handle> >;
72 Simplex_cell_maps& simplex_cell_maps = (is_boundary ? boundary_simplex_cell_maps_ : interior_simplex_cell_maps_);
74 CC_detail_list& cc_detail_list =
75 (is_boundary ? cc_boundary_detail_lists[cell_d] : cc_interior_detail_lists[cell_d]);
76 cc_detail_list.emplace_back(simplex);
79 auto map_it = simplex_cell_map.find(simplex);
80 if (map_it == simplex_cell_map.end()) {
81 hasse_cells_.push_back(
new Hasse_cell(is_boundary, cell_d));
83 simplex_cell_map.emplace(simplex, new_cell);
84 cell_simplex_map_.emplace(new_cell, simplex);
86 cc_detail_list.back().status_ = CC_detail_info::Result_type::inserted;
91 CC_detail_info& cc_info = cc_detail_list.back();
92 cc_info.trigger_ = to_string(map_it->first);
93 cc_info.status_ = CC_detail_info::Result_type::self;
95 return map_it->second;
98 void expand_level(std::size_t cell_d) {
99 bool is_manifold_with_boundary = boundary_simplex_cell_maps_.size() > 0;
100 for (
auto& sc_pair : interior_simplex_cell_maps_[cell_d - 1]) {
101 const Simplex_handle& simplex = sc_pair.first;
103 for (Simplex_handle coface : simplex.coface_range(cod_d_ + cell_d)) {
104 Hasse_cell* new_cell = insert_cell(coface, cell_d,
false);
105 new_cell->get_boundary().emplace_back(cell, 1);
109 if (is_manifold_with_boundary) {
110 for (
auto& sc_pair : boundary_simplex_cell_maps_[cell_d - 1]) {
111 const Simplex_handle& simplex = sc_pair.first;
113 if (cell_d != intr_d_)
114 for (Simplex_handle coface : simplex.coface_range(cod_d_ + cell_d + 1)) {
115 Hasse_cell* new_cell = insert_cell(coface, cell_d,
true);
116 new_cell->get_boundary().emplace_back(cell, 1);
118 auto map_it = interior_simplex_cell_maps_[cell_d].find(simplex);
119 if (map_it == interior_simplex_cell_maps_[cell_d].end())
120 std::cerr <<
"Cell_complex::expand_level error: A boundary cell does not have an interior counterpart.\n";
123 i_cell->get_boundary().emplace_back(cell, 1);
129 void construct_complex_(
const Out_simplex_map_& out_simplex_map) {
131 cc_interior_summary_lists.resize(interior_simplex_cell_maps_.size());
132 cc_interior_prejoin_lists.resize(interior_simplex_cell_maps_.size());
133 cc_interior_detail_lists.resize(interior_simplex_cell_maps_.size());
135 for (
auto& os_pair : out_simplex_map) {
136 const Simplex_handle& simplex = os_pair.first;
137 const Eigen::VectorXd& point = os_pair.second;
138 Hasse_cell* new_cell = insert_cell(simplex, 0,
false);
139 cell_point_map_.emplace(new_cell, point);
141 for (std::size_t cell_d = 1;
142 cell_d < interior_simplex_cell_maps_.size() && !interior_simplex_cell_maps_[cell_d - 1].empty(); ++cell_d) {
143 expand_level(cell_d);
147 void construct_complex_(
const Out_simplex_map_& interior_simplex_map,
const Out_simplex_map_& boundary_simplex_map) {
149 cc_interior_summary_lists.resize(interior_simplex_cell_maps_.size());
150 cc_interior_prejoin_lists.resize(interior_simplex_cell_maps_.size());
151 cc_interior_detail_lists.resize(interior_simplex_cell_maps_.size());
152 cc_boundary_summary_lists.resize(boundary_simplex_cell_maps_.size());
153 cc_boundary_prejoin_lists.resize(boundary_simplex_cell_maps_.size());
154 cc_boundary_detail_lists.resize(boundary_simplex_cell_maps_.size());
156 for (
auto& os_pair : boundary_simplex_map) {
157 const Simplex_handle& simplex = os_pair.first;
158 const Eigen::VectorXd& point = os_pair.second;
159 Hasse_cell* new_cell = insert_cell(simplex, 0,
true);
160 cell_point_map_.emplace(new_cell, point);
162 for (
auto& os_pair : interior_simplex_map) {
163 const Simplex_handle& simplex = os_pair.first;
164 const Eigen::VectorXd& point = os_pair.second;
165 Hasse_cell* new_cell = insert_cell(simplex, 0,
false);
166 cell_point_map_.emplace(new_cell, point);
169 for (
const auto& sc_pair : interior_simplex_cell_maps_[0])
170 cc_interior_summary_lists[0].push_back(CC_summary_info(sc_pair));
171 for (
const auto& sc_pair : boundary_simplex_cell_maps_[0])
172 cc_boundary_summary_lists[0].push_back(CC_summary_info(sc_pair));
175 for (std::size_t cell_d = 1;
176 cell_d < interior_simplex_cell_maps_.size() && !interior_simplex_cell_maps_[cell_d - 1].empty(); ++cell_d) {
177 expand_level(cell_d);
180 for (
const auto& sc_pair : interior_simplex_cell_maps_[cell_d])
181 cc_interior_summary_lists[cell_d].push_back(CC_summary_info(sc_pair));
182 if (cell_d < boundary_simplex_cell_maps_.size())
183 for (
const auto& sc_pair : boundary_simplex_cell_maps_[cell_d])
184 cc_boundary_summary_lists[cell_d].push_back(CC_summary_info(sc_pair));
200 interior_simplex_cell_maps_.resize(intr_d_ + 1);
201 if (!out_simplex_map.empty()) cod_d_ = out_simplex_map.begin()->first.dimension();
202 construct_complex_(out_simplex_map);
217 interior_simplex_cell_maps_.resize(limit_dimension + 1);
218 if (!out_simplex_map.empty()) cod_d_ = out_simplex_map.begin()->first.dimension();
219 construct_complex_(out_simplex_map);
234 void construct_complex(
const Out_simplex_map_& interior_simplex_map,
const Out_simplex_map_& boundary_simplex_map) {
235 interior_simplex_cell_maps_.resize(intr_d_ + 1);
236 boundary_simplex_cell_maps_.resize(intr_d_);
237 if (!interior_simplex_map.empty()) cod_d_ = interior_simplex_map.begin()->first.dimension();
238 construct_complex_(interior_simplex_map, boundary_simplex_map);
255 void construct_complex(
const Out_simplex_map_& interior_simplex_map,
const Out_simplex_map_& boundary_simplex_map,
256 std::size_t limit_dimension) {
257 interior_simplex_cell_maps_.resize(limit_dimension + 1);
258 boundary_simplex_cell_maps_.resize(limit_dimension);
259 if (!interior_simplex_map.empty()) cod_d_ = interior_simplex_map.begin()->first.dimension();
260 construct_complex_(interior_simplex_map, boundary_simplex_map);
290 return interior_simplex_cell_maps_[cell_d];
301 return boundary_simplex_cell_maps_[cell_d];
325 for (
Hasse_cell* hs_ptr : hasse_cells_) delete hs_ptr;
329 std::size_t intr_d_, cod_d_;
333 std::vector<Hasse_cell*> hasse_cells_;