@smake/eigen 1.1.0 → 1.1.1

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (431) hide show
  1. package/README.md +1 -1
  2. package/eigen/Eigen/AccelerateSupport +52 -0
  3. package/eigen/Eigen/Cholesky +18 -20
  4. package/eigen/Eigen/CholmodSupport +28 -28
  5. package/eigen/Eigen/Core +187 -120
  6. package/eigen/Eigen/Eigenvalues +16 -13
  7. package/eigen/Eigen/Geometry +18 -18
  8. package/eigen/Eigen/Householder +9 -7
  9. package/eigen/Eigen/IterativeLinearSolvers +8 -4
  10. package/eigen/Eigen/Jacobi +14 -13
  11. package/eigen/Eigen/KLUSupport +23 -21
  12. package/eigen/Eigen/LU +15 -16
  13. package/eigen/Eigen/MetisSupport +12 -12
  14. package/eigen/Eigen/OrderingMethods +54 -51
  15. package/eigen/Eigen/PaStiXSupport +23 -21
  16. package/eigen/Eigen/PardisoSupport +17 -14
  17. package/eigen/Eigen/QR +18 -20
  18. package/eigen/Eigen/QtAlignedMalloc +5 -12
  19. package/eigen/Eigen/SPQRSupport +21 -14
  20. package/eigen/Eigen/SVD +23 -17
  21. package/eigen/Eigen/Sparse +1 -2
  22. package/eigen/Eigen/SparseCholesky +18 -15
  23. package/eigen/Eigen/SparseCore +18 -17
  24. package/eigen/Eigen/SparseLU +9 -9
  25. package/eigen/Eigen/SparseQR +16 -14
  26. package/eigen/Eigen/StdDeque +5 -2
  27. package/eigen/Eigen/StdList +5 -2
  28. package/eigen/Eigen/StdVector +5 -2
  29. package/eigen/Eigen/SuperLUSupport +30 -24
  30. package/eigen/Eigen/ThreadPool +80 -0
  31. package/eigen/Eigen/UmfPackSupport +19 -17
  32. package/eigen/Eigen/Version +14 -0
  33. package/eigen/Eigen/src/AccelerateSupport/AccelerateSupport.h +423 -0
  34. package/eigen/Eigen/src/AccelerateSupport/InternalHeaderCheck.h +3 -0
  35. package/eigen/Eigen/src/Cholesky/InternalHeaderCheck.h +3 -0
  36. package/eigen/Eigen/src/Cholesky/LDLT.h +366 -405
  37. package/eigen/Eigen/src/Cholesky/LLT.h +323 -367
  38. package/eigen/Eigen/src/Cholesky/LLT_LAPACKE.h +81 -56
  39. package/eigen/Eigen/src/CholmodSupport/CholmodSupport.h +585 -529
  40. package/eigen/Eigen/src/CholmodSupport/InternalHeaderCheck.h +3 -0
  41. package/eigen/Eigen/src/Core/ArithmeticSequence.h +143 -317
  42. package/eigen/Eigen/src/Core/Array.h +329 -370
  43. package/eigen/Eigen/src/Core/ArrayBase.h +190 -203
  44. package/eigen/Eigen/src/Core/ArrayWrapper.h +126 -170
  45. package/eigen/Eigen/src/Core/Assign.h +30 -40
  46. package/eigen/Eigen/src/Core/AssignEvaluator.h +651 -604
  47. package/eigen/Eigen/src/Core/Assign_MKL.h +125 -120
  48. package/eigen/Eigen/src/Core/BandMatrix.h +267 -282
  49. package/eigen/Eigen/src/Core/Block.h +371 -390
  50. package/eigen/Eigen/src/Core/CommaInitializer.h +85 -100
  51. package/eigen/Eigen/src/Core/ConditionEstimator.h +51 -53
  52. package/eigen/Eigen/src/Core/CoreEvaluators.h +1214 -937
  53. package/eigen/Eigen/src/Core/CoreIterators.h +72 -63
  54. package/eigen/Eigen/src/Core/CwiseBinaryOp.h +112 -129
  55. package/eigen/Eigen/src/Core/CwiseNullaryOp.h +676 -702
  56. package/eigen/Eigen/src/Core/CwiseTernaryOp.h +77 -103
  57. package/eigen/Eigen/src/Core/CwiseUnaryOp.h +55 -67
  58. package/eigen/Eigen/src/Core/CwiseUnaryView.h +127 -92
  59. package/eigen/Eigen/src/Core/DenseBase.h +630 -658
  60. package/eigen/Eigen/src/Core/DenseCoeffsBase.h +511 -628
  61. package/eigen/Eigen/src/Core/DenseStorage.h +511 -590
  62. package/eigen/Eigen/src/Core/DeviceWrapper.h +153 -0
  63. package/eigen/Eigen/src/Core/Diagonal.h +168 -207
  64. package/eigen/Eigen/src/Core/DiagonalMatrix.h +346 -317
  65. package/eigen/Eigen/src/Core/DiagonalProduct.h +12 -10
  66. package/eigen/Eigen/src/Core/Dot.h +167 -217
  67. package/eigen/Eigen/src/Core/EigenBase.h +74 -85
  68. package/eigen/Eigen/src/Core/Fill.h +138 -0
  69. package/eigen/Eigen/src/Core/FindCoeff.h +464 -0
  70. package/eigen/Eigen/src/Core/ForceAlignedAccess.h +90 -113
  71. package/eigen/Eigen/src/Core/Fuzzy.h +82 -105
  72. package/eigen/Eigen/src/Core/GeneralProduct.h +315 -261
  73. package/eigen/Eigen/src/Core/GenericPacketMath.h +1182 -520
  74. package/eigen/Eigen/src/Core/GlobalFunctions.h +193 -157
  75. package/eigen/Eigen/src/Core/IO.h +131 -156
  76. package/eigen/Eigen/src/Core/IndexedView.h +209 -125
  77. package/eigen/Eigen/src/Core/InnerProduct.h +260 -0
  78. package/eigen/Eigen/src/Core/InternalHeaderCheck.h +3 -0
  79. package/eigen/Eigen/src/Core/Inverse.h +50 -59
  80. package/eigen/Eigen/src/Core/Map.h +123 -141
  81. package/eigen/Eigen/src/Core/MapBase.h +255 -282
  82. package/eigen/Eigen/src/Core/MathFunctions.h +1247 -1201
  83. package/eigen/Eigen/src/Core/MathFunctionsImpl.h +162 -99
  84. package/eigen/Eigen/src/Core/Matrix.h +463 -494
  85. package/eigen/Eigen/src/Core/MatrixBase.h +468 -470
  86. package/eigen/Eigen/src/Core/NestByValue.h +58 -52
  87. package/eigen/Eigen/src/Core/NoAlias.h +79 -86
  88. package/eigen/Eigen/src/Core/NumTraits.h +206 -206
  89. package/eigen/Eigen/src/Core/PartialReduxEvaluator.h +163 -142
  90. package/eigen/Eigen/src/Core/PermutationMatrix.h +461 -511
  91. package/eigen/Eigen/src/Core/PlainObjectBase.h +858 -972
  92. package/eigen/Eigen/src/Core/Product.h +246 -130
  93. package/eigen/Eigen/src/Core/ProductEvaluators.h +779 -671
  94. package/eigen/Eigen/src/Core/Random.h +153 -164
  95. package/eigen/Eigen/src/Core/RandomImpl.h +262 -0
  96. package/eigen/Eigen/src/Core/RealView.h +250 -0
  97. package/eigen/Eigen/src/Core/Redux.h +334 -314
  98. package/eigen/Eigen/src/Core/Ref.h +259 -257
  99. package/eigen/Eigen/src/Core/Replicate.h +92 -104
  100. package/eigen/Eigen/src/Core/Reshaped.h +215 -271
  101. package/eigen/Eigen/src/Core/ReturnByValue.h +47 -55
  102. package/eigen/Eigen/src/Core/Reverse.h +133 -148
  103. package/eigen/Eigen/src/Core/Select.h +68 -140
  104. package/eigen/Eigen/src/Core/SelfAdjointView.h +254 -290
  105. package/eigen/Eigen/src/Core/SelfCwiseBinaryOp.h +23 -20
  106. package/eigen/Eigen/src/Core/SkewSymmetricMatrix3.h +382 -0
  107. package/eigen/Eigen/src/Core/Solve.h +88 -102
  108. package/eigen/Eigen/src/Core/SolveTriangular.h +126 -124
  109. package/eigen/Eigen/src/Core/SolverBase.h +132 -133
  110. package/eigen/Eigen/src/Core/StableNorm.h +113 -147
  111. package/eigen/Eigen/src/Core/StlIterators.h +404 -248
  112. package/eigen/Eigen/src/Core/Stride.h +90 -92
  113. package/eigen/Eigen/src/Core/Swap.h +70 -39
  114. package/eigen/Eigen/src/Core/Transpose.h +258 -295
  115. package/eigen/Eigen/src/Core/Transpositions.h +270 -333
  116. package/eigen/Eigen/src/Core/TriangularMatrix.h +642 -743
  117. package/eigen/Eigen/src/Core/VectorBlock.h +59 -72
  118. package/eigen/Eigen/src/Core/VectorwiseOp.h +653 -704
  119. package/eigen/Eigen/src/Core/Visitor.h +464 -308
  120. package/eigen/Eigen/src/Core/arch/AVX/Complex.h +380 -187
  121. package/eigen/Eigen/src/Core/arch/AVX/MathFunctions.h +65 -163
  122. package/eigen/Eigen/src/Core/arch/AVX/PacketMath.h +2145 -638
  123. package/eigen/Eigen/src/Core/arch/AVX/Reductions.h +353 -0
  124. package/eigen/Eigen/src/Core/arch/AVX/TypeCasting.h +253 -60
  125. package/eigen/Eigen/src/Core/arch/AVX512/Complex.h +278 -228
  126. package/eigen/Eigen/src/Core/arch/AVX512/GemmKernel.h +1245 -0
  127. package/eigen/Eigen/src/Core/arch/AVX512/MathFunctions.h +48 -269
  128. package/eigen/Eigen/src/Core/arch/AVX512/MathFunctionsFP16.h +75 -0
  129. package/eigen/Eigen/src/Core/arch/AVX512/PacketMath.h +1597 -754
  130. package/eigen/Eigen/src/Core/arch/AVX512/PacketMathFP16.h +1413 -0
  131. package/eigen/Eigen/src/Core/arch/AVX512/Reductions.h +297 -0
  132. package/eigen/Eigen/src/Core/arch/AVX512/TrsmKernel.h +1167 -0
  133. package/eigen/Eigen/src/Core/arch/AVX512/TrsmUnrolls.inc +1219 -0
  134. package/eigen/Eigen/src/Core/arch/AVX512/TypeCasting.h +229 -41
  135. package/eigen/Eigen/src/Core/arch/AVX512/TypeCastingFP16.h +130 -0
  136. package/eigen/Eigen/src/Core/arch/AltiVec/Complex.h +420 -184
  137. package/eigen/Eigen/src/Core/arch/AltiVec/MathFunctions.h +40 -49
  138. package/eigen/Eigen/src/Core/arch/AltiVec/MatrixProduct.h +2962 -2213
  139. package/eigen/Eigen/src/Core/arch/AltiVec/MatrixProductCommon.h +196 -212
  140. package/eigen/Eigen/src/Core/arch/AltiVec/MatrixProductMMA.h +713 -441
  141. package/eigen/Eigen/src/Core/arch/AltiVec/MatrixProductMMAbfloat16.h +742 -0
  142. package/eigen/Eigen/src/Core/arch/AltiVec/MatrixVectorProduct.inc +2818 -0
  143. package/eigen/Eigen/src/Core/arch/AltiVec/PacketMath.h +2380 -1362
  144. package/eigen/Eigen/src/Core/arch/AltiVec/TypeCasting.h +153 -0
  145. package/eigen/Eigen/src/Core/arch/Default/BFloat16.h +390 -224
  146. package/eigen/Eigen/src/Core/arch/Default/ConjHelper.h +78 -67
  147. package/eigen/Eigen/src/Core/arch/Default/GenericPacketMathFunctions.h +1784 -799
  148. package/eigen/Eigen/src/Core/arch/Default/GenericPacketMathFunctionsFwd.h +167 -50
  149. package/eigen/Eigen/src/Core/arch/Default/Half.h +528 -379
  150. package/eigen/Eigen/src/Core/arch/Default/Settings.h +10 -12
  151. package/eigen/Eigen/src/Core/arch/GPU/Complex.h +244 -0
  152. package/eigen/Eigen/src/Core/arch/GPU/MathFunctions.h +41 -40
  153. package/eigen/Eigen/src/Core/arch/GPU/PacketMath.h +550 -523
  154. package/eigen/Eigen/src/Core/arch/GPU/Tuple.h +268 -0
  155. package/eigen/Eigen/src/Core/arch/GPU/TypeCasting.h +27 -30
  156. package/eigen/Eigen/src/Core/arch/HIP/hcc/math_constants.h +8 -8
  157. package/eigen/Eigen/src/Core/arch/HVX/PacketMath.h +1088 -0
  158. package/eigen/Eigen/src/Core/arch/LSX/Complex.h +520 -0
  159. package/eigen/Eigen/src/Core/arch/LSX/GeneralBlockPanelKernel.h +23 -0
  160. package/eigen/Eigen/src/Core/arch/LSX/MathFunctions.h +43 -0
  161. package/eigen/Eigen/src/Core/arch/LSX/PacketMath.h +2866 -0
  162. package/eigen/Eigen/src/Core/arch/LSX/TypeCasting.h +526 -0
  163. package/eigen/Eigen/src/Core/arch/MSA/Complex.h +54 -82
  164. package/eigen/Eigen/src/Core/arch/MSA/MathFunctions.h +84 -92
  165. package/eigen/Eigen/src/Core/arch/MSA/PacketMath.h +51 -47
  166. package/eigen/Eigen/src/Core/arch/NEON/Complex.h +454 -306
  167. package/eigen/Eigen/src/Core/arch/NEON/GeneralBlockPanelKernel.h +175 -115
  168. package/eigen/Eigen/src/Core/arch/NEON/MathFunctions.h +23 -30
  169. package/eigen/Eigen/src/Core/arch/NEON/PacketMath.h +4366 -2857
  170. package/eigen/Eigen/src/Core/arch/NEON/TypeCasting.h +616 -393
  171. package/eigen/Eigen/src/Core/arch/NEON/UnaryFunctors.h +57 -0
  172. package/eigen/Eigen/src/Core/arch/SSE/Complex.h +350 -198
  173. package/eigen/Eigen/src/Core/arch/SSE/MathFunctions.h +38 -149
  174. package/eigen/Eigen/src/Core/arch/SSE/PacketMath.h +1791 -912
  175. package/eigen/Eigen/src/Core/arch/SSE/Reductions.h +324 -0
  176. package/eigen/Eigen/src/Core/arch/SSE/TypeCasting.h +128 -40
  177. package/eigen/Eigen/src/Core/arch/SVE/MathFunctions.h +10 -6
  178. package/eigen/Eigen/src/Core/arch/SVE/PacketMath.h +156 -234
  179. package/eigen/Eigen/src/Core/arch/SVE/TypeCasting.h +6 -3
  180. package/eigen/Eigen/src/Core/arch/SYCL/InteropHeaders.h +27 -32
  181. package/eigen/Eigen/src/Core/arch/SYCL/MathFunctions.h +119 -117
  182. package/eigen/Eigen/src/Core/arch/SYCL/PacketMath.h +325 -419
  183. package/eigen/Eigen/src/Core/arch/SYCL/TypeCasting.h +15 -17
  184. package/eigen/Eigen/src/Core/arch/ZVector/Complex.h +325 -181
  185. package/eigen/Eigen/src/Core/arch/ZVector/MathFunctions.h +94 -83
  186. package/eigen/Eigen/src/Core/arch/ZVector/PacketMath.h +811 -458
  187. package/eigen/Eigen/src/Core/functors/AssignmentFunctors.h +121 -124
  188. package/eigen/Eigen/src/Core/functors/BinaryFunctors.h +576 -370
  189. package/eigen/Eigen/src/Core/functors/NullaryFunctors.h +194 -109
  190. package/eigen/Eigen/src/Core/functors/StlFunctors.h +95 -112
  191. package/eigen/Eigen/src/Core/functors/TernaryFunctors.h +34 -7
  192. package/eigen/Eigen/src/Core/functors/UnaryFunctors.h +1038 -749
  193. package/eigen/Eigen/src/Core/products/GeneralBlockPanelKernel.h +1883 -1375
  194. package/eigen/Eigen/src/Core/products/GeneralMatrixMatrix.h +312 -370
  195. package/eigen/Eigen/src/Core/products/GeneralMatrixMatrixTriangular.h +189 -176
  196. package/eigen/Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h +84 -81
  197. package/eigen/Eigen/src/Core/products/GeneralMatrixMatrix_BLAS.h +154 -73
  198. package/eigen/Eigen/src/Core/products/GeneralMatrixVector.h +292 -337
  199. package/eigen/Eigen/src/Core/products/GeneralMatrixVector_BLAS.h +80 -77
  200. package/eigen/Eigen/src/Core/products/Parallelizer.h +207 -105
  201. package/eigen/Eigen/src/Core/products/SelfadjointMatrixMatrix.h +327 -388
  202. package/eigen/Eigen/src/Core/products/SelfadjointMatrixMatrix_BLAS.h +206 -224
  203. package/eigen/Eigen/src/Core/products/SelfadjointMatrixVector.h +138 -147
  204. package/eigen/Eigen/src/Core/products/SelfadjointMatrixVector_BLAS.h +58 -61
  205. package/eigen/Eigen/src/Core/products/SelfadjointProduct.h +71 -71
  206. package/eigen/Eigen/src/Core/products/SelfadjointRank2Update.h +48 -47
  207. package/eigen/Eigen/src/Core/products/TriangularMatrixMatrix.h +294 -369
  208. package/eigen/Eigen/src/Core/products/TriangularMatrixMatrix_BLAS.h +246 -238
  209. package/eigen/Eigen/src/Core/products/TriangularMatrixVector.h +244 -247
  210. package/eigen/Eigen/src/Core/products/TriangularMatrixVector_BLAS.h +212 -192
  211. package/eigen/Eigen/src/Core/products/TriangularSolverMatrix.h +328 -277
  212. package/eigen/Eigen/src/Core/products/TriangularSolverMatrix_BLAS.h +108 -109
  213. package/eigen/Eigen/src/Core/products/TriangularSolverVector.h +68 -94
  214. package/eigen/Eigen/src/Core/util/Assert.h +158 -0
  215. package/eigen/Eigen/src/Core/util/BlasUtil.h +342 -303
  216. package/eigen/Eigen/src/Core/util/ConfigureVectorization.h +348 -317
  217. package/eigen/Eigen/src/Core/util/Constants.h +297 -262
  218. package/eigen/Eigen/src/Core/util/DisableStupidWarnings.h +130 -90
  219. package/eigen/Eigen/src/Core/util/EmulateArray.h +270 -0
  220. package/eigen/Eigen/src/Core/util/ForwardDeclarations.h +449 -247
  221. package/eigen/Eigen/src/Core/util/GpuHipCudaDefines.inc +101 -0
  222. package/eigen/Eigen/src/Core/util/GpuHipCudaUndefines.inc +45 -0
  223. package/eigen/Eigen/src/Core/util/IndexedViewHelper.h +417 -116
  224. package/eigen/Eigen/src/Core/util/IntegralConstant.h +211 -204
  225. package/eigen/Eigen/src/Core/util/MKL_support.h +39 -37
  226. package/eigen/Eigen/src/Core/util/Macros.h +655 -773
  227. package/eigen/Eigen/src/Core/util/MaxSizeVector.h +139 -0
  228. package/eigen/Eigen/src/Core/util/Memory.h +970 -748
  229. package/eigen/Eigen/src/Core/util/Meta.h +581 -633
  230. package/eigen/Eigen/src/Core/util/MoreMeta.h +638 -0
  231. package/eigen/Eigen/src/Core/util/ReenableStupidWarnings.h +32 -19
  232. package/eigen/Eigen/src/Core/util/ReshapedHelper.h +17 -17
  233. package/eigen/Eigen/src/Core/util/Serializer.h +209 -0
  234. package/eigen/Eigen/src/Core/util/StaticAssert.h +50 -166
  235. package/eigen/Eigen/src/Core/util/SymbolicIndex.h +377 -225
  236. package/eigen/Eigen/src/Core/util/XprHelper.h +784 -547
  237. package/eigen/Eigen/src/Eigenvalues/ComplexEigenSolver.h +246 -277
  238. package/eigen/Eigen/src/Eigenvalues/ComplexSchur.h +299 -319
  239. package/eigen/Eigen/src/Eigenvalues/ComplexSchur_LAPACKE.h +52 -48
  240. package/eigen/Eigen/src/Eigenvalues/EigenSolver.h +413 -456
  241. package/eigen/Eigen/src/Eigenvalues/GeneralizedEigenSolver.h +309 -325
  242. package/eigen/Eigen/src/Eigenvalues/GeneralizedSelfAdjointEigenSolver.h +157 -171
  243. package/eigen/Eigen/src/Eigenvalues/HessenbergDecomposition.h +292 -310
  244. package/eigen/Eigen/src/Eigenvalues/InternalHeaderCheck.h +3 -0
  245. package/eigen/Eigen/src/Eigenvalues/MatrixBaseEigenvalues.h +89 -105
  246. package/eigen/Eigen/src/Eigenvalues/RealQZ.h +537 -607
  247. package/eigen/Eigen/src/Eigenvalues/RealSchur.h +342 -381
  248. package/eigen/Eigen/src/Eigenvalues/RealSchur_LAPACKE.h +41 -35
  249. package/eigen/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h +541 -595
  250. package/eigen/Eigen/src/Eigenvalues/SelfAdjointEigenSolver_LAPACKE.h +47 -44
  251. package/eigen/Eigen/src/Eigenvalues/Tridiagonalization.h +430 -462
  252. package/eigen/Eigen/src/Geometry/AlignedBox.h +226 -227
  253. package/eigen/Eigen/src/Geometry/AngleAxis.h +131 -133
  254. package/eigen/Eigen/src/Geometry/EulerAngles.h +163 -74
  255. package/eigen/Eigen/src/Geometry/Homogeneous.h +285 -333
  256. package/eigen/Eigen/src/Geometry/Hyperplane.h +151 -160
  257. package/eigen/Eigen/src/Geometry/InternalHeaderCheck.h +3 -0
  258. package/eigen/Eigen/src/Geometry/OrthoMethods.h +168 -146
  259. package/eigen/Eigen/src/Geometry/ParametrizedLine.h +127 -127
  260. package/eigen/Eigen/src/Geometry/Quaternion.h +566 -506
  261. package/eigen/Eigen/src/Geometry/Rotation2D.h +107 -105
  262. package/eigen/Eigen/src/Geometry/RotationBase.h +148 -145
  263. package/eigen/Eigen/src/Geometry/Scaling.h +113 -106
  264. package/eigen/Eigen/src/Geometry/Transform.h +858 -936
  265. package/eigen/Eigen/src/Geometry/Translation.h +94 -92
  266. package/eigen/Eigen/src/Geometry/Umeyama.h +79 -84
  267. package/eigen/Eigen/src/Geometry/arch/Geometry_SIMD.h +90 -104
  268. package/eigen/Eigen/src/Householder/BlockHouseholder.h +51 -46
  269. package/eigen/Eigen/src/Householder/Householder.h +102 -124
  270. package/eigen/Eigen/src/Householder/HouseholderSequence.h +412 -453
  271. package/eigen/Eigen/src/Householder/InternalHeaderCheck.h +3 -0
  272. package/eigen/Eigen/src/IterativeLinearSolvers/BasicPreconditioners.h +149 -162
  273. package/eigen/Eigen/src/IterativeLinearSolvers/BiCGSTAB.h +124 -119
  274. package/eigen/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h +92 -104
  275. package/eigen/Eigen/src/IterativeLinearSolvers/IncompleteCholesky.h +251 -243
  276. package/eigen/Eigen/src/IterativeLinearSolvers/IncompleteLUT.h +224 -228
  277. package/eigen/Eigen/src/IterativeLinearSolvers/InternalHeaderCheck.h +3 -0
  278. package/eigen/Eigen/src/IterativeLinearSolvers/IterativeSolverBase.h +178 -227
  279. package/eigen/Eigen/src/IterativeLinearSolvers/LeastSquareConjugateGradient.h +79 -84
  280. package/eigen/Eigen/src/IterativeLinearSolvers/SolveWithGuess.h +54 -60
  281. package/eigen/Eigen/src/Jacobi/InternalHeaderCheck.h +3 -0
  282. package/eigen/Eigen/src/Jacobi/Jacobi.h +252 -308
  283. package/eigen/Eigen/src/KLUSupport/InternalHeaderCheck.h +3 -0
  284. package/eigen/Eigen/src/KLUSupport/KLUSupport.h +208 -227
  285. package/eigen/Eigen/src/LU/Determinant.h +50 -69
  286. package/eigen/Eigen/src/LU/FullPivLU.h +545 -596
  287. package/eigen/Eigen/src/LU/InternalHeaderCheck.h +3 -0
  288. package/eigen/Eigen/src/LU/InverseImpl.h +206 -285
  289. package/eigen/Eigen/src/LU/PartialPivLU.h +390 -428
  290. package/eigen/Eigen/src/LU/PartialPivLU_LAPACKE.h +54 -40
  291. package/eigen/Eigen/src/LU/arch/InverseSize4.h +72 -70
  292. package/eigen/Eigen/src/MetisSupport/InternalHeaderCheck.h +3 -0
  293. package/eigen/Eigen/src/MetisSupport/MetisSupport.h +81 -93
  294. package/eigen/Eigen/src/OrderingMethods/Amd.h +243 -265
  295. package/eigen/Eigen/src/OrderingMethods/Eigen_Colamd.h +831 -1004
  296. package/eigen/Eigen/src/OrderingMethods/InternalHeaderCheck.h +3 -0
  297. package/eigen/Eigen/src/OrderingMethods/Ordering.h +112 -119
  298. package/eigen/Eigen/src/PaStiXSupport/InternalHeaderCheck.h +3 -0
  299. package/eigen/Eigen/src/PaStiXSupport/PaStiXSupport.h +524 -570
  300. package/eigen/Eigen/src/PardisoSupport/InternalHeaderCheck.h +3 -0
  301. package/eigen/Eigen/src/PardisoSupport/PardisoSupport.h +385 -430
  302. package/eigen/Eigen/src/QR/ColPivHouseholderQR.h +479 -479
  303. package/eigen/Eigen/src/QR/ColPivHouseholderQR_LAPACKE.h +120 -56
  304. package/eigen/Eigen/src/QR/CompleteOrthogonalDecomposition.h +166 -153
  305. package/eigen/Eigen/src/QR/FullPivHouseholderQR.h +495 -475
  306. package/eigen/Eigen/src/QR/HouseholderQR.h +394 -285
  307. package/eigen/Eigen/src/QR/HouseholderQR_LAPACKE.h +32 -23
  308. package/eigen/Eigen/src/QR/InternalHeaderCheck.h +3 -0
  309. package/eigen/Eigen/src/SPQRSupport/InternalHeaderCheck.h +3 -0
  310. package/eigen/Eigen/src/SPQRSupport/SuiteSparseQRSupport.h +244 -264
  311. package/eigen/Eigen/src/SVD/BDCSVD.h +817 -713
  312. package/eigen/Eigen/src/SVD/BDCSVD_LAPACKE.h +174 -0
  313. package/eigen/Eigen/src/SVD/InternalHeaderCheck.h +3 -0
  314. package/eigen/Eigen/src/SVD/JacobiSVD.h +577 -543
  315. package/eigen/Eigen/src/SVD/JacobiSVD_LAPACKE.h +85 -49
  316. package/eigen/Eigen/src/SVD/SVDBase.h +242 -182
  317. package/eigen/Eigen/src/SVD/UpperBidiagonalization.h +200 -235
  318. package/eigen/Eigen/src/SparseCholesky/InternalHeaderCheck.h +3 -0
  319. package/eigen/Eigen/src/SparseCholesky/SimplicialCholesky.h +765 -594
  320. package/eigen/Eigen/src/SparseCholesky/SimplicialCholesky_impl.h +308 -94
  321. package/eigen/Eigen/src/SparseCore/AmbiVector.h +202 -251
  322. package/eigen/Eigen/src/SparseCore/CompressedStorage.h +184 -252
  323. package/eigen/Eigen/src/SparseCore/ConservativeSparseSparseProduct.h +134 -178
  324. package/eigen/Eigen/src/SparseCore/InternalHeaderCheck.h +3 -0
  325. package/eigen/Eigen/src/SparseCore/SparseAssign.h +149 -140
  326. package/eigen/Eigen/src/SparseCore/SparseBlock.h +403 -440
  327. package/eigen/Eigen/src/SparseCore/SparseColEtree.h +100 -112
  328. package/eigen/Eigen/src/SparseCore/SparseCompressedBase.h +525 -303
  329. package/eigen/Eigen/src/SparseCore/SparseCwiseBinaryOp.h +555 -339
  330. package/eigen/Eigen/src/SparseCore/SparseCwiseUnaryOp.h +100 -108
  331. package/eigen/Eigen/src/SparseCore/SparseDenseProduct.h +169 -197
  332. package/eigen/Eigen/src/SparseCore/SparseDiagonalProduct.h +71 -71
  333. package/eigen/Eigen/src/SparseCore/SparseDot.h +49 -47
  334. package/eigen/Eigen/src/SparseCore/SparseFuzzy.h +13 -11
  335. package/eigen/Eigen/src/SparseCore/SparseMap.h +243 -253
  336. package/eigen/Eigen/src/SparseCore/SparseMatrix.h +1603 -1245
  337. package/eigen/Eigen/src/SparseCore/SparseMatrixBase.h +403 -350
  338. package/eigen/Eigen/src/SparseCore/SparsePermutation.h +186 -115
  339. package/eigen/Eigen/src/SparseCore/SparseProduct.h +94 -97
  340. package/eigen/Eigen/src/SparseCore/SparseRedux.h +22 -24
  341. package/eigen/Eigen/src/SparseCore/SparseRef.h +268 -295
  342. package/eigen/Eigen/src/SparseCore/SparseSelfAdjointView.h +370 -416
  343. package/eigen/Eigen/src/SparseCore/SparseSolverBase.h +78 -87
  344. package/eigen/Eigen/src/SparseCore/SparseSparseProductWithPruning.h +81 -95
  345. package/eigen/Eigen/src/SparseCore/SparseTranspose.h +62 -71
  346. package/eigen/Eigen/src/SparseCore/SparseTriangularView.h +132 -144
  347. package/eigen/Eigen/src/SparseCore/SparseUtil.h +138 -115
  348. package/eigen/Eigen/src/SparseCore/SparseVector.h +426 -372
  349. package/eigen/Eigen/src/SparseCore/SparseView.h +164 -193
  350. package/eigen/Eigen/src/SparseCore/TriangularSolver.h +129 -170
  351. package/eigen/Eigen/src/SparseLU/InternalHeaderCheck.h +3 -0
  352. package/eigen/Eigen/src/SparseLU/SparseLU.h +756 -710
  353. package/eigen/Eigen/src/SparseLU/SparseLUImpl.h +61 -48
  354. package/eigen/Eigen/src/SparseLU/SparseLU_Memory.h +102 -118
  355. package/eigen/Eigen/src/SparseLU/SparseLU_Structs.h +38 -35
  356. package/eigen/Eigen/src/SparseLU/SparseLU_SupernodalMatrix.h +245 -301
  357. package/eigen/Eigen/src/SparseLU/SparseLU_Utils.h +44 -49
  358. package/eigen/Eigen/src/SparseLU/SparseLU_column_bmod.h +104 -108
  359. package/eigen/Eigen/src/SparseLU/SparseLU_column_dfs.h +89 -100
  360. package/eigen/Eigen/src/SparseLU/SparseLU_copy_to_ucol.h +57 -58
  361. package/eigen/Eigen/src/SparseLU/SparseLU_heap_relax_snode.h +43 -55
  362. package/eigen/Eigen/src/SparseLU/SparseLU_kernel_bmod.h +74 -71
  363. package/eigen/Eigen/src/SparseLU/SparseLU_panel_bmod.h +124 -132
  364. package/eigen/Eigen/src/SparseLU/SparseLU_panel_dfs.h +136 -159
  365. package/eigen/Eigen/src/SparseLU/SparseLU_pivotL.h +51 -52
  366. package/eigen/Eigen/src/SparseLU/SparseLU_pruneL.h +67 -73
  367. package/eigen/Eigen/src/SparseLU/SparseLU_relax_snode.h +24 -26
  368. package/eigen/Eigen/src/SparseQR/InternalHeaderCheck.h +3 -0
  369. package/eigen/Eigen/src/SparseQR/SparseQR.h +450 -502
  370. package/eigen/Eigen/src/StlSupport/StdDeque.h +28 -93
  371. package/eigen/Eigen/src/StlSupport/StdList.h +28 -84
  372. package/eigen/Eigen/src/StlSupport/StdVector.h +28 -108
  373. package/eigen/Eigen/src/StlSupport/details.h +48 -50
  374. package/eigen/Eigen/src/SuperLUSupport/InternalHeaderCheck.h +3 -0
  375. package/eigen/Eigen/src/SuperLUSupport/SuperLUSupport.h +634 -730
  376. package/eigen/Eigen/src/ThreadPool/Barrier.h +70 -0
  377. package/eigen/Eigen/src/ThreadPool/CoreThreadPoolDevice.h +336 -0
  378. package/eigen/Eigen/src/ThreadPool/EventCount.h +241 -0
  379. package/eigen/Eigen/src/ThreadPool/ForkJoin.h +140 -0
  380. package/eigen/Eigen/src/ThreadPool/InternalHeaderCheck.h +4 -0
  381. package/eigen/Eigen/src/ThreadPool/NonBlockingThreadPool.h +587 -0
  382. package/eigen/Eigen/src/ThreadPool/RunQueue.h +230 -0
  383. package/eigen/Eigen/src/ThreadPool/ThreadCancel.h +21 -0
  384. package/eigen/Eigen/src/ThreadPool/ThreadEnvironment.h +43 -0
  385. package/eigen/Eigen/src/ThreadPool/ThreadLocal.h +289 -0
  386. package/eigen/Eigen/src/ThreadPool/ThreadPoolInterface.h +50 -0
  387. package/eigen/Eigen/src/ThreadPool/ThreadYield.h +16 -0
  388. package/eigen/Eigen/src/UmfPackSupport/InternalHeaderCheck.h +3 -0
  389. package/eigen/Eigen/src/UmfPackSupport/UmfPackSupport.h +428 -464
  390. package/eigen/Eigen/src/misc/Image.h +41 -43
  391. package/eigen/Eigen/src/misc/InternalHeaderCheck.h +3 -0
  392. package/eigen/Eigen/src/misc/Kernel.h +39 -41
  393. package/eigen/Eigen/src/misc/RealSvd2x2.h +19 -21
  394. package/eigen/Eigen/src/misc/blas.h +83 -426
  395. package/eigen/Eigen/src/misc/lapacke.h +9972 -16179
  396. package/eigen/Eigen/src/misc/lapacke_helpers.h +163 -0
  397. package/eigen/Eigen/src/misc/lapacke_mangling.h +4 -5
  398. package/eigen/Eigen/src/plugins/ArrayCwiseBinaryOps.inc +344 -0
  399. package/eigen/Eigen/src/plugins/ArrayCwiseUnaryOps.inc +544 -0
  400. package/eigen/Eigen/src/plugins/{BlockMethods.h → BlockMethods.inc} +434 -506
  401. package/eigen/Eigen/src/plugins/CommonCwiseBinaryOps.inc +116 -0
  402. package/eigen/Eigen/src/plugins/{CommonCwiseUnaryOps.h → CommonCwiseUnaryOps.inc} +58 -68
  403. package/eigen/Eigen/src/plugins/IndexedViewMethods.inc +192 -0
  404. package/eigen/Eigen/src/plugins/InternalHeaderCheck.inc +3 -0
  405. package/eigen/Eigen/src/plugins/MatrixCwiseBinaryOps.inc +331 -0
  406. package/eigen/Eigen/src/plugins/MatrixCwiseUnaryOps.inc +118 -0
  407. package/eigen/Eigen/src/plugins/ReshapedMethods.inc +133 -0
  408. package/package.json +1 -1
  409. package/eigen/COPYING.APACHE +0 -203
  410. package/eigen/COPYING.BSD +0 -26
  411. package/eigen/COPYING.GPL +0 -674
  412. package/eigen/COPYING.LGPL +0 -502
  413. package/eigen/COPYING.MINPACK +0 -51
  414. package/eigen/COPYING.MPL2 +0 -373
  415. package/eigen/COPYING.README +0 -18
  416. package/eigen/Eigen/src/Core/BooleanRedux.h +0 -162
  417. package/eigen/Eigen/src/Core/arch/CUDA/Complex.h +0 -258
  418. package/eigen/Eigen/src/Core/arch/Default/TypeCasting.h +0 -120
  419. package/eigen/Eigen/src/Core/arch/SYCL/SyclMemoryModel.h +0 -694
  420. package/eigen/Eigen/src/Core/util/NonMPL2.h +0 -3
  421. package/eigen/Eigen/src/SparseCore/MappedSparseMatrix.h +0 -67
  422. package/eigen/Eigen/src/SparseLU/SparseLU_gemm_kernel.h +0 -280
  423. package/eigen/Eigen/src/misc/lapack.h +0 -152
  424. package/eigen/Eigen/src/plugins/ArrayCwiseBinaryOps.h +0 -358
  425. package/eigen/Eigen/src/plugins/ArrayCwiseUnaryOps.h +0 -696
  426. package/eigen/Eigen/src/plugins/CommonCwiseBinaryOps.h +0 -115
  427. package/eigen/Eigen/src/plugins/IndexedViewMethods.h +0 -262
  428. package/eigen/Eigen/src/plugins/MatrixCwiseBinaryOps.h +0 -152
  429. package/eigen/Eigen/src/plugins/MatrixCwiseUnaryOps.h +0 -95
  430. package/eigen/Eigen/src/plugins/ReshapedMethods.h +0 -149
  431. package/eigen/README.md +0 -5
@@ -1,694 +0,0 @@
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
@@ -1,3 +0,0 @@
1
- #ifdef EIGEN_MPL2_ONLY
2
- #error Including non-MPL2 code in EIGEN_MPL2_ONLY mode
3
- #endif