00001 //============================================================================= 00002 // 00003 // OpenMesh 00004 // Copyright (C) 2003 by Computer Graphics Group, RWTH Aachen 00005 // www.openmesh.org 00006 // 00007 //----------------------------------------------------------------------------- 00008 // 00009 // License 00010 // 00011 // This library is free software; you can redistribute it and/or modify it 00012 // under the terms of the GNU Library General Public License as published 00013 // by the Free Software Foundation, version 2. 00014 // 00015 // This library is distributed in the hope that it will be useful, but 00016 // WITHOUT ANY WARRANTY; without even the implied warranty of 00017 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00018 // Library General Public License for more details. 00019 // 00020 // You should have received a copy of the GNU Library General Public 00021 // License along with this library; if not, write to the Free Software 00022 // Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. 00023 // 00024 //----------------------------------------------------------------------------- 00025 // 00026 // $Revision: 1.4 $ 00027 // $Date: 2004/01/15 15:34:49 $ 00028 // 00029 //============================================================================= 00030 00035 //============================================================================= 00036 // 00037 // CLASS LoopT 00038 // 00039 //============================================================================= 00040 00041 #ifndef OPENMESH_SUBDIVIDER_UNIFORM_COMPOSITELOOPT_HH 00042 #define OPENMESH_SUBDIVIDER_UNIFORM_COMPOSITELOOPT_HH 00043 00044 00045 //== INCLUDES ================================================================= 00046 00047 #include "Composite/CompositeT.hh" 00048 #include "Composite/CompositeTraits.hh" 00049 00050 00051 //== NAMESPACE ================================================================ 00052 00053 namespace OpenMesh { // BEGIN_NS_OPENMESH 00054 namespace Subdivider { // BEGIN_NS_DECIMATER 00055 namespace Uniform { // BEGIN_NS_DECIMATER 00056 00057 00058 //== CLASS DEFINITION ========================================================= 00059 00062 template <class MeshType, class RealType=float> 00063 class CompositeLoopT : public CompositeT<MeshType, RealType> 00064 { 00065 public: 00066 00067 typedef CompositeT<MeshType, RealType> Inherited; 00068 00069 public: 00070 00071 CompositeLoopT() : Inherited() {}; 00072 CompositeLoopT(MeshType& _mesh) : Inherited(_mesh) {}; 00073 ~CompositeLoopT() {} 00074 00075 public: 00076 00077 const char *name() const { return "Uniform Composite Loop"; } 00078 00079 protected: // inherited interface 00080 00081 void apply_rules(void) { Tvv4(); VdE(); EVc(coeffs_); VdE(); EVc(coeffs_); } 00082 00083 protected: 00084 00085 #if defined(OM_CC_MIPS) 00086 typedef typename Inherited::Coeff Coeff; 00087 #endif 00088 00092 struct EVCoeff : public Coeff 00093 { 00094 EVCoeff() : Coeff() { init(50); } 00095 00096 void init(size_t _max_valence) 00097 { 00098 weights_.resize(_max_valence); 00099 std::generate(weights_.begin(), 00100 weights_.end(), compute_weight() ); 00101 } 00102 00103 double operator()(size_t _valence) { return weights_[_valence]; } 00104 00106 struct compute_weight 00107 { 00108 compute_weight() : val_(0) { } 00109 00110 double operator()(void) // Loop weights for non-boundary vertices 00111 { 00112 // 1 3 2 * pi 00113 // - * ( --- + cos ( ------- ) )² - 1.0 00114 // 2 2 valence 00115 double f1 = 1.5 + cos(2.0*M_PI/val_++); 00116 return 0.5 * f1 * f1 - 1.0; 00117 } 00118 size_t val_; 00119 00120 }; 00121 00122 std::vector<double> weights_; 00123 } coeffs_; 00124 00125 }; 00126 00127 //============================================================================= 00128 } // END_NS_UNIFORM 00129 } // END_NS_SUBDIVIDER 00130 } // END_NS_OPENMESH 00131 //============================================================================= 00132 #endif // OPENMESH_SUBDIVIDER_UNIFORM_COMPOSITELOOPT_HH defined 00133 //=============================================================================