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