tomoto 0.2.3 → 0.3.0

Sign up to get free protection for your applications and to get access to all the features.
Files changed (347) hide show
  1. checksums.yaml +4 -4
  2. data/CHANGELOG.md +6 -0
  3. data/README.md +8 -10
  4. data/ext/tomoto/extconf.rb +6 -2
  5. data/ext/tomoto/{ext.cpp → tomoto.cpp} +1 -1
  6. data/lib/tomoto/version.rb +1 -1
  7. data/lib/tomoto.rb +5 -1
  8. data/vendor/EigenRand/EigenRand/Core.h +10 -10
  9. data/vendor/EigenRand/EigenRand/Dists/Basic.h +208 -9
  10. data/vendor/EigenRand/EigenRand/Dists/Discrete.h +52 -31
  11. data/vendor/EigenRand/EigenRand/Dists/GammaPoisson.h +9 -8
  12. data/vendor/EigenRand/EigenRand/Dists/NormalExp.h +28 -21
  13. data/vendor/EigenRand/EigenRand/EigenRand +11 -6
  14. data/vendor/EigenRand/EigenRand/Macro.h +13 -7
  15. data/vendor/EigenRand/EigenRand/MorePacketMath.h +348 -740
  16. data/vendor/EigenRand/EigenRand/MvDists/Multinomial.h +5 -3
  17. data/vendor/EigenRand/EigenRand/MvDists/MvNormal.h +9 -3
  18. data/vendor/EigenRand/EigenRand/PacketFilter.h +11 -253
  19. data/vendor/EigenRand/EigenRand/PacketRandomEngine.h +21 -47
  20. data/vendor/EigenRand/EigenRand/RandUtils.h +50 -344
  21. data/vendor/EigenRand/EigenRand/arch/AVX/MorePacketMath.h +619 -0
  22. data/vendor/EigenRand/EigenRand/arch/AVX/PacketFilter.h +149 -0
  23. data/vendor/EigenRand/EigenRand/arch/AVX/RandUtils.h +228 -0
  24. data/vendor/EigenRand/EigenRand/arch/NEON/MorePacketMath.h +473 -0
  25. data/vendor/EigenRand/EigenRand/arch/NEON/PacketFilter.h +142 -0
  26. data/vendor/EigenRand/EigenRand/arch/NEON/RandUtils.h +126 -0
  27. data/vendor/EigenRand/EigenRand/arch/SSE/MorePacketMath.h +501 -0
  28. data/vendor/EigenRand/EigenRand/arch/SSE/PacketFilter.h +133 -0
  29. data/vendor/EigenRand/EigenRand/arch/SSE/RandUtils.h +120 -0
  30. data/vendor/EigenRand/EigenRand/doc.h +24 -12
  31. data/vendor/EigenRand/README.md +57 -4
  32. data/vendor/eigen/COPYING.APACHE +203 -0
  33. data/vendor/eigen/COPYING.BSD +1 -1
  34. data/vendor/eigen/COPYING.MINPACK +51 -52
  35. data/vendor/eigen/Eigen/Cholesky +0 -1
  36. data/vendor/eigen/Eigen/Core +112 -265
  37. data/vendor/eigen/Eigen/Eigenvalues +2 -3
  38. data/vendor/eigen/Eigen/Geometry +5 -8
  39. data/vendor/eigen/Eigen/Householder +0 -1
  40. data/vendor/eigen/Eigen/Jacobi +0 -1
  41. data/vendor/eigen/Eigen/KLUSupport +41 -0
  42. data/vendor/eigen/Eigen/LU +2 -5
  43. data/vendor/eigen/Eigen/OrderingMethods +0 -3
  44. data/vendor/eigen/Eigen/PaStiXSupport +1 -0
  45. data/vendor/eigen/Eigen/PardisoSupport +0 -0
  46. data/vendor/eigen/Eigen/QR +2 -3
  47. data/vendor/eigen/Eigen/QtAlignedMalloc +0 -1
  48. data/vendor/eigen/Eigen/SVD +0 -1
  49. data/vendor/eigen/Eigen/Sparse +0 -2
  50. data/vendor/eigen/Eigen/SparseCholesky +0 -8
  51. data/vendor/eigen/Eigen/SparseLU +4 -0
  52. data/vendor/eigen/Eigen/SparseQR +0 -1
  53. data/vendor/eigen/Eigen/src/Cholesky/LDLT.h +42 -27
  54. data/vendor/eigen/Eigen/src/Cholesky/LLT.h +39 -23
  55. data/vendor/eigen/Eigen/src/CholmodSupport/CholmodSupport.h +90 -47
  56. data/vendor/eigen/Eigen/src/Core/ArithmeticSequence.h +413 -0
  57. data/vendor/eigen/Eigen/src/Core/Array.h +99 -11
  58. data/vendor/eigen/Eigen/src/Core/ArrayBase.h +3 -3
  59. data/vendor/eigen/Eigen/src/Core/ArrayWrapper.h +21 -21
  60. data/vendor/eigen/Eigen/src/Core/Assign.h +1 -1
  61. data/vendor/eigen/Eigen/src/Core/AssignEvaluator.h +125 -50
  62. data/vendor/eigen/Eigen/src/Core/Assign_MKL.h +10 -10
  63. data/vendor/eigen/Eigen/src/Core/BandMatrix.h +16 -16
  64. data/vendor/eigen/Eigen/src/Core/Block.h +56 -60
  65. data/vendor/eigen/Eigen/src/Core/BooleanRedux.h +29 -31
  66. data/vendor/eigen/Eigen/src/Core/CommaInitializer.h +7 -3
  67. data/vendor/eigen/Eigen/src/Core/CoreEvaluators.h +325 -272
  68. data/vendor/eigen/Eigen/src/Core/CoreIterators.h +5 -0
  69. data/vendor/eigen/Eigen/src/Core/CwiseBinaryOp.h +21 -22
  70. data/vendor/eigen/Eigen/src/Core/CwiseNullaryOp.h +153 -18
  71. data/vendor/eigen/Eigen/src/Core/CwiseUnaryOp.h +6 -6
  72. data/vendor/eigen/Eigen/src/Core/CwiseUnaryView.h +14 -10
  73. data/vendor/eigen/Eigen/src/Core/DenseBase.h +132 -42
  74. data/vendor/eigen/Eigen/src/Core/DenseCoeffsBase.h +25 -21
  75. data/vendor/eigen/Eigen/src/Core/DenseStorage.h +153 -71
  76. data/vendor/eigen/Eigen/src/Core/Diagonal.h +21 -23
  77. data/vendor/eigen/Eigen/src/Core/DiagonalMatrix.h +50 -2
  78. data/vendor/eigen/Eigen/src/Core/DiagonalProduct.h +1 -1
  79. data/vendor/eigen/Eigen/src/Core/Dot.h +10 -10
  80. data/vendor/eigen/Eigen/src/Core/EigenBase.h +10 -9
  81. data/vendor/eigen/Eigen/src/Core/ForceAlignedAccess.h +8 -4
  82. data/vendor/eigen/Eigen/src/Core/Fuzzy.h +3 -3
  83. data/vendor/eigen/Eigen/src/Core/GeneralProduct.h +20 -10
  84. data/vendor/eigen/Eigen/src/Core/GenericPacketMath.h +599 -152
  85. data/vendor/eigen/Eigen/src/Core/GlobalFunctions.h +40 -33
  86. data/vendor/eigen/Eigen/src/Core/IO.h +40 -7
  87. data/vendor/eigen/Eigen/src/Core/IndexedView.h +237 -0
  88. data/vendor/eigen/Eigen/src/Core/Inverse.h +9 -10
  89. data/vendor/eigen/Eigen/src/Core/Map.h +7 -7
  90. data/vendor/eigen/Eigen/src/Core/MapBase.h +10 -3
  91. data/vendor/eigen/Eigen/src/Core/MathFunctions.h +767 -125
  92. data/vendor/eigen/Eigen/src/Core/MathFunctionsImpl.h +118 -19
  93. data/vendor/eigen/Eigen/src/Core/Matrix.h +131 -25
  94. data/vendor/eigen/Eigen/src/Core/MatrixBase.h +21 -3
  95. data/vendor/eigen/Eigen/src/Core/NestByValue.h +25 -50
  96. data/vendor/eigen/Eigen/src/Core/NoAlias.h +4 -3
  97. data/vendor/eigen/Eigen/src/Core/NumTraits.h +107 -20
  98. data/vendor/eigen/Eigen/src/Core/PartialReduxEvaluator.h +232 -0
  99. data/vendor/eigen/Eigen/src/Core/PermutationMatrix.h +3 -31
  100. data/vendor/eigen/Eigen/src/Core/PlainObjectBase.h +152 -59
  101. data/vendor/eigen/Eigen/src/Core/Product.h +30 -25
  102. data/vendor/eigen/Eigen/src/Core/ProductEvaluators.h +192 -125
  103. data/vendor/eigen/Eigen/src/Core/Random.h +37 -1
  104. data/vendor/eigen/Eigen/src/Core/Redux.h +180 -170
  105. data/vendor/eigen/Eigen/src/Core/Ref.h +121 -23
  106. data/vendor/eigen/Eigen/src/Core/Replicate.h +8 -8
  107. data/vendor/eigen/Eigen/src/Core/Reshaped.h +454 -0
  108. data/vendor/eigen/Eigen/src/Core/ReturnByValue.h +7 -5
  109. data/vendor/eigen/Eigen/src/Core/Reverse.h +18 -12
  110. data/vendor/eigen/Eigen/src/Core/Select.h +8 -6
  111. data/vendor/eigen/Eigen/src/Core/SelfAdjointView.h +33 -20
  112. data/vendor/eigen/Eigen/src/Core/Solve.h +14 -14
  113. data/vendor/eigen/Eigen/src/Core/SolveTriangular.h +16 -16
  114. data/vendor/eigen/Eigen/src/Core/SolverBase.h +41 -3
  115. data/vendor/eigen/Eigen/src/Core/StableNorm.h +100 -70
  116. data/vendor/eigen/Eigen/src/Core/StlIterators.h +463 -0
  117. data/vendor/eigen/Eigen/src/Core/Stride.h +9 -4
  118. data/vendor/eigen/Eigen/src/Core/Swap.h +5 -4
  119. data/vendor/eigen/Eigen/src/Core/Transpose.h +88 -27
  120. data/vendor/eigen/Eigen/src/Core/Transpositions.h +26 -47
  121. data/vendor/eigen/Eigen/src/Core/TriangularMatrix.h +93 -75
  122. data/vendor/eigen/Eigen/src/Core/VectorBlock.h +5 -5
  123. data/vendor/eigen/Eigen/src/Core/VectorwiseOp.h +159 -70
  124. data/vendor/eigen/Eigen/src/Core/Visitor.h +137 -29
  125. data/vendor/eigen/Eigen/src/Core/arch/AVX/Complex.h +50 -129
  126. data/vendor/eigen/Eigen/src/Core/arch/AVX/MathFunctions.h +126 -337
  127. data/vendor/eigen/Eigen/src/Core/arch/AVX/PacketMath.h +1092 -155
  128. data/vendor/eigen/Eigen/src/Core/arch/AVX/TypeCasting.h +65 -1
  129. data/vendor/eigen/Eigen/src/Core/arch/AVX512/Complex.h +422 -0
  130. data/vendor/eigen/Eigen/src/Core/arch/AVX512/MathFunctions.h +207 -236
  131. data/vendor/eigen/Eigen/src/Core/arch/AVX512/PacketMath.h +1482 -495
  132. data/vendor/eigen/Eigen/src/Core/arch/AVX512/TypeCasting.h +89 -0
  133. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/Complex.h +152 -165
  134. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/MathFunctions.h +19 -251
  135. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/MatrixProduct.h +2937 -0
  136. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/MatrixProductCommon.h +221 -0
  137. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/MatrixProductMMA.h +629 -0
  138. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/PacketMath.h +2042 -392
  139. data/vendor/eigen/Eigen/src/Core/arch/CUDA/Complex.h +235 -80
  140. data/vendor/eigen/Eigen/src/Core/arch/Default/BFloat16.h +700 -0
  141. data/vendor/eigen/Eigen/src/Core/arch/Default/ConjHelper.h +102 -14
  142. data/vendor/eigen/Eigen/src/Core/arch/Default/GenericPacketMathFunctions.h +1649 -0
  143. data/vendor/eigen/Eigen/src/Core/arch/Default/GenericPacketMathFunctionsFwd.h +110 -0
  144. data/vendor/eigen/Eigen/src/Core/arch/Default/Half.h +942 -0
  145. data/vendor/eigen/Eigen/src/Core/arch/Default/Settings.h +1 -1
  146. data/vendor/eigen/Eigen/src/Core/arch/Default/TypeCasting.h +120 -0
  147. data/vendor/eigen/Eigen/src/Core/arch/{CUDA → GPU}/MathFunctions.h +16 -4
  148. data/vendor/eigen/Eigen/src/Core/arch/GPU/PacketMath.h +1685 -0
  149. data/vendor/eigen/Eigen/src/Core/arch/GPU/TypeCasting.h +80 -0
  150. data/vendor/eigen/Eigen/src/Core/arch/HIP/hcc/math_constants.h +23 -0
  151. data/vendor/eigen/Eigen/src/Core/arch/MSA/Complex.h +648 -0
  152. data/vendor/eigen/Eigen/src/Core/arch/MSA/MathFunctions.h +387 -0
  153. data/vendor/eigen/Eigen/src/Core/arch/MSA/PacketMath.h +1233 -0
  154. data/vendor/eigen/Eigen/src/Core/arch/NEON/Complex.h +313 -219
  155. data/vendor/eigen/Eigen/src/Core/arch/NEON/GeneralBlockPanelKernel.h +183 -0
  156. data/vendor/eigen/Eigen/src/Core/arch/NEON/MathFunctions.h +54 -70
  157. data/vendor/eigen/Eigen/src/Core/arch/NEON/PacketMath.h +4376 -549
  158. data/vendor/eigen/Eigen/src/Core/arch/NEON/TypeCasting.h +1419 -0
  159. data/vendor/eigen/Eigen/src/Core/arch/SSE/Complex.h +59 -179
  160. data/vendor/eigen/Eigen/src/Core/arch/SSE/MathFunctions.h +65 -428
  161. data/vendor/eigen/Eigen/src/Core/arch/SSE/PacketMath.h +893 -283
  162. data/vendor/eigen/Eigen/src/Core/arch/SSE/TypeCasting.h +65 -0
  163. data/vendor/eigen/Eigen/src/Core/arch/SVE/MathFunctions.h +44 -0
  164. data/vendor/eigen/Eigen/src/Core/arch/SVE/PacketMath.h +752 -0
  165. data/vendor/eigen/Eigen/src/Core/arch/SVE/TypeCasting.h +49 -0
  166. data/vendor/eigen/Eigen/src/Core/arch/SYCL/InteropHeaders.h +232 -0
  167. data/vendor/eigen/Eigen/src/Core/arch/SYCL/MathFunctions.h +301 -0
  168. data/vendor/eigen/Eigen/src/Core/arch/SYCL/PacketMath.h +670 -0
  169. data/vendor/eigen/Eigen/src/Core/arch/SYCL/SyclMemoryModel.h +694 -0
  170. data/vendor/eigen/Eigen/src/Core/arch/SYCL/TypeCasting.h +85 -0
  171. data/vendor/eigen/Eigen/src/Core/arch/ZVector/Complex.h +212 -183
  172. data/vendor/eigen/Eigen/src/Core/arch/ZVector/MathFunctions.h +101 -5
  173. data/vendor/eigen/Eigen/src/Core/arch/ZVector/PacketMath.h +510 -395
  174. data/vendor/eigen/Eigen/src/Core/functors/AssignmentFunctors.h +11 -2
  175. data/vendor/eigen/Eigen/src/Core/functors/BinaryFunctors.h +112 -46
  176. data/vendor/eigen/Eigen/src/Core/functors/NullaryFunctors.h +31 -30
  177. data/vendor/eigen/Eigen/src/Core/functors/StlFunctors.h +32 -2
  178. data/vendor/eigen/Eigen/src/Core/functors/UnaryFunctors.h +355 -16
  179. data/vendor/eigen/Eigen/src/Core/products/GeneralBlockPanelKernel.h +1075 -586
  180. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrix.h +49 -24
  181. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrixTriangular.h +41 -35
  182. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h +6 -6
  183. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrix_BLAS.h +4 -2
  184. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixVector.h +382 -483
  185. data/vendor/eigen/Eigen/src/Core/products/Parallelizer.h +22 -5
  186. data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixMatrix.h +53 -30
  187. data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixMatrix_BLAS.h +16 -8
  188. data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixVector.h +8 -6
  189. data/vendor/eigen/Eigen/src/Core/products/SelfadjointProduct.h +4 -4
  190. data/vendor/eigen/Eigen/src/Core/products/SelfadjointRank2Update.h +5 -4
  191. data/vendor/eigen/Eigen/src/Core/products/TriangularMatrixMatrix.h +33 -27
  192. data/vendor/eigen/Eigen/src/Core/products/TriangularMatrixMatrix_BLAS.h +14 -12
  193. data/vendor/eigen/Eigen/src/Core/products/TriangularSolverMatrix.h +36 -34
  194. data/vendor/eigen/Eigen/src/Core/products/TriangularSolverMatrix_BLAS.h +8 -4
  195. data/vendor/eigen/Eigen/src/Core/products/TriangularSolverVector.h +13 -10
  196. data/vendor/eigen/Eigen/src/Core/util/BlasUtil.h +304 -119
  197. data/vendor/eigen/Eigen/src/Core/util/ConfigureVectorization.h +512 -0
  198. data/vendor/eigen/Eigen/src/Core/util/Constants.h +25 -9
  199. data/vendor/eigen/Eigen/src/Core/util/DisableStupidWarnings.h +26 -3
  200. data/vendor/eigen/Eigen/src/Core/util/ForwardDeclarations.h +29 -9
  201. data/vendor/eigen/Eigen/src/Core/util/IndexedViewHelper.h +186 -0
  202. data/vendor/eigen/Eigen/src/Core/util/IntegralConstant.h +272 -0
  203. data/vendor/eigen/Eigen/src/Core/util/MKL_support.h +8 -1
  204. data/vendor/eigen/Eigen/src/Core/util/Macros.h +709 -246
  205. data/vendor/eigen/Eigen/src/Core/util/Memory.h +222 -52
  206. data/vendor/eigen/Eigen/src/Core/util/Meta.h +355 -77
  207. data/vendor/eigen/Eigen/src/Core/util/ReenableStupidWarnings.h +5 -1
  208. data/vendor/eigen/Eigen/src/Core/util/ReshapedHelper.h +51 -0
  209. data/vendor/eigen/Eigen/src/Core/util/StaticAssert.h +8 -5
  210. data/vendor/eigen/Eigen/src/Core/util/SymbolicIndex.h +293 -0
  211. data/vendor/eigen/Eigen/src/Core/util/XprHelper.h +65 -30
  212. data/vendor/eigen/Eigen/src/Eigenvalues/ComplexEigenSolver.h +1 -1
  213. data/vendor/eigen/Eigen/src/Eigenvalues/ComplexSchur.h +7 -4
  214. data/vendor/eigen/Eigen/src/Eigenvalues/EigenSolver.h +2 -2
  215. data/vendor/eigen/Eigen/src/Eigenvalues/GeneralizedSelfAdjointEigenSolver.h +1 -1
  216. data/vendor/eigen/Eigen/src/Eigenvalues/HessenbergDecomposition.h +2 -2
  217. data/vendor/eigen/Eigen/src/Eigenvalues/MatrixBaseEigenvalues.h +2 -2
  218. data/vendor/eigen/Eigen/src/Eigenvalues/RealQZ.h +9 -6
  219. data/vendor/eigen/Eigen/src/Eigenvalues/RealSchur.h +21 -9
  220. data/vendor/eigen/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h +77 -43
  221. data/vendor/eigen/Eigen/src/Eigenvalues/Tridiagonalization.h +20 -15
  222. data/vendor/eigen/Eigen/src/Geometry/AlignedBox.h +99 -5
  223. data/vendor/eigen/Eigen/src/Geometry/AngleAxis.h +4 -4
  224. data/vendor/eigen/Eigen/src/Geometry/EulerAngles.h +3 -3
  225. data/vendor/eigen/Eigen/src/Geometry/Homogeneous.h +15 -11
  226. data/vendor/eigen/Eigen/src/Geometry/Hyperplane.h +1 -1
  227. data/vendor/eigen/Eigen/src/Geometry/OrthoMethods.h +3 -2
  228. data/vendor/eigen/Eigen/src/Geometry/ParametrizedLine.h +39 -2
  229. data/vendor/eigen/Eigen/src/Geometry/Quaternion.h +70 -14
  230. data/vendor/eigen/Eigen/src/Geometry/Rotation2D.h +3 -3
  231. data/vendor/eigen/Eigen/src/Geometry/Scaling.h +23 -5
  232. data/vendor/eigen/Eigen/src/Geometry/Transform.h +88 -67
  233. data/vendor/eigen/Eigen/src/Geometry/Translation.h +6 -12
  234. data/vendor/eigen/Eigen/src/Geometry/Umeyama.h +1 -1
  235. data/vendor/eigen/Eigen/src/Geometry/arch/Geometry_SIMD.h +168 -0
  236. data/vendor/eigen/Eigen/src/Householder/BlockHouseholder.h +9 -2
  237. data/vendor/eigen/Eigen/src/Householder/Householder.h +8 -4
  238. data/vendor/eigen/Eigen/src/Householder/HouseholderSequence.h +123 -48
  239. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/BasicPreconditioners.h +15 -15
  240. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/BiCGSTAB.h +7 -23
  241. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h +5 -22
  242. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/IncompleteCholesky.h +41 -47
  243. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/IncompleteLUT.h +51 -60
  244. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/IterativeSolverBase.h +70 -20
  245. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/LeastSquareConjugateGradient.h +2 -20
  246. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/SolveWithGuess.h +11 -9
  247. data/vendor/eigen/Eigen/src/Jacobi/Jacobi.h +31 -10
  248. data/vendor/eigen/Eigen/src/KLUSupport/KLUSupport.h +358 -0
  249. data/vendor/eigen/Eigen/src/LU/Determinant.h +35 -19
  250. data/vendor/eigen/Eigen/src/LU/FullPivLU.h +29 -43
  251. data/vendor/eigen/Eigen/src/LU/InverseImpl.h +25 -8
  252. data/vendor/eigen/Eigen/src/LU/PartialPivLU.h +71 -58
  253. data/vendor/eigen/Eigen/src/LU/arch/InverseSize4.h +351 -0
  254. data/vendor/eigen/Eigen/src/OrderingMethods/Amd.h +7 -17
  255. data/vendor/eigen/Eigen/src/OrderingMethods/Eigen_Colamd.h +297 -277
  256. data/vendor/eigen/Eigen/src/OrderingMethods/Ordering.h +6 -10
  257. data/vendor/eigen/Eigen/src/PaStiXSupport/PaStiXSupport.h +1 -1
  258. data/vendor/eigen/Eigen/src/PardisoSupport/PardisoSupport.h +9 -7
  259. data/vendor/eigen/Eigen/src/QR/ColPivHouseholderQR.h +41 -20
  260. data/vendor/eigen/Eigen/src/QR/CompleteOrthogonalDecomposition.h +100 -27
  261. data/vendor/eigen/Eigen/src/QR/FullPivHouseholderQR.h +59 -22
  262. data/vendor/eigen/Eigen/src/QR/HouseholderQR.h +48 -23
  263. data/vendor/eigen/Eigen/src/SPQRSupport/SuiteSparseQRSupport.h +25 -3
  264. data/vendor/eigen/Eigen/src/SVD/BDCSVD.h +183 -63
  265. data/vendor/eigen/Eigen/src/SVD/JacobiSVD.h +22 -14
  266. data/vendor/eigen/Eigen/src/SVD/SVDBase.h +83 -22
  267. data/vendor/eigen/Eigen/src/SVD/UpperBidiagonalization.h +3 -3
  268. data/vendor/eigen/Eigen/src/SparseCholesky/SimplicialCholesky.h +17 -9
  269. data/vendor/eigen/Eigen/src/SparseCholesky/SimplicialCholesky_impl.h +12 -37
  270. data/vendor/eigen/Eigen/src/SparseCore/AmbiVector.h +3 -2
  271. data/vendor/eigen/Eigen/src/SparseCore/CompressedStorage.h +16 -0
  272. data/vendor/eigen/Eigen/src/SparseCore/ConservativeSparseSparseProduct.h +6 -6
  273. data/vendor/eigen/Eigen/src/SparseCore/SparseAssign.h +81 -27
  274. data/vendor/eigen/Eigen/src/SparseCore/SparseBlock.h +25 -57
  275. data/vendor/eigen/Eigen/src/SparseCore/SparseCompressedBase.h +40 -11
  276. data/vendor/eigen/Eigen/src/SparseCore/SparseCwiseBinaryOp.h +11 -15
  277. data/vendor/eigen/Eigen/src/SparseCore/SparseCwiseUnaryOp.h +4 -2
  278. data/vendor/eigen/Eigen/src/SparseCore/SparseDenseProduct.h +30 -8
  279. data/vendor/eigen/Eigen/src/SparseCore/SparseMatrix.h +126 -11
  280. data/vendor/eigen/Eigen/src/SparseCore/SparseMatrixBase.h +5 -12
  281. data/vendor/eigen/Eigen/src/SparseCore/SparseProduct.h +13 -1
  282. data/vendor/eigen/Eigen/src/SparseCore/SparseRef.h +7 -7
  283. data/vendor/eigen/Eigen/src/SparseCore/SparseSelfAdjointView.h +5 -2
  284. data/vendor/eigen/Eigen/src/SparseCore/SparseUtil.h +8 -0
  285. data/vendor/eigen/Eigen/src/SparseCore/SparseVector.h +1 -1
  286. data/vendor/eigen/Eigen/src/SparseCore/SparseView.h +1 -0
  287. data/vendor/eigen/Eigen/src/SparseLU/SparseLU.h +162 -12
  288. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_Memory.h +1 -1
  289. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_SupernodalMatrix.h +76 -2
  290. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_column_dfs.h +2 -2
  291. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_gemm_kernel.h +1 -1
  292. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_panel_bmod.h +1 -1
  293. data/vendor/eigen/Eigen/src/SparseQR/SparseQR.h +19 -6
  294. data/vendor/eigen/Eigen/src/StlSupport/StdDeque.h +2 -12
  295. data/vendor/eigen/Eigen/src/StlSupport/StdList.h +2 -2
  296. data/vendor/eigen/Eigen/src/StlSupport/StdVector.h +2 -2
  297. data/vendor/eigen/Eigen/src/SuperLUSupport/SuperLUSupport.h +6 -8
  298. data/vendor/eigen/Eigen/src/UmfPackSupport/UmfPackSupport.h +175 -39
  299. data/vendor/eigen/Eigen/src/misc/lapacke.h +5 -4
  300. data/vendor/eigen/Eigen/src/plugins/ArrayCwiseBinaryOps.h +28 -2
  301. data/vendor/eigen/Eigen/src/plugins/ArrayCwiseUnaryOps.h +155 -11
  302. data/vendor/eigen/Eigen/src/plugins/BlockMethods.h +626 -242
  303. data/vendor/eigen/Eigen/src/plugins/CommonCwiseUnaryOps.h +14 -0
  304. data/vendor/eigen/Eigen/src/plugins/IndexedViewMethods.h +262 -0
  305. data/vendor/eigen/Eigen/src/plugins/MatrixCwiseBinaryOps.h +4 -4
  306. data/vendor/eigen/Eigen/src/plugins/MatrixCwiseUnaryOps.h +10 -0
  307. data/vendor/eigen/Eigen/src/plugins/ReshapedMethods.h +149 -0
  308. data/vendor/eigen/README.md +2 -0
  309. data/vendor/eigen/bench/btl/README +1 -1
  310. data/vendor/eigen/bench/tensors/README +6 -7
  311. data/vendor/eigen/ci/README.md +56 -0
  312. data/vendor/eigen/demos/mix_eigen_and_c/README +1 -1
  313. data/vendor/eigen/unsupported/Eigen/CXX11/src/Tensor/README.md +213 -158
  314. data/vendor/eigen/unsupported/README.txt +1 -1
  315. data/vendor/tomotopy/README.kr.rst +21 -0
  316. data/vendor/tomotopy/README.rst +20 -0
  317. data/vendor/tomotopy/src/Labeling/FoRelevance.cpp +2 -2
  318. data/vendor/tomotopy/src/Labeling/Phraser.hpp +1 -1
  319. data/vendor/tomotopy/src/TopicModel/CTModel.hpp +2 -1
  320. data/vendor/tomotopy/src/TopicModel/DMRModel.hpp +2 -1
  321. data/vendor/tomotopy/src/TopicModel/DTModel.hpp +1 -1
  322. data/vendor/tomotopy/src/TopicModel/GDMRModel.hpp +2 -2
  323. data/vendor/tomotopy/src/TopicModel/HDP.h +1 -0
  324. data/vendor/tomotopy/src/TopicModel/HDPModel.hpp +53 -2
  325. data/vendor/tomotopy/src/TopicModel/HLDAModel.hpp +1 -1
  326. data/vendor/tomotopy/src/TopicModel/HPAModel.hpp +1 -0
  327. data/vendor/tomotopy/src/TopicModel/LDACVB0Model.hpp +2 -2
  328. data/vendor/tomotopy/src/TopicModel/LDAModel.hpp +16 -5
  329. data/vendor/tomotopy/src/TopicModel/LLDAModel.hpp +1 -0
  330. data/vendor/tomotopy/src/TopicModel/MGLDAModel.hpp +1 -0
  331. data/vendor/tomotopy/src/TopicModel/PLDAModel.hpp +1 -0
  332. data/vendor/tomotopy/src/TopicModel/PT.h +3 -1
  333. data/vendor/tomotopy/src/TopicModel/PTModel.hpp +31 -1
  334. data/vendor/tomotopy/src/TopicModel/SLDAModel.hpp +2 -2
  335. data/vendor/tomotopy/src/TopicModel/TopicModel.hpp +7 -5
  336. data/vendor/tomotopy/src/Utils/EigenAddonOps.hpp +36 -1
  337. data/vendor/tomotopy/src/Utils/exception.h +6 -0
  338. data/vendor/tomotopy/src/Utils/sample.hpp +14 -12
  339. data/vendor/tomotopy/src/Utils/sse_gamma.h +0 -3
  340. metadata +60 -14
  341. data/vendor/eigen/Eigen/CMakeLists.txt +0 -19
  342. data/vendor/eigen/Eigen/src/Core/arch/CUDA/Half.h +0 -674
  343. data/vendor/eigen/Eigen/src/Core/arch/CUDA/PacketMath.h +0 -333
  344. data/vendor/eigen/Eigen/src/Core/arch/CUDA/PacketMathHalf.h +0 -1124
  345. data/vendor/eigen/Eigen/src/Core/arch/CUDA/TypeCasting.h +0 -212
  346. data/vendor/eigen/Eigen/src/Geometry/arch/Geometry_SSE.h +0 -161
  347. data/vendor/eigen/Eigen/src/LU/arch/Inverse_SSE.h +0 -338
