Developer Documentation
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
smooth_mesh.hh
1 /* ========================================================================= *
2  * *
3  * OpenMesh *
4  * Copyright (c) 2001-2015, 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  * $Revision$ *
45  * $Date$ *
46  * *
47 \*===========================================================================*/
48 
49 #ifndef SMOOTH_MESH_HH
50 #define SMOOTH_MESH_HH
51 
52 //== INCLUDES =================================================================
53 
54 #include <OpenMesh/Core/Utils/Property.hh>
55 
56 //== NAMESPACE ================================================================
57 
58 namespace OpenMesh { //BEGIN_NS_OPENMESH
59 
60 template <class _Mesh, class _PropertyHandle>
61 void smooth_mesh_property(unsigned int _n_iters, _Mesh& _m, _PropertyHandle _pph)
62 {
63  typedef typename _PropertyHandle::Value Value;
64 
65  std::vector<Value> temp_values(_m.n_vertices());
66 
67  for (unsigned int i=0; i < _n_iters; ++i)
68  {
69  for ( typename _Mesh::ConstVertexIter cv_it = _m.vertices_begin();
70  cv_it != _m.vertices_end(); ++cv_it)
71  {
72  unsigned int valence = 0;
73 
74  Value& temp_value = temp_values[cv_it->idx()];
75 
76  temp_value.vectorize(0);
77 
78  for ( typename _Mesh::ConstVertexVertexIter cvv_it = _m.cvv_iter(cv_it);
79  cvv_it; ++cvv_it)
80  {
81  temp_value += _m.property(_pph,cvv_it);
82  ++valence;
83  }
84  if (valence > 0)
85  {//guard against isolated vertices
86  temp_value *= (typename Value::value_type)(1.0 / valence);
87  }
88  else
89  {
90  temp_value = _m.property(_pph, cv_it);
91  }
92  }
93 
94  for ( typename _Mesh::ConstVertexIter cv_it = _m.vertices_begin();
95  cv_it != _m.vertices_end(); ++cv_it)
96  {
97  _m.property(_pph,cv_it) = temp_values[cv_it->idx()];
98  }
99  }
100 }
101 
102 template <class _Mesh>
103 void smooth_mesh(_Mesh& _m, uint _n_iters)
104 {
105  smooth_mesh_property(_n_iters, _m, _m.points_pph());
106 }
107 
108 };//namespace OpenMesh
109 
110 #endif//SMOOTH_MESH_HH