extern/bullet/LinearMath/SimdScalar.h

   1 /*
   2 Copyright (c) 2003-2006 Gino van den Bergen / Erwin Coumans  http://continuousphysics.com/Bullet/
   3
   4 This software is provided 'as-is', without any express or implied warranty.
   5 In no event will the authors be held liable for any damages arising from the use of this software.
   6 Permission is granted to anyone to use this software for any purpose,
   7 including commercial applications, and to alter it and redistribute it freely,
   8 subject to the following restrictions:
   9
  10 1. The origin of this software must not be misrepresented; you must not claim that you wrote the original software. If you use this software in a product, an acknowledgment in the product documentation would be appreciated but is not required.
  11 2. Altered source versions must be plainly marked as such, and must not be misrepresented as being the original software.
  12 3. This notice may not be removed or altered from any source distribution.
  13 */
  14
  15
  16
  17 #ifndef SIMD___SCALAR_H
  18 #define SIMD___SCALAR_H
  19
  20 #include <math.h>
  21 #undef max
  22
  23
  24
  25 #include <cstdlib>
  26 #include <cfloat>
  27 #include <float.h>
  28
  29 #ifdef WIN32
  30 #pragma warning(disable:4530)
  31 #pragma warning(disable:4996)
  32 #ifdef __MINGW32__
  33 #define SIMD_FORCE_INLINE inline
  34 #else
  35 #define SIMD_FORCE_INLINE __forceinline
  36 #endif //__MINGW32__
  37
  38 //#define ATTRIBUTE_ALIGNED16(a) __declspec(align(16)) a
  39 #define ATTRIBUTE_ALIGNED16(a) a
  40
  41 #include <assert.h>
  42 #define ASSERT assert
  43 #else
  44 #define SIMD_FORCE_INLINE inline
  45 #define ATTRIBUTE_ALIGNED16(a) a
  46 #ifndef assert
  47 #include <assert.h>
  48 #endif
  49
  50
  51 #define ASSERT assert
  52
  53 #endif
  54
  55
  56
  57 typedef float    SimdScalar;
  58
  59 #if defined (__sun) || defined (__sun__) || defined (__sparc) || defined (__APPLE__)
  60 //use double float precision operation on those platforms for Blender
  61
  62 SIMD_FORCE_INLINE SimdScalar SimdSqrt(SimdScalar x) { return sqrt(x); }
  63 SIMD_FORCE_INLINE SimdScalar SimdFabs(SimdScalar x) { return fabs(x); }
  64 SIMD_FORCE_INLINE SimdScalar SimdCos(SimdScalar x) { return cos(x); }
  65 SIMD_FORCE_INLINE SimdScalar SimdSin(SimdScalar x) { return sin(x); }
  66 SIMD_FORCE_INLINE SimdScalar SimdTan(SimdScalar x) { return tan(x); }
  67 SIMD_FORCE_INLINE SimdScalar SimdAcos(SimdScalar x) { return acos(x); }
  68 SIMD_FORCE_INLINE SimdScalar SimdAsin(SimdScalar x) { return asin(x); }
  69 SIMD_FORCE_INLINE SimdScalar SimdAtan(SimdScalar x) { return atan(x); }
  70 SIMD_FORCE_INLINE SimdScalar SimdAtan2(SimdScalar x, SimdScalar y) { return atan2(x, y); }
  71 SIMD_FORCE_INLINE SimdScalar SimdExp(SimdScalar x) { return exp(x); }
  72 SIMD_FORCE_INLINE SimdScalar SimdLog(SimdScalar x) { return log(x); }
  73 SIMD_FORCE_INLINE SimdScalar SimdPow(SimdScalar x,SimdScalar y) { return pow(x,y); }
  74
  75 #else
  76
  77 SIMD_FORCE_INLINE SimdScalar SimdSqrt(SimdScalar x) { return sqrtf(x); }
  78 SIMD_FORCE_INLINE SimdScalar SimdFabs(SimdScalar x) { return fabsf(x); }
  79 SIMD_FORCE_INLINE SimdScalar SimdCos(SimdScalar x) { return cosf(x); }
  80 SIMD_FORCE_INLINE SimdScalar SimdSin(SimdScalar x) { return sinf(x); }
  81 SIMD_FORCE_INLINE SimdScalar SimdTan(SimdScalar x) { return tanf(x); }
  82 SIMD_FORCE_INLINE SimdScalar SimdAcos(SimdScalar x) { return acosf(x); }
  83 SIMD_FORCE_INLINE SimdScalar SimdAsin(SimdScalar x) { return asinf(x); }
  84 SIMD_FORCE_INLINE SimdScalar SimdAtan(SimdScalar x) { return atanf(x); }
  85 SIMD_FORCE_INLINE SimdScalar SimdAtan2(SimdScalar x, SimdScalar y) { return atan2f(x, y); }
  86 SIMD_FORCE_INLINE SimdScalar SimdExp(SimdScalar x) { return expf(x); }
  87 SIMD_FORCE_INLINE SimdScalar SimdLog(SimdScalar x) { return logf(x); }
  88 SIMD_FORCE_INLINE SimdScalar SimdPow(SimdScalar x,SimdScalar y) { return powf(x,y); }
  89
  90 #endif
  91
  92
  93 const SimdScalar  SIMD_2_PI         = 6.283185307179586232f;
  94 const SimdScalar  SIMD_PI           = SIMD_2_PI * SimdScalar(0.5f);
  95 const SimdScalar  SIMD_HALF_PI           = SIMD_2_PI * SimdScalar(0.25f);
  96 const SimdScalar  SIMD_RADS_PER_DEG = SIMD_2_PI / SimdScalar(360.0f);
  97 const SimdScalar  SIMD_DEGS_PER_RAD = SimdScalar(360.0f) / SIMD_2_PI;
  98 const SimdScalar  SIMD_EPSILON      = FLT_EPSILON;
  99 const SimdScalar  SIMD_INFINITY     = FLT_MAX;
 100
 101 SIMD_FORCE_INLINE bool      SimdFuzzyZero(SimdScalar x) { return SimdFabs(x) < SIMD_EPSILON; }
 102
 103 SIMD_FORCE_INLINE bool  SimdEqual(SimdScalar a, SimdScalar eps) {
 104         return (((a) <= eps) && !((a) < -eps));
 105 }
 106 SIMD_FORCE_INLINE bool  SimdGreaterEqual (SimdScalar a, SimdScalar eps) {
 107         return (!((a) <= eps));
 108 }
 109
 110 /*SIMD_FORCE_INLINE SimdScalar SimdCos(SimdScalar x) { return cosf(x); }
 111 SIMD_FORCE_INLINE SimdScalar SimdSin(SimdScalar x) { return sinf(x); }
 112 SIMD_FORCE_INLINE SimdScalar SimdTan(SimdScalar x) { return tanf(x); }
 113 SIMD_FORCE_INLINE SimdScalar SimdAcos(SimdScalar x) { return acosf(x); }
 114 SIMD_FORCE_INLINE SimdScalar SimdAsin(SimdScalar x) { return asinf(x); }
 115 SIMD_FORCE_INLINE SimdScalar SimdAtan(SimdScalar x) { return atanf(x); }
 116 SIMD_FORCE_INLINE SimdScalar SimdAtan2(SimdScalar x, SimdScalar y) { return atan2f(x, y); }
 117 */
 118
 119 SIMD_FORCE_INLINE int       SimdSign(SimdScalar x) {
 120     return x < 0.0f ? -1 : x > 0.0f ? 1 : 0;
 121 }
 122
 123 SIMD_FORCE_INLINE SimdScalar SimdRadians(SimdScalar x) { return x * SIMD_RADS_PER_DEG; }
 124 SIMD_FORCE_INLINE SimdScalar SimdDegrees(SimdScalar x) { return x * SIMD_DEGS_PER_RAD; }
 125
 126
 127
 128 #endif //SIMD___SCALAR_H