umappp 0.1.2

Sign up to get free protection for your applications and to get access to all the features.
Files changed (395) hide show
  1. checksums.yaml +7 -0
  2. data/LICENSE.txt +25 -0
  3. data/README.md +110 -0
  4. data/ext/umappp/extconf.rb +25 -0
  5. data/ext/umappp/numo.hpp +867 -0
  6. data/ext/umappp/umappp.cpp +225 -0
  7. data/lib/umappp/version.rb +5 -0
  8. data/lib/umappp.rb +41 -0
  9. data/vendor/Eigen/Cholesky +45 -0
  10. data/vendor/Eigen/CholmodSupport +48 -0
  11. data/vendor/Eigen/Core +384 -0
  12. data/vendor/Eigen/Dense +7 -0
  13. data/vendor/Eigen/Eigen +2 -0
  14. data/vendor/Eigen/Eigenvalues +60 -0
  15. data/vendor/Eigen/Geometry +59 -0
  16. data/vendor/Eigen/Householder +29 -0
  17. data/vendor/Eigen/IterativeLinearSolvers +48 -0
  18. data/vendor/Eigen/Jacobi +32 -0
  19. data/vendor/Eigen/KLUSupport +41 -0
  20. data/vendor/Eigen/LU +47 -0
  21. data/vendor/Eigen/MetisSupport +35 -0
  22. data/vendor/Eigen/OrderingMethods +70 -0
  23. data/vendor/Eigen/PaStiXSupport +49 -0
  24. data/vendor/Eigen/PardisoSupport +35 -0
  25. data/vendor/Eigen/QR +50 -0
  26. data/vendor/Eigen/QtAlignedMalloc +39 -0
  27. data/vendor/Eigen/SPQRSupport +34 -0
  28. data/vendor/Eigen/SVD +50 -0
  29. data/vendor/Eigen/Sparse +34 -0
  30. data/vendor/Eigen/SparseCholesky +37 -0
  31. data/vendor/Eigen/SparseCore +69 -0
  32. data/vendor/Eigen/SparseLU +50 -0
  33. data/vendor/Eigen/SparseQR +36 -0
  34. data/vendor/Eigen/StdDeque +27 -0
  35. data/vendor/Eigen/StdList +26 -0
  36. data/vendor/Eigen/StdVector +27 -0
  37. data/vendor/Eigen/SuperLUSupport +64 -0
  38. data/vendor/Eigen/UmfPackSupport +40 -0
  39. data/vendor/Eigen/src/Cholesky/LDLT.h +688 -0
  40. data/vendor/Eigen/src/Cholesky/LLT.h +558 -0
  41. data/vendor/Eigen/src/Cholesky/LLT_LAPACKE.h +99 -0
  42. data/vendor/Eigen/src/CholmodSupport/CholmodSupport.h +682 -0
  43. data/vendor/Eigen/src/Core/ArithmeticSequence.h +413 -0
  44. data/vendor/Eigen/src/Core/Array.h +417 -0
  45. data/vendor/Eigen/src/Core/ArrayBase.h +226 -0
  46. data/vendor/Eigen/src/Core/ArrayWrapper.h +209 -0
  47. data/vendor/Eigen/src/Core/Assign.h +90 -0
  48. data/vendor/Eigen/src/Core/AssignEvaluator.h +1010 -0
  49. data/vendor/Eigen/src/Core/Assign_MKL.h +178 -0
  50. data/vendor/Eigen/src/Core/BandMatrix.h +353 -0
  51. data/vendor/Eigen/src/Core/Block.h +448 -0
  52. data/vendor/Eigen/src/Core/BooleanRedux.h +162 -0
  53. data/vendor/Eigen/src/Core/CommaInitializer.h +164 -0
  54. data/vendor/Eigen/src/Core/ConditionEstimator.h +175 -0
  55. data/vendor/Eigen/src/Core/CoreEvaluators.h +1741 -0
  56. data/vendor/Eigen/src/Core/CoreIterators.h +132 -0
  57. data/vendor/Eigen/src/Core/CwiseBinaryOp.h +183 -0
  58. data/vendor/Eigen/src/Core/CwiseNullaryOp.h +1001 -0
  59. data/vendor/Eigen/src/Core/CwiseTernaryOp.h +197 -0
  60. data/vendor/Eigen/src/Core/CwiseUnaryOp.h +103 -0
  61. data/vendor/Eigen/src/Core/CwiseUnaryView.h +132 -0
  62. data/vendor/Eigen/src/Core/DenseBase.h +701 -0
  63. data/vendor/Eigen/src/Core/DenseCoeffsBase.h +685 -0
  64. data/vendor/Eigen/src/Core/DenseStorage.h +652 -0
  65. data/vendor/Eigen/src/Core/Diagonal.h +258 -0
  66. data/vendor/Eigen/src/Core/DiagonalMatrix.h +391 -0
  67. data/vendor/Eigen/src/Core/DiagonalProduct.h +28 -0
  68. data/vendor/Eigen/src/Core/Dot.h +318 -0
  69. data/vendor/Eigen/src/Core/EigenBase.h +160 -0
  70. data/vendor/Eigen/src/Core/ForceAlignedAccess.h +150 -0
  71. data/vendor/Eigen/src/Core/Fuzzy.h +155 -0
  72. data/vendor/Eigen/src/Core/GeneralProduct.h +465 -0
  73. data/vendor/Eigen/src/Core/GenericPacketMath.h +1040 -0
  74. data/vendor/Eigen/src/Core/GlobalFunctions.h +194 -0
  75. data/vendor/Eigen/src/Core/IO.h +258 -0
  76. data/vendor/Eigen/src/Core/IndexedView.h +237 -0
  77. data/vendor/Eigen/src/Core/Inverse.h +117 -0
  78. data/vendor/Eigen/src/Core/Map.h +171 -0
  79. data/vendor/Eigen/src/Core/MapBase.h +310 -0
  80. data/vendor/Eigen/src/Core/MathFunctions.h +2057 -0
  81. data/vendor/Eigen/src/Core/MathFunctionsImpl.h +200 -0
  82. data/vendor/Eigen/src/Core/Matrix.h +565 -0
  83. data/vendor/Eigen/src/Core/MatrixBase.h +547 -0
  84. data/vendor/Eigen/src/Core/NestByValue.h +85 -0
  85. data/vendor/Eigen/src/Core/NoAlias.h +109 -0
  86. data/vendor/Eigen/src/Core/NumTraits.h +335 -0
  87. data/vendor/Eigen/src/Core/PartialReduxEvaluator.h +232 -0
  88. data/vendor/Eigen/src/Core/PermutationMatrix.h +605 -0
  89. data/vendor/Eigen/src/Core/PlainObjectBase.h +1128 -0
  90. data/vendor/Eigen/src/Core/Product.h +191 -0
  91. data/vendor/Eigen/src/Core/ProductEvaluators.h +1179 -0
  92. data/vendor/Eigen/src/Core/Random.h +218 -0
  93. data/vendor/Eigen/src/Core/Redux.h +515 -0
  94. data/vendor/Eigen/src/Core/Ref.h +381 -0
  95. data/vendor/Eigen/src/Core/Replicate.h +142 -0
  96. data/vendor/Eigen/src/Core/Reshaped.h +454 -0
  97. data/vendor/Eigen/src/Core/ReturnByValue.h +119 -0
  98. data/vendor/Eigen/src/Core/Reverse.h +217 -0
  99. data/vendor/Eigen/src/Core/Select.h +164 -0
  100. data/vendor/Eigen/src/Core/SelfAdjointView.h +365 -0
  101. data/vendor/Eigen/src/Core/SelfCwiseBinaryOp.h +47 -0
  102. data/vendor/Eigen/src/Core/Solve.h +188 -0
  103. data/vendor/Eigen/src/Core/SolveTriangular.h +235 -0
  104. data/vendor/Eigen/src/Core/SolverBase.h +168 -0
  105. data/vendor/Eigen/src/Core/StableNorm.h +251 -0
  106. data/vendor/Eigen/src/Core/StlIterators.h +463 -0
  107. data/vendor/Eigen/src/Core/Stride.h +116 -0
  108. data/vendor/Eigen/src/Core/Swap.h +68 -0
  109. data/vendor/Eigen/src/Core/Transpose.h +464 -0
  110. data/vendor/Eigen/src/Core/Transpositions.h +386 -0
  111. data/vendor/Eigen/src/Core/TriangularMatrix.h +1001 -0
  112. data/vendor/Eigen/src/Core/VectorBlock.h +96 -0
  113. data/vendor/Eigen/src/Core/VectorwiseOp.h +784 -0
  114. data/vendor/Eigen/src/Core/Visitor.h +381 -0
  115. data/vendor/Eigen/src/Core/arch/AVX/Complex.h +372 -0
  116. data/vendor/Eigen/src/Core/arch/AVX/MathFunctions.h +228 -0
  117. data/vendor/Eigen/src/Core/arch/AVX/PacketMath.h +1574 -0
  118. data/vendor/Eigen/src/Core/arch/AVX/TypeCasting.h +115 -0
  119. data/vendor/Eigen/src/Core/arch/AVX512/Complex.h +422 -0
  120. data/vendor/Eigen/src/Core/arch/AVX512/MathFunctions.h +362 -0
  121. data/vendor/Eigen/src/Core/arch/AVX512/PacketMath.h +2303 -0
  122. data/vendor/Eigen/src/Core/arch/AVX512/TypeCasting.h +89 -0
  123. data/vendor/Eigen/src/Core/arch/AltiVec/Complex.h +417 -0
  124. data/vendor/Eigen/src/Core/arch/AltiVec/MathFunctions.h +90 -0
  125. data/vendor/Eigen/src/Core/arch/AltiVec/MatrixProduct.h +2937 -0
  126. data/vendor/Eigen/src/Core/arch/AltiVec/MatrixProductCommon.h +221 -0
  127. data/vendor/Eigen/src/Core/arch/AltiVec/MatrixProductMMA.h +629 -0
  128. data/vendor/Eigen/src/Core/arch/AltiVec/PacketMath.h +2711 -0
  129. data/vendor/Eigen/src/Core/arch/CUDA/Complex.h +258 -0
  130. data/vendor/Eigen/src/Core/arch/Default/BFloat16.h +700 -0
  131. data/vendor/Eigen/src/Core/arch/Default/ConjHelper.h +117 -0
  132. data/vendor/Eigen/src/Core/arch/Default/GenericPacketMathFunctions.h +1649 -0
  133. data/vendor/Eigen/src/Core/arch/Default/GenericPacketMathFunctionsFwd.h +110 -0
  134. data/vendor/Eigen/src/Core/arch/Default/Half.h +942 -0
  135. data/vendor/Eigen/src/Core/arch/Default/Settings.h +49 -0
  136. data/vendor/Eigen/src/Core/arch/Default/TypeCasting.h +120 -0
  137. data/vendor/Eigen/src/Core/arch/GPU/MathFunctions.h +103 -0
  138. data/vendor/Eigen/src/Core/arch/GPU/PacketMath.h +1685 -0
  139. data/vendor/Eigen/src/Core/arch/GPU/TypeCasting.h +80 -0
  140. data/vendor/Eigen/src/Core/arch/HIP/hcc/math_constants.h +23 -0
  141. data/vendor/Eigen/src/Core/arch/MSA/Complex.h +648 -0
  142. data/vendor/Eigen/src/Core/arch/MSA/MathFunctions.h +387 -0
  143. data/vendor/Eigen/src/Core/arch/MSA/PacketMath.h +1233 -0
  144. data/vendor/Eigen/src/Core/arch/NEON/Complex.h +584 -0
  145. data/vendor/Eigen/src/Core/arch/NEON/GeneralBlockPanelKernel.h +183 -0
  146. data/vendor/Eigen/src/Core/arch/NEON/MathFunctions.h +75 -0
  147. data/vendor/Eigen/src/Core/arch/NEON/PacketMath.h +4587 -0
  148. data/vendor/Eigen/src/Core/arch/NEON/TypeCasting.h +1419 -0
  149. data/vendor/Eigen/src/Core/arch/SSE/Complex.h +351 -0
  150. data/vendor/Eigen/src/Core/arch/SSE/MathFunctions.h +199 -0
  151. data/vendor/Eigen/src/Core/arch/SSE/PacketMath.h +1505 -0
  152. data/vendor/Eigen/src/Core/arch/SSE/TypeCasting.h +142 -0
  153. data/vendor/Eigen/src/Core/arch/SVE/MathFunctions.h +44 -0
  154. data/vendor/Eigen/src/Core/arch/SVE/PacketMath.h +752 -0
  155. data/vendor/Eigen/src/Core/arch/SVE/TypeCasting.h +49 -0
  156. data/vendor/Eigen/src/Core/arch/SYCL/InteropHeaders.h +232 -0
  157. data/vendor/Eigen/src/Core/arch/SYCL/MathFunctions.h +301 -0
  158. data/vendor/Eigen/src/Core/arch/SYCL/PacketMath.h +670 -0
  159. data/vendor/Eigen/src/Core/arch/SYCL/SyclMemoryModel.h +694 -0
  160. data/vendor/Eigen/src/Core/arch/SYCL/TypeCasting.h +85 -0
  161. data/vendor/Eigen/src/Core/arch/ZVector/Complex.h +426 -0
  162. data/vendor/Eigen/src/Core/arch/ZVector/MathFunctions.h +233 -0
  163. data/vendor/Eigen/src/Core/arch/ZVector/PacketMath.h +1060 -0
  164. data/vendor/Eigen/src/Core/functors/AssignmentFunctors.h +177 -0
  165. data/vendor/Eigen/src/Core/functors/BinaryFunctors.h +541 -0
  166. data/vendor/Eigen/src/Core/functors/NullaryFunctors.h +189 -0
  167. data/vendor/Eigen/src/Core/functors/StlFunctors.h +166 -0
  168. data/vendor/Eigen/src/Core/functors/TernaryFunctors.h +25 -0
  169. data/vendor/Eigen/src/Core/functors/UnaryFunctors.h +1131 -0
  170. data/vendor/Eigen/src/Core/products/GeneralBlockPanelKernel.h +2645 -0
  171. data/vendor/Eigen/src/Core/products/GeneralMatrixMatrix.h +517 -0
  172. data/vendor/Eigen/src/Core/products/GeneralMatrixMatrixTriangular.h +317 -0
  173. data/vendor/Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h +145 -0
  174. data/vendor/Eigen/src/Core/products/GeneralMatrixMatrix_BLAS.h +124 -0
  175. data/vendor/Eigen/src/Core/products/GeneralMatrixVector.h +518 -0
  176. data/vendor/Eigen/src/Core/products/GeneralMatrixVector_BLAS.h +136 -0
  177. data/vendor/Eigen/src/Core/products/Parallelizer.h +180 -0
  178. data/vendor/Eigen/src/Core/products/SelfadjointMatrixMatrix.h +544 -0
  179. data/vendor/Eigen/src/Core/products/SelfadjointMatrixMatrix_BLAS.h +295 -0
  180. data/vendor/Eigen/src/Core/products/SelfadjointMatrixVector.h +262 -0
  181. data/vendor/Eigen/src/Core/products/SelfadjointMatrixVector_BLAS.h +118 -0
  182. data/vendor/Eigen/src/Core/products/SelfadjointProduct.h +133 -0
  183. data/vendor/Eigen/src/Core/products/SelfadjointRank2Update.h +94 -0
  184. data/vendor/Eigen/src/Core/products/TriangularMatrixMatrix.h +472 -0
  185. data/vendor/Eigen/src/Core/products/TriangularMatrixMatrix_BLAS.h +317 -0
  186. data/vendor/Eigen/src/Core/products/TriangularMatrixVector.h +350 -0
  187. data/vendor/Eigen/src/Core/products/TriangularMatrixVector_BLAS.h +255 -0
  188. data/vendor/Eigen/src/Core/products/TriangularSolverMatrix.h +337 -0
  189. data/vendor/Eigen/src/Core/products/TriangularSolverMatrix_BLAS.h +167 -0
  190. data/vendor/Eigen/src/Core/products/TriangularSolverVector.h +148 -0
  191. data/vendor/Eigen/src/Core/util/BlasUtil.h +583 -0
  192. data/vendor/Eigen/src/Core/util/ConfigureVectorization.h +512 -0
  193. data/vendor/Eigen/src/Core/util/Constants.h +563 -0
  194. data/vendor/Eigen/src/Core/util/DisableStupidWarnings.h +106 -0
  195. data/vendor/Eigen/src/Core/util/ForwardDeclarations.h +322 -0
  196. data/vendor/Eigen/src/Core/util/IndexedViewHelper.h +186 -0
  197. data/vendor/Eigen/src/Core/util/IntegralConstant.h +272 -0
  198. data/vendor/Eigen/src/Core/util/MKL_support.h +137 -0
  199. data/vendor/Eigen/src/Core/util/Macros.h +1464 -0
  200. data/vendor/Eigen/src/Core/util/Memory.h +1163 -0
  201. data/vendor/Eigen/src/Core/util/Meta.h +812 -0
  202. data/vendor/Eigen/src/Core/util/NonMPL2.h +3 -0
  203. data/vendor/Eigen/src/Core/util/ReenableStupidWarnings.h +31 -0
  204. data/vendor/Eigen/src/Core/util/ReshapedHelper.h +51 -0
  205. data/vendor/Eigen/src/Core/util/StaticAssert.h +221 -0
  206. data/vendor/Eigen/src/Core/util/SymbolicIndex.h +293 -0
  207. data/vendor/Eigen/src/Core/util/XprHelper.h +856 -0
  208. data/vendor/Eigen/src/Eigenvalues/ComplexEigenSolver.h +346 -0
  209. data/vendor/Eigen/src/Eigenvalues/ComplexSchur.h +462 -0
  210. data/vendor/Eigen/src/Eigenvalues/ComplexSchur_LAPACKE.h +91 -0
  211. data/vendor/Eigen/src/Eigenvalues/EigenSolver.h +622 -0
  212. data/vendor/Eigen/src/Eigenvalues/GeneralizedEigenSolver.h +418 -0
  213. data/vendor/Eigen/src/Eigenvalues/GeneralizedSelfAdjointEigenSolver.h +226 -0
  214. data/vendor/Eigen/src/Eigenvalues/HessenbergDecomposition.h +374 -0
  215. data/vendor/Eigen/src/Eigenvalues/MatrixBaseEigenvalues.h +158 -0
  216. data/vendor/Eigen/src/Eigenvalues/RealQZ.h +657 -0
  217. data/vendor/Eigen/src/Eigenvalues/RealSchur.h +558 -0
  218. data/vendor/Eigen/src/Eigenvalues/RealSchur_LAPACKE.h +77 -0
  219. data/vendor/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h +904 -0
  220. data/vendor/Eigen/src/Eigenvalues/SelfAdjointEigenSolver_LAPACKE.h +87 -0
  221. data/vendor/Eigen/src/Eigenvalues/Tridiagonalization.h +561 -0
  222. data/vendor/Eigen/src/Geometry/AlignedBox.h +486 -0
  223. data/vendor/Eigen/src/Geometry/AngleAxis.h +247 -0
  224. data/vendor/Eigen/src/Geometry/EulerAngles.h +114 -0
  225. data/vendor/Eigen/src/Geometry/Homogeneous.h +501 -0
  226. data/vendor/Eigen/src/Geometry/Hyperplane.h +282 -0
  227. data/vendor/Eigen/src/Geometry/OrthoMethods.h +235 -0
  228. data/vendor/Eigen/src/Geometry/ParametrizedLine.h +232 -0
  229. data/vendor/Eigen/src/Geometry/Quaternion.h +870 -0
  230. data/vendor/Eigen/src/Geometry/Rotation2D.h +199 -0
  231. data/vendor/Eigen/src/Geometry/RotationBase.h +206 -0
  232. data/vendor/Eigen/src/Geometry/Scaling.h +188 -0
  233. data/vendor/Eigen/src/Geometry/Transform.h +1563 -0
  234. data/vendor/Eigen/src/Geometry/Translation.h +202 -0
  235. data/vendor/Eigen/src/Geometry/Umeyama.h +166 -0
  236. data/vendor/Eigen/src/Geometry/arch/Geometry_SIMD.h +168 -0
  237. data/vendor/Eigen/src/Householder/BlockHouseholder.h +110 -0
  238. data/vendor/Eigen/src/Householder/Householder.h +176 -0
  239. data/vendor/Eigen/src/Householder/HouseholderSequence.h +545 -0
  240. data/vendor/Eigen/src/IterativeLinearSolvers/BasicPreconditioners.h +226 -0
  241. data/vendor/Eigen/src/IterativeLinearSolvers/BiCGSTAB.h +212 -0
  242. data/vendor/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h +229 -0
  243. data/vendor/Eigen/src/IterativeLinearSolvers/IncompleteCholesky.h +394 -0
  244. data/vendor/Eigen/src/IterativeLinearSolvers/IncompleteLUT.h +453 -0
  245. data/vendor/Eigen/src/IterativeLinearSolvers/IterativeSolverBase.h +444 -0
  246. data/vendor/Eigen/src/IterativeLinearSolvers/LeastSquareConjugateGradient.h +198 -0
  247. data/vendor/Eigen/src/IterativeLinearSolvers/SolveWithGuess.h +117 -0
  248. data/vendor/Eigen/src/Jacobi/Jacobi.h +483 -0
  249. data/vendor/Eigen/src/KLUSupport/KLUSupport.h +358 -0
  250. data/vendor/Eigen/src/LU/Determinant.h +117 -0
  251. data/vendor/Eigen/src/LU/FullPivLU.h +877 -0
  252. data/vendor/Eigen/src/LU/InverseImpl.h +432 -0
  253. data/vendor/Eigen/src/LU/PartialPivLU.h +624 -0
  254. data/vendor/Eigen/src/LU/PartialPivLU_LAPACKE.h +83 -0
  255. data/vendor/Eigen/src/LU/arch/InverseSize4.h +351 -0
  256. data/vendor/Eigen/src/MetisSupport/MetisSupport.h +137 -0
  257. data/vendor/Eigen/src/OrderingMethods/Amd.h +435 -0
  258. data/vendor/Eigen/src/OrderingMethods/Eigen_Colamd.h +1863 -0
  259. data/vendor/Eigen/src/OrderingMethods/Ordering.h +153 -0
  260. data/vendor/Eigen/src/PaStiXSupport/PaStiXSupport.h +678 -0
  261. data/vendor/Eigen/src/PardisoSupport/PardisoSupport.h +545 -0
  262. data/vendor/Eigen/src/QR/ColPivHouseholderQR.h +674 -0
  263. data/vendor/Eigen/src/QR/ColPivHouseholderQR_LAPACKE.h +97 -0
  264. data/vendor/Eigen/src/QR/CompleteOrthogonalDecomposition.h +635 -0
  265. data/vendor/Eigen/src/QR/FullPivHouseholderQR.h +713 -0
  266. data/vendor/Eigen/src/QR/HouseholderQR.h +434 -0
  267. data/vendor/Eigen/src/QR/HouseholderQR_LAPACKE.h +68 -0
  268. data/vendor/Eigen/src/SPQRSupport/SuiteSparseQRSupport.h +335 -0
  269. data/vendor/Eigen/src/SVD/BDCSVD.h +1366 -0
  270. data/vendor/Eigen/src/SVD/JacobiSVD.h +812 -0
  271. data/vendor/Eigen/src/SVD/JacobiSVD_LAPACKE.h +91 -0
  272. data/vendor/Eigen/src/SVD/SVDBase.h +376 -0
  273. data/vendor/Eigen/src/SVD/UpperBidiagonalization.h +414 -0
  274. data/vendor/Eigen/src/SparseCholesky/SimplicialCholesky.h +697 -0
  275. data/vendor/Eigen/src/SparseCholesky/SimplicialCholesky_impl.h +174 -0
  276. data/vendor/Eigen/src/SparseCore/AmbiVector.h +378 -0
  277. data/vendor/Eigen/src/SparseCore/CompressedStorage.h +274 -0
  278. data/vendor/Eigen/src/SparseCore/ConservativeSparseSparseProduct.h +352 -0
  279. data/vendor/Eigen/src/SparseCore/MappedSparseMatrix.h +67 -0
  280. data/vendor/Eigen/src/SparseCore/SparseAssign.h +270 -0
  281. data/vendor/Eigen/src/SparseCore/SparseBlock.h +571 -0
  282. data/vendor/Eigen/src/SparseCore/SparseColEtree.h +206 -0
  283. data/vendor/Eigen/src/SparseCore/SparseCompressedBase.h +370 -0
  284. data/vendor/Eigen/src/SparseCore/SparseCwiseBinaryOp.h +722 -0
  285. data/vendor/Eigen/src/SparseCore/SparseCwiseUnaryOp.h +150 -0
  286. data/vendor/Eigen/src/SparseCore/SparseDenseProduct.h +342 -0
  287. data/vendor/Eigen/src/SparseCore/SparseDiagonalProduct.h +138 -0
  288. data/vendor/Eigen/src/SparseCore/SparseDot.h +98 -0
  289. data/vendor/Eigen/src/SparseCore/SparseFuzzy.h +29 -0
  290. data/vendor/Eigen/src/SparseCore/SparseMap.h +305 -0
  291. data/vendor/Eigen/src/SparseCore/SparseMatrix.h +1518 -0
  292. data/vendor/Eigen/src/SparseCore/SparseMatrixBase.h +398 -0
  293. data/vendor/Eigen/src/SparseCore/SparsePermutation.h +178 -0
  294. data/vendor/Eigen/src/SparseCore/SparseProduct.h +181 -0
  295. data/vendor/Eigen/src/SparseCore/SparseRedux.h +49 -0
  296. data/vendor/Eigen/src/SparseCore/SparseRef.h +397 -0
  297. data/vendor/Eigen/src/SparseCore/SparseSelfAdjointView.h +659 -0
  298. data/vendor/Eigen/src/SparseCore/SparseSolverBase.h +124 -0
  299. data/vendor/Eigen/src/SparseCore/SparseSparseProductWithPruning.h +198 -0
  300. data/vendor/Eigen/src/SparseCore/SparseTranspose.h +92 -0
  301. data/vendor/Eigen/src/SparseCore/SparseTriangularView.h +189 -0
  302. data/vendor/Eigen/src/SparseCore/SparseUtil.h +186 -0
  303. data/vendor/Eigen/src/SparseCore/SparseVector.h +478 -0
  304. data/vendor/Eigen/src/SparseCore/SparseView.h +254 -0
  305. data/vendor/Eigen/src/SparseCore/TriangularSolver.h +315 -0
  306. data/vendor/Eigen/src/SparseLU/SparseLU.h +923 -0
  307. data/vendor/Eigen/src/SparseLU/SparseLUImpl.h +66 -0
  308. data/vendor/Eigen/src/SparseLU/SparseLU_Memory.h +226 -0
  309. data/vendor/Eigen/src/SparseLU/SparseLU_Structs.h +110 -0
  310. data/vendor/Eigen/src/SparseLU/SparseLU_SupernodalMatrix.h +375 -0
  311. data/vendor/Eigen/src/SparseLU/SparseLU_Utils.h +80 -0
  312. data/vendor/Eigen/src/SparseLU/SparseLU_column_bmod.h +181 -0
  313. data/vendor/Eigen/src/SparseLU/SparseLU_column_dfs.h +179 -0
  314. data/vendor/Eigen/src/SparseLU/SparseLU_copy_to_ucol.h +107 -0
  315. data/vendor/Eigen/src/SparseLU/SparseLU_gemm_kernel.h +280 -0
  316. data/vendor/Eigen/src/SparseLU/SparseLU_heap_relax_snode.h +126 -0
  317. data/vendor/Eigen/src/SparseLU/SparseLU_kernel_bmod.h +130 -0
  318. data/vendor/Eigen/src/SparseLU/SparseLU_panel_bmod.h +223 -0
  319. data/vendor/Eigen/src/SparseLU/SparseLU_panel_dfs.h +258 -0
  320. data/vendor/Eigen/src/SparseLU/SparseLU_pivotL.h +137 -0
  321. data/vendor/Eigen/src/SparseLU/SparseLU_pruneL.h +136 -0
  322. data/vendor/Eigen/src/SparseLU/SparseLU_relax_snode.h +83 -0
  323. data/vendor/Eigen/src/SparseQR/SparseQR.h +758 -0
  324. data/vendor/Eigen/src/StlSupport/StdDeque.h +116 -0
  325. data/vendor/Eigen/src/StlSupport/StdList.h +106 -0
  326. data/vendor/Eigen/src/StlSupport/StdVector.h +131 -0
  327. data/vendor/Eigen/src/StlSupport/details.h +84 -0
  328. data/vendor/Eigen/src/SuperLUSupport/SuperLUSupport.h +1025 -0
  329. data/vendor/Eigen/src/UmfPackSupport/UmfPackSupport.h +642 -0
  330. data/vendor/Eigen/src/misc/Image.h +82 -0
  331. data/vendor/Eigen/src/misc/Kernel.h +79 -0
  332. data/vendor/Eigen/src/misc/RealSvd2x2.h +55 -0
  333. data/vendor/Eigen/src/misc/blas.h +440 -0
  334. data/vendor/Eigen/src/misc/lapack.h +152 -0
  335. data/vendor/Eigen/src/misc/lapacke.h +16292 -0
  336. data/vendor/Eigen/src/misc/lapacke_mangling.h +17 -0
  337. data/vendor/Eigen/src/plugins/ArrayCwiseBinaryOps.h +358 -0
  338. data/vendor/Eigen/src/plugins/ArrayCwiseUnaryOps.h +696 -0
  339. data/vendor/Eigen/src/plugins/BlockMethods.h +1442 -0
  340. data/vendor/Eigen/src/plugins/CommonCwiseBinaryOps.h +115 -0
  341. data/vendor/Eigen/src/plugins/CommonCwiseUnaryOps.h +177 -0
  342. data/vendor/Eigen/src/plugins/IndexedViewMethods.h +262 -0
  343. data/vendor/Eigen/src/plugins/MatrixCwiseBinaryOps.h +152 -0
  344. data/vendor/Eigen/src/plugins/MatrixCwiseUnaryOps.h +95 -0
  345. data/vendor/Eigen/src/plugins/ReshapedMethods.h +149 -0
  346. data/vendor/aarand/aarand.hpp +114 -0
  347. data/vendor/annoy/annoylib.h +1495 -0
  348. data/vendor/annoy/kissrandom.h +120 -0
  349. data/vendor/annoy/mman.h +242 -0
  350. data/vendor/hnswlib/bruteforce.h +152 -0
  351. data/vendor/hnswlib/hnswalg.h +1192 -0
  352. data/vendor/hnswlib/hnswlib.h +108 -0
  353. data/vendor/hnswlib/space_ip.h +282 -0
  354. data/vendor/hnswlib/space_l2.h +281 -0
  355. data/vendor/hnswlib/visited_list_pool.h +79 -0
  356. data/vendor/irlba/irlba.hpp +575 -0
  357. data/vendor/irlba/lanczos.hpp +212 -0
  358. data/vendor/irlba/parallel.hpp +474 -0
  359. data/vendor/irlba/utils.hpp +224 -0
  360. data/vendor/irlba/wrappers.hpp +228 -0
  361. data/vendor/kmeans/Base.hpp +75 -0
  362. data/vendor/kmeans/Details.hpp +79 -0
  363. data/vendor/kmeans/HartiganWong.hpp +492 -0
  364. data/vendor/kmeans/InitializeKmeansPP.hpp +144 -0
  365. data/vendor/kmeans/InitializeNone.hpp +44 -0
  366. data/vendor/kmeans/InitializePCAPartition.hpp +309 -0
  367. data/vendor/kmeans/InitializeRandom.hpp +91 -0
  368. data/vendor/kmeans/Kmeans.hpp +161 -0
  369. data/vendor/kmeans/Lloyd.hpp +134 -0
  370. data/vendor/kmeans/MiniBatch.hpp +269 -0
  371. data/vendor/kmeans/QuickSearch.hpp +179 -0
  372. data/vendor/kmeans/compute_centroids.hpp +32 -0
  373. data/vendor/kmeans/compute_wcss.hpp +27 -0
  374. data/vendor/kmeans/is_edge_case.hpp +42 -0
  375. data/vendor/kmeans/random.hpp +55 -0
  376. data/vendor/knncolle/Annoy/Annoy.hpp +193 -0
  377. data/vendor/knncolle/BruteForce/BruteForce.hpp +120 -0
  378. data/vendor/knncolle/Hnsw/Hnsw.hpp +225 -0
  379. data/vendor/knncolle/Kmknn/Kmknn.hpp +286 -0
  380. data/vendor/knncolle/VpTree/VpTree.hpp +256 -0
  381. data/vendor/knncolle/knncolle.hpp +34 -0
  382. data/vendor/knncolle/utils/Base.hpp +100 -0
  383. data/vendor/knncolle/utils/NeighborQueue.hpp +94 -0
  384. data/vendor/knncolle/utils/distances.hpp +98 -0
  385. data/vendor/knncolle/utils/find_nearest_neighbors.hpp +112 -0
  386. data/vendor/powerit/PowerIterations.hpp +157 -0
  387. data/vendor/umappp/NeighborList.hpp +37 -0
  388. data/vendor/umappp/Umap.hpp +662 -0
  389. data/vendor/umappp/combine_neighbor_sets.hpp +95 -0
  390. data/vendor/umappp/find_ab.hpp +157 -0
  391. data/vendor/umappp/neighbor_similarities.hpp +136 -0
  392. data/vendor/umappp/optimize_layout.hpp +285 -0
  393. data/vendor/umappp/spectral_init.hpp +181 -0
  394. data/vendor/umappp/umappp.hpp +13 -0
  395. metadata +465 -0
