OpenMesh
OpenMesh/Tools/Subdivider/Uniform/CompositeLoopT.hh
Go to the documentation of this file.
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 
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 //=============================================================================