OpenMesh: OpenMesh/Tools/Utils/MeshCheckerT_impl.hh Source File

| OpenMesh

MeshCheckerT_impl.hh
1 /* ========================================================================= *
2 * *
3 * OpenMesh *
4 * Copyright (c) 2001-2022, RWTH-Aachen University *
5 * Department of Computer Graphics and Multimedia *
6 * All rights reserved. *
7 * www.openmesh.org *
8 * *
9 *---------------------------------------------------------------------------*
10 * This file is part of OpenMesh. *
11 *---------------------------------------------------------------------------*
12 * *
13 * Redistribution and use in source and binary forms, with or without *
14 * modification, are permitted provided that the following conditions *
15 * are met: *
16 * *
17 * 1. Redistributions of source code must retain the above copyright notice, *
18 * this list of conditions and the following disclaimer. *
19 * *
20 * 2. Redistributions in binary form must reproduce the above copyright *
21 * notice, this list of conditions and the following disclaimer in the *
22 * documentation and/or other materials provided with the distribution. *
23 * *
24 * 3. Neither the name of the copyright holder nor the names of its *
25 * contributors may be used to endorse or promote products derived from *
26 * this software without specific prior written permission. *
27 * *
28 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS *
29 * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED *
30 * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A *
31 * PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER *
32 * OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, *
33 * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, *
34 * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR *
35 * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF *
36 * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING *
37 * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS *
38 * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. *
39 * *
40 * ========================================================================= */
41
42
43
44
45 #define OPENMESH_MESHCHECKER_C
46
47
48 //== INCLUDES =================================================================
49
50
51 #include <OpenMesh/Tools/Utils/MeshCheckerT.hh>
52
53
54 //== NAMESPACES ==============================================================
55
56
57 namespace OpenMesh {
58 namespace Utils {
59
60 //== IMPLEMENTATION ==========================================================
61
62
63 template <class Mesh>
64 bool
65 MeshCheckerT<Mesh>::
66 check(unsigned int _targets, std::ostream& _os)
67 {
68 bool ok(true);
69
70
71
72 //--- vertex checks ---
73
74 if (_targets & CHECK_VERTICES)
75 {
76 typename Mesh::ConstVertexIter v_it(mesh_.vertices_begin()), v_end(mesh_.vertices_end());
77 typename Mesh::VertexHandle vh;
78 typename Mesh::ConstVertexVertexCWIter vv_it;
79 typename Mesh::HalfedgeHandle heh;
80 unsigned int count;
81 const unsigned int max_valence(10000);
82
83
84 for (; v_it != v_end; ++v_it)
85 {
86 if (!is_deleted(*v_it))
87 {
88 vh = *v_it;
89
90
91 /* The outgoing halfedge of a boundary vertex has to be a boundary halfedge */
92 heh = mesh_.halfedge_handle(vh);
93 if (heh.is_valid() && !mesh_.is_boundary(heh))
94 {
95 for (typename Mesh::ConstVertexOHalfedgeIter vh_it(mesh_, vh);
96 vh_it.is_valid(); ++vh_it)
97 {
98 if (mesh_.is_boundary(*vh_it))
99 {
100 _os << "MeshChecker: vertex " << vh
101 << ": outgoing halfedge not on boundary error\n";
102 ok = false;
103 }
104 }
105 }
106
107
108
109 // outgoing halfedge has to refer back to vertex
110 if (mesh_.halfedge_handle(vh).is_valid() &&
111 mesh_.from_vertex_handle(mesh_.halfedge_handle(vh)) != vh)
112 {
113 _os << "MeshChecker: vertex " << vh
114 << ": outgoing halfedge does not reference vertex\n";
115 ok = false;
116 }
117
118
119 // check whether circulators are still in order
120 vv_it = mesh_.cvv_cwiter(vh);
121 for (count=0; vv_it.is_valid() && (count < max_valence); ++vv_it, ++count) {};
122 if (count == max_valence)
123 {
124 _os << "MeshChecker: vertex " << vh
125 << ": ++circulator problem, one ring corrupt\n";
126 ok = false;
127 }
128 vv_it = mesh_.cvv_cwiter(vh);
129 for (count=0; vv_it.is_valid() && (count < max_valence); --vv_it, ++count) {};
130 if (count == max_valence)
131 {
132 _os << "MeshChecker: vertex " << vh
133 << ": --circulator problem, one ring corrupt\n";
134 ok = false;
135 }
136 }
137 }
138 }
139
140
141
142 //--- halfedge checks ---
143
144 if (_targets & CHECK_EDGES)
145 {
146 typename Mesh::ConstHalfedgeIter h_it(mesh_.halfedges_begin()),
147 h_end(mesh_.halfedges_end());
148 typename Mesh::HalfedgeHandle hh, hstart, hhh;
149 size_t count, n_halfedges = 2*mesh_.n_edges();
150
151 for (; h_it != h_end; ++h_it)
152 {
153 if (!is_deleted(mesh_.edge_handle(*h_it)))
154 {
155 hh = *h_it;
156
157
158 // degenerated halfedge ?
159 if (mesh_.from_vertex_handle(hh) == mesh_.to_vertex_handle(hh))
160 {
161 _os << "MeshChecker: halfedge " << hh
162 << ": to-vertex == from-vertex\n";
163 ok = false;
164 }
165
166
167 // next <-> prev check
168 if (mesh_.next_halfedge_handle(mesh_.prev_halfedge_handle(hh)) != hh)
169 {
170 _os << "MeshChecker: halfedge " << hh
171 << ": prev->next != this\n";
172 ok = false;
173 }
174
175
176 // halfedges should form a cycle
177 count=0; hstart=hhh=hh;
178 do
179 {
180 hhh = mesh_.next_halfedge_handle(hhh);
181 ++count;
182 } while (hhh != hstart && count < n_halfedges);
183
184 if (count == n_halfedges)
185 {
186 _os << "MeshChecker: halfedges starting from " << hh
187 << " do not form a cycle\n";
188 ok = false;
189 }
190 }
191 }
192 }
193
194
195
196 //--- face checks ---
197
198 if (_targets & CHECK_FACES)
199 {
200 typename Mesh::ConstFaceIter f_it(mesh_.faces_begin()),
201 f_end(mesh_.faces_end());
202 typename Mesh::FaceHandle fh;
203 typename Mesh::ConstFaceHalfedgeIter fh_it;
204
205 for (; f_it != f_end; ++f_it)
206 {
207 if (!is_deleted(*f_it))
208 {
209 fh = *f_it;
210
211 for (fh_it=mesh_.cfh_iter(fh); fh_it.is_valid(); ++fh_it)
212 {
213 if (mesh_.face_handle(*fh_it) != fh)
214 {
215 _os << "MeshChecker: face " << fh
216 << ": its halfedge does not reference face\n";
217 ok = false;
218 }
219 }
220 }
221 }
222 }
223
224
225
226 return ok;
227 }
228
229
230 //=============================================================================
231 } // naespace Utils
232 } // namespace OpenMesh
233 //=============================================================================
OpenMesh
Contains all the mesh ingredients like the polygonal mesh, the triangle mesh, different mesh kernels ...
Definition: MeshItems.hh:59
OpenMesh::PolyMeshT::ConstFaceIter
Kernel::ConstFaceIter ConstFaceIter
Scalar type.
Definition: PolyMeshT.hh:151
OpenMesh::PolyMeshT::VertexHandle
Kernel::VertexHandle VertexHandle
Handle for referencing the corresponding item.
Definition: PolyMeshT.hh:136
OpenMesh::PolyMeshT::ConstFaceHalfedgeIter
Kernel::ConstFaceHalfedgeIter ConstFaceHalfedgeIter
Circulator.
Definition: PolyMeshT.hh:178
OpenMesh::PolyMeshT::ConstHalfedgeIter
Kernel::ConstHalfedgeIter ConstHalfedgeIter
Scalar type.
Definition: PolyMeshT.hh:149
OpenMesh::PolyMeshT::FaceHandle
Kernel::FaceHandle FaceHandle
Scalar type.
Definition: PolyMeshT.hh:139
OpenMesh::PolyMeshT::HalfedgeHandle
Kernel::HalfedgeHandle HalfedgeHandle
Scalar type.
Definition: PolyMeshT.hh:137
OpenMesh::PolyMeshT::ConstVertexIter
Kernel::ConstVertexIter ConstVertexIter
Scalar type.
Definition: PolyMeshT.hh:148
OpenMesh::PolyMeshT::ConstVertexOHalfedgeIter
Kernel::ConstVertexOHalfedgeIter ConstVertexOHalfedgeIter
Circulator.
Definition: PolyMeshT.hh:173
OpenMesh::Utils::MeshCheckerT::check
bool check(unsigned int _targets=CHECK_ALL, std::ostream &_os=omerr())
check it, return true iff ok
Definition: MeshCheckerT_impl.hh:66