@@ -0,0 +1,512 @@
1
+ // This file is part of Eigen, a lightweight C++ template library
2
+ // for linear algebra.
3
+ //
4
+ // Copyright (C) 2008-2018 Gael Guennebaud <gael.guennebaud@inria.fr>
5
+ // Copyright (C) 2020, Arm Limited and Contributors
6
+ //
7
+ // This Source Code Form is subject to the terms of the Mozilla
8
+ // Public License v. 2.0. If a copy of the MPL was not distributed
9
+ // with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
10
+
11
+ #ifndef EIGEN_CONFIGURE_VECTORIZATION_H
12
+ #define EIGEN_CONFIGURE_VECTORIZATION_H
13
+
14
+ //------------------------------------------------------------------------------------------
15
+ // Static and dynamic alignment control
16
+ //
17
+ // The main purpose of this section is to define EIGEN_MAX_ALIGN_BYTES and EIGEN_MAX_STATIC_ALIGN_BYTES
18
+ // as the maximal boundary in bytes on which dynamically and statically allocated data may be alignment respectively.
19
+ // The values of EIGEN_MAX_ALIGN_BYTES and EIGEN_MAX_STATIC_ALIGN_BYTES can be specified by the user. If not,
20
+ // a default value is automatically computed based on architecture, compiler, and OS.
21
+ //
22
+ // This section also defines macros EIGEN_ALIGN_TO_BOUNDARY(N) and the shortcuts EIGEN_ALIGN{8,16,32,_MAX}
23
+ // to be used to declare statically aligned buffers.
24
+ //------------------------------------------------------------------------------------------
25
+
26
+
27
+ /* EIGEN_ALIGN_TO_BOUNDARY(n) forces data to be n-byte aligned. This is used to satisfy SIMD requirements.
28
+ * However, we do that EVEN if vectorization (EIGEN_VECTORIZE) is disabled,
29
+ * so that vectorization doesn't affect binary compatibility.
30
+ *
31
+ * If we made alignment depend on whether or not EIGEN_VECTORIZE is defined, it would be impossible to link
32
+ * vectorized and non-vectorized code.
33
+ *
34
+ * FIXME: this code can be cleaned up once we switch to proper C++11 only.
35
+ */
36
+ #if (defined EIGEN_CUDACC)
37
+ #define EIGEN_ALIGN_TO_BOUNDARY(n) __align__(n)
38
+ #define EIGEN_ALIGNOF(x) __alignof(x)
39
+ #elif EIGEN_HAS_ALIGNAS
40
+ #define EIGEN_ALIGN_TO_BOUNDARY(n) alignas(n)
41
+ #define EIGEN_ALIGNOF(x) alignof(x)
42
+ #elif EIGEN_COMP_GNUC || EIGEN_COMP_PGI || EIGEN_COMP_IBM || EIGEN_COMP_ARM
43
+ #define EIGEN_ALIGN_TO_BOUNDARY(n) __attribute__((aligned(n)))
44
+ #define EIGEN_ALIGNOF(x) __alignof(x)
45
+ #elif EIGEN_COMP_MSVC
46
+ #define EIGEN_ALIGN_TO_BOUNDARY(n) __declspec(align(n))
47
+ #define EIGEN_ALIGNOF(x) __alignof(x)
48
+ #elif EIGEN_COMP_SUNCC
49
+ // FIXME not sure about this one:
50
+ #define EIGEN_ALIGN_TO_BOUNDARY(n) __attribute__((aligned(n)))
51
+ #define EIGEN_ALIGNOF(x) __alignof(x)
52
+ #else
53
+ #error Please tell me what is the equivalent of alignas(n) and alignof(x) for your compiler
54
+ #endif
55
+
56
+ // If the user explicitly disable vectorization, then we also disable alignment
57
+ #if defined(EIGEN_DONT_VECTORIZE)
58
+ #if defined(EIGEN_GPUCC)
59
+ // GPU code is always vectorized and requires memory alignment for
60
+ // statically allocated buffers.
61
+ #define EIGEN_IDEAL_MAX_ALIGN_BYTES 16
62
+ #else
63
+ #define EIGEN_IDEAL_MAX_ALIGN_BYTES 0
64
+ #endif
65
+ #elif defined(__AVX512F__)
66
+ // 64 bytes static alignment is preferred only if really required
67
+ #define EIGEN_IDEAL_MAX_ALIGN_BYTES 64
68
+ #elif defined(__AVX__)
69
+ // 32 bytes static alignment is preferred only if really required
70
+ #define EIGEN_IDEAL_MAX_ALIGN_BYTES 32
71
+ #else
72
+ #define EIGEN_IDEAL_MAX_ALIGN_BYTES 16
73
+ #endif
74
+
75
+
76
+ // EIGEN_MIN_ALIGN_BYTES defines the minimal value for which the notion of explicit alignment makes sense
77
+ #define EIGEN_MIN_ALIGN_BYTES 16
78
+
79
+ // Defined the boundary (in bytes) on which the data needs to be aligned. Note
80
+ // that unless EIGEN_ALIGN is defined and not equal to 0, the data may not be
81
+ // aligned at all regardless of the value of this #define.
82
+
83
+ #if (defined(EIGEN_DONT_ALIGN_STATICALLY) || defined(EIGEN_DONT_ALIGN)) && defined(EIGEN_MAX_STATIC_ALIGN_BYTES) && EIGEN_MAX_STATIC_ALIGN_BYTES>0
84
+ #error EIGEN_MAX_STATIC_ALIGN_BYTES and EIGEN_DONT_ALIGN[_STATICALLY] are both defined with EIGEN_MAX_STATIC_ALIGN_BYTES!=0. Use EIGEN_MAX_STATIC_ALIGN_BYTES=0 as a synonym of EIGEN_DONT_ALIGN_STATICALLY.
85
+ #endif
86
+
87
+ // EIGEN_DONT_ALIGN_STATICALLY and EIGEN_DONT_ALIGN are deprecated
88
+ // They imply EIGEN_MAX_STATIC_ALIGN_BYTES=0
89
+ #if defined(EIGEN_DONT_ALIGN_STATICALLY) || defined(EIGEN_DONT_ALIGN)
90
+ #ifdef EIGEN_MAX_STATIC_ALIGN_BYTES
91
+ #undef EIGEN_MAX_STATIC_ALIGN_BYTES
92
+ #endif
93
+ #define EIGEN_MAX_STATIC_ALIGN_BYTES 0
94
+ #endif
95
+
96
+ #ifndef EIGEN_MAX_STATIC_ALIGN_BYTES
97
+
98
+ // Try to automatically guess what is the best default value for EIGEN_MAX_STATIC_ALIGN_BYTES
99
+
100
+ // 16 byte alignment is only useful for vectorization. Since it affects the ABI, we need to enable
101
+ // 16 byte alignment on all platforms where vectorization might be enabled. In theory we could always
102
+ // enable alignment, but it can be a cause of problems on some platforms, so we just disable it in
103
+ // certain common platform (compiler+architecture combinations) to avoid these problems.
104
+ // Only static alignment is really problematic (relies on nonstandard compiler extensions),
105
+ // try to keep heap alignment even when we have to disable static alignment.
106
+ #if EIGEN_COMP_GNUC && !(EIGEN_ARCH_i386_OR_x86_64 || EIGEN_ARCH_ARM_OR_ARM64 || EIGEN_ARCH_PPC || EIGEN_ARCH_IA64 || EIGEN_ARCH_MIPS)
107
+ #define EIGEN_GCC_AND_ARCH_DOESNT_WANT_STACK_ALIGNMENT 1
108
+ #elif EIGEN_ARCH_ARM_OR_ARM64 && EIGEN_COMP_GNUC_STRICT && EIGEN_GNUC_AT_MOST(4, 6)
109
+ // Old versions of GCC on ARM, at least 4.4, were once seen to have buggy static alignment support.
110
+ // Not sure which version fixed it, hopefully it doesn't affect 4.7, which is still somewhat in use.
111
+ // 4.8 and newer seem definitely unaffected.
112
+ #define EIGEN_GCC_AND_ARCH_DOESNT_WANT_STACK_ALIGNMENT 1
113
+ #else
114
+ #define EIGEN_GCC_AND_ARCH_DOESNT_WANT_STACK_ALIGNMENT 0
115
+ #endif
116
+
117
+ // static alignment is completely disabled with GCC 3, Sun Studio, and QCC/QNX
118
+ #if !EIGEN_GCC_AND_ARCH_DOESNT_WANT_STACK_ALIGNMENT \
119
+ && !EIGEN_GCC3_OR_OLDER \
120
+ && !EIGEN_COMP_SUNCC \
121
+ && !EIGEN_OS_QNX
122
+ #define EIGEN_ARCH_WANTS_STACK_ALIGNMENT 1
123
+ #else
124
+ #define EIGEN_ARCH_WANTS_STACK_ALIGNMENT 0
125
+ #endif
126
+
127
+ #if EIGEN_ARCH_WANTS_STACK_ALIGNMENT
128
+ #define EIGEN_MAX_STATIC_ALIGN_BYTES EIGEN_IDEAL_MAX_ALIGN_BYTES
129
+ #else
130
+ #define EIGEN_MAX_STATIC_ALIGN_BYTES 0
131
+ #endif
132
+
133
+ #endif
134
+
135
+ // If EIGEN_MAX_ALIGN_BYTES is defined, then it is considered as an upper bound for EIGEN_MAX_STATIC_ALIGN_BYTES
136
+ #if defined(EIGEN_MAX_ALIGN_BYTES) && EIGEN_MAX_ALIGN_BYTES<EIGEN_MAX_STATIC_ALIGN_BYTES
137
+ #undef EIGEN_MAX_STATIC_ALIGN_BYTES
138
+ #define EIGEN_MAX_STATIC_ALIGN_BYTES EIGEN_MAX_ALIGN_BYTES
139
+ #endif
140
+
141
+ #if EIGEN_MAX_STATIC_ALIGN_BYTES==0 && !defined(EIGEN_DISABLE_UNALIGNED_ARRAY_ASSERT)
142
+ #define EIGEN_DISABLE_UNALIGNED_ARRAY_ASSERT
143
+ #endif
144
+
145
+ // At this stage, EIGEN_MAX_STATIC_ALIGN_BYTES>0 is the true test whether we want to align arrays on the stack or not.
146
+ // It takes into account both the user choice to explicitly enable/disable alignment (by setting EIGEN_MAX_STATIC_ALIGN_BYTES)
147
+ // and the architecture config (EIGEN_ARCH_WANTS_STACK_ALIGNMENT).
148
+ // Henceforth, only EIGEN_MAX_STATIC_ALIGN_BYTES should be used.
149
+
150
+
151
+ // Shortcuts to EIGEN_ALIGN_TO_BOUNDARY
152
+ #define EIGEN_ALIGN8 EIGEN_ALIGN_TO_BOUNDARY(8)
153
+ #define EIGEN_ALIGN16 EIGEN_ALIGN_TO_BOUNDARY(16)
154
+ #define EIGEN_ALIGN32 EIGEN_ALIGN_TO_BOUNDARY(32)
155
+ #define EIGEN_ALIGN64 EIGEN_ALIGN_TO_BOUNDARY(64)
156
+ #if EIGEN_MAX_STATIC_ALIGN_BYTES>0
157
+ #define EIGEN_ALIGN_MAX EIGEN_ALIGN_TO_BOUNDARY(EIGEN_MAX_STATIC_ALIGN_BYTES)
158
+ #else
159
+ #define EIGEN_ALIGN_MAX
160
+ #endif
161
+
162
+
163
+ // Dynamic alignment control
164
+
165
+ #if defined(EIGEN_DONT_ALIGN) && defined(EIGEN_MAX_ALIGN_BYTES) && EIGEN_MAX_ALIGN_BYTES>0
166
+ #error EIGEN_MAX_ALIGN_BYTES and EIGEN_DONT_ALIGN are both defined with EIGEN_MAX_ALIGN_BYTES!=0. Use EIGEN_MAX_ALIGN_BYTES=0 as a synonym of EIGEN_DONT_ALIGN.
167
+ #endif
168
+
169
+ #ifdef EIGEN_DONT_ALIGN
170
+ #ifdef EIGEN_MAX_ALIGN_BYTES
171
+ #undef EIGEN_MAX_ALIGN_BYTES
172
+ #endif
173
+ #define EIGEN_MAX_ALIGN_BYTES 0
174
+ #elif !defined(EIGEN_MAX_ALIGN_BYTES)
175
+ #define EIGEN_MAX_ALIGN_BYTES EIGEN_IDEAL_MAX_ALIGN_BYTES
176
+ #endif
177
+
178
+ #if EIGEN_IDEAL_MAX_ALIGN_BYTES > EIGEN_MAX_ALIGN_BYTES
179
+ #define EIGEN_DEFAULT_ALIGN_BYTES EIGEN_IDEAL_MAX_ALIGN_BYTES
180
+ #else
181
+ #define EIGEN_DEFAULT_ALIGN_BYTES EIGEN_MAX_ALIGN_BYTES
182
+ #endif
183
+
184
+
185
+ #ifndef EIGEN_UNALIGNED_VECTORIZE
186
+ #define EIGEN_UNALIGNED_VECTORIZE 1
187
+ #endif
188
+
189
+ //----------------------------------------------------------------------
190
+
191
+ // if alignment is disabled, then disable vectorization. Note: EIGEN_MAX_ALIGN_BYTES is the proper check, it takes into
192
+ // account both the user's will (EIGEN_MAX_ALIGN_BYTES,EIGEN_DONT_ALIGN) and our own platform checks
193
+ #if EIGEN_MAX_ALIGN_BYTES==0
194
+ #ifndef EIGEN_DONT_VECTORIZE
195
+ #define EIGEN_DONT_VECTORIZE
196
+ #endif
197
+ #endif
198
+
199
+
200
+ // The following (except #include <malloc.h> and _M_IX86_FP ??) can likely be
201
+ // removed as gcc 4.1 and msvc 2008 are not supported anyways.
202
+ #if EIGEN_COMP_MSVC
203
+ #include <malloc.h> // for _aligned_malloc -- need it regardless of whether vectorization is enabled
204
+ #if (EIGEN_COMP_MSVC >= 1500) // 2008 or later
205
+ // a user reported that in 64-bit mode, MSVC doesn't care to define _M_IX86_FP.
206
+ #if (defined(_M_IX86_FP) && (_M_IX86_FP >= 2)) || EIGEN_ARCH_x86_64
207
+ #define EIGEN_SSE2_ON_MSVC_2008_OR_LATER
208
+ #endif
209
+ #endif
210
+ #else
211
+ #if (defined __SSE2__) && ( (!EIGEN_COMP_GNUC) || EIGEN_COMP_ICC || EIGEN_GNUC_AT_LEAST(4,2) )
212
+ #define EIGEN_SSE2_ON_NON_MSVC_BUT_NOT_OLD_GCC
213
+ #endif
214
+ #endif
215
+
216
+ #if !(defined(EIGEN_DONT_VECTORIZE) || defined(EIGEN_GPUCC))
217
+
218
+ #if defined (EIGEN_SSE2_ON_NON_MSVC_BUT_NOT_OLD_GCC) || defined(EIGEN_SSE2_ON_MSVC_2008_OR_LATER)
219
+
220
+ // Defines symbols for compile-time detection of which instructions are
221
+ // used.
222
+ // EIGEN_VECTORIZE_YY is defined if and only if the instruction set YY is used
223
+ #define EIGEN_VECTORIZE
224
+ #define EIGEN_VECTORIZE_SSE
225
+ #define EIGEN_VECTORIZE_SSE2
226
+
227
+ // Detect sse3/ssse3/sse4:
228
+ // gcc and icc defines __SSE3__, ...
229
+ // there is no way to know about this on msvc. You can define EIGEN_VECTORIZE_SSE* if you
230
+ // want to force the use of those instructions with msvc.
231
+ #ifdef __SSE3__
232
+ #define EIGEN_VECTORIZE_SSE3
233
+ #endif
234
+ #ifdef __SSSE3__
235
+ #define EIGEN_VECTORIZE_SSSE3
236
+ #endif
237
+ #ifdef __SSE4_1__
238
+ #define EIGEN_VECTORIZE_SSE4_1
239
+ #endif
240
+ #ifdef __SSE4_2__
241
+ #define EIGEN_VECTORIZE_SSE4_2
242
+ #endif
243
+ #ifdef __AVX__
244
+ #ifndef EIGEN_USE_SYCL
245
+ #define EIGEN_VECTORIZE_AVX
246
+ #endif
247
+ #define EIGEN_VECTORIZE_SSE3
248
+ #define EIGEN_VECTORIZE_SSSE3
249
+ #define EIGEN_VECTORIZE_SSE4_1
250
+ #define EIGEN_VECTORIZE_SSE4_2
251
+ #endif
252
+ #ifdef __AVX2__
253
+ #ifndef EIGEN_USE_SYCL
254
+ #define EIGEN_VECTORIZE_AVX2
255
+ #define EIGEN_VECTORIZE_AVX
256
+ #endif
257
+ #define EIGEN_VECTORIZE_SSE3
258
+ #define EIGEN_VECTORIZE_SSSE3
259
+ #define EIGEN_VECTORIZE_SSE4_1
260
+ #define EIGEN_VECTORIZE_SSE4_2
261
+ #endif
262
+ #if defined(__FMA__) || (EIGEN_COMP_MSVC && defined(__AVX2__))
263
+ // MSVC does not expose a switch dedicated for FMA
264
+ // For MSVC, AVX2 => FMA
265
+ #define EIGEN_VECTORIZE_FMA
266
+ #endif
267
+ #if defined(__AVX512F__)
268
+ #ifndef EIGEN_VECTORIZE_FMA
269
+ #if EIGEN_COMP_GNUC
270
+ #error Please add -mfma to your compiler flags: compiling with -mavx512f alone without SSE/AVX FMA is not supported (bug 1638).
271
+ #else
272
+ #error Please enable FMA in your compiler flags (e.g. -mfma): compiling with AVX512 alone without SSE/AVX FMA is not supported (bug 1638).
273
+ #endif
274
+ #endif
275
+ #ifndef EIGEN_USE_SYCL
276
+ #define EIGEN_VECTORIZE_AVX512
277
+ #define EIGEN_VECTORIZE_AVX2
278
+ #define EIGEN_VECTORIZE_AVX
279
+ #endif
280
+ #define EIGEN_VECTORIZE_FMA
281
+ #define EIGEN_VECTORIZE_SSE3
282
+ #define EIGEN_VECTORIZE_SSSE3
283
+ #define EIGEN_VECTORIZE_SSE4_1
284
+ #define EIGEN_VECTORIZE_SSE4_2
285
+ #ifndef EIGEN_USE_SYCL
286
+ #ifdef __AVX512DQ__
287
+ #define EIGEN_VECTORIZE_AVX512DQ
288
+ #endif
289
+ #ifdef __AVX512ER__
290
+ #define EIGEN_VECTORIZE_AVX512ER
291
+ #endif
292
+ #ifdef __AVX512BF16__
293
+ #define EIGEN_VECTORIZE_AVX512BF16
294
+ #endif
295
+ #endif
296
+ #endif
297
+
298
+ // Disable AVX support on broken xcode versions
299
+ #if defined(__apple_build_version__) && (__apple_build_version__ == 11000033 ) && ( __MAC_OS_X_VERSION_MIN_REQUIRED == 101500 )
300
+ // A nasty bug in the clang compiler shipped with xcode in a common compilation situation
301
+ // when XCode 11.0 and Mac deployment target macOS 10.15 is https://trac.macports.org/ticket/58776#no1
302
+ #ifdef EIGEN_VECTORIZE_AVX
303
+ #undef EIGEN_VECTORIZE_AVX
304
+ #warning "Disabling AVX support: clang compiler shipped with XCode 11.[012] generates broken assembly with -macosx-version-min=10.15 and AVX enabled. "
305
+ #ifdef EIGEN_VECTORIZE_AVX2
306
+ #undef EIGEN_VECTORIZE_AVX2
307
+ #endif
308
+ #ifdef EIGEN_VECTORIZE_FMA
309
+ #undef EIGEN_VECTORIZE_FMA
310
+ #endif
311
+ #ifdef EIGEN_VECTORIZE_AVX512
312
+ #undef EIGEN_VECTORIZE_AVX512
313
+ #endif
314
+ #ifdef EIGEN_VECTORIZE_AVX512DQ
315
+ #undef EIGEN_VECTORIZE_AVX512DQ
316
+ #endif
317
+ #ifdef EIGEN_VECTORIZE_AVX512ER
318
+ #undef EIGEN_VECTORIZE_AVX512ER
319
+ #endif
320
+ #endif
321
+ // NOTE: Confirmed test failures in XCode 11.0, and XCode 11.2 with -macosx-version-min=10.15 and AVX
322
+ // NOTE using -macosx-version-min=10.15 with Xcode 11.0 results in runtime segmentation faults in many tests, 11.2 produce core dumps in 3 tests
323
+ // NOTE using -macosx-version-min=10.14 produces functioning and passing tests in all cases
324
+ // NOTE __clang_version__ "11.0.0 (clang-1100.0.33.8)" XCode 11.0 <- Produces many segfault and core dumping tests
325
+ // with -macosx-version-min=10.15 and AVX
326
+ // NOTE __clang_version__ "11.0.0 (clang-1100.0.33.12)" XCode 11.2 <- Produces 3 core dumping tests with
327
+ // -macosx-version-min=10.15 and AVX
328
+ #endif
329
+
330
+ // include files
331
+
332
+ // This extern "C" works around a MINGW-w64 compilation issue
333
+ // https://sourceforge.net/tracker/index.php?func=detail&aid=3018394&group_id=202880&atid=983354
334
+ // In essence, intrin.h is included by windows.h and also declares intrinsics (just as emmintrin.h etc. below do).
335
+ // However, intrin.h uses an extern "C" declaration, and g++ thus complains of duplicate declarations
336
+ // with conflicting linkage. The linkage for intrinsics doesn't matter, but at that stage the compiler doesn't know;
337
+ // so, to avoid compile errors when windows.h is included after Eigen/Core, ensure intrinsics are extern "C" here too.
338
+ // notice that since these are C headers, the extern "C" is theoretically needed anyways.
339
+ extern "C" {
340
+ // In theory we should only include immintrin.h and not the other *mmintrin.h header files directly.
341
+ // Doing so triggers some issues with ICC. However old gcc versions seems to not have this file, thus:
342
+ #if EIGEN_COMP_ICC >= 1110
343
+ #include <immintrin.h>
344
+ #else
345
+ #include <mmintrin.h>
346
+ #include <emmintrin.h>
347
+ #include <xmmintrin.h>
348
+ #ifdef EIGEN_VECTORIZE_SSE3
349
+ #include <pmmintrin.h>
350
+ #endif
351
+ #ifdef EIGEN_VECTORIZE_SSSE3
352
+ #include <tmmintrin.h>
353
+ #endif
354
+ #ifdef EIGEN_VECTORIZE_SSE4_1
355
+ #include <smmintrin.h>
356
+ #endif
357
+ #ifdef EIGEN_VECTORIZE_SSE4_2
358
+ #include <nmmintrin.h>
359
+ #endif
360
+ #if defined(EIGEN_VECTORIZE_AVX) || defined(EIGEN_VECTORIZE_AVX512)
361
+ #include <immintrin.h>
362
+ #endif
363
+ #endif
364
+ } // end extern "C"
365
+
366
+ #elif defined __VSX__
367
+
368
+ #define EIGEN_VECTORIZE
369
+ #define EIGEN_VECTORIZE_VSX
370
+ #include <altivec.h>
371
+ // We need to #undef all these ugly tokens defined in <altivec.h>
372
+ // => use __vector instead of vector
373
+ #undef bool
374
+ #undef vector
375
+ #undef pixel
376
+
377
+ #elif defined __ALTIVEC__
378
+
379
+ #define EIGEN_VECTORIZE
380
+ #define EIGEN_VECTORIZE_ALTIVEC
381
+ #include <altivec.h>
382
+ // We need to #undef all these ugly tokens defined in <altivec.h>
383
+ // => use __vector instead of vector
384
+ #undef bool
385
+ #undef vector
386
+ #undef pixel
387
+
388
+ #elif ((defined __ARM_NEON) || (defined __ARM_NEON__)) && !(defined EIGEN_ARM64_USE_SVE)
389
+
390
+ #define EIGEN_VECTORIZE
391
+ #define EIGEN_VECTORIZE_NEON
392
+ #include <arm_neon.h>
393
+
394
+ // We currently require SVE to be enabled explicitly via EIGEN_ARM64_USE_SVE and
395
+ // will not select the backend automatically
396
+ #elif (defined __ARM_FEATURE_SVE) && (defined EIGEN_ARM64_USE_SVE)
397
+
398
+ #define EIGEN_VECTORIZE
399
+ #define EIGEN_VECTORIZE_SVE
400
+ #include <arm_sve.h>
401
+
402
+ // Since we depend on knowing SVE vector lengths at compile-time, we need
403
+ // to ensure a fixed lengths is set
404
+ #if defined __ARM_FEATURE_SVE_BITS
405
+ #define EIGEN_ARM64_SVE_VL __ARM_FEATURE_SVE_BITS
406
+ #else
407
+ #error "Eigen requires a fixed SVE lector length but EIGEN_ARM64_SVE_VL is not set."
408
+ #endif
409
+
410
+ #elif (defined __s390x__ && defined __VEC__)
411
+
412
+ #define EIGEN_VECTORIZE
413
+ #define EIGEN_VECTORIZE_ZVECTOR
414
+ #include <vecintrin.h>
415
+
416
+ #elif defined __mips_msa
417
+
418
+ // Limit MSA optimizations to little-endian CPUs for now.
419
+ // TODO: Perhaps, eventually support MSA optimizations on big-endian CPUs?
420
+ #if defined(__BYTE_ORDER__) && (__BYTE_ORDER__ == __ORDER_LITTLE_ENDIAN__)
421
+ #if defined(__LP64__)
422
+ #define EIGEN_MIPS_64
423
+ #else
424
+ #define EIGEN_MIPS_32
425
+ #endif
426
+ #define EIGEN_VECTORIZE
427
+ #define EIGEN_VECTORIZE_MSA
428
+ #include <msa.h>
429
+ #endif
430
+
431
+ #endif
432
+ #endif
433
+
434
+ // Following the Arm ACLE arm_neon.h should also include arm_fp16.h but not all
435
+ // compilers seem to follow this. We therefore include it explicitly.
436
+ // See also: https://bugs.llvm.org/show_bug.cgi?id=47955
437
+ #if defined(EIGEN_HAS_ARM64_FP16_SCALAR_ARITHMETIC)
438
+ #include <arm_fp16.h>
439
+ #endif
440
+
441
+ #if defined(__F16C__) && (!defined(EIGEN_GPUCC) && (!defined(EIGEN_COMP_CLANG) || EIGEN_COMP_CLANG>=380))
442
+ // We can use the optimized fp16 to float and float to fp16 conversion routines
443
+ #define EIGEN_HAS_FP16_C
444
+
445
+ #if defined(EIGEN_COMP_CLANG)
446
+ // Workaround for clang: The FP16C intrinsics for clang are included by
447
+ // immintrin.h, as opposed to emmintrin.h as suggested by Intel:
448
+ // https://software.intel.com/sites/landingpage/IntrinsicsGuide/#othertechs=FP16C&expand=1711
449
+ #include <immintrin.h>
450
+ #endif
451
+ #endif
452
+
453
+ #if defined EIGEN_CUDACC
454
+ #define EIGEN_VECTORIZE_GPU
455
+ #include <vector_types.h>
456
+ #if EIGEN_CUDA_SDK_VER >= 70500
457
+ #define EIGEN_HAS_CUDA_FP16
458
+ #endif
459
+ #endif
460
+
461
+ #if defined(EIGEN_HAS_CUDA_FP16)
462
+ #include <cuda_runtime_api.h>
463
+ #include <cuda_fp16.h>
464
+ #endif
465
+
466
+ #if defined(EIGEN_HIPCC)
467
+ #define EIGEN_VECTORIZE_GPU
468
+ #include <hip/hip_vector_types.h>
469
+ #define EIGEN_HAS_HIP_FP16
470
+ #include <hip/hip_fp16.h>
471
+ #endif
472
+
473
+
474
+ /** \brief Namespace containing all symbols from the %Eigen library. */
475
+ namespace Eigen {
476
+
477
+ inline static const char *SimdInstructionSetsInUse(void) {
478
+ #if defined(EIGEN_VECTORIZE_AVX512)
479
+ return "AVX512, FMA, AVX2, AVX, SSE, SSE2, SSE3, SSSE3, SSE4.1, SSE4.2";
480
+ #elif defined(EIGEN_VECTORIZE_AVX)
481
+ return "AVX SSE, SSE2, SSE3, SSSE3, SSE4.1, SSE4.2";
482
+ #elif defined(EIGEN_VECTORIZE_SSE4_2)
483
+ return "SSE, SSE2, SSE3, SSSE3, SSE4.1, SSE4.2";
484
+ #elif defined(EIGEN_VECTORIZE_SSE4_1)
485
+ return "SSE, SSE2, SSE3, SSSE3, SSE4.1";
486
+ #elif defined(EIGEN_VECTORIZE_SSSE3)
487
+ return "SSE, SSE2, SSE3, SSSE3";
488
+ #elif defined(EIGEN_VECTORIZE_SSE3)
489
+ return "SSE, SSE2, SSE3";
490
+ #elif defined(EIGEN_VECTORIZE_SSE2)
491
+ return "SSE, SSE2";
492
+ #elif defined(EIGEN_VECTORIZE_ALTIVEC)
493
+ return "AltiVec";
494
+ #elif defined(EIGEN_VECTORIZE_VSX)
495
+ return "VSX";
496
+ #elif defined(EIGEN_VECTORIZE_NEON)
497
+ return "ARM NEON";
498
+ #elif defined(EIGEN_VECTORIZE_SVE)
499
+ return "ARM SVE";
500
+ #elif defined(EIGEN_VECTORIZE_ZVECTOR)
501
+ return "S390X ZVECTOR";
502
+ #elif defined(EIGEN_VECTORIZE_MSA)
503
+ return "MIPS MSA";
504
+ #else
505
+ return "None";
506
+ #endif
507
+ }
508
+
509
+ } // end namespace Eigen
510
+
511
+
512
+ #endif // EIGEN_CONFIGURE_VECTORIZATION_H
@@ -3,6 +3,7 @@
3
3
  //