@@ -0,0 +1,694 @@
1
+ /***************************************************************************
2
+ * Copyright (C) 2017 Codeplay Software Limited
3
+ * This Source Code Form is subject to the terms of the Mozilla
4
+ * Public License v. 2.0. If a copy of the MPL was not distributed
5
+ * with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
6
+ *
7
+ *
8
+ * SyclMemoryModel.h
9
+ *
10
+ * Description:
11
+ * Interface for SYCL buffers to behave as a non-dereferenceable pointer
12
+ * Interface for Placeholder accessor to behave as a pointer on both host
13
+ * and device
14
+ *
15
+ * Authors:
16
+ *
17
+ * Ruyman Reyes Codeplay Software Ltd.
18
+ * Mehdi Goli Codeplay Software Ltd.
19
+ * Vanya Yaneva Codeplay Software Ltd.
20
+ *
21
+ **************************************************************************/
22
+
23
+ #if defined(EIGEN_USE_SYCL) && \
24
+ !defined(EIGEN_CXX11_TENSOR_TENSOR_SYCL_STORAGE_MEMORY_H)
25
+ #define EIGEN_CXX11_TENSOR_TENSOR_SYCL_STORAGE_MEMORY_H
26
+
27
+ #include <CL/sycl.hpp>
28
+ #ifdef EIGEN_EXCEPTIONS
29
+ #include <stdexcept>
30
+ #endif
31
+ #include <cstddef>
32
+ #include <queue>
33
+ #include <set>
34
+ #include <unordered_map>
35
+
36
+ namespace Eigen {
37
+ namespace TensorSycl {
38
+ namespace internal {
39
+
40
+ using sycl_acc_target = cl::sycl::access::target;
41
+ using sycl_acc_mode = cl::sycl::access::mode;
42
+
43
+ /**
44
+ * Default values for template arguments
45
+ */
46
+ using buffer_data_type_t = uint8_t;
47
+ const sycl_acc_target default_acc_target = sycl_acc_target::global_buffer;
48
+ const sycl_acc_mode default_acc_mode = sycl_acc_mode::read_write;
49
+
50
+ /**
51
+ * PointerMapper
52
+ * Associates fake pointers with buffers.
53
+ *
54
+ */
55
+ class PointerMapper {
56
+ public:
57
+ using base_ptr_t = std::intptr_t;
58
+
59
+ /* Structure of a virtual pointer
60
+ *
61
+ * |================================================|
62
+ * | POINTER ADDRESS |
63
+ * |================================================|
64
+ */
65
+ struct virtual_pointer_t {
66
+ /* Type for the pointers
67
+ */
68
+ base_ptr_t m_contents;
69
+
70
+ /** Conversions from virtual_pointer_t to
71
+ * void * should just reinterpret_cast the integer number
72
+ */
73
+ operator void *() const { return reinterpret_cast<void *>(m_contents); }
74
+
75
+ /**
76
+ * Convert back to the integer number.
77
+ */
78
+ operator base_ptr_t() const { return m_contents; }
79
+
80
+ /**
81
+ * Add a certain value to the pointer to create a
82
+ * new pointer to that offset
83
+ */
84
+ virtual_pointer_t operator+(size_t off) { return m_contents + off; }
85
+
86
+ /* Numerical order for sorting pointers in containers. */
87
+ bool operator<(virtual_pointer_t rhs) const {
88
+ return (static_cast<base_ptr_t>(m_contents) <
89
+ static_cast<base_ptr_t>(rhs.m_contents));
90
+ }
91
+
92
+ bool operator>(virtual_pointer_t rhs) const {
93
+ return (static_cast<base_ptr_t>(m_contents) >
94
+ static_cast<base_ptr_t>(rhs.m_contents));
95
+ }
96
+
97
+ /**
98
+ * Numerical order for sorting pointers in containers
99
+ */
100
+ bool operator==(virtual_pointer_t rhs) const {
101
+ return (static_cast<base_ptr_t>(m_contents) ==
102
+ static_cast<base_ptr_t>(rhs.m_contents));
103
+ }
104
+
105
+ /**
106
+ * Simple forward to the equality overload.
107
+ */
108
+ bool operator!=(virtual_pointer_t rhs) const {
109
+ return !(this->operator==(rhs));
110
+ }
111
+
112
+ /**
113
+ * Converts a void * into a virtual pointer structure.
114
+ * Note that this will only work if the void * was
115
+ * already a virtual_pointer_t, but we have no way of
116
+ * checking
117
+ */
118
+ virtual_pointer_t(const void *ptr)
119
+ : m_contents(reinterpret_cast<base_ptr_t>(ptr)){};
120
+
121
+ /**
122
+ * Creates a virtual_pointer_t from the given integer
123
+ * number
124
+ */
125
+ virtual_pointer_t(base_ptr_t u) : m_contents(u){};
126
+ };
127
+
128
+ /* Definition of a null pointer
129
+ */
130
+ const virtual_pointer_t null_virtual_ptr = nullptr;
131
+
132
+ /**
133
+ * Whether if a pointer is null or not.
134
+ * A pointer is nullptr if the value is of null_virtual_ptr
135
+ */
136
+ static inline bool is_nullptr(virtual_pointer_t ptr) {
137
+ return (static_cast<void *>(ptr) == nullptr);
138
+ }
139
+
140
+ /* basic type for all buffers
141
+ */
142
+ using buffer_t = cl::sycl::buffer_mem;
143
+
144
+ /**
145
+ * Node that stores information about a device allocation.
146
+ * Nodes are sorted by size to organise a free list of nodes
147
+ * that can be recovered.
148
+ */
149
+ struct pMapNode_t {
150
+ buffer_t m_buffer;
151
+ size_t m_size;
152
+ bool m_free;
153
+
154
+ pMapNode_t(buffer_t b, size_t size, bool f)
155
+ : m_buffer{b}, m_size{size}, m_free{f} {
156
+ m_buffer.set_final_data(nullptr);
157
+ }
158
+
159
+ bool operator<=(const pMapNode_t &rhs) { return (m_size <= rhs.m_size); }
160
+ };
161
+
162
+ /** Storage of the pointer / buffer tree
163
+ */
164
+ using pointerMap_t = std::map<virtual_pointer_t, pMapNode_t>;
165
+
166
+ /**
167
+ * Obtain the insertion point in the pointer map for
168
+ * a pointer of the given size.
169
+ * \param requiredSize Size attemted to reclaim
170
+ */
171
+ typename pointerMap_t::iterator get_insertion_point(size_t requiredSize) {
172
+ typename pointerMap_t::iterator retVal;
173
+ bool reuse = false;
174
+ if (!m_freeList.empty()) {
175
+ // try to re-use an existing block
176
+ for (auto freeElem : m_freeList) {
177
+ if (freeElem->second.m_size >= requiredSize) {
178
+ retVal = freeElem;
179
+ reuse = true;
180
+ // Element is not going to be free anymore
181
+ m_freeList.erase(freeElem);
182
+ break;
183
+ }
184
+ }
185
+ }
186
+ if (!reuse) {
187
+ retVal = std::prev(m_pointerMap.end());
188
+ }
189
+ return retVal;
190
+ }
191
+
192
+ /**
193
+ * Returns an iterator to the node that stores the information
194
+ * of the given virtual pointer from the given pointer map structure.
195
+ * If pointer is not found, throws std::out_of_range.
196
+ * If the pointer map structure is empty, throws std::out_of_range
197
+ *
198
+ * \param pMap the pointerMap_t structure storing all the pointers
199
+ * \param virtual_pointer_ptr The virtual pointer to obtain the node of
200
+ * \throws std::out:of_range if the pointer is not found or pMap is empty
201
+ */
202
+ typename pointerMap_t::iterator get_node(const virtual_pointer_t ptr) {
203
+ if (this->count() == 0) {
204
+ m_pointerMap.clear();
205
+ EIGEN_THROW_X(std::out_of_range("There are no pointers allocated\n"));
206
+
207
+ }
208
+ if (is_nullptr(ptr)) {
209
+ m_pointerMap.clear();
210
+ EIGEN_THROW_X(std::out_of_range("Cannot access null pointer\n"));
211
+ }
212
+ // The previous element to the lower bound is the node that
213
+ // holds this memory address
214
+ auto node = m_pointerMap.lower_bound(ptr);
215
+ // If the value of the pointer is not the one of the node
216
+ // then we return the previous one
217
+ if (node == std::end(m_pointerMap)) {
218
+ --node;
219
+ } else if (node->first != ptr) {
220
+ if (node == std::begin(m_pointerMap)) {
221
+ m_pointerMap.clear();
222
+ EIGEN_THROW_X(
223
+ std::out_of_range("The pointer is not registered in the map\n"));
224
+
225
+ }
226
+ --node;
227
+ }
228
+
229
+ return node;
230
+ }
231
+
232
+ /* get_buffer.
233
+ * Returns a buffer from the map using the pointer address
234
+ */
235
+ template <typename buffer_data_type = buffer_data_type_t>
236
+ cl::sycl::buffer<buffer_data_type, 1> get_buffer(
237
+ const virtual_pointer_t ptr) {
238
+ using sycl_buffer_t = cl::sycl::buffer<buffer_data_type, 1>;
239
+
240
+ // get_node() returns a `buffer_mem`, so we need to cast it to a `buffer<>`.
241
+ // We can do this without the `buffer_mem` being a pointer, as we
242
+ // only declare member variables in the base class (`buffer_mem`) and not in
243
+ // the child class (`buffer<>).
244
+ auto node = get_node(ptr);
245
+ eigen_assert(node->first == ptr || node->first < ptr);
246
+ eigen_assert(ptr < static_cast<virtual_pointer_t>(node->second.m_size +
247
+ node->first));
248
+ return *(static_cast<sycl_buffer_t *>(&node->second.m_buffer));
249
+ }
250
+
251
+ /**
252
+ * @brief Returns an accessor to the buffer of the given virtual pointer
253
+ * @param accessMode
254
+ * @param accessTarget
255
+ * @param ptr The virtual pointer
256
+ */
257
+ template <sycl_acc_mode access_mode = default_acc_mode,
258
+ sycl_acc_target access_target = default_acc_target,
259
+ typename buffer_data_type = buffer_data_type_t>
260
+ cl::sycl::accessor<buffer_data_type, 1, access_mode, access_target>
261
+ get_access(const virtual_pointer_t ptr) {
262
+ auto buf = get_buffer<buffer_data_type>(ptr);
263
+ return buf.template get_access<access_mode, access_target>();
264
+ }
265
+
266
+ /**
267
+ * @brief Returns an accessor to the buffer of the given virtual pointer
268
+ * in the given command group scope
269
+ * @param accessMode
270
+ * @param accessTarget
271
+ * @param ptr The virtual pointer
272
+ * @param cgh Reference to the command group scope
273
+ */
274
+ template <sycl_acc_mode access_mode = default_acc_mode,
275
+ sycl_acc_target access_target = default_acc_target,
276
+ typename buffer_data_type = buffer_data_type_t>
277
+ cl::sycl::accessor<buffer_data_type, 1, access_mode, access_target>
278
+ get_access(const virtual_pointer_t ptr, cl::sycl::handler &cgh) {
279
+ auto buf = get_buffer<buffer_data_type>(ptr);
280
+ return buf.template get_access<access_mode, access_target>(cgh);
281
+ }
282
+
283
+ /*
284
+ * Returns the offset from the base address of this pointer.
285
+ */
286
+ inline std::ptrdiff_t get_offset(const virtual_pointer_t ptr) {
287
+ // The previous element to the lower bound is the node that
288
+ // holds this memory address
289
+ auto node = get_node(ptr);
290
+ auto start = node->first;
291
+ eigen_assert(start == ptr || start < ptr);
292
+ eigen_assert(ptr < start + node->second.m_size);
293
+ return (ptr - start);
294
+ }
295
+
296
+ /*
297
+ * Returns the number of elements by which the given pointer is offset from
298
+ * the base address.
299
+ */
300
+ template <typename buffer_data_type>
301
+ inline size_t get_element_offset(const virtual_pointer_t ptr) {
302
+ return get_offset(ptr) / sizeof(buffer_data_type);
303
+ }
304
+
305
+ /**
306
+ * Constructs the PointerMapper structure.
307
+ */
308
+ PointerMapper(base_ptr_t baseAddress = 4096)
309
+ : m_pointerMap{}, m_freeList{}, m_baseAddress{baseAddress} {
310
+ if (m_baseAddress == 0) {
311
+ EIGEN_THROW_X(std::invalid_argument("Base address cannot be zero\n"));
312
+ }
313
+ };
314
+
315
+ /**
316
+ * PointerMapper cannot be copied or moved
317
+ */
318
+ PointerMapper(const PointerMapper &) = delete;
319
+
320
+ /**
321
+ * Empty the pointer list
322
+ */
323
+ inline void clear() {
324
+ m_freeList.clear();
325
+ m_pointerMap.clear();
326
+ }
327
+
328
+ /* add_pointer.
329
+ * Adds an existing pointer to the map and returns the virtual pointer id.
330
+ */
331
+ inline virtual_pointer_t add_pointer(const buffer_t &b) {
332
+ return add_pointer_impl(b);
333
+ }
334
+
335
+ /* add_pointer.
336
+ * Adds a pointer to the map and returns the virtual pointer id.
337
+ */
338
+ inline virtual_pointer_t add_pointer(buffer_t &&b) {
339
+ return add_pointer_impl(b);
340
+ }
341
+
342
+ /**
343
+ * @brief Fuses the given node with the previous nodes in the
344
+ * pointer map if they are free
345
+ *
346
+ * @param node A reference to the free node to be fused
347
+ */
348
+ void fuse_forward(typename pointerMap_t::iterator &node) {
349
+ while (node != std::prev(m_pointerMap.end())) {
350
+ // if following node is free
351
+ // remove it and extend the current node with its size
352
+ auto fwd_node = std::next(node);
353
+ if (!fwd_node->second.m_free) {
354
+ break;
355
+ }
356
+ auto fwd_size = fwd_node->second.m_size;
357
+ m_freeList.erase(fwd_node);
358
+ m_pointerMap.erase(fwd_node);
359
+
360
+ node->second.m_size += fwd_size;
361
+ }
362
+ }
363
+
364
+ /**
365
+ * @brief Fuses the given node with the following nodes in the
366
+ * pointer map if they are free
367
+ *
368
+ * @param node A reference to the free node to be fused
369
+ */
370
+ void fuse_backward(typename pointerMap_t::iterator &node) {
371
+ while (node != m_pointerMap.begin()) {
372
+ // if previous node is free, extend it
373
+ // with the size of the current one
374
+ auto prev_node = std::prev(node);
375
+ if (!prev_node->second.m_free) {
376
+ break;
377
+ }
378
+ prev_node->second.m_size += node->second.m_size;
379
+
380
+ // remove the current node
381
+ m_freeList.erase(node);
382
+ m_pointerMap.erase(node);
383
+
384
+ // point to the previous node
385
+ node = prev_node;
386
+ }
387
+ }
388
+
389
+ /* remove_pointer.
390
+ * Removes the given pointer from the map.
391
+ * The pointer is allowed to be reused only if ReUse if true.
392
+ */
393
+ template <bool ReUse = true>
394
+ void remove_pointer(const virtual_pointer_t ptr) {
395
+ if (is_nullptr(ptr)) {
396
+ return;
397
+ }
398
+ auto node = this->get_node(ptr);
399
+
400
+ node->second.m_free = true;
401
+ m_freeList.emplace(node);
402
+
403
+ // Fuse the node
404
+ // with free nodes before and after it
405
+ fuse_forward(node);
406
+ fuse_backward(node);
407
+
408
+ // If after fusing the node is the last one
409
+ // simply remove it (since it is free)
410
+ if (node == std::prev(m_pointerMap.end())) {
411
+ m_freeList.erase(node);
412
+ m_pointerMap.erase(node);
413
+ }
414
+ }
415
+
416
+ /* count.
417
+ * Return the number of active pointers (i.e, pointers that
418
+ * have been malloc but not freed).
419
+ */
420
+ size_t count() const { return (m_pointerMap.size() - m_freeList.size()); }
421
+
422
+ private:
423
+ /* add_pointer_impl.
424
+ * Adds a pointer to the map and returns the virtual pointer id.
425
+ * BufferT is either a const buffer_t& or a buffer_t&&.
426
+ */
427
+ template <class BufferT>
428
+ virtual_pointer_t add_pointer_impl(BufferT b) {
429
+ virtual_pointer_t retVal = nullptr;
430
+ size_t bufSize = b.get_count();
431
+ pMapNode_t p{b, bufSize, false};
432
+ // If this is the first pointer:
433
+ if (m_pointerMap.empty()) {
434
+ virtual_pointer_t initialVal{m_baseAddress};
435
+ m_pointerMap.emplace(initialVal, p);
436
+ return initialVal;
437
+ }
438
+
439
+ auto lastElemIter = get_insertion_point(bufSize);
440
+ // We are recovering an existing free node
441
+ if (lastElemIter->second.m_free) {
442
+ lastElemIter->second.m_buffer = b;
443
+ lastElemIter->second.m_free = false;
444
+
445
+ // If the recovered node is bigger than the inserted one
446
+ // add a new free node with the remaining space
447
+ if (lastElemIter->second.m_size > bufSize) {
448
+ // create a new node with the remaining space
449
+ auto remainingSize = lastElemIter->second.m_size - bufSize;
450
+ pMapNode_t p2{b, remainingSize, true};
451
+
452
+ // update size of the current node
453
+ lastElemIter->second.m_size = bufSize;
454
+
455
+ // add the new free node
456
+ auto newFreePtr = lastElemIter->first + bufSize;
457
+ auto freeNode = m_pointerMap.emplace(newFreePtr, p2).first;
458
+ m_freeList.emplace(freeNode);
459
+ }
460
+
461
+ retVal = lastElemIter->first;
462
+ } else {
463
+ size_t lastSize = lastElemIter->second.m_size;
464
+ retVal = lastElemIter->first + lastSize;
465
+ m_pointerMap.emplace(retVal, p);
466
+ }
467
+ return retVal;
468
+ }
469
+
470
+ /**
471
+ * Compare two iterators to pointer map entries according to
472
+ * the size of the allocation on the device.
473
+ */
474
+ struct SortBySize {
475
+ bool operator()(typename pointerMap_t::iterator a,
476
+ typename pointerMap_t::iterator b) const {
477
+ return ((a->first < b->first) && (a->second <= b->second)) ||
478
+ ((a->first < b->first) && (b->second <= a->second));
479
+ }
480
+ };
481
+
482
+ /* Maps the pointer addresses to buffer and size pairs.
483
+ */
484
+ pointerMap_t m_pointerMap;
485
+
486
+ /* List of free nodes available for re-using
487
+ */
488
+ std::set<typename pointerMap_t::iterator, SortBySize> m_freeList;
489
+
490
+ /* Base address used when issuing the first virtual pointer, allows users
491
+ * to specify alignment. Cannot be zero. */
492
+ std::intptr_t m_baseAddress;
493
+ };
494
+
495
+ /* remove_pointer.
496
+ * Removes the given pointer from the map.
497
+ * The pointer is allowed to be reused only if ReUse if true.
498
+ */
499
+ template <>
500
+ inline void PointerMapper::remove_pointer<false>(const virtual_pointer_t ptr) {
501
+ if (is_nullptr(ptr)) {
502
+ return;
503
+ }
504
+ m_pointerMap.erase(this->get_node(ptr));
505
+ }
506
+
507
+ /**
508
+ * Malloc-like interface to the pointer-mapper.
509
+ * Given a size, creates a byte-typed buffer and returns a
510
+ * fake pointer to keep track of it.
511
+ * \param size Size in bytes of the desired allocation
512
+ * \throw cl::sycl::exception if error while creating the buffer
513
+ */
514
+ inline void *SYCLmalloc(size_t size, PointerMapper &pMap) {
515
+ if (size == 0) {
516
+ return nullptr;
517
+ }
518
+ // Create a generic buffer of the given size
519
+ using buffer_t = cl::sycl::buffer<buffer_data_type_t, 1>;
520
+ auto thePointer = pMap.add_pointer(buffer_t(cl::sycl::range<1>{size}));
521
+ // Store the buffer on the global list
522
+ return static_cast<void *>(thePointer);
523
+ }
524
+
525
+ /**
526
+ * Free-like interface to the pointer mapper.
527
+ * Given a fake-pointer created with the virtual-pointer malloc,
528
+ * destroys the buffer and remove it from the list.
529
+ * If ReUse is false, the pointer is not added to the freeList,
530
+ * it should be false only for sub-buffers.
531
+ */
532
+ template <bool ReUse = true, typename PointerMapper>
533
+ inline void SYCLfree(void *ptr, PointerMapper &pMap) {
534
+ pMap.template remove_pointer<ReUse>(ptr);
535
+ }
536
+
537
+ /**
538
+ * Clear all the memory allocated by SYCL.
539
+ */
540
+ template <typename PointerMapper>
541
+ inline void SYCLfreeAll(PointerMapper &pMap) {
542
+ pMap.clear();
543
+ }
544
+
545
+ template <cl::sycl::access::mode AcMd, typename T>
546
+ struct RangeAccess {
547
+ static const auto global_access = cl::sycl::access::target::global_buffer;
548
+ static const auto is_place_holder = cl::sycl::access::placeholder::true_t;
549
+ typedef T scalar_t;
550
+ typedef scalar_t &ref_t;
551
+ typedef typename cl::sycl::global_ptr<scalar_t>::pointer_t ptr_t;
552
+
553
+ // the accessor type does not necessarily the same as T
554
+ typedef cl::sycl::accessor<scalar_t, 1, AcMd, global_access, is_place_holder>
555
+ accessor;
556
+
557
+ typedef RangeAccess<AcMd, T> self_t;
558
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE RangeAccess(accessor access,
559
+ size_t offset,
560
+ std::intptr_t virtual_ptr)
561
+ : access_(access), offset_(offset), virtual_ptr_(virtual_ptr) {}
562
+
563
+ RangeAccess(cl::sycl::buffer<scalar_t, 1> buff =
564
+ cl::sycl::buffer<scalar_t, 1>(cl::sycl::range<1>(1)))
565
+ : access_{accessor{buff}}, offset_(0), virtual_ptr_(-1) {}
566
+
567
+ // This should be only used for null constructor on the host side
568
+ RangeAccess(std::nullptr_t) : RangeAccess() {}
569
+ // This template parameter must be removed and scalar_t should be replaced
570
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE ptr_t get_pointer() const {
571
+ return (access_.get_pointer().get() + offset_);
572
+ }
573
+ template <typename Index>
574
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE self_t &operator+=(Index offset) {
575
+ offset_ += (offset);
576
+ return *this;
577
+ }
578
+ template <typename Index>
579
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE self_t operator+(Index offset) const {
580
+ return self_t(access_, offset_ + offset, virtual_ptr_);
581
+ }
582
+ template <typename Index>
583
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE self_t operator-(Index offset) const {
584
+ return self_t(access_, offset_ - offset, virtual_ptr_);
585
+ }
586
+ template <typename Index>
587
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE self_t &operator-=(Index offset) {
588
+ offset_ -= offset;
589
+ return *this;
590
+ }
591
+
592
+ // THIS IS FOR NULL COMPARISON ONLY
593
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE friend bool operator==(
594
+ const RangeAccess &lhs, std::nullptr_t) {
595
+ return ((lhs.virtual_ptr_ == -1));
596
+ }
597
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE friend bool operator!=(
598
+ const RangeAccess &lhs, std::nullptr_t i) {
599
+ return !(lhs == i);
600
+ }
601
+
602
+ // THIS IS FOR NULL COMPARISON ONLY
603
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE friend bool operator==(
604
+ std::nullptr_t, const RangeAccess &rhs) {
605
+ return ((rhs.virtual_ptr_ == -1));
606
+ }
607
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE friend bool operator!=(
608
+ std::nullptr_t i, const RangeAccess &rhs) {
609
+ return !(i == rhs);
610
+ }
611
+ // Prefix operator (Increment and return value)
612
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE self_t &operator++() {
613
+ offset_++;
614
+ return (*this);
615
+ }
616
+
617
+ // Postfix operator (Return value and increment)
618
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE self_t operator++(int i) {
619
+ EIGEN_UNUSED_VARIABLE(i);
620
+ self_t temp_iterator(*this);
621
+ offset_++;
622
+ return temp_iterator;
623
+ }
624
+
625
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE std::ptrdiff_t get_size() const {
626
+ return (access_.get_count() - offset_);
627
+ }
628
+
629
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE std::ptrdiff_t get_offset() const {
630
+ return offset_;
631
+ }
632
+
633
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE void set_offset(std::ptrdiff_t offset) {
634
+ offset_ = offset;
635
+ }
636
+
637
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE ref_t operator*() const {
638
+ return *get_pointer();
639
+ }
640
+
641
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE ref_t operator*() {
642
+ return *get_pointer();
643
+ }
644
+
645
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE ptr_t operator->() = delete;
646
+
647
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE ref_t operator[](int x) {
648
+ return *(get_pointer() + x);
649
+ }
650
+
651
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE ref_t operator[](int x) const {
652
+ return *(get_pointer() + x);
653
+ }
654
+
655
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE scalar_t *get_virtual_pointer() const {
656
+ return reinterpret_cast<scalar_t *>(virtual_ptr_ +
657
+ (offset_ * sizeof(scalar_t)));
658
+ }
659
+
660
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE explicit operator bool() const {
661
+ return (virtual_ptr_ != -1);
662
+ }
663
+
664
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE operator RangeAccess<AcMd, const T>() {
665
+ return RangeAccess<AcMd, const T>(access_, offset_, virtual_ptr_);
666
+ }
667
+
668
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
669
+ operator RangeAccess<AcMd, const T>() const {
670
+ return RangeAccess<AcMd, const T>(access_, offset_, virtual_ptr_);
671
+ }
672
+ // binding placeholder accessors to a command group handler for SYCL
673
+ EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE void bind(
674
+ cl::sycl::handler &cgh) const {
675
+ cgh.require(access_);
676
+ }
677
+
678
+ private:
679
+ accessor access_;
680
+ size_t offset_;
681
+ std::intptr_t virtual_ptr_; // the location of the buffer in the map
682
+ };
683
+
684
+ template <cl::sycl::access::mode AcMd, typename T>
685
+ struct RangeAccess<AcMd, const T> : RangeAccess<AcMd, T> {
686
+ typedef RangeAccess<AcMd, T> Base;
687
+ using Base::Base;
688
+ };
689
+
690
+ } // namespace internal
691
+ } // namespace TensorSycl
692
+ } // namespace Eigen
693
+
694
+ #endif // EIGEN_CXX11_TENSOR_TENSOR_SYCL_STORAGE_MEMORY_H