OpenMesh: OpenMesh/Core/Geometry/VectorT.hh Source File

| OpenMesh

VectorT.hh
1 /* ========================================================================= *
2 * *
3 * OpenMesh *
4 * Copyright (c) 2001-2022, RWTH-Aachen University *
5 * Department of Computer Graphics and Multimedia *
6 * All rights reserved. *
7 * www.openmesh.org *
8 * *
9 *---------------------------------------------------------------------------*
10 * This file is part of OpenMesh. *
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
45
46 //=============================================================================
47 //
48 // CLASS VectorT
49 //
50 //=============================================================================
51
52 // Don't parse this header file with doxygen since
53 // for some reason (obviously due to a bug in doxygen,
54 // bugreport: https://bugzilla.gnome.org/show_bug.cgi?id=629182)
55 // macro expansion and preprocessor defines
56 // don't work properly.
57
58 #if ((defined(_MSC_VER) && (_MSC_VER >= 1900)) || __cplusplus > 199711L || defined(__GXX_EXPERIMENTAL_CXX0X__)) && !defined(OPENMESH_VECTOR_LEGACY)
59 #include "Vector11T.hh"
60 #else
61 #ifndef DOXYGEN
62
63 #ifndef OPENMESH_VECTOR_HH
64 #define OPENMESH_VECTOR_HH
65
66
67 //== INCLUDES =================================================================
68
69 #include <OpenMesh/Core/System/config.h>
70 #include <ostream>
71 #include <cmath>
72 #include <cassert>
73 #include <cstring>
74
75 #if defined(__GNUC__) && defined(__SSE__)
76 #include <xmmintrin.h>
77 #endif
78
79 //== NAMESPACES ===============================================================
80
81
82 namespace OpenMesh {
83
84
85 //== CLASS DEFINITION =========================================================
86
87
97 template<typename Scalar, int N> class VectorDataT {
98 public:
99 Scalar values_[N];
100 };
101
102
103 #if defined(__GNUC__) && defined(__SSE__)
104
106 template<> class VectorDataT<float, 4> {
107 public:
108 union {
109 __m128 m128;
110 float values_[4];
111 };
112 };
113
114 #endif
115
116
117
118
119 //== CLASS DEFINITION =========================================================
120
121
122 #define DIM N
123 #define TEMPLATE_HEADER template <typename Scalar, int N>
124 #define CLASSNAME VectorT
125 #define DERIVED VectorDataT<Scalar,N>
126 #define unroll(expr) for (int i=0; i<N; ++i) expr(i)
127
133 #include "VectorT_inc.hh"
134
135 #undef DIM
136 #undef TEMPLATE_HEADER
137 #undef CLASSNAME
138 #undef DERIVED
139 #undef unroll
140
141
142
143
144 //== PARTIAL TEMPLATE SPECIALIZATIONS =========================================
145 #if OM_PARTIAL_SPECIALIZATION
146
147
148 #define TEMPLATE_HEADER template <typename Scalar>
149 #define CLASSNAME VectorT<Scalar,DIM>
150 #define DERIVED VectorDataT<Scalar,DIM>
151
152
153 #define DIM 2
154 #define unroll(expr) expr(0) expr(1)
155 #define unroll_comb(expr, op) expr(0) op expr(1)
156 #define unroll_csv(expr) expr(0), expr(1)
157 #include "VectorT_inc.hh"
158 #undef DIM
159 #undef unroll
160 #undef unroll_comb
161 #undef unroll_csv
162
163
164 #define DIM 3
165 #define unroll(expr) expr(0) expr(1) expr(2)
166 #define unroll_comb(expr, op) expr(0) op expr(1) op expr(2)
167 #define unroll_csv(expr) expr(0), expr(1), expr(2)
168 #include "VectorT_inc.hh"
169 #undef DIM
170 #undef unroll
171 #undef unroll_comb
172 #undef unroll_csv
173
174
175 #define DIM 4
176 #define unroll(expr) expr(0) expr(1) expr(2) expr(3)
177 #define unroll_comb(expr, op) expr(0) op expr(1) op expr(2) op expr(3)
178 #define unroll_csv(expr) expr(0), expr(1), expr(2), expr(3)
179 #include "VectorT_inc.hh"
180 #undef DIM
181 #undef unroll
182 #undef unroll_comb
183 #undef unroll_csv
184
185 #define DIM 5
186 #define unroll(expr) expr(0) expr(1) expr(2) expr(3) expr(4)
187 #define unroll_comb(expr, op) expr(0) op expr(1) op expr(2) op expr(3) op expr(4)
188 #define unroll_csv(expr) expr(0), expr(1), expr(2), expr(3), expr(4)
189 #include "VectorT_inc.hh"
190 #undef DIM
191 #undef unroll
192 #undef unroll_comb
193 #undef unroll_csv
194
195 #define DIM 6
196 #define unroll(expr) expr(0) expr(1) expr(2) expr(3) expr(4) expr(5)
197 #define unroll_comb(expr, op) expr(0) op expr(1) op expr(2) op expr(3) op expr(4) op expr(5)
198 #define unroll_csv(expr) expr(0), expr(1), expr(2), expr(3), expr(4), expr(5)
199 #include "VectorT_inc.hh"
200 #undef DIM
201 #undef unroll
202 #undef unroll_comb
203 #undef unroll_csv
204
205
206 #undef TEMPLATE_HEADER
207 #undef CLASSNAME
208 #undef DERIVED
209
210
211
212
213 //== FULL TEMPLATE SPECIALIZATIONS ============================================
214 #else
215
217 template<>
218 inline VectorT<float,3>
219 VectorT<float,3>::operator%(const VectorT<float,3>& _rhs) const
220 {
221 return
222 VectorT<float,3>(values_[1]*_rhs.values_[2]-values_[2]*_rhs.values_[1],
223 values_[2]*_rhs.values_[0]-values_[0]*_rhs.values_[2],
224 values_[0]*_rhs.values_[1]-values_[1]*_rhs.values_[0]);
225 }
226
227
229 template<>
230 inline VectorT<double,3>
231 VectorT<double,3>::operator%(const VectorT<double,3>& _rhs) const
232 {
233 return
234 VectorT<double,3>(values_[1]*_rhs.values_[2]-values_[2]*_rhs.values_[1],
235 values_[2]*_rhs.values_[0]-values_[0]*_rhs.values_[2],
236 values_[0]*_rhs.values_[1]-values_[1]*_rhs.values_[0]);
237 }
238
239 #endif
240
241
242
243 //== GLOBAL FUNCTIONS =========================================================
244
245
248 template<typename Scalar1, typename Scalar2,int N>
249 inline VectorT<Scalar1,N> operator*(Scalar2 _s, const VectorT<Scalar1,N>& _v) {
250 return _v*_s;
251 }
252
253
256 template<typename Scalar, int N>
257 inline Scalar
258 dot(const VectorT<Scalar,N>& _v1, const VectorT<Scalar,N>& _v2) {
259 return (_v1 | _v2);
260 }
261
262
265 template<typename Scalar, int N>
266 inline VectorT<Scalar,N>
267 cross(const VectorT<Scalar,N>& _v1, const VectorT<Scalar,N>& _v2) {
268 return (_v1 % _v2);
269 }
270
271
274 template<typename Scalar, int DIM>
275 Scalar norm(const VectorT<Scalar, DIM>& _v) {
276 return _v.norm();
277 }
278
279
282 template<typename Scalar, int DIM>
283 Scalar sqrnorm(const VectorT<Scalar, DIM>& _v) {
284 return _v.sqrnorm();
285 }
286
287
290 template<typename Scalar, int DIM, typename OtherScalar>
291 VectorT<Scalar, DIM>& vectorize(VectorT<Scalar, DIM>& _v, OtherScalar const& _val) {
292 return _v.vectorize(_val);
293 }
294
295
298 template<typename Scalar, int DIM>
299 VectorT<Scalar, DIM>& normalize(VectorT<Scalar, DIM>& _v) {
300 return _v.normalize();
301 }
302
303
306 template<typename Scalar, int DIM>
307 VectorT<Scalar, DIM>& maximize(VectorT<Scalar, DIM>& _v1, VectorT<Scalar, DIM>& _v2) {
308 return _v1.maximize(_v2);
309 }
310
311
314 template<typename Scalar, int DIM>
315 VectorT<Scalar, DIM>& minimize(VectorT<Scalar, DIM>& _v1, VectorT<Scalar, DIM>& _v2) {
316 return _v1.minimize(_v2);
317 }
318
319
322 template<typename Scalar, int DIM>
323 VectorT<Scalar, DIM> max(VectorT<Scalar, DIM>& _v1, VectorT<Scalar, DIM>& _v2) {
324 return VectorT<Scalar, DIM>(_v1).maximize(_v2);
325 }
326
327
330 template<typename Scalar, int DIM>
331 VectorT<Scalar, DIM> min(VectorT<Scalar, DIM>& _v1, VectorT<Scalar, DIM>& _v2) {
332 return VectorT<Scalar, DIM>(_v1).minimize(_v2);
333 }
334
335
336 //== TYPEDEFS =================================================================
337
339 typedef VectorT<signed char,1> Vec1c;
341 typedef VectorT<unsigned char,1> Vec1uc;
343 typedef VectorT<signed short int,1> Vec1s;
345 typedef VectorT<unsigned short int,1> Vec1us;
347 typedef VectorT<signed int,1> Vec1i;
349 typedef VectorT<unsigned int,1> Vec1ui;
351 typedef VectorT<float,1> Vec1f;
353 typedef VectorT<double,1> Vec1d;
354
356 typedef VectorT<signed char,2> Vec2c;
358 typedef VectorT<unsigned char,2> Vec2uc;
360 typedef VectorT<signed short int,2> Vec2s;
362 typedef VectorT<unsigned short int,2> Vec2us;
364 typedef VectorT<signed int,2> Vec2i;
366 typedef VectorT<unsigned int,2> Vec2ui;
368 typedef VectorT<float,2> Vec2f;
370 typedef VectorT<double,2> Vec2d;
371
373 typedef VectorT<signed char,3> Vec3c;
375 typedef VectorT<unsigned char,3> Vec3uc;
377 typedef VectorT<signed short int,3> Vec3s;
379 typedef VectorT<unsigned short int,3> Vec3us;
381 typedef VectorT<signed int,3> Vec3i;
383 typedef VectorT<unsigned int,3> Vec3ui;
385 typedef VectorT<float,3> Vec3f;
387 typedef VectorT<double,3> Vec3d;
389 typedef VectorT<bool,3> Vec3b;
390
392 typedef VectorT<signed char,4> Vec4c;
394 typedef VectorT<unsigned char,4> Vec4uc;
396 typedef VectorT<signed short int,4> Vec4s;
398 typedef VectorT<unsigned short int,4> Vec4us;
400 typedef VectorT<signed int,4> Vec4i;
402 typedef VectorT<unsigned int,4> Vec4ui;
404 typedef VectorT<float,4> Vec4f;
406 typedef VectorT<double,4> Vec4d;
407
409 typedef VectorT<signed char, 5> Vec5c;
411 typedef VectorT<unsigned char, 5> Vec5uc;
413 typedef VectorT<signed short int, 5> Vec5s;
415 typedef VectorT<unsigned short int, 5> Vec5us;
417 typedef VectorT<signed int, 5> Vec5i;
419 typedef VectorT<unsigned int, 5> Vec5ui;
421 typedef VectorT<float, 5> Vec5f;
423 typedef VectorT<double, 5> Vec5d;
424
426 typedef VectorT<signed char,6> Vec6c;
428 typedef VectorT<unsigned char,6> Vec6uc;
430 typedef VectorT<signed short int,6> Vec6s;
432 typedef VectorT<unsigned short int,6> Vec6us;
434 typedef VectorT<signed int,6> Vec6i;
436 typedef VectorT<unsigned int,6> Vec6ui;
438 typedef VectorT<float,6> Vec6f;
440 typedef VectorT<double,6> Vec6d;
441
442
443 //=============================================================================
444 } // namespace OpenMesh
445 //=============================================================================
446
447
448 #endif // OPENMESH_VECTOR_HH defined
449 //=============================================================================
450 #endif // DOXYGEN
451 #endif // C++11
OpenMesh
Contains all the mesh ingredients like the polygonal mesh, the triangle mesh, different mesh kernels ...
Definition: MeshItems.hh:59
OpenMesh::Vec3c
VectorT< signed char, 3 > Vec3c
3-byte signed vector
Definition: Vector11T.hh:838
OpenMesh::Vec3d
VectorT< double, 3 > Vec3d
3-double vector
Definition: Vector11T.hh:852
OpenMesh::Vec1i
VectorT< signed int, 1 > Vec1i
1-int signed vector
Definition: Vector11T.hh:812
OpenMesh::Vec5uc
VectorT< unsigned char, 5 > Vec5uc
5-byte unsigned vector
Definition: Vector11T.hh:876
OpenMesh::Vec2ui
VectorT< unsigned int, 2 > Vec2ui
2-int unsigned vector
Definition: Vector11T.hh:831
OpenMesh::Vec1d
VectorT< double, 1 > Vec1d
1-double vector
Definition: Vector11T.hh:818
OpenMesh::Vec4d
VectorT< double, 4 > Vec4d
4-double vector
Definition: Vector11T.hh:871
OpenMesh::Vec2f
VectorT< float, 2 > Vec2f
2-float vector
Definition: Vector11T.hh:833
OpenMesh::Vec3s
VectorT< signed short int, 3 > Vec3s
3-short signed vector
Definition: Vector11T.hh:842
OpenMesh::Vec2d
VectorT< double, 2 > Vec2d
2-double vector
Definition: Vector11T.hh:835
OpenMesh::Vec5i
VectorT< signed int, 5 > Vec5i
5-int signed vector
Definition: Vector11T.hh:882
OpenMesh::Vec4c
VectorT< signed char, 4 > Vec4c
4-byte signed vector
Definition: Vector11T.hh:857
OpenMesh::Vec5c
VectorT< signed char, 5 > Vec5c
5-byte signed vector
Definition: Vector11T.hh:874
OpenMesh::Vec5ui
VectorT< unsigned int, 5 > Vec5ui
5-int unsigned vector
Definition: Vector11T.hh:884
OpenMesh::Vec5s
VectorT< signed short int, 5 > Vec5s
5-short signed vector
Definition: Vector11T.hh:878
OpenMesh::Vec4ui
VectorT< unsigned int, 4 > Vec4ui
4-int unsigned vector
Definition: Vector11T.hh:867
OpenMesh::Vec3us
VectorT< unsigned short int, 3 > Vec3us
3-short unsigned vector
Definition: Vector11T.hh:844
OpenMesh::Vec3uc
VectorT< unsigned char, 3 > Vec3uc
3-byte unsigned vector
Definition: Vector11T.hh:840
OpenMesh::Vec2s
VectorT< signed short int, 2 > Vec2s
2-short signed vector
Definition: Vector11T.hh:825
OpenMesh::Vec2i
VectorT< signed int, 2 > Vec2i
2-int signed vector
Definition: Vector11T.hh:829
OpenMesh::operator*
auto operator*(const OtherScalar &_s, const VectorT< Scalar, DIM > &rhs) -> decltype(rhs.operator*(_s))
Component wise multiplication from the left.
Definition: Vector11T.hh:692
OpenMesh::Vec2c
VectorT< signed char, 2 > Vec2c
2-byte signed vector
Definition: Vector11T.hh:821
OpenMesh::Vec2us
VectorT< unsigned short int, 2 > Vec2us
2-short unsigned vector
Definition: Vector11T.hh:827
OpenMesh::Vec6c
VectorT< signed char, 6 > Vec6c
6-byte signed vector
Definition: Vector11T.hh:891
OpenMesh::Vec6i
VectorT< signed int, 6 > Vec6i
6-int signed vector
Definition: Vector11T.hh:899
OpenMesh::Vec4i
VectorT< signed int, 4 > Vec4i
4-int signed vector
Definition: Vector11T.hh:865
OpenMesh::Vec4s
VectorT< signed short int, 4 > Vec4s
4-short signed vector
Definition: Vector11T.hh:861
OpenMesh::Vec6uc
VectorT< unsigned char, 6 > Vec6uc
6-byte unsigned vector
Definition: Vector11T.hh:893
OpenMesh::Vec6f
VectorT< float, 6 > Vec6f
6-float vector
Definition: Vector11T.hh:903
OpenMesh::Vec5f
VectorT< float, 5 > Vec5f
5-float vector
Definition: Vector11T.hh:886
OpenMesh::Vec2uc
VectorT< unsigned char, 2 > Vec2uc
2-byte unsigned vector
Definition: Vector11T.hh:823
OpenMesh::Vec6ui
VectorT< unsigned int, 6 > Vec6ui
6-int unsigned vector
Definition: Vector11T.hh:901
OpenMesh::Vec4f
VectorT< float, 4 > Vec4f
4-float vector
Definition: Vector11T.hh:869
OpenMesh::Vec1uc
VectorT< unsigned char, 1 > Vec1uc
1-byte unsigned vector
Definition: Vector11T.hh:806
OpenMesh::Vec4uc
VectorT< unsigned char, 4 > Vec4uc
4-byte unsigned vector
Definition: Vector11T.hh:859
OpenMesh::Vec1us
VectorT< unsigned short int, 1 > Vec1us
1-short unsigned vector
Definition: Vector11T.hh:810
OpenMesh::Vec1f
VectorT< float, 1 > Vec1f
1-float vector
Definition: Vector11T.hh:816
OpenMesh::Vec6d
VectorT< double, 6 > Vec6d
6-double vector
Definition: Vector11T.hh:905
OpenMesh::Vec1s
VectorT< signed short int, 1 > Vec1s
1-short signed vector
Definition: Vector11T.hh:808
OpenMesh::Vec3b
VectorT< bool, 3 > Vec3b
3-bool vector
Definition: Vector11T.hh:854
OpenMesh::Vec1ui
VectorT< unsigned int, 1 > Vec1ui
1-int unsigned vector
Definition: Vector11T.hh:814
OpenMesh::Vec1c
VectorT< signed char, 1 > Vec1c
1-byte signed vector
Definition: Vector11T.hh:804
OpenMesh::Vec5d
VectorT< double, 5 > Vec5d
5-double vector
Definition: Vector11T.hh:888
OpenMesh::Vec4us
VectorT< unsigned short int, 4 > Vec4us
4-short unsigned vector
Definition: Vector11T.hh:863
OpenMesh::Vec3ui
VectorT< unsigned int, 3 > Vec3ui
3-int unsigned vector
Definition: Vector11T.hh:848
OpenMesh::Vec6us
VectorT< unsigned short int, 6 > Vec6us
6-short unsigned vector
Definition: Vector11T.hh:897
OpenMesh::Vec5us
VectorT< unsigned short int, 5 > Vec5us
5-short unsigned vector
Definition: Vector11T.hh:880
OpenMesh::Vec3f
VectorT< float, 3 > Vec3f
3-float vector
Definition: Vector11T.hh:850
OpenMesh::Vec6s
VectorT< signed short int, 6 > Vec6s
6-short signed vector
Definition: Vector11T.hh:895
OpenMesh::Vec3i
VectorT< signed int, 3 > Vec3i
3-int signed vector
Definition: Vector11T.hh:846
OpenMesh::VectorT::normalize
VectorT< Scalar, DIM > & normalize(VectorT< Scalar, DIM > &_v)
non-member normalize
Definition: Vector11T.hh:768
OpenMesh::VectorT::dot
Scalar dot(const VectorT< Scalar, DIM > &_v1, const VectorT< Scalar, DIM > &_v2)
symmetric version of the dot product
Definition: Vector11T.hh:724
OpenMesh::VectorT::minimize
VectorT< Scalar, DIM > & minimize(VectorT< Scalar, DIM > &_v1, VectorT< Scalar, DIM > &_v2)
non-member minimize
Definition: Vector11T.hh:782
OpenMesh::VectorT::vectorize
VectorT< Scalar, DIM > & vectorize(VectorT< Scalar, DIM > &_v, OtherScalar const &_val)
non-member vectorize
Definition: Vector11T.hh:761
OpenMesh::VectorT::operator%
auto operator%(const VectorT< OtherScalar, DIM > &_rhs) const -> typename std::enable_if< DIM==3, VectorT< decltype((*this)[0] *_rhs[0] -(*this)[0] *_rhs[0]), DIM >>::type
cross product: only defined for Vec3* as specialization
Definition: Vector11T.hh:374
OpenMesh::VectorT::min
VectorT< Scalar, DIM > min(const VectorT< Scalar, DIM > &_v1, const VectorT< Scalar, DIM > &_v2)
non-member min
Definition: Vector11T.hh:796
OpenMesh::VectorT::sqrnorm
Scalar sqrnorm(const VectorT< Scalar, DIM > &_v)
non-member sqrnorm
Definition: Vector11T.hh:755
OpenMesh::VectorT::cross
auto cross(const VectorT< LScalar, DIM > &_v1, const VectorT< RScalar, DIM > &_v2) -> decltype(_v1 % _v2)
symmetric version of the cross product
Definition: Vector11T.hh:732
OpenMesh::VectorT::maximize
VectorT< Scalar, DIM > & maximize(VectorT< Scalar, DIM > &_v1, VectorT< Scalar, DIM > &_v2)
non-member maximize
Definition: Vector11T.hh:775
OpenMesh::VectorT::max
VectorT< Scalar, DIM > max(const VectorT< Scalar, DIM > &_v1, const VectorT< Scalar, DIM > &_v2)
non-member max
Definition: Vector11T.hh:789
OpenMesh::VectorT::norm
Scalar norm(const VectorT< Scalar, DIM > &_v)
non-member norm
Definition: Vector11T.hh:748
VectorT
Definition: VectorT_inc.hh:68
VectorT::minimize
vector_type & minimize(const vector_type &_rhs)
minimize values: same as *this = min(*this, _rhs), but faster
Definition: VectorT_inc.hh:547
VectorT::norm
Scalar norm() const
compute euclidean norm
Definition: VectorT_inc.hh:410
VectorT::operator%
VectorT< Scalar, 3 > operator%(const VectorT< Scalar, 3 > &_rhs) const
cross product: only defined for Vec3* as specialization
VectorT::sqrnorm
Scalar sqrnorm() const
compute squared euclidean norm
Definition: VectorT_inc.hh:414
VectorT::maximize
vector_type & maximize(const vector_type &_rhs)
maximize values: same as *this = max(*this, _rhs), but faster
Definition: VectorT_inc.hh:564
VectorT::vectorize
vector_type & vectorize(const Scalar &_s)
store the same value in each component (e.g. to clear all entries)
Definition: VectorT_inc.hh:605
VectorT::normalize
vector_type & normalize()
normalize vector, return normalized vector
Definition: VectorT_inc.hh:432