4
4
  // Copyright (C) 2008-2015 Gael Guennebaud <gael.guennebaud@inria.fr>
5
5
  // Copyright (C) 2007-2009 Benoit Jacob <jacob.benoit.1@gmail.com>
6
+ // Copyright (C) 2020, Arm Limited and Contributors
6
7
  //
7
8
  // This Source Code Form is subject to the terms of the Mozilla
8
9
  // Public License v. 2.0. If a copy of the MPL was not distributed
@@ -25,6 +26,10 @@ const int Dynamic = -1;
25
26
  */
26
27
  const int DynamicIndex = 0xffffff;
27
28
 
29
+ /** This value means that the increment to go from one value to another in a sequence is not constant for each step.
30
+ */
31
+ const int UndefinedIncr = 0xfffffe;
32
+
28
33
  /** This value means +Infinity; it is currently used only as the p parameter to MatrixBase::lpNorm<int>().
29
34
  * The value Infinity there means the L-infinity norm.
30
35
  */
@@ -152,7 +157,7 @@ const unsigned int DirectAccessBit = 0x40;
152
157
  /** \deprecated \ingroup flags
153
158
  *
154
159
  * means the first coefficient packet is guaranteed to be aligned.
155
- * An expression cannot has the AlignedBit without the PacketAccessBit flag.
160
+ * An expression cannot have the AlignedBit without the PacketAccessBit flag.
156
161
  * In other words, this means we are allow to perform an aligned packet access to the first element regardless
157
162
  * of the expression kind:
158
163
  * \code
@@ -250,12 +255,6 @@ enum AlignmentType {
250
255
  #endif
251
256
  };
252
257
 
253
- /** \ingroup enums
254
- * Enum used by DenseBase::corner() in Eigen2 compatibility mode. */
255
- // FIXME after the corner() API change, this was not needed anymore, except by AlignedBox
256
- // TODO: find out what to do with that. Adapt the AlignedBox API ?
257
- enum CornerType { TopLeft, TopRight, BottomLeft, BottomRight };
258
-
259
258
  /** \ingroup enums
260
259
  * Enum containing possible values for the \p Direction parameter of
261
260
  * Reverse, PartialReduxExpr and VectorwiseOp. */
