OpenMesh
OpenMesh/Tools/Smoother/smooth_mesh.hh
00001 /*===========================================================================*\
00002  *                                                                           *
00003  *                               OpenMesh                                    *
00004  *      Copyright (C) 2001-2011 by Computer Graphics Group, RWTH Aachen      *
00005  *                           www.openmesh.org                                *
00006  *                                                                           *
00007  *---------------------------------------------------------------------------* 
00008  *  This file is part of OpenMesh.                                           *
00009  *                                                                           *
00010  *  OpenMesh is free software: you can redistribute it and/or modify         * 
00011  *  it under the terms of the GNU Lesser General Public License as           *
00012  *  published by the Free Software Foundation, either version 3 of           *
00013  *  the License, or (at your option) any later version with the              *
00014  *  following exceptions:                                                    *
00015  *                                                                           *
00016  *  If other files instantiate templates or use macros                       *
00017  *  or inline functions from this file, or you compile this file and         *
00018  *  link it with other files to produce an executable, this file does        *
00019  *  not by itself cause the resulting executable to be covered by the        *
00020  *  GNU Lesser General Public License. This exception does not however       *
00021  *  invalidate any other reasons why the executable file might be            *
00022  *  covered by the GNU Lesser General Public License.                        *
00023  *                                                                           *
00024  *  OpenMesh is distributed in the hope that it will be useful,              *
00025  *  but WITHOUT ANY WARRANTY; without even the implied warranty of           *
00026  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the            *
00027  *  GNU Lesser General Public License for more details.                      *
00028  *                                                                           *
00029  *  You should have received a copy of the GNU LesserGeneral Public          *
00030  *  License along with OpenMesh.  If not,                                    *
00031  *  see <http://www.gnu.org/licenses/>.                                      *
00032  *                                                                           *
00033 \*===========================================================================*/ 
00034 
00035 /*===========================================================================*\
00036  *                                                                           *             
00037  *   $Revision: 362 $                                                         *
00038  *   $Date: 2011-01-26 10:21:12 +0100 (Mi, 26 Jan 2011) $                   *
00039  *                                                                           *
00040 \*===========================================================================*/
00041 
00042 #ifndef SMOOTH_MESH_HH
00043 #define SMOOTH_MESH_HH
00044 
00045 //== INCLUDES =================================================================
00046 
00047 #include <OpenMesh/Core/Utils/Property.hh>
00048 
00049 //== NAMESPACE ================================================================
00050 
00051 namespace OpenMesh { //BEGIN_NS_OPENMESH
00052 
00053 template <class _Mesh, class _PropertyHandle>
00054 void smooth_mesh_property(unsigned int _n_iters, _Mesh& _m, _PropertyHandle _pph)
00055 {
00056   typedef typename _PropertyHandle::Value   Value;
00057 
00058   std::vector<Value> temp_values(_m.n_vertices());
00059 
00060   for (unsigned int i=0; i < _n_iters; ++i)
00061   {
00062     for ( typename _Mesh::ConstVertexIter cv_it = _m.vertices_begin();
00063           cv_it != _m.vertices_end(); ++cv_it)
00064     {
00065       unsigned int valence = 0;
00066 
00067       Value& temp_value = temp_values[cv_it.handle().idx()];
00068 
00069       temp_value.vectorize(0);
00070 
00071       for ( typename _Mesh::ConstVertexVertexIter cvv_it = _m.cvv_iter(cv_it);
00072             cvv_it; ++cvv_it)
00073       {
00074         temp_value += _m.property(_pph,cvv_it);
00075         ++valence;
00076       }
00077       if (valence > 0)
00078       {//guard against isolated vertices
00079         temp_value *= (typename Value::value_type)(1.0 / valence);
00080       }
00081       else
00082       {
00083         temp_value = _m.property(_pph, cv_it);
00084       }
00085     }
00086 
00087     for ( typename _Mesh::ConstVertexIter cv_it = _m.vertices_begin();
00088           cv_it != _m.vertices_end(); ++cv_it)
00089     {
00090       _m.property(_pph,cv_it) = temp_values[cv_it.handle().idx()];
00091     }
00092   }
00093 }
00094 
00095 template <class _Mesh>
00096 void smooth_mesh(_Mesh& _m, uint _n_iters)
00097 {
00098   smooth_mesh_property(_n_iters, _m, _m.points_pph());
00099 }
00100 
00101 };//namespace OpenMesh
00102 
00103 #endif//SMOOTH_MESH_HH