1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
|
/**************************************************************************
* *
* Regina - A Normal Surface Theory Calculator *
* Python Interface *
* *
* Copyright (c) 1999-2025, Ben Burton *
* For further details contact Ben Burton (bab@debian.org). *
* *
* This program is free software; you can redistribute it and/or *
* modify it under the terms of the GNU General Public License as *
* published by the Free Software Foundation; either version 2 of the *
* License, or (at your option) any later version. *
* *
* As an exception, when this program is distributed through (i) the *
* App Store by Apple Inc.; (ii) the Mac App Store by Apple Inc.; or *
* (iii) Google Play by Google Inc., then that store may impose any *
* digital rights management, device limits and/or redistribution *
* restrictions that are required by its terms of service. *
* *
* This program is distributed in the hope that it will be useful, but *
* WITHOUT ANY WARRANTY; without even the implied warranty of *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU *
* General Public License for more details. *
* *
* You should have received a copy of the GNU General Public License *
* along with this program. If not, see <https://www.gnu.org/licenses/>. *
* *
**************************************************************************/
#include <pybind11/pybind11.h>
#include <pybind11/stl.h>
#include "triangulation/generic.h"
#include "../helpers.h"
#include "../generic/facehelper.h"
#include "../docstrings/triangulation/generic/boundarycomponent.h"
#include "../docstrings/triangulation/detail/boundarycomponent.h"
using regina::BoundaryComponent;
using regina::python::invalidFaceDimension;
template <int dim>
void addBoundaryComponent(pybind11::module_& m, pybind11::module_& internal,
const char* name) {
// In higher dimensions:
// - we do not store lower-dimensional faces;
// - we do not recognise ideal or invalid vertices;
// - we can still triangulate a real boundary component.
// We use the global scope here because all of BoundaryComponent's members
// are inherited, and so BoundaryComponent's own docstring namespace
// does not exist.
RDOC_SCOPE_BEGIN_MAIN
RDOC_SCOPE_BASE(detail::BoundaryComponentBase)
auto c = pybind11::class_<BoundaryComponent<dim>>(m, name,
rdoc::BoundaryComponent)
.def("index", &BoundaryComponent<dim>::index, rbase::index)
.def("size", &BoundaryComponent<dim>::size, rbase::size)
.def("countRidges", &BoundaryComponent<dim>::countRidges,
rbase::countRidges)
.def("countFaces", [](const BoundaryComponent<dim>& b, int subdim) {
if (subdim == dim - 1)
return b.template countFaces<dim - 1>();
else if (subdim == dim - 2)
return b.template countFaces<dim - 2>();
else {
invalidFaceDimension("countFaces", dim - 2, dim - 1);
// This throws, but the compiler wants us to return a value.
return (size_t)0;
}
}, pybind11::arg("subdim"), rbase::countFaces)
.def("facets", &BoundaryComponent<dim>::facets, rbase::facets)
.def("faces", [](const BoundaryComponent<dim>& b, int subdim) {
if (subdim != dim - 1)
invalidFaceDimension("faces", dim - 1, dim - 1);
return b.template faces<dim - 1>();
}, pybind11::arg("subdim"), rbase::faces)
.def("facet", &BoundaryComponent<dim>::facet,
pybind11::return_value_policy::reference, rbase::facet)
.def("face", [](const BoundaryComponent<dim>& b, int subdim,
size_t index) {
if (subdim != dim - 1)
invalidFaceDimension("face", dim - 1, dim - 1);
return b.template face<dim - 1>(index);
}, pybind11::return_value_policy::reference,
pybind11::arg("subdim"), pybind11::arg("index"), rbase::face)
.def("component", &BoundaryComponent<dim>::component,
pybind11::return_value_policy::reference, rbase::component)
.def("triangulation", &BoundaryComponent<dim>::triangulation,
rbase::triangulation)
.def("build", [](const BoundaryComponent<dim>& b) {
// Return a clone of the resulting triangulation.
// This is because Python cannot enforce the constness of
// the reference that would normally be returned.
return new regina::Triangulation<dim-1>(b.build());
}, rbase::build)
.def("isReal", &BoundaryComponent<dim>::isReal, rbase::isReal)
.def("isIdeal", &BoundaryComponent<dim>::isIdeal, rbase::isIdeal)
.def("isInvalidVertex", &BoundaryComponent<dim>::isInvalidVertex,
rbase::isInvalidVertex)
.def("isOrientable", &BoundaryComponent<dim>::isOrientable,
rbase::isOrientable)
.def_readonly_static("dimension", &BoundaryComponent<dim>::dimension)
.def_readonly_static("allFaces", &BoundaryComponent<dim>::allFaces)
.def_readonly_static("allowVertex",
&BoundaryComponent<dim>::allowVertex)
.def_readonly_static("canBuild", &BoundaryComponent<dim>::canBuild)
;
if constexpr (dim == 5) {
c.def("countPentachora", &BoundaryComponent<dim>::countPentachora,
rbase::countPentachora);
c.def("pentachoron", &BoundaryComponent<dim>::pentachoron,
pybind11::return_value_policy::reference, rbase::pentachoron);
c.def("pentachora", &BoundaryComponent<dim>::pentachora,
rbase::pentachora);
}
if constexpr (dim == 6) {
c.def("countPentachora", &BoundaryComponent<dim>::countPentachora,
rbase::countPentachora);
}
regina::python::add_output(c);
regina::python::add_eq_operators(c);
RDOC_SCOPE_END
regina::python::addListView<
decltype(std::declval<BoundaryComponent<dim>>().facets())>(internal,
(std::string(name) + "_facets").c_str());
}
|
