Developer Documentation
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
DiffGeoT.hh
1 /*===========================================================================*\
2 * *
3 * OpenFlipper *
4  * Copyright (c) 2001-2015, RWTH-Aachen University *
5  * Department of Computer Graphics and Multimedia *
6  * All rights reserved. *
7  * www.openflipper.org *
8  * *
9  *---------------------------------------------------------------------------*
10  * This file is part of OpenFlipper. *
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 * $Revision$ *
45 * $LastChangedBy$ *
46 * $Date$ *
47 * *
48 \*===========================================================================*/
49 
50 //=============================================================================
51 //
52 // CLASS DiffGeoT
53 //
54 //=============================================================================
55 
56 
57 #ifndef DIFFGEO_HH
58 #define DIFFGEO_HH
59 
60 
61 //== INCLUDES =================================================================
62 
63 #include <OpenMesh/Core/Utils/Property.hh>
64 
65 
66 //== NAMESPACES ===============================================================
67 
68 namespace Remeshing {
69 
70 //== CLASS DEFINITION =========================================================
71 
72 
73 template <class Mesh>
74 class DiffGeoT
75 {
76 public:
77 
78  typedef typename Mesh::Scalar Scalar;
79  typedef typename Mesh::VertexHandle VertexHandle;
80 
81 
82  DiffGeoT(Mesh& _mesh);
83  ~DiffGeoT();
84 
85 
86  void compute(unsigned int _post_smoothing_iters=0);
87 
88  void compute_edge_weights();
89  void compute_area();
90  void compute_gauss_curvature();
91  void compute_mean_curvature();
92  void post_smoothing(unsigned int _iters);
93 
94 
95  Scalar compute_area(VertexHandle _vh) const;
96 
97 
98  Scalar area(VertexHandle _vh) const {
99  return mesh_.property(area_, _vh);
100  }
101 
102  Scalar gauss_curvature(VertexHandle _vh) const {
103  return mesh_.property(gauss_curvature_, _vh);
104  }
105 
106  Scalar mean_curvature(VertexHandle _vh) const {
107  return mesh_.property(mean_curvature_, _vh);
108  }
109 
110  Scalar min_curvature(VertexHandle _vh) const {
111  const Scalar zero(0.0);
112  Scalar H = mean_curvature(_vh);
113  Scalar K = gauss_curvature(_vh);
114  return H - sqrt(std::max(zero, H*H-K));
115  }
116 
117  Scalar max_curvature(VertexHandle _vh) const {
118  const Scalar zero(0.0);
119  Scalar H = mean_curvature(_vh);
120  Scalar K = gauss_curvature(_vh);
121  return H + sqrt(std::max(zero, H*H-K));
122  }
123 
124 
125 private:
126 
127  Mesh& mesh_;
128 
130  OpenMesh::VPropHandleT<Scalar> gauss_curvature_;
131  OpenMesh::VPropHandleT<Scalar> mean_curvature_;
132  OpenMesh::EPropHandleT<Scalar> edge_weight_;
133 
134  bool weights_computed_, area_computed_;
135 };
136 
137 
138 //=============================================================================
139 } // namespace Remeshing
140 //=============================================================================
141 #if defined(INCLUDE_TEMPLATES) && !defined(DIFFGEOT_C)
142 #define DIFFGEO_TEMPLATES
143 #include "DiffGeoT.cc"
144 #endif
145 //=============================================================================
146 #endif // DIFFGEO_HH defined
147 //=============================================================================
148 
Kernel::Scalar Scalar
Scalar type.
Definition: PolyMeshT.hh:113
Kernel::VertexHandle VertexHandle
Handle for referencing the corresponding item.
Definition: PolyMeshT.hh:139