2020-12-19 13:50:20 +00:00
|
|
|
// Copyright 2009-2020 Intel Corporation
|
|
|
|
// SPDX-License-Identifier: Apache-2.0
|
|
|
|
|
|
|
|
#pragma once
|
|
|
|
|
|
|
|
#include "../sys/platform.h"
|
|
|
|
|
|
|
|
#include <limits>
|
|
|
|
|
|
|
|
#define _USE_MATH_DEFINES
|
|
|
|
#include <math.h> // using cmath causes issues under Windows
|
|
|
|
#include <cfloat>
|
|
|
|
#include <climits>
|
|
|
|
|
2021-05-04 09:07:12 +00:00
|
|
|
// Math constants may not be defined in libcxx + mingw + strict C++ standard
|
|
|
|
#if defined(__MINGW32__)
|
|
|
|
|
|
|
|
// TODO(LTE): use constexpr
|
|
|
|
#ifndef M_PI
|
|
|
|
#define M_PI 3.14159265358979323846
|
|
|
|
#endif
|
|
|
|
#ifndef M_1_PI
|
|
|
|
#define M_1_PI 0.31830988618379067154
|
|
|
|
#endif
|
|
|
|
|
|
|
|
#endif // __MINGW32__
|
|
|
|
|
2020-12-19 13:50:20 +00:00
|
|
|
namespace embree
|
|
|
|
{
|
|
|
|
static MAYBE_UNUSED const float one_over_255 = 1.0f/255.0f;
|
|
|
|
static MAYBE_UNUSED const float min_rcp_input = 1E-18f; // for abs(x) >= min_rcp_input the newton raphson rcp calculation does not fail
|
|
|
|
|
|
|
|
/* we consider floating point numbers in that range as valid input numbers */
|
|
|
|
static MAYBE_UNUSED float FLT_LARGE = 1.844E18f;
|
|
|
|
|
|
|
|
struct TrueTy {
|
|
|
|
__forceinline operator bool( ) const { return true; }
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED TrueTy True;
|
|
|
|
|
|
|
|
struct FalseTy {
|
|
|
|
__forceinline operator bool( ) const { return false; }
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED FalseTy False;
|
|
|
|
|
|
|
|
struct ZeroTy
|
|
|
|
{
|
|
|
|
__forceinline operator double ( ) const { return 0; }
|
|
|
|
__forceinline operator float ( ) const { return 0; }
|
|
|
|
__forceinline operator long long( ) const { return 0; }
|
|
|
|
__forceinline operator unsigned long long( ) const { return 0; }
|
|
|
|
__forceinline operator long ( ) const { return 0; }
|
|
|
|
__forceinline operator unsigned long ( ) const { return 0; }
|
|
|
|
__forceinline operator int ( ) const { return 0; }
|
|
|
|
__forceinline operator unsigned int ( ) const { return 0; }
|
|
|
|
__forceinline operator short ( ) const { return 0; }
|
|
|
|
__forceinline operator unsigned short ( ) const { return 0; }
|
2021-05-04 09:07:12 +00:00
|
|
|
__forceinline operator int8_t ( ) const { return 0; }
|
|
|
|
__forceinline operator uint8_t ( ) const { return 0; }
|
2020-12-19 13:50:20 +00:00
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED ZeroTy zero;
|
|
|
|
|
|
|
|
struct OneTy
|
|
|
|
{
|
|
|
|
__forceinline operator double ( ) const { return 1; }
|
|
|
|
__forceinline operator float ( ) const { return 1; }
|
|
|
|
__forceinline operator long long( ) const { return 1; }
|
|
|
|
__forceinline operator unsigned long long( ) const { return 1; }
|
|
|
|
__forceinline operator long ( ) const { return 1; }
|
|
|
|
__forceinline operator unsigned long ( ) const { return 1; }
|
|
|
|
__forceinline operator int ( ) const { return 1; }
|
|
|
|
__forceinline operator unsigned int ( ) const { return 1; }
|
|
|
|
__forceinline operator short ( ) const { return 1; }
|
|
|
|
__forceinline operator unsigned short ( ) const { return 1; }
|
2021-05-04 09:07:12 +00:00
|
|
|
__forceinline operator int8_t ( ) const { return 1; }
|
|
|
|
__forceinline operator uint8_t ( ) const { return 1; }
|
2020-12-19 13:50:20 +00:00
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED OneTy one;
|
|
|
|
|
|
|
|
struct NegInfTy
|
|
|
|
{
|
|
|
|
__forceinline operator double ( ) const { return -std::numeric_limits<double>::infinity(); }
|
|
|
|
__forceinline operator float ( ) const { return -std::numeric_limits<float>::infinity(); }
|
|
|
|
__forceinline operator long long( ) const { return std::numeric_limits<long long>::min(); }
|
|
|
|
__forceinline operator unsigned long long( ) const { return std::numeric_limits<unsigned long long>::min(); }
|
|
|
|
__forceinline operator long ( ) const { return std::numeric_limits<long>::min(); }
|
|
|
|
__forceinline operator unsigned long ( ) const { return std::numeric_limits<unsigned long>::min(); }
|
|
|
|
__forceinline operator int ( ) const { return std::numeric_limits<int>::min(); }
|
|
|
|
__forceinline operator unsigned int ( ) const { return std::numeric_limits<unsigned int>::min(); }
|
|
|
|
__forceinline operator short ( ) const { return std::numeric_limits<short>::min(); }
|
|
|
|
__forceinline operator unsigned short ( ) const { return std::numeric_limits<unsigned short>::min(); }
|
2021-05-04 09:07:12 +00:00
|
|
|
__forceinline operator int8_t ( ) const { return std::numeric_limits<int8_t>::min(); }
|
|
|
|
__forceinline operator uint8_t ( ) const { return std::numeric_limits<uint8_t>::min(); }
|
2020-12-19 13:50:20 +00:00
|
|
|
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED NegInfTy neg_inf;
|
|
|
|
|
|
|
|
struct PosInfTy
|
|
|
|
{
|
|
|
|
__forceinline operator double ( ) const { return std::numeric_limits<double>::infinity(); }
|
|
|
|
__forceinline operator float ( ) const { return std::numeric_limits<float>::infinity(); }
|
|
|
|
__forceinline operator long long( ) const { return std::numeric_limits<long long>::max(); }
|
|
|
|
__forceinline operator unsigned long long( ) const { return std::numeric_limits<unsigned long long>::max(); }
|
|
|
|
__forceinline operator long ( ) const { return std::numeric_limits<long>::max(); }
|
|
|
|
__forceinline operator unsigned long ( ) const { return std::numeric_limits<unsigned long>::max(); }
|
|
|
|
__forceinline operator int ( ) const { return std::numeric_limits<int>::max(); }
|
|
|
|
__forceinline operator unsigned int ( ) const { return std::numeric_limits<unsigned int>::max(); }
|
|
|
|
__forceinline operator short ( ) const { return std::numeric_limits<short>::max(); }
|
|
|
|
__forceinline operator unsigned short ( ) const { return std::numeric_limits<unsigned short>::max(); }
|
2021-05-04 09:07:12 +00:00
|
|
|
__forceinline operator int8_t ( ) const { return std::numeric_limits<int8_t>::max(); }
|
|
|
|
__forceinline operator uint8_t ( ) const { return std::numeric_limits<uint8_t>::max(); }
|
2020-12-19 13:50:20 +00:00
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED PosInfTy inf;
|
|
|
|
extern MAYBE_UNUSED PosInfTy pos_inf;
|
|
|
|
|
|
|
|
struct NaNTy
|
|
|
|
{
|
|
|
|
__forceinline operator double( ) const { return std::numeric_limits<double>::quiet_NaN(); }
|
|
|
|
__forceinline operator float ( ) const { return std::numeric_limits<float>::quiet_NaN(); }
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED NaNTy nan;
|
|
|
|
|
|
|
|
struct UlpTy
|
|
|
|
{
|
|
|
|
__forceinline operator double( ) const { return std::numeric_limits<double>::epsilon(); }
|
|
|
|
__forceinline operator float ( ) const { return std::numeric_limits<float>::epsilon(); }
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED UlpTy ulp;
|
|
|
|
|
|
|
|
struct PiTy
|
|
|
|
{
|
|
|
|
__forceinline operator double( ) const { return double(M_PI); }
|
|
|
|
__forceinline operator float ( ) const { return float(M_PI); }
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED PiTy pi;
|
|
|
|
|
|
|
|
struct OneOverPiTy
|
|
|
|
{
|
|
|
|
__forceinline operator double( ) const { return double(M_1_PI); }
|
|
|
|
__forceinline operator float ( ) const { return float(M_1_PI); }
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED OneOverPiTy one_over_pi;
|
|
|
|
|
|
|
|
struct TwoPiTy
|
|
|
|
{
|
|
|
|
__forceinline operator double( ) const { return double(2.0*M_PI); }
|
|
|
|
__forceinline operator float ( ) const { return float(2.0*M_PI); }
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED TwoPiTy two_pi;
|
|
|
|
|
|
|
|
struct OneOverTwoPiTy
|
|
|
|
{
|
|
|
|
__forceinline operator double( ) const { return double(0.5*M_1_PI); }
|
|
|
|
__forceinline operator float ( ) const { return float(0.5*M_1_PI); }
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED OneOverTwoPiTy one_over_two_pi;
|
|
|
|
|
|
|
|
struct FourPiTy
|
|
|
|
{
|
|
|
|
__forceinline operator double( ) const { return double(4.0*M_PI); }
|
|
|
|
__forceinline operator float ( ) const { return float(4.0*M_PI); }
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED FourPiTy four_pi;
|
|
|
|
|
|
|
|
struct OneOverFourPiTy
|
|
|
|
{
|
|
|
|
__forceinline operator double( ) const { return double(0.25*M_1_PI); }
|
|
|
|
__forceinline operator float ( ) const { return float(0.25*M_1_PI); }
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED OneOverFourPiTy one_over_four_pi;
|
|
|
|
|
|
|
|
struct StepTy {
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED StepTy step;
|
|
|
|
|
|
|
|
struct ReverseStepTy {
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED ReverseStepTy reverse_step;
|
|
|
|
|
|
|
|
struct EmptyTy {
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED EmptyTy empty;
|
|
|
|
|
|
|
|
struct FullTy {
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED FullTy full;
|
|
|
|
|
|
|
|
struct UndefinedTy {
|
|
|
|
};
|
|
|
|
|
|
|
|
extern MAYBE_UNUSED UndefinedTy undefined;
|
2021-05-04 09:07:12 +00:00
|
|
|
|
|
|
|
#if defined(__aarch64__)
|
|
|
|
extern const uint32x4_t movemask_mask;
|
|
|
|
extern const uint32x4_t vzero;
|
|
|
|
extern const uint32x4_t v0x80000000;
|
|
|
|
extern const uint32x4_t v0x7fffffff;
|
|
|
|
extern const uint32x4_t v000F;
|
|
|
|
extern const uint32x4_t v00F0;
|
|
|
|
extern const uint32x4_t v00FF;
|
|
|
|
extern const uint32x4_t v0F00;
|
|
|
|
extern const uint32x4_t v0F0F;
|
|
|
|
extern const uint32x4_t v0FF0;
|
|
|
|
extern const uint32x4_t v0FFF;
|
|
|
|
extern const uint32x4_t vF000;
|
|
|
|
extern const uint32x4_t vF00F;
|
|
|
|
extern const uint32x4_t vF0F0;
|
|
|
|
extern const uint32x4_t vF0FF;
|
|
|
|
extern const uint32x4_t vFF00;
|
|
|
|
extern const uint32x4_t vFF0F;
|
|
|
|
extern const uint32x4_t vFFF0;
|
|
|
|
extern const uint32x4_t vFFFF;
|
|
|
|
extern const uint8x16_t v0022;
|
|
|
|
extern const uint8x16_t v1133;
|
|
|
|
extern const uint8x16_t v0101;
|
|
|
|
extern const float32x4_t vOne;
|
|
|
|
extern const float32x4_t vmOne;
|
|
|
|
extern const float32x4_t vInf;
|
|
|
|
extern const float32x4_t vmInf;
|
|
|
|
#endif
|
2020-12-19 13:50:20 +00:00
|
|
|
}
|