Macros.h revision 7faaa9f3f0df9d23790277834d426c3d992ac3ba
1// This file is part of Eigen, a lightweight C++ template library 2// for linear algebra. 3// 4// Copyright (C) 2008-2010 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_MACROS_H 12#define EIGEN_MACROS_H 13 14#define EIGEN_WORLD_VERSION 3 15#define EIGEN_MAJOR_VERSION 2 16#define EIGEN_MINOR_VERSION 2 17 18#define EIGEN_VERSION_AT_LEAST(x,y,z) (EIGEN_WORLD_VERSION>x || (EIGEN_WORLD_VERSION>=x && \ 19 (EIGEN_MAJOR_VERSION>y || (EIGEN_MAJOR_VERSION>=y && \ 20 EIGEN_MINOR_VERSION>=z)))) 21#ifdef __GNUC__ 22 #define EIGEN_GNUC_AT_LEAST(x,y) ((__GNUC__==x && __GNUC_MINOR__>=y) || __GNUC__>x) 23#else 24 #define EIGEN_GNUC_AT_LEAST(x,y) 0 25#endif 26 27#ifdef __GNUC__ 28 #define EIGEN_GNUC_AT_MOST(x,y) ((__GNUC__==x && __GNUC_MINOR__<=y) || __GNUC__<x) 29#else 30 #define EIGEN_GNUC_AT_MOST(x,y) 0 31#endif 32 33#if EIGEN_GNUC_AT_MOST(4,3) && !defined(__clang__) 34 // see bug 89 35 #define EIGEN_SAFE_TO_USE_STANDARD_ASSERT_MACRO 0 36#else 37 #define EIGEN_SAFE_TO_USE_STANDARD_ASSERT_MACRO 1 38#endif 39 40#if defined(__GNUC__) && (__GNUC__ <= 3) 41#define EIGEN_GCC3_OR_OLDER 1 42#else 43#define EIGEN_GCC3_OR_OLDER 0 44#endif 45 46// 16 byte alignment is only useful for vectorization. Since it affects the ABI, we need to enable 47// 16 byte alignment on all platforms where vectorization might be enabled. In theory we could always 48// enable alignment, but it can be a cause of problems on some platforms, so we just disable it in 49// certain common platform (compiler+architecture combinations) to avoid these problems. 50// Only static alignment is really problematic (relies on nonstandard compiler extensions that don't 51// work everywhere, for example don't work on GCC/ARM), try to keep heap alignment even 52// when we have to disable static alignment. 53#if defined(__GNUC__) && !(defined(__i386__) || defined(__x86_64__) || defined(__powerpc__) || defined(__ppc__) || defined(__ia64__)) 54#define EIGEN_GCC_AND_ARCH_DOESNT_WANT_STACK_ALIGNMENT 1 55#else 56#define EIGEN_GCC_AND_ARCH_DOESNT_WANT_STACK_ALIGNMENT 0 57#endif 58 59// static alignment is completely disabled with GCC 3, Sun Studio, and QCC/QNX 60#if !EIGEN_GCC_AND_ARCH_DOESNT_WANT_STACK_ALIGNMENT \ 61 && !EIGEN_GCC3_OR_OLDER \ 62 && !defined(__SUNPRO_CC) \ 63 && !defined(__QNXNTO__) 64 #define EIGEN_ARCH_WANTS_STACK_ALIGNMENT 1 65#else 66 #define EIGEN_ARCH_WANTS_STACK_ALIGNMENT 0 67#endif 68 69#ifdef EIGEN_DONT_ALIGN 70 #ifndef EIGEN_DONT_ALIGN_STATICALLY 71 #define EIGEN_DONT_ALIGN_STATICALLY 72 #endif 73 #define EIGEN_ALIGN 0 74#else 75 #define EIGEN_ALIGN 1 76#endif 77 78// EIGEN_ALIGN_STATICALLY is the true test whether we want to align arrays on the stack or not. It takes into account both the user choice to explicitly disable 79// alignment (EIGEN_DONT_ALIGN_STATICALLY) and the architecture config (EIGEN_ARCH_WANTS_STACK_ALIGNMENT). Henceforth, only EIGEN_ALIGN_STATICALLY should be used. 80#if EIGEN_ARCH_WANTS_STACK_ALIGNMENT && !defined(EIGEN_DONT_ALIGN_STATICALLY) 81 #define EIGEN_ALIGN_STATICALLY 1 82#else 83 #define EIGEN_ALIGN_STATICALLY 0 84 #ifndef EIGEN_DISABLE_UNALIGNED_ARRAY_ASSERT 85 #define EIGEN_DISABLE_UNALIGNED_ARRAY_ASSERT 86 #endif 87#endif 88 89#ifdef EIGEN_DEFAULT_TO_ROW_MAJOR 90#define EIGEN_DEFAULT_MATRIX_STORAGE_ORDER_OPTION RowMajor 91#else 92#define EIGEN_DEFAULT_MATRIX_STORAGE_ORDER_OPTION ColMajor 93#endif 94 95#ifndef EIGEN_DEFAULT_DENSE_INDEX_TYPE 96#define EIGEN_DEFAULT_DENSE_INDEX_TYPE std::ptrdiff_t 97#endif 98 99/** Allows to disable some optimizations which might affect the accuracy of the result. 100 * Such optimization are enabled by default, and set EIGEN_FAST_MATH to 0 to disable them. 101 * They currently include: 102 * - single precision Cwise::sin() and Cwise::cos() when SSE vectorization is enabled. 103 */ 104#ifndef EIGEN_FAST_MATH 105#define EIGEN_FAST_MATH 1 106#endif 107 108#define EIGEN_DEBUG_VAR(x) std::cerr << #x << " = " << x << std::endl; 109 110// concatenate two tokens 111#define EIGEN_CAT2(a,b) a ## b 112#define EIGEN_CAT(a,b) EIGEN_CAT2(a,b) 113 114// convert a token to a string 115#define EIGEN_MAKESTRING2(a) #a 116#define EIGEN_MAKESTRING(a) EIGEN_MAKESTRING2(a) 117 118// EIGEN_STRONG_INLINE is a stronger version of the inline, using __forceinline on MSVC, 119// but it still doesn't use GCC's always_inline. This is useful in (common) situations where MSVC needs forceinline 120// but GCC is still doing fine with just inline. 121#if (defined _MSC_VER) || (defined __INTEL_COMPILER) 122#define EIGEN_STRONG_INLINE __forceinline 123#else 124#define EIGEN_STRONG_INLINE inline 125#endif 126 127// EIGEN_ALWAYS_INLINE is the stronget, it has the effect of making the function inline and adding every possible 128// attribute to maximize inlining. This should only be used when really necessary: in particular, 129// it uses __attribute__((always_inline)) on GCC, which most of the time is useless and can severely harm compile times. 130// FIXME with the always_inline attribute, 131// gcc 3.4.x reports the following compilation error: 132// Eval.h:91: sorry, unimplemented: inlining failed in call to 'const Eigen::Eval<Derived> Eigen::MatrixBase<Scalar, Derived>::eval() const' 133// : function body not available 134#if EIGEN_GNUC_AT_LEAST(4,0) 135#define EIGEN_ALWAYS_INLINE __attribute__((always_inline)) inline 136#else 137#define EIGEN_ALWAYS_INLINE EIGEN_STRONG_INLINE 138#endif 139 140#if (defined __GNUC__) 141#define EIGEN_DONT_INLINE __attribute__((noinline)) 142#elif (defined _MSC_VER) 143#define EIGEN_DONT_INLINE __declspec(noinline) 144#else 145#define EIGEN_DONT_INLINE 146#endif 147 148#if (defined __GNUC__) 149#define EIGEN_PERMISSIVE_EXPR __extension__ 150#else 151#define EIGEN_PERMISSIVE_EXPR 152#endif 153 154// this macro allows to get rid of linking errors about multiply defined functions. 155// - static is not very good because it prevents definitions from different object files to be merged. 156// So static causes the resulting linked executable to be bloated with multiple copies of the same function. 157// - inline is not perfect either as it unwantedly hints the compiler toward inlining the function. 158#define EIGEN_DECLARE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS 159#define EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS inline 160 161#ifdef NDEBUG 162# ifndef EIGEN_NO_DEBUG 163# define EIGEN_NO_DEBUG 164# endif 165#endif 166 167// eigen_plain_assert is where we implement the workaround for the assert() bug in GCC <= 4.3, see bug 89 168#ifdef EIGEN_NO_DEBUG 169 #define eigen_plain_assert(x) 170#else 171 #if EIGEN_SAFE_TO_USE_STANDARD_ASSERT_MACRO 172 namespace Eigen { 173 namespace internal { 174 inline bool copy_bool(bool b) { return b; } 175 } 176 } 177 #define eigen_plain_assert(x) assert(x) 178 #else 179 // work around bug 89 180 #include <cstdlib> // for abort 181 #include <iostream> // for std::cerr 182 183 namespace Eigen { 184 namespace internal { 185 // trivial function copying a bool. Must be EIGEN_DONT_INLINE, so we implement it after including Eigen headers. 186 // see bug 89. 187 namespace { 188 EIGEN_DONT_INLINE bool copy_bool(bool b) { return b; } 189 } 190 inline void assert_fail(const char *condition, const char *function, const char *file, int line) 191 { 192 std::cerr << "assertion failed: " << condition << " in function " << function << " at " << file << ":" << line << std::endl; 193 abort(); 194 } 195 } 196 } 197 #define eigen_plain_assert(x) \ 198 do { \ 199 if(!Eigen::internal::copy_bool(x)) \ 200 Eigen::internal::assert_fail(EIGEN_MAKESTRING(x), __PRETTY_FUNCTION__, __FILE__, __LINE__); \ 201 } while(false) 202 #endif 203#endif 204 205// eigen_assert can be overridden 206#ifndef eigen_assert 207#define eigen_assert(x) eigen_plain_assert(x) 208#endif 209 210#ifdef EIGEN_INTERNAL_DEBUGGING 211#define eigen_internal_assert(x) eigen_assert(x) 212#else 213#define eigen_internal_assert(x) 214#endif 215 216#ifdef EIGEN_NO_DEBUG 217#define EIGEN_ONLY_USED_FOR_DEBUG(x) (void)x 218#else 219#define EIGEN_ONLY_USED_FOR_DEBUG(x) 220#endif 221 222#ifndef EIGEN_NO_DEPRECATED_WARNING 223 #if (defined __GNUC__) 224 #define EIGEN_DEPRECATED __attribute__((deprecated)) 225 #elif (defined _MSC_VER) 226 #define EIGEN_DEPRECATED __declspec(deprecated) 227 #else 228 #define EIGEN_DEPRECATED 229 #endif 230#else 231 #define EIGEN_DEPRECATED 232#endif 233 234#if (defined __GNUC__) 235#define EIGEN_UNUSED __attribute__((unused)) 236#else 237#define EIGEN_UNUSED 238#endif 239 240// Suppresses 'unused variable' warnings. 241namespace Eigen { 242 namespace internal { 243 template<typename T> void ignore_unused_variable(const T&) {} 244 } 245} 246#define EIGEN_UNUSED_VARIABLE(var) Eigen::internal::ignore_unused_variable(var); 247 248#if !defined(EIGEN_ASM_COMMENT) 249 #if (defined __GNUC__) && ( defined(__i386__) || defined(__x86_64__) ) 250 #define EIGEN_ASM_COMMENT(X) asm("#" X) 251 #else 252 #define EIGEN_ASM_COMMENT(X) 253 #endif 254#endif 255 256/* EIGEN_ALIGN_TO_BOUNDARY(n) forces data to be n-byte aligned. This is used to satisfy SIMD requirements. 257 * However, we do that EVEN if vectorization (EIGEN_VECTORIZE) is disabled, 258 * so that vectorization doesn't affect binary compatibility. 259 * 260 * If we made alignment depend on whether or not EIGEN_VECTORIZE is defined, it would be impossible to link 261 * vectorized and non-vectorized code. 262 */ 263#if (defined __GNUC__) || (defined __PGI) || (defined __IBMCPP__) || (defined __ARMCC_VERSION) 264 #define EIGEN_ALIGN_TO_BOUNDARY(n) __attribute__((aligned(n))) 265#elif (defined _MSC_VER) 266 #define EIGEN_ALIGN_TO_BOUNDARY(n) __declspec(align(n)) 267#elif (defined __SUNPRO_CC) 268 // FIXME not sure about this one: 269 #define EIGEN_ALIGN_TO_BOUNDARY(n) __attribute__((aligned(n))) 270#else 271 #error Please tell me what is the equivalent of __attribute__((aligned(n))) for your compiler 272#endif 273 274#define EIGEN_ALIGN16 EIGEN_ALIGN_TO_BOUNDARY(16) 275 276#if EIGEN_ALIGN_STATICALLY 277#define EIGEN_USER_ALIGN_TO_BOUNDARY(n) EIGEN_ALIGN_TO_BOUNDARY(n) 278#define EIGEN_USER_ALIGN16 EIGEN_ALIGN16 279#else 280#define EIGEN_USER_ALIGN_TO_BOUNDARY(n) 281#define EIGEN_USER_ALIGN16 282#endif 283 284#ifdef EIGEN_DONT_USE_RESTRICT_KEYWORD 285 #define EIGEN_RESTRICT 286#endif 287#ifndef EIGEN_RESTRICT 288 #define EIGEN_RESTRICT __restrict 289#endif 290 291#ifndef EIGEN_STACK_ALLOCATION_LIMIT 292// 131072 == 128 KB 293#define EIGEN_STACK_ALLOCATION_LIMIT 131072 294#endif 295 296#ifndef EIGEN_DEFAULT_IO_FORMAT 297#ifdef EIGEN_MAKING_DOCS 298// format used in Eigen's documentation 299// needed to define it here as escaping characters in CMake add_definition's argument seems very problematic. 300#define EIGEN_DEFAULT_IO_FORMAT Eigen::IOFormat(3, 0, " ", "\n", "", "") 301#else 302#define EIGEN_DEFAULT_IO_FORMAT Eigen::IOFormat() 303#endif 304#endif 305 306// just an empty macro ! 307#define EIGEN_EMPTY 308 309#if defined(_MSC_VER) && (!defined(__INTEL_COMPILER)) 310#define EIGEN_INHERIT_ASSIGNMENT_EQUAL_OPERATOR(Derived) \ 311 using Base::operator =; 312#elif defined(__clang__) // workaround clang bug (see http://forum.kde.org/viewtopic.php?f=74&t=102653) 313#define EIGEN_INHERIT_ASSIGNMENT_EQUAL_OPERATOR(Derived) \ 314 using Base::operator =; \ 315 EIGEN_STRONG_INLINE Derived& operator=(const Derived& other) { Base::operator=(other); return *this; } \ 316 template <typename OtherDerived> \ 317 EIGEN_STRONG_INLINE Derived& operator=(const DenseBase<OtherDerived>& other) { Base::operator=(other.derived()); return *this; } 318#else 319#define EIGEN_INHERIT_ASSIGNMENT_EQUAL_OPERATOR(Derived) \ 320 using Base::operator =; \ 321 EIGEN_STRONG_INLINE Derived& operator=(const Derived& other) \ 322 { \ 323 Base::operator=(other); \ 324 return *this; \ 325 } 326#endif 327 328#define EIGEN_INHERIT_ASSIGNMENT_OPERATORS(Derived) \ 329 EIGEN_INHERIT_ASSIGNMENT_EQUAL_OPERATOR(Derived) 330 331/** 332* Just a side note. Commenting within defines works only by documenting 333* behind the object (via '!<'). Comments cannot be multi-line and thus 334* we have these extra long lines. What is confusing doxygen over here is 335* that we use '\' and basically have a bunch of typedefs with their 336* documentation in a single line. 337**/ 338 339#define EIGEN_GENERIC_PUBLIC_INTERFACE(Derived) \ 340 typedef typename Eigen::internal::traits<Derived>::Scalar Scalar; /*!< \brief Numeric type, e.g. float, double, int or std::complex<float>. */ \ 341 typedef typename Eigen::NumTraits<Scalar>::Real RealScalar; /*!< \brief The underlying numeric type for composed scalar types. \details In cases where Scalar is e.g. std::complex<T>, T were corresponding to RealScalar. */ \ 342 typedef typename Base::CoeffReturnType CoeffReturnType; /*!< \brief The return type for coefficient access. \details Depending on whether the object allows direct coefficient access (e.g. for a MatrixXd), this type is either 'const Scalar&' or simply 'Scalar' for objects that do not allow direct coefficient access. */ \ 343 typedef typename Eigen::internal::nested<Derived>::type Nested; \ 344 typedef typename Eigen::internal::traits<Derived>::StorageKind StorageKind; \ 345 typedef typename Eigen::internal::traits<Derived>::Index Index; \ 346 enum { RowsAtCompileTime = Eigen::internal::traits<Derived>::RowsAtCompileTime, \ 347 ColsAtCompileTime = Eigen::internal::traits<Derived>::ColsAtCompileTime, \ 348 Flags = Eigen::internal::traits<Derived>::Flags, \ 349 CoeffReadCost = Eigen::internal::traits<Derived>::CoeffReadCost, \ 350 SizeAtCompileTime = Base::SizeAtCompileTime, \ 351 MaxSizeAtCompileTime = Base::MaxSizeAtCompileTime, \ 352 IsVectorAtCompileTime = Base::IsVectorAtCompileTime }; 353 354 355#define EIGEN_DENSE_PUBLIC_INTERFACE(Derived) \ 356 typedef typename Eigen::internal::traits<Derived>::Scalar Scalar; /*!< \brief Numeric type, e.g. float, double, int or std::complex<float>. */ \ 357 typedef typename Eigen::NumTraits<Scalar>::Real RealScalar; /*!< \brief The underlying numeric type for composed scalar types. \details In cases where Scalar is e.g. std::complex<T>, T were corresponding to RealScalar. */ \ 358 typedef typename Base::PacketScalar PacketScalar; \ 359 typedef typename Base::CoeffReturnType CoeffReturnType; /*!< \brief The return type for coefficient access. \details Depending on whether the object allows direct coefficient access (e.g. for a MatrixXd), this type is either 'const Scalar&' or simply 'Scalar' for objects that do not allow direct coefficient access. */ \ 360 typedef typename Eigen::internal::nested<Derived>::type Nested; \ 361 typedef typename Eigen::internal::traits<Derived>::StorageKind StorageKind; \ 362 typedef typename Eigen::internal::traits<Derived>::Index Index; \ 363 enum { RowsAtCompileTime = Eigen::internal::traits<Derived>::RowsAtCompileTime, \ 364 ColsAtCompileTime = Eigen::internal::traits<Derived>::ColsAtCompileTime, \ 365 MaxRowsAtCompileTime = Eigen::internal::traits<Derived>::MaxRowsAtCompileTime, \ 366 MaxColsAtCompileTime = Eigen::internal::traits<Derived>::MaxColsAtCompileTime, \ 367 Flags = Eigen::internal::traits<Derived>::Flags, \ 368 CoeffReadCost = Eigen::internal::traits<Derived>::CoeffReadCost, \ 369 SizeAtCompileTime = Base::SizeAtCompileTime, \ 370 MaxSizeAtCompileTime = Base::MaxSizeAtCompileTime, \ 371 IsVectorAtCompileTime = Base::IsVectorAtCompileTime }; \ 372 using Base::derived; \ 373 using Base::const_cast_derived; 374 375 376#define EIGEN_PLAIN_ENUM_MIN(a,b) (((int)a <= (int)b) ? (int)a : (int)b) 377#define EIGEN_PLAIN_ENUM_MAX(a,b) (((int)a >= (int)b) ? (int)a : (int)b) 378 379// EIGEN_SIZE_MIN_PREFER_DYNAMIC gives the min between compile-time sizes. 0 has absolute priority, followed by 1, 380// followed by Dynamic, followed by other finite values. The reason for giving Dynamic the priority over 381// finite values is that min(3, Dynamic) should be Dynamic, since that could be anything between 0 and 3. 382#define EIGEN_SIZE_MIN_PREFER_DYNAMIC(a,b) (((int)a == 0 || (int)b == 0) ? 0 \ 383 : ((int)a == 1 || (int)b == 1) ? 1 \ 384 : ((int)a == Dynamic || (int)b == Dynamic) ? Dynamic \ 385 : ((int)a <= (int)b) ? (int)a : (int)b) 386 387// EIGEN_SIZE_MIN_PREFER_FIXED is a variant of EIGEN_SIZE_MIN_PREFER_DYNAMIC comparing MaxSizes. The difference is that finite values 388// now have priority over Dynamic, so that min(3, Dynamic) gives 3. Indeed, whatever the actual value is 389// (between 0 and 3), it is not more than 3. 390#define EIGEN_SIZE_MIN_PREFER_FIXED(a,b) (((int)a == 0 || (int)b == 0) ? 0 \ 391 : ((int)a == 1 || (int)b == 1) ? 1 \ 392 : ((int)a == Dynamic && (int)b == Dynamic) ? Dynamic \ 393 : ((int)a == Dynamic) ? (int)b \ 394 : ((int)b == Dynamic) ? (int)a \ 395 : ((int)a <= (int)b) ? (int)a : (int)b) 396 397// see EIGEN_SIZE_MIN_PREFER_DYNAMIC. No need for a separate variant for MaxSizes here. 398#define EIGEN_SIZE_MAX(a,b) (((int)a == Dynamic || (int)b == Dynamic) ? Dynamic \ 399 : ((int)a >= (int)b) ? (int)a : (int)b) 400 401#define EIGEN_LOGICAL_XOR(a,b) (((a) || (b)) && !((a) && (b))) 402 403#define EIGEN_IMPLIES(a,b) (!(a) || (b)) 404 405#define EIGEN_MAKE_CWISE_BINARY_OP(METHOD,FUNCTOR) \ 406 template<typename OtherDerived> \ 407 EIGEN_STRONG_INLINE const CwiseBinaryOp<FUNCTOR<Scalar>, const Derived, const OtherDerived> \ 408 (METHOD)(const EIGEN_CURRENT_STORAGE_BASE_CLASS<OtherDerived> &other) const \ 409 { \ 410 return CwiseBinaryOp<FUNCTOR<Scalar>, const Derived, const OtherDerived>(derived(), other.derived()); \ 411 } 412 413// the expression type of a cwise product 414#define EIGEN_CWISE_PRODUCT_RETURN_TYPE(LHS,RHS) \ 415 CwiseBinaryOp< \ 416 internal::scalar_product_op< \ 417 typename internal::traits<LHS>::Scalar, \ 418 typename internal::traits<RHS>::Scalar \ 419 >, \ 420 const LHS, \ 421 const RHS \ 422 > 423 424#endif // EIGEN_MACROS_H 425