@@ -330,9 +329,20 @@ enum StorageOptions {
330
329
  * Enum for specifying whether to apply or solve on the left or right. */
331
330
  enum SideType {
332
331
  /** Apply transformation on the left. */
333
- OnTheLeft = 1,
332
+ OnTheLeft = 1,
334
333
  /** Apply transformation on the right. */
335
- OnTheRight = 2
334
+ OnTheRight = 2
335
+ };
336
+
337
+ /** \ingroup enums
338
+ * Enum for specifying NaN-propagation behavior, e.g. for coeff-wise min/max. */
339
+ enum NaNPropagationOptions {
340
+ /** Implementation defined behavior if NaNs are present. */
341
+ PropagateFast = 0,
342
+ /** Always propagate NaNs. */
343
+ PropagateNaN,
344
+ /** Always propagate not-NaNs. */
345
+ PropagateNumbers
336
346
  };
337
347
 
338
348
  /* the following used to be written as:
@@ -464,6 +474,8 @@ namespace Architecture
464
474
  AltiVec = 0x2,
465
475
  VSX = 0x3,
466
476
  NEON = 0x4,
477
+ MSA = 0x5,
478
+ SVE = 0x6,
467
479
  #if defined EIGEN_VECTORIZE_SSE
468
480
  Target = SSE
469
481
  #elif defined EIGEN_VECTORIZE_ALTIVEC
@@ -472,6 +484,10 @@ namespace Architecture
472
484
  Target = VSX
473
485
  #elif defined EIGEN_VECTORIZE_NEON
474
486
  Target = NEON
487
+ #elif defined EIGEN_VECTORIZE_SVE
488
+ Target = SVE
489
+ #elif defined EIGEN_VECTORIZE_MSA
490
+ Target = MSA
475
491
  #else
476
492
  Target = Generic
477
493
  #endif
@@ -4,7 +4,6 @@
4
4
  #ifdef _MSC_VER
5
5
  // 4100 - unreferenced formal parameter (occurred e.g. in aligned_allocator::destroy(pointer p))
6
6
  // 4101 - unreferenced local variable
7
- // 4127 - conditional expression is constant
8
7
  // 4181 - qualifier applied to reference type ignored
9
8
  // 4211 - nonstandard extension used : redefined extern to static
10
9
  // 4244 - 'argument' : conversion from 'type1' to 'type2', possible loss of data
@@ -20,7 +19,7 @@
20
19
  #ifndef EIGEN_PERMANENTLY_DISABLE_STUPID_WARNINGS
21
20
  #pragma warning( push )
22
21
  #endif
23
- #pragma warning( disable : 4100 4101 4127 4181 4211 4244 4273 4324 4503 4512 4522 4700 4714 4717 4800)
22
+ #pragma warning( disable : 4100 4101 4181 4211 4244 4273 4324 4503 4512 4522 4700 4714 4717 4800)
24
23
 
25
24
  #elif defined __INTEL_COMPILER
26
25
  // 2196 - routine is both "inline" and "noinline" ("noinline" assumed)
@@ -42,6 +41,17 @@
42
41
  #pragma clang diagnostic push
43
42
  #endif
44
43
  #pragma clang diagnostic ignored "-Wconstant-logical-operand"
44
+ #if __clang_major__ >= 3 && __clang_minor__ >= 5
45
+ #pragma clang diagnostic ignored "-Wabsolute-value"
46
+ #endif
47
+ #if __clang_major__ >= 10
48
+ #pragma clang diagnostic ignored "-Wimplicit-int-float-conversion"
49
+ #endif
50
+ #if ( defined(__ALTIVEC__) || defined(__VSX__) ) && __cplusplus < 201103L
51
+ // warning: generic selections are a C11-specific feature
52
+ // ignoring warnings thrown at vec_ctf in Altivec/PacketMath.h
53
+ #pragma clang diagnostic ignored "-Wc11-extensions"
54
+ #endif
45
55
 
46
56
  #elif defined __GNUC__
47
57
 
@@ -57,10 +67,14 @@
57
67
  #if __GNUC__>=6
58
68
  #pragma GCC diagnostic ignored "-Wignored-attributes"
59
69
  #endif
60
-
70
+ #if __GNUC__==7
71
+ // See: https://gcc.gnu.org/bugzilla/show_bug.cgi?id=89325
72
+ #pragma GCC diagnostic ignored "-Wattributes"
73
+ #endif
61
74
  #endif
62
75
 
63
76
  #if defined __NVCC__
77
+ #pragma diag_suppress boolean_controlling_expr_is_constant
64
78
  // Disable the "statement is unreachable" message
65
79
  #pragma diag_suppress code_is_unreachable
66
80
  // Disable the "dynamic initialization in unreachable code" message
@@ -78,6 +92,15 @@
78
92
  #pragma diag_suppress 2671
79
93
  #pragma diag_suppress 2735
80
94
  #pragma diag_suppress 2737
95
+ #pragma diag_suppress 2739
81
96
  #endif
82
97
 
98
+ #else
99
+ // warnings already disabled:
100
+ # ifndef EIGEN_WARNINGS_DISABLED_2
101
+ # define EIGEN_WARNINGS_DISABLED_2
102
+ # elif defined(EIGEN_INTERNAL_DEBUGGING)
103
+ # error "Do not include \"DisableStupidWarnings.h\" recursively more than twice!"
104
+ # endif
105
+
83
106
  #endif // not EIGEN_WARNINGS_DISABLED