vec3ia_sycl.h 9.2 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178
  1. // Copyright 2009-2021 Intel Corporation
  2. // SPDX-License-Identifier: Apache-2.0
  3. #pragma once
  4. #include "../sys/alloc.h"
  5. #include "emath.h"
  6. #include "../simd/sse.h"
  7. namespace embree
  8. {
  9. ////////////////////////////////////////////////////////////////////////////////
  10. /// SSE Vec3ia Type
  11. ////////////////////////////////////////////////////////////////////////////////
  12. struct __aligned(16) Vec3ia
  13. {
  14. ALIGNED_STRUCT_(16);
  15. struct { int x,y,z; };
  16. typedef int Scalar;
  17. enum { N = 3 };
  18. ////////////////////////////////////////////////////////////////////////////////
  19. /// Constructors, Assignment & Cast Operators
  20. ////////////////////////////////////////////////////////////////////////////////
  21. __forceinline Vec3ia( ) {}
  22. //__forceinline Vec3ia( const __m128i a ) : m128(a) {}
  23. __forceinline Vec3ia( const Vec3ia& other ) : x(other.x), y(other.y), z(other.z) {}
  24. __forceinline Vec3ia& operator =(const Vec3ia& other) { x = other.x; y = other.y; z = other.z; return *this; }
  25. __forceinline explicit Vec3ia( const int a ) : x(a), y(a), z(a) {}
  26. __forceinline Vec3ia( const int x, const int y, const int z) : x(x), y(y), z(z) {}
  27. //__forceinline explicit Vec3ia( const __m128 a ) : m128(_mm_cvtps_epi32(a)) {}
  28. __forceinline explicit Vec3ia(const vint4& a) : x(a[0]), y(a[1]), z(a[2]) {}
  29. __forceinline explicit Vec3ia( const Vec3fa& a );
  30. //__forceinline operator const __m128i&() const { return m128; }
  31. //__forceinline operator __m128i&() { return m128; }
  32. __forceinline operator vint4() const { return vint4(x,y,z,z); }
  33. ////////////////////////////////////////////////////////////////////////////////
  34. /// Constants
  35. ////////////////////////////////////////////////////////////////////////////////
  36. __forceinline Vec3ia( ZeroTy ) : x(0), y(0), z(0) {}
  37. __forceinline Vec3ia( OneTy ) : x(1), y(1), z(1) {}
  38. __forceinline Vec3ia( PosInfTy ) : x(0x7FFFFFFF), y(0x7FFFFFFF), z(0x7FFFFFFF) {}
  39. __forceinline Vec3ia( NegInfTy ) : x(0x80000000), y(0x80000000), z(0x80000000) {}
  40. ////////////////////////////////////////////////////////////////////////////////
  41. /// Array Access
  42. ////////////////////////////////////////////////////////////////////////////////
  43. __forceinline const int& operator []( const size_t index ) const { assert(index < 3); return (&x)[index]; }
  44. __forceinline int& operator []( const size_t index ) { assert(index < 3); return (&x)[index]; }
  45. };
  46. ////////////////////////////////////////////////////////////////////////////////
  47. /// Unary Operators
  48. ////////////////////////////////////////////////////////////////////////////////
  49. __forceinline Vec3ia operator +( const Vec3ia& a ) { return Vec3ia(+a.x,+a.y,+a.z); }
  50. __forceinline Vec3ia operator -( const Vec3ia& a ) { return Vec3ia(-a.x,-a.y,-a.z); }
  51. __forceinline Vec3ia abs ( const Vec3ia& a ) { return Vec3ia(sycl::abs(a.x),sycl::abs(a.y),sycl::abs(a.z)); }
  52. ////////////////////////////////////////////////////////////////////////////////
  53. /// Binary Operators
  54. ////////////////////////////////////////////////////////////////////////////////
  55. __forceinline Vec3ia operator +( const Vec3ia& a, const Vec3ia& b ) { return Vec3ia(a.x+b.x, a.y+b.y, a.z+b.z); }
  56. __forceinline Vec3ia operator +( const Vec3ia& a, const int b ) { return a+Vec3ia(b); }
  57. __forceinline Vec3ia operator +( const int a, const Vec3ia& b ) { return Vec3ia(a)+b; }
  58. __forceinline Vec3ia operator -( const Vec3ia& a, const Vec3ia& b ) { return Vec3ia(a.x-b.x, a.y-b.y, a.z-b.z); }
  59. __forceinline Vec3ia operator -( const Vec3ia& a, const int b ) { return a-Vec3ia(b); }
  60. __forceinline Vec3ia operator -( const int a, const Vec3ia& b ) { return Vec3ia(a)-b; }
  61. __forceinline Vec3ia operator *( const Vec3ia& a, const Vec3ia& b ) { return Vec3ia(a.x*b.x, a.y*b.y, a.z*b.z); }
  62. __forceinline Vec3ia operator *( const Vec3ia& a, const int b ) { return a * Vec3ia(b); }
  63. __forceinline Vec3ia operator *( const int a, const Vec3ia& b ) { return Vec3ia(a) * b; }
  64. __forceinline Vec3ia operator &( const Vec3ia& a, const Vec3ia& b ) { return Vec3ia(a.x&b.x, a.y&b.y, a.z&b.z); }
  65. __forceinline Vec3ia operator &( const Vec3ia& a, const int b ) { return a & Vec3ia(b); }
  66. __forceinline Vec3ia operator &( const int a, const Vec3ia& b ) { return Vec3ia(a) & b; }
  67. __forceinline Vec3ia operator |( const Vec3ia& a, const Vec3ia& b ) { return Vec3ia(a.x|b.x, a.y|b.y, a.z|b.z); }
  68. __forceinline Vec3ia operator |( const Vec3ia& a, const int b ) { return a | Vec3ia(b); }
  69. __forceinline Vec3ia operator |( const int a, const Vec3ia& b ) { return Vec3ia(a) | b; }
  70. __forceinline Vec3ia operator ^( const Vec3ia& a, const Vec3ia& b ) { return Vec3ia(a.x^b.x, a.y^b.y, a.z^b.z); }
  71. __forceinline Vec3ia operator ^( const Vec3ia& a, const int b ) { return a ^ Vec3ia(b); }
  72. __forceinline Vec3ia operator ^( const int a, const Vec3ia& b ) { return Vec3ia(a) ^ b; }
  73. __forceinline Vec3ia operator <<( const Vec3ia& a, const int n ) { return Vec3ia(a.x<<n, a.y<<n, a.z<<n); }
  74. __forceinline Vec3ia operator >>( const Vec3ia& a, const int n ) { return Vec3ia(a.x>>n, a.y>>n, a.z>>n); }
  75. __forceinline Vec3ia sll ( const Vec3ia& a, const int b ) { return Vec3ia(a.x<<b, a.y<<b, a.z<<b); }
  76. __forceinline Vec3ia sra ( const Vec3ia& a, const int b ) { return Vec3ia(a.x>>b, a.y>>b, a.z>>b); }
  77. __forceinline Vec3ia srl ( const Vec3ia& a, const int b ) { return Vec3ia(unsigned(a.x)>>b, unsigned(a.y)>>b, unsigned(a.z)>>b); }
  78. ////////////////////////////////////////////////////////////////////////////////
  79. /// Assignment Operators
  80. ////////////////////////////////////////////////////////////////////////////////
  81. __forceinline Vec3ia& operator +=( Vec3ia& a, const Vec3ia& b ) { return a = a + b; }
  82. __forceinline Vec3ia& operator +=( Vec3ia& a, const int& b ) { return a = a + b; }
  83. __forceinline Vec3ia& operator -=( Vec3ia& a, const Vec3ia& b ) { return a = a - b; }
  84. __forceinline Vec3ia& operator -=( Vec3ia& a, const int& b ) { return a = a - b; }
  85. __forceinline Vec3ia& operator *=( Vec3ia& a, const Vec3ia& b ) { return a = a * b; }
  86. __forceinline Vec3ia& operator *=( Vec3ia& a, const int& b ) { return a = a * b; }
  87. __forceinline Vec3ia& operator &=( Vec3ia& a, const Vec3ia& b ) { return a = a & b; }
  88. __forceinline Vec3ia& operator &=( Vec3ia& a, const int& b ) { return a = a & b; }
  89. __forceinline Vec3ia& operator |=( Vec3ia& a, const Vec3ia& b ) { return a = a | b; }
  90. __forceinline Vec3ia& operator |=( Vec3ia& a, const int& b ) { return a = a | b; }
  91. __forceinline Vec3ia& operator <<=( Vec3ia& a, const int& b ) { return a = a << b; }
  92. __forceinline Vec3ia& operator >>=( Vec3ia& a, const int& b ) { return a = a >> b; }
  93. ////////////////////////////////////////////////////////////////////////////////
  94. /// Reductions
  95. ////////////////////////////////////////////////////////////////////////////////
  96. __forceinline int reduce_add(const Vec3ia& v) { return v.x+v.y+v.z; }
  97. __forceinline int reduce_mul(const Vec3ia& v) { return v.x*v.y*v.z; }
  98. __forceinline int reduce_min(const Vec3ia& v) { return sycl::min(sycl::min(v.x,v.y),v.z); }
  99. __forceinline int reduce_max(const Vec3ia& v) { return sycl::max(sycl::max(v.x,v.y),v.z); }
  100. ////////////////////////////////////////////////////////////////////////////////
  101. /// Comparison Operators
  102. ////////////////////////////////////////////////////////////////////////////////
  103. __forceinline bool operator ==( const Vec3ia& a, const Vec3ia& b ) { return a.x == b.x & a.y == b.y & a.z == b.z; }
  104. __forceinline bool operator !=( const Vec3ia& a, const Vec3ia& b ) { return a.x != b.x & a.y != b.y & a.z != b.z; }
  105. /*
  106. __forceinline bool operator < ( const Vec3ia& a, const Vec3ia& b ) {
  107. if (a.x != b.x) return a.x < b.x;
  108. if (a.y != b.y) return a.y < b.y;
  109. if (a.z != b.z) return a.z < b.z;
  110. return false;
  111. }
  112. */
  113. __forceinline Vec3ba eq_mask( const Vec3ia& a, const Vec3ia& b ) { return Vec3ba(a.x == b.x, a.y == b.y, a.z == b.z); }
  114. __forceinline Vec3ba lt_mask( const Vec3ia& a, const Vec3ia& b ) { return Vec3ba(a.x < b.x, a.y < b.y, a.z < b.z); }
  115. __forceinline Vec3ba gt_mask( const Vec3ia& a, const Vec3ia& b ) { return Vec3ba(a.x > b.x, a.y > b.y, a.z > b.z); }
  116. ////////////////////////////////////////////////////////////////////////////////
  117. /// Select
  118. ////////////////////////////////////////////////////////////////////////////////
  119. __forceinline Vec3ia select( const Vec3ba& m, const Vec3ia& t, const Vec3ia& f ) {
  120. const int x = m.x ? t.x : f.x;
  121. const int y = m.y ? t.y : f.y;
  122. const int z = m.z ? t.z : f.z;
  123. return Vec3ia(x,y,z);
  124. }
  125. __forceinline Vec3ia min( const Vec3ia& a, const Vec3ia& b ) { return Vec3ia(sycl::min(a.x,b.x), sycl::min(a.y,b.y), sycl::min(a.z,b.z)); }
  126. __forceinline Vec3ia max( const Vec3ia& a, const Vec3ia& b ) { return Vec3ia(sycl::max(a.x,b.x), sycl::max(a.y,b.y), sycl::max(a.z,b.z)); }
  127. ////////////////////////////////////////////////////////////////////////////////
  128. /// Output Operators
  129. ////////////////////////////////////////////////////////////////////////////////
  130. inline embree_ostream operator<<(embree_ostream cout, const Vec3ia& a) {
  131. return cout;
  132. }
  133. }