Please, help us to better know about our user community by answering the following short survey: https://forms.gle/wpyrxWi18ox9Z5ae9
Eigen  3.4.0
 
Loading...
Searching...
No Matches
Meta.h
1// This file is part of Eigen, a lightweight C++ template library
2// for linear algebra.
3//
4// Copyright (C) 2008-2015 Gael Guennebaud <gael.guennebaud@inria.fr>
5// Copyright (C) 2006-2008 Benoit Jacob <jacob.benoit.1@gmail.com>
6//
7// This Source Code Form is subject to the terms of the Mozilla
8// Public License v. 2.0. If a copy of the MPL was not distributed
9// with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
10
11#ifndef EIGEN_META_H
12#define EIGEN_META_H
13
14#if defined(EIGEN_GPU_COMPILE_PHASE)
15
16 #include <cfloat>
17
18 #if defined(EIGEN_CUDA_ARCH)
19 #include <math_constants.h>
20 #endif
21
22 #if defined(EIGEN_HIP_DEVICE_COMPILE)
23 #include "Eigen/src/Core/arch/HIP/hcc/math_constants.h"
24 #endif
25
26#endif
27
28// Recent versions of ICC require <cstdint> for pointer types below.
29#define EIGEN_ICC_NEEDS_CSTDINT (EIGEN_COMP_ICC>=1600 && EIGEN_COMP_CXXVER >= 11)
30
31// Define portable (u)int{32,64} types
32#if EIGEN_HAS_CXX11 || EIGEN_ICC_NEEDS_CSTDINT
33#include <cstdint>
34namespace Eigen {
35namespace numext {
36typedef std::uint8_t uint8_t;
37typedef std::int8_t int8_t;
38typedef std::uint16_t uint16_t;
39typedef std::int16_t int16_t;
40typedef std::uint32_t uint32_t;
41typedef std::int32_t int32_t;
42typedef std::uint64_t uint64_t;
43typedef std::int64_t int64_t;
44}
45}
46#else
47// Without c++11, all compilers able to compile Eigen also
48// provide the C99 stdint.h header file.
49#include <stdint.h>
50namespace Eigen {
51namespace numext {
52typedef ::uint8_t uint8_t;
53typedef ::int8_t int8_t;
54typedef ::uint16_t uint16_t;
55typedef ::int16_t int16_t;
56typedef ::uint32_t uint32_t;
57typedef ::int32_t int32_t;
58typedef ::uint64_t uint64_t;
59typedef ::int64_t int64_t;
60}
61}
62#endif
63
64namespace Eigen {
65
66typedef EIGEN_DEFAULT_DENSE_INDEX_TYPE DenseIndex;
67
74typedef EIGEN_DEFAULT_DENSE_INDEX_TYPE Index;
75
76namespace internal {
77
85// Only recent versions of ICC complain about using ptrdiff_t to hold pointers,
86// and older versions do not provide *intptr_t types.
87#if EIGEN_ICC_NEEDS_CSTDINT
88typedef std::intptr_t IntPtr;
89typedef std::uintptr_t UIntPtr;
90#else
91typedef std::ptrdiff_t IntPtr;
92typedef std::size_t UIntPtr;
93#endif
94#undef EIGEN_ICC_NEEDS_CSTDINT
95
96struct true_type { enum { value = 1 }; };
97struct false_type { enum { value = 0 }; };
98
99template<bool Condition>
100struct bool_constant;
101
102template<>
103struct bool_constant<true> : true_type {};
104
105template<>
106struct bool_constant<false> : false_type {};
107
108template<bool Condition, typename Then, typename Else>
109struct conditional { typedef Then type; };
110
111template<typename Then, typename Else>
112struct conditional <false, Then, Else> { typedef Else type; };
113
114template<typename T> struct remove_reference { typedef T type; };
115template<typename T> struct remove_reference<T&> { typedef T type; };
116
117template<typename T> struct remove_pointer { typedef T type; };
118template<typename T> struct remove_pointer<T*> { typedef T type; };
119template<typename T> struct remove_pointer<T*const> { typedef T type; };
120
121template <class T> struct remove_const { typedef T type; };
122template <class T> struct remove_const<const T> { typedef T type; };
123template <class T> struct remove_const<const T[]> { typedef T type[]; };
124template <class T, unsigned int Size> struct remove_const<const T[Size]> { typedef T type[Size]; };
125
126template<typename T> struct remove_all { typedef T type; };
127template<typename T> struct remove_all<const T> { typedef typename remove_all<T>::type type; };
128template<typename T> struct remove_all<T const&> { typedef typename remove_all<T>::type type; };
129template<typename T> struct remove_all<T&> { typedef typename remove_all<T>::type type; };
130template<typename T> struct remove_all<T const*> { typedef typename remove_all<T>::type type; };
131template<typename T> struct remove_all<T*> { typedef typename remove_all<T>::type type; };
132
133template<typename T> struct is_arithmetic { enum { value = false }; };
134template<> struct is_arithmetic<float> { enum { value = true }; };
135template<> struct is_arithmetic<double> { enum { value = true }; };
136template<> struct is_arithmetic<long double> { enum { value = true }; };
137template<> struct is_arithmetic<bool> { enum { value = true }; };
138template<> struct is_arithmetic<char> { enum { value = true }; };
139template<> struct is_arithmetic<signed char> { enum { value = true }; };
140template<> struct is_arithmetic<unsigned char> { enum { value = true }; };
141template<> struct is_arithmetic<signed short> { enum { value = true }; };
142template<> struct is_arithmetic<unsigned short>{ enum { value = true }; };
143template<> struct is_arithmetic<signed int> { enum { value = true }; };
144template<> struct is_arithmetic<unsigned int> { enum { value = true }; };
145template<> struct is_arithmetic<signed long> { enum { value = true }; };
146template<> struct is_arithmetic<unsigned long> { enum { value = true }; };
147
148template<typename T, typename U> struct is_same { enum { value = 0 }; };
149template<typename T> struct is_same<T,T> { enum { value = 1 }; };
150
151template< class T >
152struct is_void : is_same<void, typename remove_const<T>::type> {};
153
154#if EIGEN_HAS_CXX11
155template<> struct is_arithmetic<signed long long> { enum { value = true }; };
156template<> struct is_arithmetic<unsigned long long> { enum { value = true }; };
157using std::is_integral;
158#else
159template<typename T> struct is_integral { enum { value = false }; };
160template<> struct is_integral<bool> { enum { value = true }; };
161template<> struct is_integral<char> { enum { value = true }; };
162template<> struct is_integral<signed char> { enum { value = true }; };
163template<> struct is_integral<unsigned char> { enum { value = true }; };
164template<> struct is_integral<signed short> { enum { value = true }; };
165template<> struct is_integral<unsigned short> { enum { value = true }; };
166template<> struct is_integral<signed int> { enum { value = true }; };
167template<> struct is_integral<unsigned int> { enum { value = true }; };
168template<> struct is_integral<signed long> { enum { value = true }; };
169template<> struct is_integral<unsigned long> { enum { value = true }; };
170#if EIGEN_COMP_MSVC
171template<> struct is_integral<signed __int64> { enum { value = true }; };
172template<> struct is_integral<unsigned __int64> { enum { value = true }; };
173#endif
174#endif
175
176#if EIGEN_HAS_CXX11
177using std::make_unsigned;
178#else
179// TODO: Possibly improve this implementation of make_unsigned.
180// It is currently used only by
181// template<typename Scalar> struct random_default_impl<Scalar, false, true>.
182template<typename> struct make_unsigned;
183template<> struct make_unsigned<char> { typedef unsigned char type; };
184template<> struct make_unsigned<signed char> { typedef unsigned char type; };
185template<> struct make_unsigned<unsigned char> { typedef unsigned char type; };
186template<> struct make_unsigned<signed short> { typedef unsigned short type; };
187template<> struct make_unsigned<unsigned short> { typedef unsigned short type; };
188template<> struct make_unsigned<signed int> { typedef unsigned int type; };
189template<> struct make_unsigned<unsigned int> { typedef unsigned int type; };
190template<> struct make_unsigned<signed long> { typedef unsigned long type; };
191template<> struct make_unsigned<unsigned long> { typedef unsigned long type; };
192#if EIGEN_COMP_MSVC
193template<> struct make_unsigned<signed __int64> { typedef unsigned __int64 type; };
194template<> struct make_unsigned<unsigned __int64> { typedef unsigned __int64 type; };
195#endif
196
197// Some platforms define int64_t as `long long` even for C++03, where
198// `long long` is not guaranteed by the standard. In this case we are missing
199// the definition for make_unsigned. If we just define it, we run into issues
200// where `long long` doesn't exist in some compilers for C++03. We therefore add
201// the specialization for these platforms only.
202#if EIGEN_OS_MAC || EIGEN_COMP_MINGW
203template<> struct make_unsigned<unsigned long long> { typedef unsigned long long type; };
204template<> struct make_unsigned<long long> { typedef unsigned long long type; };
205#endif
206#endif
207
208template <typename T> struct add_const { typedef const T type; };
209template <typename T> struct add_const<T&> { typedef T& type; };
210
211template <typename T> struct is_const { enum { value = 0 }; };
212template <typename T> struct is_const<T const> { enum { value = 1 }; };
213
214template<typename T> struct add_const_on_value_type { typedef const T type; };
215template<typename T> struct add_const_on_value_type<T&> { typedef T const& type; };
216template<typename T> struct add_const_on_value_type<T*> { typedef T const* type; };
217template<typename T> struct add_const_on_value_type<T* const> { typedef T const* const type; };
218template<typename T> struct add_const_on_value_type<T const* const> { typedef T const* const type; };
219
220#if EIGEN_HAS_CXX11
221
222using std::is_convertible;
223
224#else
225
226template<typename From, typename To>
227struct is_convertible_impl
228{
229private:
230 struct any_conversion
231 {
232 template <typename T> any_conversion(const volatile T&);
233 template <typename T> any_conversion(T&);
234 };
235 struct yes {int a[1];};
236 struct no {int a[2];};
237
238 template<typename T>
239 static yes test(T, int);
240
241 template<typename T>
242 static no test(any_conversion, ...);
243
244public:
245 static typename internal::remove_reference<From>::type* ms_from;
246#ifdef __INTEL_COMPILER
247 #pragma warning push
248 #pragma warning ( disable : 2259 )
249#endif
250 enum { value = sizeof(test<To>(*ms_from, 0))==sizeof(yes) };
251#ifdef __INTEL_COMPILER
252 #pragma warning pop
253#endif
254};
255
256template<typename From, typename To>
257struct is_convertible
258{
259 enum { value = is_convertible_impl<From,To>::value };
260};
261
262template<typename T>
263struct is_convertible<T,T&> { enum { value = false }; };
264
265template<typename T>
266struct is_convertible<const T,const T&> { enum { value = true }; };
267
268#endif
269
273template<bool Condition, typename T=void> struct enable_if;
274
275template<typename T> struct enable_if<true,T>
276{ typedef T type; };
277
278#if defined(EIGEN_GPU_COMPILE_PHASE) && !EIGEN_HAS_CXX11
279#if !defined(__FLT_EPSILON__)
280#define __FLT_EPSILON__ FLT_EPSILON
281#define __DBL_EPSILON__ DBL_EPSILON
282#endif
283
284namespace device {
285
286template<typename T> struct numeric_limits
287{
288 EIGEN_DEVICE_FUNC
289 static EIGEN_CONSTEXPR T epsilon() { return 0; }
290 static T (max)() { assert(false && "Highest not supported for this type"); }
291 static T (min)() { assert(false && "Lowest not supported for this type"); }
292 static T infinity() { assert(false && "Infinity not supported for this type"); }
293 static T quiet_NaN() { assert(false && "quiet_NaN not supported for this type"); }
294};
295template<> struct numeric_limits<float>
296{
297 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
298 static float epsilon() { return __FLT_EPSILON__; }
299 EIGEN_DEVICE_FUNC
300 static float (max)() {
301 #if defined(EIGEN_CUDA_ARCH)
302 return CUDART_MAX_NORMAL_F;
303 #else
304 return HIPRT_MAX_NORMAL_F;
305 #endif
306 }
307 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
308 static float (min)() { return FLT_MIN; }
309 EIGEN_DEVICE_FUNC
310 static float infinity() {
311 #if defined(EIGEN_CUDA_ARCH)
312 return CUDART_INF_F;
313 #else
314 return HIPRT_INF_F;
315 #endif
316 }
317 EIGEN_DEVICE_FUNC
318 static float quiet_NaN() {
319 #if defined(EIGEN_CUDA_ARCH)
320 return CUDART_NAN_F;
321 #else
322 return HIPRT_NAN_F;
323 #endif
324 }
325};
326template<> struct numeric_limits<double>
327{
328 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
329 static double epsilon() { return __DBL_EPSILON__; }
330 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
331 static double (max)() { return DBL_MAX; }
332 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
333 static double (min)() { return DBL_MIN; }
334 EIGEN_DEVICE_FUNC
335 static double infinity() {
336 #if defined(EIGEN_CUDA_ARCH)
337 return CUDART_INF;
338 #else
339 return HIPRT_INF;
340 #endif
341 }
342 EIGEN_DEVICE_FUNC
343 static double quiet_NaN() {
344 #if defined(EIGEN_CUDA_ARCH)
345 return CUDART_NAN;
346 #else
347 return HIPRT_NAN;
348 #endif
349 }
350};
351template<> struct numeric_limits<int>
352{
353 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
354 static int epsilon() { return 0; }
355 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
356 static int (max)() { return INT_MAX; }
357 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
358 static int (min)() { return INT_MIN; }
359};
360template<> struct numeric_limits<unsigned int>
361{
362 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
363 static unsigned int epsilon() { return 0; }
364 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
365 static unsigned int (max)() { return UINT_MAX; }
366 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
367 static unsigned int (min)() { return 0; }
368};
369template<> struct numeric_limits<long>
370{
371 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
372 static long epsilon() { return 0; }
373 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
374 static long (max)() { return LONG_MAX; }
375 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
376 static long (min)() { return LONG_MIN; }
377};
378template<> struct numeric_limits<unsigned long>
379{
380 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
381 static unsigned long epsilon() { return 0; }
382 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
383 static unsigned long (max)() { return ULONG_MAX; }
384 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
385 static unsigned long (min)() { return 0; }
386};
387template<> struct numeric_limits<long long>
388{
389 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
390 static long long epsilon() { return 0; }
391 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
392 static long long (max)() { return LLONG_MAX; }
393 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
394 static long long (min)() { return LLONG_MIN; }
395};
396template<> struct numeric_limits<unsigned long long>
397{
398 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
399 static unsigned long long epsilon() { return 0; }
400 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
401 static unsigned long long (max)() { return ULLONG_MAX; }
402 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
403 static unsigned long long (min)() { return 0; }
404};
405template<> struct numeric_limits<bool>
406{
407 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
408 static bool epsilon() { return false; }
409 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
410 static bool (max)() { return true; }
411 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
412 static bool (min)() { return false; }
413};
414
415}
416
417#endif // defined(EIGEN_GPU_COMPILE_PHASE) && !EIGEN_HAS_CXX11
418
422class noncopyable
423{
424 EIGEN_DEVICE_FUNC noncopyable(const noncopyable&);
425 EIGEN_DEVICE_FUNC const noncopyable& operator=(const noncopyable&);
426protected:
427 EIGEN_DEVICE_FUNC noncopyable() {}
428 EIGEN_DEVICE_FUNC ~noncopyable() {}
429};
430
445template<typename T, typename EnableIf = void> struct array_size {
446 enum { value = Dynamic };
447};
448
449template<typename T> struct array_size<T,typename internal::enable_if<((T::SizeAtCompileTime&0)==0)>::type> {
450 enum { value = T::SizeAtCompileTime };
451};
452
453template<typename T, int N> struct array_size<const T (&)[N]> {
454 enum { value = N };
455};
456template<typename T, int N> struct array_size<T (&)[N]> {
457 enum { value = N };
458};
459
460#if EIGEN_HAS_CXX11
461template<typename T, std::size_t N> struct array_size<const std::array<T,N> > {
462 enum { value = N };
463};
464template<typename T, std::size_t N> struct array_size<std::array<T,N> > {
465 enum { value = N };
466};
467#endif
468
478template<typename T>
479EIGEN_CONSTEXPR Index size(const T& x) { return x.size(); }
480
481template<typename T,std::size_t N>
482EIGEN_CONSTEXPR Index size(const T (&) [N]) { return N; }
483
499#if EIGEN_HAS_STD_INVOKE_RESULT
500template<typename T> struct result_of;
501
502template<typename F, typename... ArgTypes>
503struct result_of<F(ArgTypes...)> {
504 typedef typename std::invoke_result<F, ArgTypes...>::type type1;
505 typedef typename remove_all<type1>::type type;
506};
507#elif EIGEN_HAS_STD_RESULT_OF
508template<typename T> struct result_of {
509 typedef typename std::result_of<T>::type type1;
510 typedef typename remove_all<type1>::type type;
511};
512#else
513template<typename T> struct result_of { };
514
515struct has_none {int a[1];};
516struct has_std_result_type {int a[2];};
517struct has_tr1_result {int a[3];};
518
519template<typename Func, int SizeOf>
520struct nullary_result_of_select {};
521
522template<typename Func>
523struct nullary_result_of_select<Func, sizeof(has_std_result_type)> {typedef typename Func::result_type type;};
524
525template<typename Func>
526struct nullary_result_of_select<Func, sizeof(has_tr1_result)> {typedef typename Func::template result<Func()>::type type;};
527
528template<typename Func>
529struct result_of<Func()> {
530 template<typename T>
531 static has_std_result_type testFunctor(T const *, typename T::result_type const * = 0);
532 template<typename T>
533 static has_tr1_result testFunctor(T const *, typename T::template result<T()>::type const * = 0);
534 static has_none testFunctor(...);
535
536 // note that the following indirection is needed for gcc-3.3
537 enum {FunctorType = sizeof(testFunctor(static_cast<Func*>(0)))};
538 typedef typename nullary_result_of_select<Func, FunctorType>::type type;
539};
540
541template<typename Func, typename ArgType, int SizeOf=sizeof(has_none)>
542struct unary_result_of_select {typedef typename internal::remove_all<ArgType>::type type;};
543
544template<typename Func, typename ArgType>
545struct unary_result_of_select<Func, ArgType, sizeof(has_std_result_type)> {typedef typename Func::result_type type;};
546
547template<typename Func, typename ArgType>
548struct unary_result_of_select<Func, ArgType, sizeof(has_tr1_result)> {typedef typename Func::template result<Func(ArgType)>::type type;};
549
550template<typename Func, typename ArgType>
551struct result_of<Func(ArgType)> {
552 template<typename T>
553 static has_std_result_type testFunctor(T const *, typename T::result_type const * = 0);
554 template<typename T>
555 static has_tr1_result testFunctor(T const *, typename T::template result<T(ArgType)>::type const * = 0);
556 static has_none testFunctor(...);
557
558 // note that the following indirection is needed for gcc-3.3
559 enum {FunctorType = sizeof(testFunctor(static_cast<Func*>(0)))};
560 typedef typename unary_result_of_select<Func, ArgType, FunctorType>::type type;
561};
562
563template<typename Func, typename ArgType0, typename ArgType1, int SizeOf=sizeof(has_none)>
564struct binary_result_of_select {typedef typename internal::remove_all<ArgType0>::type type;};
565
566template<typename Func, typename ArgType0, typename ArgType1>
567struct binary_result_of_select<Func, ArgType0, ArgType1, sizeof(has_std_result_type)>
568{typedef typename Func::result_type type;};
569
570template<typename Func, typename ArgType0, typename ArgType1>
571struct binary_result_of_select<Func, ArgType0, ArgType1, sizeof(has_tr1_result)>
572{typedef typename Func::template result<Func(ArgType0,ArgType1)>::type type;};
573
574template<typename Func, typename ArgType0, typename ArgType1>
575struct result_of<Func(ArgType0,ArgType1)> {
576 template<typename T>
577 static has_std_result_type testFunctor(T const *, typename T::result_type const * = 0);
578 template<typename T>
579 static has_tr1_result testFunctor(T const *, typename T::template result<T(ArgType0,ArgType1)>::type const * = 0);
580 static has_none testFunctor(...);
581
582 // note that the following indirection is needed for gcc-3.3
583 enum {FunctorType = sizeof(testFunctor(static_cast<Func*>(0)))};
584 typedef typename binary_result_of_select<Func, ArgType0, ArgType1, FunctorType>::type type;
585};
586
587template<typename Func, typename ArgType0, typename ArgType1, typename ArgType2, int SizeOf=sizeof(has_none)>
588struct ternary_result_of_select {typedef typename internal::remove_all<ArgType0>::type type;};
589
590template<typename Func, typename ArgType0, typename ArgType1, typename ArgType2>
591struct ternary_result_of_select<Func, ArgType0, ArgType1, ArgType2, sizeof(has_std_result_type)>
592{typedef typename Func::result_type type;};
593
594template<typename Func, typename ArgType0, typename ArgType1, typename ArgType2>
595struct ternary_result_of_select<Func, ArgType0, ArgType1, ArgType2, sizeof(has_tr1_result)>
596{typedef typename Func::template result<Func(ArgType0,ArgType1,ArgType2)>::type type;};
597
598template<typename Func, typename ArgType0, typename ArgType1, typename ArgType2>
599struct result_of<Func(ArgType0,ArgType1,ArgType2)> {
600 template<typename T>
601 static has_std_result_type testFunctor(T const *, typename T::result_type const * = 0);
602 template<typename T>
603 static has_tr1_result testFunctor(T const *, typename T::template result<T(ArgType0,ArgType1,ArgType2)>::type const * = 0);
604 static has_none testFunctor(...);
605
606 // note that the following indirection is needed for gcc-3.3
607 enum {FunctorType = sizeof(testFunctor(static_cast<Func*>(0)))};
608 typedef typename ternary_result_of_select<Func, ArgType0, ArgType1, ArgType2, FunctorType>::type type;
609};
610
611#endif
612
613#if EIGEN_HAS_STD_INVOKE_RESULT
614template<typename F, typename... ArgTypes>
615struct invoke_result {
616 typedef typename std::invoke_result<F, ArgTypes...>::type type1;
617 typedef typename remove_all<type1>::type type;
618};
619#elif EIGEN_HAS_CXX11
620template<typename F, typename... ArgTypes>
621struct invoke_result {
622 typedef typename result_of<F(ArgTypes...)>::type type1;
623 typedef typename remove_all<type1>::type type;
624};
625#else
626template<typename F, typename ArgType0 = void, typename ArgType1 = void, typename ArgType2 = void>
627struct invoke_result {
628 typedef typename result_of<F(ArgType0, ArgType1, ArgType2)>::type type1;
629 typedef typename remove_all<type1>::type type;
630};
631
632template<typename F>
633struct invoke_result<F, void, void, void> {
634 typedef typename result_of<F()>::type type1;
635 typedef typename remove_all<type1>::type type;
636};
637
638template<typename F, typename ArgType0>
639struct invoke_result<F, ArgType0, void, void> {
640 typedef typename result_of<F(ArgType0)>::type type1;
641 typedef typename remove_all<type1>::type type;
642};
643
644template<typename F, typename ArgType0, typename ArgType1>
645struct invoke_result<F, ArgType0, ArgType1, void> {
646 typedef typename result_of<F(ArgType0, ArgType1)>::type type1;
647 typedef typename remove_all<type1>::type type;
648};
649#endif
650
651struct meta_yes { char a[1]; };
652struct meta_no { char a[2]; };
653
654// Check whether T::ReturnType does exist
655template <typename T>
656struct has_ReturnType
657{
658 template <typename C> static meta_yes testFunctor(C const *, typename C::ReturnType const * = 0);
659 template <typename C> static meta_no testFunctor(...);
660
661 enum { value = sizeof(testFunctor<T>(static_cast<T*>(0))) == sizeof(meta_yes) };
662};
663
664template<typename T> const T* return_ptr();
665
666template <typename T, typename IndexType=Index>
667struct has_nullary_operator
668{
669 template <typename C> static meta_yes testFunctor(C const *,typename enable_if<(sizeof(return_ptr<C>()->operator()())>0)>::type * = 0);
670 static meta_no testFunctor(...);
671
672 enum { value = sizeof(testFunctor(static_cast<T*>(0))) == sizeof(meta_yes) };
673};
674
675template <typename T, typename IndexType=Index>
676struct has_unary_operator
677{
678 template <typename C> static meta_yes testFunctor(C const *,typename enable_if<(sizeof(return_ptr<C>()->operator()(IndexType(0)))>0)>::type * = 0);
679 static meta_no testFunctor(...);
680
681 enum { value = sizeof(testFunctor(static_cast<T*>(0))) == sizeof(meta_yes) };
682};
683
684template <typename T, typename IndexType=Index>
685struct has_binary_operator
686{
687 template <typename C> static meta_yes testFunctor(C const *,typename enable_if<(sizeof(return_ptr<C>()->operator()(IndexType(0),IndexType(0)))>0)>::type * = 0);
688 static meta_no testFunctor(...);
689
690 enum { value = sizeof(testFunctor(static_cast<T*>(0))) == sizeof(meta_yes) };
691};
692
696template<int Y,
697 int InfX = 0,
698 int SupX = ((Y==1) ? 1 : Y/2),
699 bool Done = ((SupX-InfX)<=1 ? true : ((SupX*SupX <= Y) && ((SupX+1)*(SupX+1) > Y))) >
700 // use ?: instead of || just to shut up a stupid gcc 4.3 warning
701class meta_sqrt
702{
703 enum {
704 MidX = (InfX+SupX)/2,
705 TakeInf = MidX*MidX > Y ? 1 : 0,
706 NewInf = int(TakeInf) ? InfX : int(MidX),
707 NewSup = int(TakeInf) ? int(MidX) : SupX
708 };
709 public:
710 enum { ret = meta_sqrt<Y,NewInf,NewSup>::ret };
711};
712
713template<int Y, int InfX, int SupX>
714class meta_sqrt<Y, InfX, SupX, true> { public: enum { ret = (SupX*SupX <= Y) ? SupX : InfX }; };
715
716
720template<int A, int B, int K=1, bool Done = ((A*K)%B)==0, bool Big=(A>=B)>
721struct meta_least_common_multiple
722{
723 enum { ret = meta_least_common_multiple<A,B,K+1>::ret };
724};
725template<int A, int B, int K, bool Done>
726struct meta_least_common_multiple<A,B,K,Done,false>
727{
728 enum { ret = meta_least_common_multiple<B,A,K>::ret };
729};
730template<int A, int B, int K>
731struct meta_least_common_multiple<A,B,K,true,true>
732{
733 enum { ret = A*K };
734};
735
736
738template<typename T, typename U> struct scalar_product_traits
739{
740 enum { Defined = 0 };
741};
742
743// FIXME quick workaround around current limitation of result_of
744// template<typename Scalar, typename ArgType0, typename ArgType1>
745// struct result_of<scalar_product_op<Scalar>(ArgType0,ArgType1)> {
746// typedef typename scalar_product_traits<typename remove_all<ArgType0>::type, typename remove_all<ArgType1>::type>::ReturnType type;
747// };
748
752template<unsigned Len, unsigned Align>
753struct aligned_storage {
754 struct type {
755 EIGEN_ALIGN_TO_BOUNDARY(Align) unsigned char data[Len];
756 };
757};
758
759} // end namespace internal
760
761namespace numext {
762
763#if defined(EIGEN_GPU_COMPILE_PHASE)
764template<typename T> EIGEN_DEVICE_FUNC void swap(T &a, T &b) { T tmp = b; b = a; a = tmp; }
765#else
766template<typename T> EIGEN_STRONG_INLINE void swap(T &a, T &b) { std::swap(a,b); }
767#endif
768
769#if defined(EIGEN_GPU_COMPILE_PHASE) && !EIGEN_HAS_CXX11
770using internal::device::numeric_limits;
771#else
772using std::numeric_limits;
773#endif
774
775// Integer division with rounding up.
776// T is assumed to be an integer type with a>=0, and b>0
777template<typename T>
778EIGEN_DEVICE_FUNC
779T div_ceil(const T &a, const T &b)
780{
781 return (a+b-1) / b;
782}
783
784// The aim of the following functions is to bypass -Wfloat-equal warnings
785// when we really want a strict equality comparison on floating points.
786template<typename X, typename Y> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC
787bool equal_strict(const X& x,const Y& y) { return x == y; }
788
789#if !defined(EIGEN_GPU_COMPILE_PHASE) || (!defined(EIGEN_CUDA_ARCH) && defined(EIGEN_CONSTEXPR_ARE_DEVICE_FUNC))
790template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC
791bool equal_strict(const float& x,const float& y) { return std::equal_to<float>()(x,y); }
792
793template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC
794bool equal_strict(const double& x,const double& y) { return std::equal_to<double>()(x,y); }
795#endif
796
797template<typename X, typename Y> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC
798bool not_equal_strict(const X& x,const Y& y) { return x != y; }
799
800#if !defined(EIGEN_GPU_COMPILE_PHASE) || (!defined(EIGEN_CUDA_ARCH) && defined(EIGEN_CONSTEXPR_ARE_DEVICE_FUNC))
801template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC
802bool not_equal_strict(const float& x,const float& y) { return std::not_equal_to<float>()(x,y); }
803
804template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC
805bool not_equal_strict(const double& x,const double& y) { return std::not_equal_to<double>()(x,y); }
806#endif
807
808} // end namespace numext
809
810} // end namespace Eigen
811
812#endif // EIGEN_META_H
Namespace containing all symbols from the Eigen library.
Definition: Core:141
EIGEN_DEFAULT_DENSE_INDEX_TYPE Index
The Index type as used for the API.
Definition: Meta.h:74
const int Dynamic
Definition: Constants.h:22