00001 /*===========================================================================*\ 00002 * * 00003 * OpenMesh * 00004 * Copyright (C) 2001-2010 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: 325 $ * 00038 * $Date: 2010-06-17 12:45:58 +0200 (Do, 17 Jun 2010) $ * 00039 * * 00040 \*===========================================================================*/ 00041 00046 //============================================================================= 00047 // 00048 // CLASS LoopT 00049 // 00050 //============================================================================= 00051 00052 #ifndef OPENMESH_SUBDIVIDER_UNIFORM_COMPOSITELOOPT_HH 00053 #define OPENMESH_SUBDIVIDER_UNIFORM_COMPOSITELOOPT_HH 00054 00055 00056 //== INCLUDES ================================================================= 00057 00058 #include "Composite/CompositeT.hh" 00059 #include "Composite/CompositeTraits.hh" 00060 00061 00062 //== NAMESPACE ================================================================ 00063 00064 namespace OpenMesh { // BEGIN_NS_OPENMESH 00065 namespace Subdivider { // BEGIN_NS_DECIMATER 00066 namespace Uniform { // BEGIN_NS_DECIMATER 00067 00068 00069 //== CLASS DEFINITION ========================================================= 00070 00073 template <class MeshType, class RealType=float> 00074 class CompositeLoopT : public CompositeT<MeshType, RealType> 00075 { 00076 public: 00077 00078 typedef CompositeT<MeshType, RealType> Inherited; 00079 00080 public: 00081 00082 CompositeLoopT() : Inherited() {}; 00083 CompositeLoopT(MeshType& _mesh) : Inherited(_mesh) {}; 00084 ~CompositeLoopT() {} 00085 00086 public: 00087 00088 const char *name() const { return "Uniform Composite Loop"; } 00089 00090 protected: // inherited interface 00091 00092 void apply_rules(void) 00093 { 00094 Inherited::Tvv4(); 00095 Inherited::VdE(); 00096 Inherited::EVc(coeffs_); 00097 Inherited::VdE(); 00098 Inherited::EVc(coeffs_); 00099 } 00100 00101 protected: 00102 00103 typedef typename Inherited::Coeff Coeff; 00104 00105 00109 struct EVCoeff : public Coeff 00110 { 00111 EVCoeff() : Coeff() { init(50); } 00112 00113 void init(size_t _max_valence) 00114 { 00115 weights_.resize(_max_valence); 00116 std::generate(weights_.begin(), 00117 weights_.end(), compute_weight() ); 00118 } 00119 00120 double operator()(size_t _valence) { return weights_[_valence]; } 00121 00123 struct compute_weight 00124 { 00125 compute_weight() : val_(0) { } 00126 00127 double operator()(void) // Loop weights for non-boundary vertices 00128 { 00129 // 1 3 2 * pi 00130 // - * ( --- + cos ( ------- ) )² - 1.0 00131 // 2 2 valence 00132 double f1 = 1.5 + cos(2.0*M_PI/val_++); 00133 return 0.5 * f1 * f1 - 1.0; 00134 } 00135 size_t val_; 00136 00137 }; 00138 00139 std::vector<double> weights_; 00140 } coeffs_; 00141 00142 }; 00143 00144 //============================================================================= 00145 } // END_NS_UNIFORM 00146 } // END_NS_SUBDIVIDER 00147 } // END_NS_OPENMESH 00148 //============================================================================= 00149 #endif // OPENMESH_SUBDIVIDER_UNIFORM_COMPOSITELOOPT_HH defined 00150 //=============================================================================