songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

mismatched volumes error detected

This commit is contained in:
Dmitry Anisimov 2020-11-24 11:48:18 +01:00
parent df3703729d
commit 1c81d34b5d
1 changed files with 23 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

32
Kinetic_shape_reconstruction/include/CGAL/KSR_3/Data_structure.h
View File

@ -2359,10 +2359,10 @@ public:
// std::to_string(volume_index) + "-" +
// std::to_string(pface.first) + "-" +
// std::to_string(pface.second) << std::endl;
dump_pface(pface, "bnd-pface-" +
std::to_string(volume_index) + "-" +
std::to_string(pface.first) + "-" +
std::to_string(pface.second));
// dump_pface(pface, "bnd-pface-" +
// std::to_string(volume_index) + "-" +
// std::to_string(pface.first) + "-" +
// std::to_string(pface.second));
std::vector<PFace> nfaces;
for (const auto pedge : pedges_of_pface(pface)) {
@ -2409,17 +2409,29 @@ public:
std::set<int> neighbors;
std::vector<PFace> nfaces;
for (const auto pedge : pedges_of_pface(pface)) {
const auto pedges = pedges_of_pface(pface);
for (const auto pedge : pedges) {
CGAL_assertion(has_iedge(pedge));
incident_faces(this->iedge(pedge), nfaces);
std::cout << "pedge: " << segment_3(pedge) << std::endl;
std::cout << "nfaces size: " << nfaces.size() << std::endl;
for (const auto& nface : nfaces) {
if (nface == pface) continue;
if (nface.first >= 6) continue; // TODO: skip interior pfaces (see below)
if (nface == pface) {
std::cout << "skip pface" << "; ";
continue;
}
// TODO: skip interior pfaces (see below)
if (nface.first >= 6) {
std::cout << "skip interior" << "; ";
continue;
}
const auto& npair = map_volumes.at(nface);
CGAL_assertion(npair.first != -1);
neighbors.insert(npair.first);
std::cout << npair.first << "; ";
// std::cout << "dumping nface: " <<
// std::to_string(nface.first) + "-" + std::to_string(nface.second) << std::endl;
@ -2428,7 +2440,9 @@ public:
if (npair.second == -1) continue;
neighbors.insert(npair.second);
std::cout << npair.second << "; ";
}
std::cout << std::endl << std::endl;
}
if (neighbors.size() == 2) { // interior pfaces, which have boundary pfaces as neighbors
@ -2444,11 +2458,11 @@ public:
} else { // pure interior pfaces
std::cout << "found volumes: ";
for (const int neighbor : neighbors)
std::cout << neighbor << " ";
std::cout << std::endl << "NOT FOUND size: " << neighbors.size() << std::endl;
std::cout << std::endl << "neighbors size: " << neighbors.size() << std::endl;
CGAL_assertion_msg(false, "TODO: HANDLE PURE INTERIOR FACE!");
}
return false;
}