51 #define OPENMESH_MULTIPLE_CHOICE_DECIMATER_DECIMATERT_CC
58 #if defined(OM_CC_MIPS)
65 # include <OpenMesh/Core/Utils/RandomNumberGenerator.hh>
78 mesh_(_mesh), randomSamples_(10) {
81 mesh_.request_vertex_status();
82 mesh_.request_halfedge_status();
83 mesh_.request_edge_status();
84 mesh_.request_face_status();
93 mesh_.release_vertex_status();
94 mesh_.release_edge_status();
95 mesh_.release_halfedge_status();
96 mesh_.release_face_status();
104 if (!this->is_initialized())
107 unsigned int n_collapses(0);
109 bool collapsesUnchanged =
false;
111 unsigned int oldCollapses = 0;
114 unsigned int noCollapses = 0;
120 const bool update_normals = mesh_.has_face_normals();
122 while ( n_collapses < _n_collapses) {
124 if (noCollapses > 20) {
125 omlog() <<
"[McDecimater] : no collapses performed in over 20 iterations in a row\n";
132 double bestEnergy = FLT_MAX;
133 double energy = FLT_MAX;
136 for (
int i = 0; i < (int)randomSamples_; ++i) {
140 tmpHandle =
typename Mesh::HalfedgeHandle(
int(randGen.getRand() * double(mesh_.n_halfedges() - 1.0)) );
142 tmpHandle =
typename Mesh::HalfedgeHandle( (
double(rand()) /
double(RAND_MAX) ) *
double(mesh_.n_halfedges()-1) );
146 if ( ! mesh_.status(tmpHandle).deleted() ) {
151 if (this->is_collapse_legal(ci)) {
152 energy = this->collapse_priority(ci);
156 if ( energy < bestEnergy ) {
158 bestHandle = tmpHandle;
169 if ( bestEnergy != FLT_MAX ) {
175 if (!this->is_collapse_legal(ci))
179 this->preprocess_collapse(ci);
182 mesh_.collapse(bestHandle);
186 oldCollapses = n_collapses;
188 collapsesUnchanged =
false;
194 for (; vf_it.is_valid(); ++vf_it)
195 if (!mesh_.status(*vf_it).deleted())
196 mesh_.set_normal(*vf_it, mesh_.calc_face_normal(*vf_it));
200 this->postprocess_collapse(ci);
203 if (!this->notify_observer(n_collapses))
207 if (oldCollapses == n_collapses) {
208 if (collapsesUnchanged ==
false) {
210 collapsesUnchanged =
true;
228 if (!this->is_initialized())
232 if ( (_nv == 0) && (_nf == 1) )
233 return decimate_constraints_only(1.0);
235 size_t nv = mesh_.n_vertices();
236 size_t nf = mesh_.n_faces();
237 unsigned int n_collapses(0);
240 bool collapsesUnchanged =
false;
242 unsigned int oldCollapses = 0;
245 unsigned int noCollapses = 0;
251 const bool update_normals = mesh_.has_face_normals();
253 while ((_nv < nv) && (_nf < nf)) {
255 if (noCollapses > 20) {
256 omlog() <<
"[McDecimater] : no collapses performed in over 20 iterations in a row\n";
263 double bestEnergy = FLT_MAX;
264 double energy = FLT_MAX;
267 for (
int i = 0; i < (int) randomSamples_; ++i) {
271 tmpHandle =
typename Mesh::HalfedgeHandle(
int(randGen.getRand() * double(mesh_.n_halfedges() - 1.0)) );
273 tmpHandle =
typename Mesh::HalfedgeHandle( (
double(rand()) /
double(RAND_MAX) ) *
double(mesh_.n_halfedges() - 1));
277 if (!mesh_.status(tmpHandle).deleted()) {
282 if (this->is_collapse_legal(ci)) {
283 energy = this->collapse_priority(ci);
287 if (energy < bestEnergy) {
289 bestHandle = tmpHandle;
300 if ( bestEnergy != FLT_MAX ) {
306 if (!this->is_collapse_legal(ci))
316 if (mesh_.is_boundary(ci.
v0v1) || mesh_.is_boundary(ci.
v1v0))
322 this->preprocess_collapse(ci);
325 mesh_.collapse(bestHandle);
329 oldCollapses = n_collapses;
331 collapsesUnchanged =
false;
337 for (; vf_it.is_valid(); ++vf_it)
338 if (!mesh_.status(*vf_it).deleted())
339 mesh_.set_normal(*vf_it, mesh_.calc_face_normal(*vf_it));
343 this->postprocess_collapse(ci);
346 if (!this->notify_observer(n_collapses))
350 if (oldCollapses == n_collapses) {
351 if (collapsesUnchanged ==
false) {
353 collapsesUnchanged =
true;
371 if (!this->is_initialized())
375 this->set_error_tolerance_factor(_factor);
377 unsigned int n_collapses(0);
378 size_t nv = mesh_.n_vertices();
379 size_t nf = mesh_.n_faces();
381 bool lastCollapseIllegal =
false;
383 unsigned int illegalCollapses = 0;
385 bool collapsesUnchanged =
false;
387 unsigned int oldCollapses = 0;
390 unsigned int noCollapses = 0;
392 double energy = FLT_MAX;
393 double bestEnergy = FLT_MAX;
399 while ((noCollapses <= 50) && (illegalCollapses <= 50) && (nv > 0) && (nf > 1)) {
404 bestEnergy = FLT_MAX;
408 const double randomNormalizer = (1.0 / RAND_MAX) * (mesh_.n_halfedges() - 1);
412 for (
int i = 0; i < (int) randomSamples_; ++i) {
416 tmpHandle =
typename Mesh::HalfedgeHandle(
int(randGen.getRand() * double(mesh_.n_halfedges() - 1.0)) );
422 if (!mesh_.status(mesh_.edge_handle(tmpHandle)).deleted()) {
427 if (this->is_collapse_legal(ci)) {
429 energy = this->collapse_priority(ci);
432 if (lastCollapseIllegal) {
435 illegalCollapses = 1;
436 lastCollapseIllegal =
true;
440 illegalCollapses = 0;
441 lastCollapseIllegal =
false;
444 if (energy < bestEnergy) {
446 bestHandle = tmpHandle;
460 if ( bestEnergy != FLT_MAX ) {
466 if (!this->is_collapse_legal(ci))
476 if (mesh_.is_boundary(ci.
v0v1) || mesh_.is_boundary(ci.
v1v0))
482 this->preprocess_collapse(ci);
485 mesh_.collapse(bestHandle);
489 oldCollapses = n_collapses;
491 collapsesUnchanged =
false;
496 for (; vf_it.is_valid(); ++vf_it)
497 if (!mesh_.status(*vf_it).deleted())
498 mesh_.set_normal(*vf_it, mesh_.calc_face_normal(*vf_it));
501 this->postprocess_collapse(ci);
504 if (!this->notify_observer(n_collapses))
508 if (oldCollapses == n_collapses) {
509 if (collapsesUnchanged ==
false) {
511 collapsesUnchanged =
true;
521 this->set_error_tolerance_factor(1.0);