tomoto 0.2.2 → 0.3.0

Sign up to get free protection for your applications and to get access to all the features.
Files changed (369) hide show
  1. checksums.yaml +4 -4
  2. data/CHANGELOG.md +10 -0
  3. data/README.md +8 -10
  4. data/ext/tomoto/ct.cpp +11 -11
  5. data/ext/tomoto/dmr.cpp +14 -13
  6. data/ext/tomoto/dt.cpp +14 -14
  7. data/ext/tomoto/extconf.rb +7 -5
  8. data/ext/tomoto/gdmr.cpp +7 -7
  9. data/ext/tomoto/hdp.cpp +9 -9
  10. data/ext/tomoto/hlda.cpp +13 -13
  11. data/ext/tomoto/hpa.cpp +5 -5
  12. data/ext/tomoto/lda.cpp +42 -39
  13. data/ext/tomoto/llda.cpp +6 -6
  14. data/ext/tomoto/mglda.cpp +15 -15
  15. data/ext/tomoto/pa.cpp +6 -6
  16. data/ext/tomoto/plda.cpp +6 -6
  17. data/ext/tomoto/slda.cpp +8 -8
  18. data/ext/tomoto/{ext.cpp → tomoto.cpp} +8 -8
  19. data/ext/tomoto/utils.h +16 -70
  20. data/lib/tomoto/version.rb +1 -1
  21. data/lib/tomoto.rb +5 -1
  22. data/vendor/EigenRand/EigenRand/Core.h +10 -10
  23. data/vendor/EigenRand/EigenRand/Dists/Basic.h +208 -9
  24. data/vendor/EigenRand/EigenRand/Dists/Discrete.h +52 -31
  25. data/vendor/EigenRand/EigenRand/Dists/GammaPoisson.h +9 -8
  26. data/vendor/EigenRand/EigenRand/Dists/NormalExp.h +28 -21
  27. data/vendor/EigenRand/EigenRand/EigenRand +11 -6
  28. data/vendor/EigenRand/EigenRand/Macro.h +13 -7
  29. data/vendor/EigenRand/EigenRand/MorePacketMath.h +348 -740
  30. data/vendor/EigenRand/EigenRand/MvDists/Multinomial.h +5 -3
  31. data/vendor/EigenRand/EigenRand/MvDists/MvNormal.h +9 -3
  32. data/vendor/EigenRand/EigenRand/PacketFilter.h +11 -253
  33. data/vendor/EigenRand/EigenRand/PacketRandomEngine.h +21 -47
  34. data/vendor/EigenRand/EigenRand/RandUtils.h +50 -344
  35. data/vendor/EigenRand/EigenRand/arch/AVX/MorePacketMath.h +619 -0
  36. data/vendor/EigenRand/EigenRand/arch/AVX/PacketFilter.h +149 -0
  37. data/vendor/EigenRand/EigenRand/arch/AVX/RandUtils.h +228 -0
  38. data/vendor/EigenRand/EigenRand/arch/NEON/MorePacketMath.h +473 -0
  39. data/vendor/EigenRand/EigenRand/arch/NEON/PacketFilter.h +142 -0
  40. data/vendor/EigenRand/EigenRand/arch/NEON/RandUtils.h +126 -0
  41. data/vendor/EigenRand/EigenRand/arch/SSE/MorePacketMath.h +501 -0
  42. data/vendor/EigenRand/EigenRand/arch/SSE/PacketFilter.h +133 -0
  43. data/vendor/EigenRand/EigenRand/arch/SSE/RandUtils.h +120 -0
  44. data/vendor/EigenRand/EigenRand/doc.h +24 -12
  45. data/vendor/EigenRand/README.md +57 -4
  46. data/vendor/eigen/COPYING.APACHE +203 -0
  47. data/vendor/eigen/COPYING.BSD +1 -1
  48. data/vendor/eigen/COPYING.MINPACK +51 -52
  49. data/vendor/eigen/Eigen/Cholesky +0 -1
  50. data/vendor/eigen/Eigen/Core +112 -265
  51. data/vendor/eigen/Eigen/Eigenvalues +2 -3
  52. data/vendor/eigen/Eigen/Geometry +5 -8
  53. data/vendor/eigen/Eigen/Householder +0 -1
  54. data/vendor/eigen/Eigen/Jacobi +0 -1
  55. data/vendor/eigen/Eigen/KLUSupport +41 -0
  56. data/vendor/eigen/Eigen/LU +2 -5
  57. data/vendor/eigen/Eigen/OrderingMethods +0 -3
  58. data/vendor/eigen/Eigen/PaStiXSupport +1 -0
  59. data/vendor/eigen/Eigen/PardisoSupport +0 -0
  60. data/vendor/eigen/Eigen/QR +2 -3
  61. data/vendor/eigen/Eigen/QtAlignedMalloc +0 -1
  62. data/vendor/eigen/Eigen/SVD +0 -1
  63. data/vendor/eigen/Eigen/Sparse +0 -2
  64. data/vendor/eigen/Eigen/SparseCholesky +0 -8
  65. data/vendor/eigen/Eigen/SparseLU +4 -0
  66. data/vendor/eigen/Eigen/SparseQR +0 -1
  67. data/vendor/eigen/Eigen/src/Cholesky/LDLT.h +42 -27
  68. data/vendor/eigen/Eigen/src/Cholesky/LLT.h +39 -23
  69. data/vendor/eigen/Eigen/src/CholmodSupport/CholmodSupport.h +90 -47
  70. data/vendor/eigen/Eigen/src/Core/ArithmeticSequence.h +413 -0
  71. data/vendor/eigen/Eigen/src/Core/Array.h +99 -11
  72. data/vendor/eigen/Eigen/src/Core/ArrayBase.h +3 -3
  73. data/vendor/eigen/Eigen/src/Core/ArrayWrapper.h +21 -21
  74. data/vendor/eigen/Eigen/src/Core/Assign.h +1 -1
  75. data/vendor/eigen/Eigen/src/Core/AssignEvaluator.h +125 -50
  76. data/vendor/eigen/Eigen/src/Core/Assign_MKL.h +10 -10
  77. data/vendor/eigen/Eigen/src/Core/BandMatrix.h +16 -16
  78. data/vendor/eigen/Eigen/src/Core/Block.h +56 -60
  79. data/vendor/eigen/Eigen/src/Core/BooleanRedux.h +29 -31
  80. data/vendor/eigen/Eigen/src/Core/CommaInitializer.h +7 -3
  81. data/vendor/eigen/Eigen/src/Core/CoreEvaluators.h +325 -272
  82. data/vendor/eigen/Eigen/src/Core/CoreIterators.h +5 -0
  83. data/vendor/eigen/Eigen/src/Core/CwiseBinaryOp.h +21 -22
  84. data/vendor/eigen/Eigen/src/Core/CwiseNullaryOp.h +153 -18
  85. data/vendor/eigen/Eigen/src/Core/CwiseUnaryOp.h +6 -6
  86. data/vendor/eigen/Eigen/src/Core/CwiseUnaryView.h +14 -10
  87. data/vendor/eigen/Eigen/src/Core/DenseBase.h +132 -42
  88. data/vendor/eigen/Eigen/src/Core/DenseCoeffsBase.h +25 -21
  89. data/vendor/eigen/Eigen/src/Core/DenseStorage.h +153 -71
  90. data/vendor/eigen/Eigen/src/Core/Diagonal.h +21 -23
  91. data/vendor/eigen/Eigen/src/Core/DiagonalMatrix.h +50 -2
  92. data/vendor/eigen/Eigen/src/Core/DiagonalProduct.h +1 -1
  93. data/vendor/eigen/Eigen/src/Core/Dot.h +10 -10
  94. data/vendor/eigen/Eigen/src/Core/EigenBase.h +10 -9
  95. data/vendor/eigen/Eigen/src/Core/ForceAlignedAccess.h +8 -4
  96. data/vendor/eigen/Eigen/src/Core/Fuzzy.h +3 -3
  97. data/vendor/eigen/Eigen/src/Core/GeneralProduct.h +20 -10
  98. data/vendor/eigen/Eigen/src/Core/GenericPacketMath.h +599 -152
  99. data/vendor/eigen/Eigen/src/Core/GlobalFunctions.h +40 -33
  100. data/vendor/eigen/Eigen/src/Core/IO.h +40 -7
  101. data/vendor/eigen/Eigen/src/Core/IndexedView.h +237 -0
  102. data/vendor/eigen/Eigen/src/Core/Inverse.h +9 -10
  103. data/vendor/eigen/Eigen/src/Core/Map.h +7 -7
  104. data/vendor/eigen/Eigen/src/Core/MapBase.h +10 -3
  105. data/vendor/eigen/Eigen/src/Core/MathFunctions.h +767 -125
  106. data/vendor/eigen/Eigen/src/Core/MathFunctionsImpl.h +118 -19
  107. data/vendor/eigen/Eigen/src/Core/Matrix.h +131 -25
  108. data/vendor/eigen/Eigen/src/Core/MatrixBase.h +21 -3
  109. data/vendor/eigen/Eigen/src/Core/NestByValue.h +25 -50
  110. data/vendor/eigen/Eigen/src/Core/NoAlias.h +4 -3
  111. data/vendor/eigen/Eigen/src/Core/NumTraits.h +107 -20
  112. data/vendor/eigen/Eigen/src/Core/PartialReduxEvaluator.h +232 -0
  113. data/vendor/eigen/Eigen/src/Core/PermutationMatrix.h +3 -31
  114. data/vendor/eigen/Eigen/src/Core/PlainObjectBase.h +152 -59
  115. data/vendor/eigen/Eigen/src/Core/Product.h +30 -25
  116. data/vendor/eigen/Eigen/src/Core/ProductEvaluators.h +192 -125
  117. data/vendor/eigen/Eigen/src/Core/Random.h +37 -1
  118. data/vendor/eigen/Eigen/src/Core/Redux.h +180 -170
  119. data/vendor/eigen/Eigen/src/Core/Ref.h +121 -23
  120. data/vendor/eigen/Eigen/src/Core/Replicate.h +8 -8
  121. data/vendor/eigen/Eigen/src/Core/Reshaped.h +454 -0
  122. data/vendor/eigen/Eigen/src/Core/ReturnByValue.h +7 -5
  123. data/vendor/eigen/Eigen/src/Core/Reverse.h +18 -12
  124. data/vendor/eigen/Eigen/src/Core/Select.h +8 -6
  125. data/vendor/eigen/Eigen/src/Core/SelfAdjointView.h +33 -20
  126. data/vendor/eigen/Eigen/src/Core/Solve.h +14 -14
  127. data/vendor/eigen/Eigen/src/Core/SolveTriangular.h +16 -16
  128. data/vendor/eigen/Eigen/src/Core/SolverBase.h +41 -3
  129. data/vendor/eigen/Eigen/src/Core/StableNorm.h +100 -70
  130. data/vendor/eigen/Eigen/src/Core/StlIterators.h +463 -0
  131. data/vendor/eigen/Eigen/src/Core/Stride.h +9 -4
  132. data/vendor/eigen/Eigen/src/Core/Swap.h +5 -4
  133. data/vendor/eigen/Eigen/src/Core/Transpose.h +88 -27
  134. data/vendor/eigen/Eigen/src/Core/Transpositions.h +26 -47
  135. data/vendor/eigen/Eigen/src/Core/TriangularMatrix.h +93 -75
  136. data/vendor/eigen/Eigen/src/Core/VectorBlock.h +5 -5
  137. data/vendor/eigen/Eigen/src/Core/VectorwiseOp.h +159 -70
  138. data/vendor/eigen/Eigen/src/Core/Visitor.h +137 -29
  139. data/vendor/eigen/Eigen/src/Core/arch/AVX/Complex.h +50 -129
  140. data/vendor/eigen/Eigen/src/Core/arch/AVX/MathFunctions.h +126 -337
  141. data/vendor/eigen/Eigen/src/Core/arch/AVX/PacketMath.h +1092 -155
  142. data/vendor/eigen/Eigen/src/Core/arch/AVX/TypeCasting.h +65 -1
  143. data/vendor/eigen/Eigen/src/Core/arch/AVX512/Complex.h +422 -0
  144. data/vendor/eigen/Eigen/src/Core/arch/AVX512/MathFunctions.h +207 -236
  145. data/vendor/eigen/Eigen/src/Core/arch/AVX512/PacketMath.h +1482 -495
  146. data/vendor/eigen/Eigen/src/Core/arch/AVX512/TypeCasting.h +89 -0
  147. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/Complex.h +152 -165
  148. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/MathFunctions.h +19 -251
  149. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/MatrixProduct.h +2937 -0
  150. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/MatrixProductCommon.h +221 -0
  151. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/MatrixProductMMA.h +629 -0
  152. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/PacketMath.h +2042 -392
  153. data/vendor/eigen/Eigen/src/Core/arch/CUDA/Complex.h +235 -80
  154. data/vendor/eigen/Eigen/src/Core/arch/Default/BFloat16.h +700 -0
  155. data/vendor/eigen/Eigen/src/Core/arch/Default/ConjHelper.h +102 -14
  156. data/vendor/eigen/Eigen/src/Core/arch/Default/GenericPacketMathFunctions.h +1649 -0
  157. data/vendor/eigen/Eigen/src/Core/arch/Default/GenericPacketMathFunctionsFwd.h +110 -0
  158. data/vendor/eigen/Eigen/src/Core/arch/Default/Half.h +942 -0
  159. data/vendor/eigen/Eigen/src/Core/arch/Default/Settings.h +1 -1
  160. data/vendor/eigen/Eigen/src/Core/arch/Default/TypeCasting.h +120 -0
  161. data/vendor/eigen/Eigen/src/Core/arch/{CUDA → GPU}/MathFunctions.h +16 -4
  162. data/vendor/eigen/Eigen/src/Core/arch/GPU/PacketMath.h +1685 -0
  163. data/vendor/eigen/Eigen/src/Core/arch/GPU/TypeCasting.h +80 -0
  164. data/vendor/eigen/Eigen/src/Core/arch/HIP/hcc/math_constants.h +23 -0
  165. data/vendor/eigen/Eigen/src/Core/arch/MSA/Complex.h +648 -0
  166. data/vendor/eigen/Eigen/src/Core/arch/MSA/MathFunctions.h +387 -0
  167. data/vendor/eigen/Eigen/src/Core/arch/MSA/PacketMath.h +1233 -0
  168. data/vendor/eigen/Eigen/src/Core/arch/NEON/Complex.h +313 -219
  169. data/vendor/eigen/Eigen/src/Core/arch/NEON/GeneralBlockPanelKernel.h +183 -0
  170. data/vendor/eigen/Eigen/src/Core/arch/NEON/MathFunctions.h +54 -70
  171. data/vendor/eigen/Eigen/src/Core/arch/NEON/PacketMath.h +4376 -549
  172. data/vendor/eigen/Eigen/src/Core/arch/NEON/TypeCasting.h +1419 -0
  173. data/vendor/eigen/Eigen/src/Core/arch/SSE/Complex.h +59 -179
  174. data/vendor/eigen/Eigen/src/Core/arch/SSE/MathFunctions.h +65 -428
  175. data/vendor/eigen/Eigen/src/Core/arch/SSE/PacketMath.h +893 -283
  176. data/vendor/eigen/Eigen/src/Core/arch/SSE/TypeCasting.h +65 -0
  177. data/vendor/eigen/Eigen/src/Core/arch/SVE/MathFunctions.h +44 -0
  178. data/vendor/eigen/Eigen/src/Core/arch/SVE/PacketMath.h +752 -0
  179. data/vendor/eigen/Eigen/src/Core/arch/SVE/TypeCasting.h +49 -0
  180. data/vendor/eigen/Eigen/src/Core/arch/SYCL/InteropHeaders.h +232 -0
  181. data/vendor/eigen/Eigen/src/Core/arch/SYCL/MathFunctions.h +301 -0
  182. data/vendor/eigen/Eigen/src/Core/arch/SYCL/PacketMath.h +670 -0
  183. data/vendor/eigen/Eigen/src/Core/arch/SYCL/SyclMemoryModel.h +694 -0
  184. data/vendor/eigen/Eigen/src/Core/arch/SYCL/TypeCasting.h +85 -0
  185. data/vendor/eigen/Eigen/src/Core/arch/ZVector/Complex.h +212 -183
  186. data/vendor/eigen/Eigen/src/Core/arch/ZVector/MathFunctions.h +101 -5
  187. data/vendor/eigen/Eigen/src/Core/arch/ZVector/PacketMath.h +510 -395
  188. data/vendor/eigen/Eigen/src/Core/functors/AssignmentFunctors.h +11 -2
  189. data/vendor/eigen/Eigen/src/Core/functors/BinaryFunctors.h +112 -46
  190. data/vendor/eigen/Eigen/src/Core/functors/NullaryFunctors.h +31 -30
  191. data/vendor/eigen/Eigen/src/Core/functors/StlFunctors.h +32 -2
  192. data/vendor/eigen/Eigen/src/Core/functors/UnaryFunctors.h +355 -16
  193. data/vendor/eigen/Eigen/src/Core/products/GeneralBlockPanelKernel.h +1075 -586
  194. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrix.h +49 -24
  195. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrixTriangular.h +41 -35
  196. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h +6 -6
  197. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrix_BLAS.h +4 -2
  198. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixVector.h +382 -483
  199. data/vendor/eigen/Eigen/src/Core/products/Parallelizer.h +22 -5
  200. data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixMatrix.h +53 -30
  201. data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixMatrix_BLAS.h +16 -8
  202. data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixVector.h +8 -6
  203. data/vendor/eigen/Eigen/src/Core/products/SelfadjointProduct.h +4 -4
  204. data/vendor/eigen/Eigen/src/Core/products/SelfadjointRank2Update.h +5 -4
  205. data/vendor/eigen/Eigen/src/Core/products/TriangularMatrixMatrix.h +33 -27
  206. data/vendor/eigen/Eigen/src/Core/products/TriangularMatrixMatrix_BLAS.h +14 -12
  207. data/vendor/eigen/Eigen/src/Core/products/TriangularSolverMatrix.h +36 -34
  208. data/vendor/eigen/Eigen/src/Core/products/TriangularSolverMatrix_BLAS.h +8 -4
  209. data/vendor/eigen/Eigen/src/Core/products/TriangularSolverVector.h +13 -10
  210. data/vendor/eigen/Eigen/src/Core/util/BlasUtil.h +304 -119
  211. data/vendor/eigen/Eigen/src/Core/util/ConfigureVectorization.h +512 -0
  212. data/vendor/eigen/Eigen/src/Core/util/Constants.h +25 -9
  213. data/vendor/eigen/Eigen/src/Core/util/DisableStupidWarnings.h +26 -3
  214. data/vendor/eigen/Eigen/src/Core/util/ForwardDeclarations.h +29 -9
  215. data/vendor/eigen/Eigen/src/Core/util/IndexedViewHelper.h +186 -0
  216. data/vendor/eigen/Eigen/src/Core/util/IntegralConstant.h +272 -0
  217. data/vendor/eigen/Eigen/src/Core/util/MKL_support.h +8 -1
  218. data/vendor/eigen/Eigen/src/Core/util/Macros.h +709 -246
  219. data/vendor/eigen/Eigen/src/Core/util/Memory.h +222 -52
  220. data/vendor/eigen/Eigen/src/Core/util/Meta.h +355 -77
  221. data/vendor/eigen/Eigen/src/Core/util/ReenableStupidWarnings.h +5 -1
  222. data/vendor/eigen/Eigen/src/Core/util/ReshapedHelper.h +51 -0
  223. data/vendor/eigen/Eigen/src/Core/util/StaticAssert.h +8 -5
  224. data/vendor/eigen/Eigen/src/Core/util/SymbolicIndex.h +293 -0
  225. data/vendor/eigen/Eigen/src/Core/util/XprHelper.h +65 -30
  226. data/vendor/eigen/Eigen/src/Eigenvalues/ComplexEigenSolver.h +1 -1
  227. data/vendor/eigen/Eigen/src/Eigenvalues/ComplexSchur.h +7 -4
  228. data/vendor/eigen/Eigen/src/Eigenvalues/EigenSolver.h +2 -2
  229. data/vendor/eigen/Eigen/src/Eigenvalues/GeneralizedSelfAdjointEigenSolver.h +1 -1
  230. data/vendor/eigen/Eigen/src/Eigenvalues/HessenbergDecomposition.h +2 -2
  231. data/vendor/eigen/Eigen/src/Eigenvalues/MatrixBaseEigenvalues.h +2 -2
  232. data/vendor/eigen/Eigen/src/Eigenvalues/RealQZ.h +9 -6
  233. data/vendor/eigen/Eigen/src/Eigenvalues/RealSchur.h +21 -9
  234. data/vendor/eigen/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h +77 -43
  235. data/vendor/eigen/Eigen/src/Eigenvalues/Tridiagonalization.h +20 -15
  236. data/vendor/eigen/Eigen/src/Geometry/AlignedBox.h +99 -5
  237. data/vendor/eigen/Eigen/src/Geometry/AngleAxis.h +4 -4
  238. data/vendor/eigen/Eigen/src/Geometry/EulerAngles.h +3 -3
  239. data/vendor/eigen/Eigen/src/Geometry/Homogeneous.h +15 -11
  240. data/vendor/eigen/Eigen/src/Geometry/Hyperplane.h +1 -1
  241. data/vendor/eigen/Eigen/src/Geometry/OrthoMethods.h +3 -2
  242. data/vendor/eigen/Eigen/src/Geometry/ParametrizedLine.h +39 -2
  243. data/vendor/eigen/Eigen/src/Geometry/Quaternion.h +70 -14
  244. data/vendor/eigen/Eigen/src/Geometry/Rotation2D.h +3 -3
  245. data/vendor/eigen/Eigen/src/Geometry/Scaling.h +23 -5
  246. data/vendor/eigen/Eigen/src/Geometry/Transform.h +88 -67
  247. data/vendor/eigen/Eigen/src/Geometry/Translation.h +6 -12
  248. data/vendor/eigen/Eigen/src/Geometry/Umeyama.h +1 -1
  249. data/vendor/eigen/Eigen/src/Geometry/arch/Geometry_SIMD.h +168 -0
  250. data/vendor/eigen/Eigen/src/Householder/BlockHouseholder.h +9 -2
  251. data/vendor/eigen/Eigen/src/Householder/Householder.h +8 -4
  252. data/vendor/eigen/Eigen/src/Householder/HouseholderSequence.h +123 -48
  253. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/BasicPreconditioners.h +15 -15
  254. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/BiCGSTAB.h +7 -23
  255. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h +5 -22
  256. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/IncompleteCholesky.h +41 -47
  257. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/IncompleteLUT.h +51 -60
  258. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/IterativeSolverBase.h +70 -20
  259. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/LeastSquareConjugateGradient.h +2 -20
  260. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/SolveWithGuess.h +11 -9
  261. data/vendor/eigen/Eigen/src/Jacobi/Jacobi.h +31 -10
  262. data/vendor/eigen/Eigen/src/KLUSupport/KLUSupport.h +358 -0
  263. data/vendor/eigen/Eigen/src/LU/Determinant.h +35 -19
  264. data/vendor/eigen/Eigen/src/LU/FullPivLU.h +29 -43
  265. data/vendor/eigen/Eigen/src/LU/InverseImpl.h +25 -8
  266. data/vendor/eigen/Eigen/src/LU/PartialPivLU.h +71 -58
  267. data/vendor/eigen/Eigen/src/LU/arch/InverseSize4.h +351 -0
  268. data/vendor/eigen/Eigen/src/OrderingMethods/Amd.h +7 -17
  269. data/vendor/eigen/Eigen/src/OrderingMethods/Eigen_Colamd.h +297 -277
  270. data/vendor/eigen/Eigen/src/OrderingMethods/Ordering.h +6 -10
  271. data/vendor/eigen/Eigen/src/PaStiXSupport/PaStiXSupport.h +1 -1
  272. data/vendor/eigen/Eigen/src/PardisoSupport/PardisoSupport.h +9 -7
  273. data/vendor/eigen/Eigen/src/QR/ColPivHouseholderQR.h +41 -20
  274. data/vendor/eigen/Eigen/src/QR/CompleteOrthogonalDecomposition.h +100 -27
  275. data/vendor/eigen/Eigen/src/QR/FullPivHouseholderQR.h +59 -22
  276. data/vendor/eigen/Eigen/src/QR/HouseholderQR.h +48 -23
  277. data/vendor/eigen/Eigen/src/SPQRSupport/SuiteSparseQRSupport.h +25 -3
  278. data/vendor/eigen/Eigen/src/SVD/BDCSVD.h +183 -63
  279. data/vendor/eigen/Eigen/src/SVD/JacobiSVD.h +22 -14
  280. data/vendor/eigen/Eigen/src/SVD/SVDBase.h +83 -22
  281. data/vendor/eigen/Eigen/src/SVD/UpperBidiagonalization.h +3 -3
  282. data/vendor/eigen/Eigen/src/SparseCholesky/SimplicialCholesky.h +17 -9
  283. data/vendor/eigen/Eigen/src/SparseCholesky/SimplicialCholesky_impl.h +12 -37
  284. data/vendor/eigen/Eigen/src/SparseCore/AmbiVector.h +3 -2
  285. data/vendor/eigen/Eigen/src/SparseCore/CompressedStorage.h +16 -0
  286. data/vendor/eigen/Eigen/src/SparseCore/ConservativeSparseSparseProduct.h +6 -6
  287. data/vendor/eigen/Eigen/src/SparseCore/SparseAssign.h +81 -27
  288. data/vendor/eigen/Eigen/src/SparseCore/SparseBlock.h +25 -57
  289. data/vendor/eigen/Eigen/src/SparseCore/SparseCompressedBase.h +40 -11
  290. data/vendor/eigen/Eigen/src/SparseCore/SparseCwiseBinaryOp.h +11 -15
  291. data/vendor/eigen/Eigen/src/SparseCore/SparseCwiseUnaryOp.h +4 -2
  292. data/vendor/eigen/Eigen/src/SparseCore/SparseDenseProduct.h +30 -8
  293. data/vendor/eigen/Eigen/src/SparseCore/SparseMatrix.h +126 -11
  294. data/vendor/eigen/Eigen/src/SparseCore/SparseMatrixBase.h +5 -12
  295. data/vendor/eigen/Eigen/src/SparseCore/SparseProduct.h +13 -1
  296. data/vendor/eigen/Eigen/src/SparseCore/SparseRef.h +7 -7
  297. data/vendor/eigen/Eigen/src/SparseCore/SparseSelfAdjointView.h +5 -2
  298. data/vendor/eigen/Eigen/src/SparseCore/SparseUtil.h +8 -0
  299. data/vendor/eigen/Eigen/src/SparseCore/SparseVector.h +1 -1
  300. data/vendor/eigen/Eigen/src/SparseCore/SparseView.h +1 -0
  301. data/vendor/eigen/Eigen/src/SparseLU/SparseLU.h +162 -12
  302. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_Memory.h +1 -1
  303. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_SupernodalMatrix.h +76 -2
  304. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_column_dfs.h +2 -2
  305. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_gemm_kernel.h +1 -1
  306. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_panel_bmod.h +1 -1
  307. data/vendor/eigen/Eigen/src/SparseQR/SparseQR.h +19 -6
  308. data/vendor/eigen/Eigen/src/StlSupport/StdDeque.h +2 -12
  309. data/vendor/eigen/Eigen/src/StlSupport/StdList.h +2 -2
  310. data/vendor/eigen/Eigen/src/StlSupport/StdVector.h +2 -2
  311. data/vendor/eigen/Eigen/src/SuperLUSupport/SuperLUSupport.h +6 -8
  312. data/vendor/eigen/Eigen/src/UmfPackSupport/UmfPackSupport.h +175 -39
  313. data/vendor/eigen/Eigen/src/misc/lapacke.h +5 -4
  314. data/vendor/eigen/Eigen/src/plugins/ArrayCwiseBinaryOps.h +28 -2
  315. data/vendor/eigen/Eigen/src/plugins/ArrayCwiseUnaryOps.h +155 -11
  316. data/vendor/eigen/Eigen/src/plugins/BlockMethods.h +626 -242
  317. data/vendor/eigen/Eigen/src/plugins/CommonCwiseUnaryOps.h +14 -0
  318. data/vendor/eigen/Eigen/src/plugins/IndexedViewMethods.h +262 -0
  319. data/vendor/eigen/Eigen/src/plugins/MatrixCwiseBinaryOps.h +4 -4
  320. data/vendor/eigen/Eigen/src/plugins/MatrixCwiseUnaryOps.h +10 -0
  321. data/vendor/eigen/Eigen/src/plugins/ReshapedMethods.h +149 -0
  322. data/vendor/eigen/README.md +2 -0
  323. data/vendor/eigen/bench/btl/README +1 -1
  324. data/vendor/eigen/bench/tensors/README +6 -7
  325. data/vendor/eigen/ci/README.md +56 -0
  326. data/vendor/eigen/demos/mix_eigen_and_c/README +1 -1
  327. data/vendor/eigen/unsupported/Eigen/CXX11/src/Tensor/README.md +213 -158
  328. data/vendor/eigen/unsupported/README.txt +1 -1
  329. data/vendor/tomotopy/README.kr.rst +78 -0
  330. data/vendor/tomotopy/README.rst +75 -0
  331. data/vendor/tomotopy/src/Labeling/FoRelevance.cpp +2 -2
  332. data/vendor/tomotopy/src/Labeling/Phraser.hpp +4 -4
  333. data/vendor/tomotopy/src/TopicModel/CTModel.hpp +7 -3
  334. data/vendor/tomotopy/src/TopicModel/DMRModel.hpp +7 -3
  335. data/vendor/tomotopy/src/TopicModel/DTModel.hpp +6 -3
  336. data/vendor/tomotopy/src/TopicModel/GDMRModel.hpp +2 -2
  337. data/vendor/tomotopy/src/TopicModel/HDP.h +1 -0
  338. data/vendor/tomotopy/src/TopicModel/HDPModel.hpp +57 -6
  339. data/vendor/tomotopy/src/TopicModel/HLDAModel.hpp +6 -3
  340. data/vendor/tomotopy/src/TopicModel/HPAModel.hpp +3 -2
  341. data/vendor/tomotopy/src/TopicModel/LDA.h +3 -3
  342. data/vendor/tomotopy/src/TopicModel/LDACVB0Model.hpp +5 -5
  343. data/vendor/tomotopy/src/TopicModel/LDAModel.hpp +50 -19
  344. data/vendor/tomotopy/src/TopicModel/LLDAModel.hpp +6 -2
  345. data/vendor/tomotopy/src/TopicModel/MGLDAModel.hpp +3 -2
  346. data/vendor/tomotopy/src/TopicModel/PAModel.hpp +1 -1
  347. data/vendor/tomotopy/src/TopicModel/PLDAModel.hpp +6 -2
  348. data/vendor/tomotopy/src/TopicModel/PT.h +3 -1
  349. data/vendor/tomotopy/src/TopicModel/PTModel.hpp +36 -3
  350. data/vendor/tomotopy/src/TopicModel/SLDAModel.hpp +6 -3
  351. data/vendor/tomotopy/src/TopicModel/TopicModel.hpp +55 -26
  352. data/vendor/tomotopy/src/Utils/AliasMethod.hpp +5 -4
  353. data/vendor/tomotopy/src/Utils/Dictionary.h +2 -2
  354. data/vendor/tomotopy/src/Utils/EigenAddonOps.hpp +36 -1
  355. data/vendor/tomotopy/src/Utils/MultiNormalDistribution.hpp +1 -1
  356. data/vendor/tomotopy/src/Utils/TruncMultiNormal.hpp +1 -1
  357. data/vendor/tomotopy/src/Utils/exception.h +6 -0
  358. data/vendor/tomotopy/src/Utils/math.h +2 -2
  359. data/vendor/tomotopy/src/Utils/sample.hpp +14 -12
  360. data/vendor/tomotopy/src/Utils/serializer.hpp +30 -5
  361. data/vendor/tomotopy/src/Utils/sse_gamma.h +0 -3
  362. metadata +64 -18
  363. data/vendor/eigen/Eigen/CMakeLists.txt +0 -19
  364. data/vendor/eigen/Eigen/src/Core/arch/CUDA/Half.h +0 -674
  365. data/vendor/eigen/Eigen/src/Core/arch/CUDA/PacketMath.h +0 -333
  366. data/vendor/eigen/Eigen/src/Core/arch/CUDA/PacketMathHalf.h +0 -1124
  367. data/vendor/eigen/Eigen/src/Core/arch/CUDA/TypeCasting.h +0 -212
  368. data/vendor/eigen/Eigen/src/Geometry/arch/Geometry_SSE.h +0 -161
  369. data/vendor/eigen/Eigen/src/LU/arch/Inverse_SSE.h +0 -338
data/ext/tomoto/lda.cpp CHANGED
@@ -3,9 +3,7 @@
3
3
 
4
4
  #include <LDA.h>
5
5
 
6
- #include <rice/Class.hpp>
7
- #include <rice/Hash.hpp>
8
- #include <rice/Module.hpp>
6
+ #include <rice/rice.hpp>
9
7
 
10
8
  #include "utils.h"
11
9
 
@@ -22,7 +20,7 @@ void init_lda(Rice::Module& m) {
22
20
  Rice::define_class_under<DocumentObject>(m, "Document")
23
21
  .define_method(
24
22
  "topics",
25
- *[](DocumentObject& self) {
23
+ [](DocumentObject& self) {
26
24
  Rice::Hash res;
27
25
  auto topics = self.tm->getTopicsByDoc(self.doc);
28
26
  for (size_t i = 0; i < topics.size(); i++) {
@@ -32,9 +30,9 @@ void init_lda(Rice::Module& m) {
32
30
  });
33
31
 
34
32
  Rice::define_class_under<tomoto::ILDAModel>(m, "LDA")
35
- .define_singleton_method(
33
+ .define_singleton_function(
36
34
  "_new",
37
- *[](size_t tw, size_t k, tomoto::Float alpha, tomoto::Float eta, size_t seed) {
35
+ [](size_t tw, size_t k, tomoto::Float alpha, tomoto::Float eta, size_t seed) {
38
36
  tomoto::LDAArgs args;
39
37
  args.k = k;
40
38
  args.alpha = {alpha};
@@ -43,10 +41,10 @@ void init_lda(Rice::Module& m) {
43
41
  args.seed = seed;
44
42
  }
45
43
  return tomoto::ILDAModel::create((tomoto::TermWeight)tw, args);
46
- })
44
+ }, Rice::Return().takeOwnership())
47
45
  .define_method(
48
46
  "_add_doc",
49
- *[](tomoto::ILDAModel& self, std::vector<std::string> words) {
47
+ [](tomoto::ILDAModel& self, std::vector<std::string> words) {
50
48
  return self.addDoc(buildDoc(words));
51
49
  })
52
50
  .define_method(
@@ -75,7 +73,7 @@ void init_lda(Rice::Module& m) {
75
73
  })
76
74
  .define_method(
77
75
  "alpha",
78
- *[](tomoto::ILDAModel& self) {
76
+ [](tomoto::ILDAModel& self) {
79
77
  Array res;
80
78
  for (size_t i = 0; i < self.getK(); i++) {
81
79
  res.push(self.getAlpha(i));
@@ -84,18 +82,18 @@ void init_lda(Rice::Module& m) {
84
82
  })
85
83
  .define_method(
86
84
  "burn_in",
87
- *[](tomoto::ILDAModel& self) {
85
+ [](tomoto::ILDAModel& self) {
88
86
  return self.getBurnInIteration();
89
87
  })
90
88
  .define_method(
91
89
  "burn_in=",
92
- *[](tomoto::ILDAModel& self, size_t iteration) {
90
+ [](tomoto::ILDAModel& self, size_t iteration) {
93
91
  self.setBurnInIteration(iteration);
94
92
  return iteration;
95
93
  })
96
94
  .define_method(
97
95
  "_count_by_topics",
98
- *[](tomoto::ILDAModel& self) {
96
+ [](tomoto::ILDAModel& self) {
99
97
  Array res;
100
98
  for (auto const& v : self.getCountByTopic()) {
101
99
  res.push(v);
@@ -104,32 +102,33 @@ void init_lda(Rice::Module& m) {
104
102
  })
105
103
  .define_method(
106
104
  "docs",
107
- *[](tomoto::ILDAModel& self) {
105
+ [](tomoto::ILDAModel& self) {
108
106
  Array res;
109
107
  auto n = self.getNumDocs();
110
108
  for (size_t i = 0; i < n; i++) {
111
- res.push(DocumentObject(self.getDoc(i), &self));
109
+ auto v = DocumentObject(self.getDoc(i), &self);
110
+ res.push(Object(Rice::detail::To_Ruby<DocumentObject>().convert(v)));
112
111
  }
113
112
  return res;
114
113
  })
115
114
  .define_method(
116
115
  "eta",
117
- *[](tomoto::ILDAModel& self) {
116
+ [](tomoto::ILDAModel& self) {
118
117
  return self.getEta();
119
118
  })
120
119
  .define_method(
121
120
  "global_step",
122
- *[](tomoto::ILDAModel& self) {
121
+ [](tomoto::ILDAModel& self) {
123
122
  return self.getGlobalStep();
124
123
  })
125
124
  .define_method(
126
125
  "k",
127
- *[](tomoto::ILDAModel& self) {
126
+ [](tomoto::ILDAModel& self) {
128
127
  return self.getK();
129
128
  })
130
129
  .define_method(
131
130
  "_load",
132
- *[](tomoto::ILDAModel& self, const char* filename) {
131
+ [](tomoto::ILDAModel& self, const char* filename) {
133
132
  std::ifstream str{ filename, std::ios_base::binary };
134
133
  if (!str) throw std::runtime_error{ std::string("cannot open file '") + filename + std::string("'") };
135
134
  std::vector<uint8_t> extra_data;
@@ -137,48 +136,48 @@ void init_lda(Rice::Module& m) {
137
136
  })
138
137
  .define_method(
139
138
  "ll_per_word",
140
- *[](tomoto::ILDAModel& self) {
139
+ [](tomoto::ILDAModel& self) {
141
140
  return self.getLLPerWord();
142
141
  })
143
142
  .define_method(
144
143
  "num_docs",
145
- *[](tomoto::ILDAModel& self) {
144
+ [](tomoto::ILDAModel& self) {
146
145
  return self.getNumDocs();
147
146
  })
148
147
  .define_method(
149
148
  "num_vocabs",
150
- *[](tomoto::ILDAModel& self) {
149
+ [](tomoto::ILDAModel& self) {
151
150
  return self.getV();
152
151
  })
153
152
  .define_method(
154
153
  "num_words",
155
- *[](tomoto::ILDAModel& self) {
154
+ [](tomoto::ILDAModel& self) {
156
155
  return self.getN();
157
156
  })
158
157
  .define_method(
159
158
  "optim_interval",
160
- *[](tomoto::ILDAModel& self) {
159
+ [](tomoto::ILDAModel& self) {
161
160
  return self.getOptimInterval();
162
161
  })
163
162
  .define_method(
164
163
  "optim_interval=",
165
- *[](tomoto::ILDAModel& self, size_t value) {
164
+ [](tomoto::ILDAModel& self, size_t value) {
166
165
  self.setOptimInterval(value);
167
166
  return value;
168
167
  })
169
168
  .define_method(
170
169
  "perplexity",
171
- *[](tomoto::ILDAModel& self) {
170
+ [](tomoto::ILDAModel& self) {
172
171
  return self.getPerplexity();
173
172
  })
174
173
  .define_method(
175
174
  "_prepare",
176
- *[](tomoto::ILDAModel& self, size_t minCnt, size_t minDf, size_t rmTop) {
175
+ [](tomoto::ILDAModel& self, size_t minCnt, size_t minDf, size_t rmTop) {
177
176
  self.prepare(true, minCnt, minDf, rmTop);
178
177
  })
179
178
  .define_method(
180
179
  "_removed_top_words",
181
- *[](tomoto::ILDAModel& self, size_t rmTop) {
180
+ [](tomoto::ILDAModel& self, size_t rmTop) {
182
181
  Array res;
183
182
  auto dict = self.getVocabDict();
184
183
  size_t size = dict.size();
@@ -189,14 +188,14 @@ void init_lda(Rice::Module& m) {
189
188
  })
190
189
  .define_method(
191
190
  "_save",
192
- *[](tomoto::ILDAModel& self, const char* filename, bool full) {
191
+ [](tomoto::ILDAModel& self, const char* filename, bool full) {
193
192
  std::ofstream str{ filename, std::ios_base::binary };
194
193
  std::vector<uint8_t> extra_data;
195
194
  self.saveModel(str, full, &extra_data);
196
195
  })
197
196
  .define_method(
198
197
  "_topic_words",
199
- *[](tomoto::ILDAModel& self, size_t topicId, size_t topN) {
198
+ [](tomoto::ILDAModel& self, size_t topicId, size_t topN) {
200
199
  Rice::Hash res;
201
200
  for (auto const& v : self.getWordsByTopicSorted(topicId, topN)) {
202
201
  res[v.first] = v.second;
@@ -205,17 +204,17 @@ void init_lda(Rice::Module& m) {
205
204
  })
206
205
  .define_method(
207
206
  "_train",
208
- *[](tomoto::ILDAModel& self, size_t iteration, size_t workers, size_t ps) {
207
+ [](tomoto::ILDAModel& self, size_t iteration, size_t workers, size_t ps) {
209
208
  self.train(iteration, workers, (tomoto::ParallelScheme)ps);
210
209
  })
211
210
  .define_method(
212
211
  "_tw",
213
- *[](tomoto::ILDAModel& self) {
212
+ [](tomoto::ILDAModel& self) {
214
213
  return (int)self.getTermWeight();
215
214
  })
216
215
  .define_method(
217
216
  "used_vocab_df",
218
- *[](tomoto::ILDAModel& self) {
217
+ [](tomoto::ILDAModel& self) {
219
218
  auto vocab = self.getVocabDf();
220
219
  Array res;
221
220
  for (size_t i = 0; i < self.getV(); i++) {
@@ -225,7 +224,7 @@ void init_lda(Rice::Module& m) {
225
224
  })
226
225
  .define_method(
227
226
  "used_vocab_freq",
228
- *[](tomoto::ILDAModel& self) {
227
+ [](tomoto::ILDAModel& self) {
229
228
  auto vocab = self.getVocabCf();
230
229
  Array res;
231
230
  for (size_t i = 0; i < self.getV(); i++) {
@@ -235,18 +234,20 @@ void init_lda(Rice::Module& m) {
235
234
  })
236
235
  .define_method(
237
236
  "used_vocabs",
238
- *[](tomoto::ILDAModel& self) {
237
+ [](tomoto::ILDAModel& self) {
239
238
  auto dict = self.getVocabDict();
240
239
  Array res;
241
240
  auto utf8 = Rice::Class(rb_cEncoding).call("const_get", "UTF_8");
242
241
  for (size_t i = 0; i < self.getV(); i++) {
243
- res.push(to_ruby<std::string>(dict.toWord(i)).call("force_encoding", utf8));
242
+ VALUE value = Rice::detail::To_Ruby<std::string>().convert(dict.toWord(i));
243
+ Object obj(value);
244
+ res.push(obj.call("force_encoding", utf8));
244
245
  }
245
246
  return res;
246
247
  })
247
248
  .define_method(
248
249
  "vocab_df",
249
- *[](tomoto::ILDAModel& self) {
250
+ [](tomoto::ILDAModel& self) {
250
251
  auto vocab = self.getVocabDf();
251
252
  Array res;
252
253
  for (size_t i = 0; i < vocab.size(); i++) {
@@ -256,7 +257,7 @@ void init_lda(Rice::Module& m) {
256
257
  })
257
258
  .define_method(
258
259
  "vocab_freq",
259
- *[](tomoto::ILDAModel& self) {
260
+ [](tomoto::ILDAModel& self) {
260
261
  auto vocab = self.getVocabCf();
261
262
  Array res;
262
263
  for (size_t i = 0; i < vocab.size(); i++) {
@@ -266,12 +267,14 @@ void init_lda(Rice::Module& m) {
266
267
  })
267
268
  .define_method(
268
269
  "vocabs",
269
- *[](tomoto::ILDAModel& self) {
270
+ [](tomoto::ILDAModel& self) {
270
271
  auto dict = self.getVocabDict();
271
272
  Array res;
272
273
  auto utf8 = Rice::Class(rb_cEncoding).call("const_get", "UTF_8");
273
274
  for (size_t i = 0; i < dict.size(); i++) {
274
- res.push(to_ruby<std::string>(dict.toWord(i)).call("force_encoding", utf8));
275
+ VALUE value = Rice::detail::To_Ruby<std::string>().convert(dict.toWord(i));
276
+ Object obj(value);
277
+ res.push(obj.call("force_encoding", utf8));
275
278
  }
276
279
  return res;
277
280
  });
data/ext/tomoto/llda.cpp CHANGED
@@ -1,14 +1,14 @@
1
1
  #include <LLDA.h>
2
2
 
3
- #include <rice/Module.hpp>
3
+ #include <rice/rice.hpp>
4
4
 
5
5
  #include "utils.h"
6
6
 
7
7
  void init_llda(Rice::Module& m) {
8
8
  Rice::define_class_under<tomoto::ILLDAModel, tomoto::ILDAModel>(m, "LLDA")
9
- .define_singleton_method(
9
+ .define_singleton_function(
10
10
  "_new",
11
- *[](size_t tw, size_t k, tomoto::Float alpha, tomoto::Float eta, size_t seed) {
11
+ [](size_t tw, size_t k, tomoto::Float alpha, tomoto::Float eta, size_t seed) {
12
12
  tomoto::LDAArgs args;
13
13
  args.k = k;
14
14
  args.alpha = {alpha};
@@ -17,17 +17,17 @@ void init_llda(Rice::Module& m) {
17
17
  args.seed = seed;
18
18
  }
19
19
  return tomoto::ILLDAModel::create((tomoto::TermWeight)tw, args);
20
- })
20
+ }, Rice::Return().takeOwnership())
21
21
  .define_method(
22
22
  "_add_doc",
23
- *[](tomoto::ILLDAModel& self, std::vector<std::string> words, std::vector<std::string> labels) {
23
+ [](tomoto::ILLDAModel& self, std::vector<std::string> words, std::vector<std::string> labels) {
24
24
  auto doc = buildDoc(words);
25
25
  doc.misc["labels"] = labels;
26
26
  return self.addDoc(doc);
27
27
  })
28
28
  .define_method(
29
29
  "topics_per_label",
30
- *[](tomoto::ILLDAModel& self) {
30
+ [](tomoto::ILLDAModel& self) {
31
31
  return self.getNumTopicsPerLabel();
32
32
  });
33
33
  }
data/ext/tomoto/mglda.cpp CHANGED
@@ -1,14 +1,14 @@
1
1
  #include <MGLDA.h>
2
2
 
3
- #include <rice/Module.hpp>
3
+ #include <rice/rice.hpp>
4
4
 
5
5
  #include "utils.h"
6
6
 
7
7
  void init_mglda(Rice::Module& m) {
8
8
  Rice::define_class_under<tomoto::IMGLDAModel, tomoto::ILDAModel>(m, "MGLDA")
9
- .define_singleton_method(
9
+ .define_singleton_function(
10
10
  "_new",
11
- *[](size_t tw, size_t k_g, size_t k_l, size_t t, tomoto::Float alpha_g, tomoto::Float alpha_l, tomoto::Float alpha_mg, tomoto::Float alpha_ml, tomoto::Float eta_g) {
11
+ [](size_t tw, size_t k_g, size_t k_l, size_t t, tomoto::Float alpha_g, tomoto::Float alpha_l, tomoto::Float alpha_mg, tomoto::Float alpha_ml, tomoto::Float eta_g) {
12
12
  tomoto::MGLDAArgs args;
13
13
  args.k = k_g;
14
14
  args.kL = k_l;
@@ -20,62 +20,62 @@ void init_mglda(Rice::Module& m) {
20
20
  args.eta = eta_g;
21
21
  // TODO more args
22
22
  return tomoto::IMGLDAModel::create((tomoto::TermWeight)tw, args);
23
- })
23
+ }, Rice::Return().takeOwnership())
24
24
  .define_method(
25
25
  "_add_doc",
26
- *[](tomoto::IMGLDAModel& self, std::vector<std::string> words, std::string delimiter) {
26
+ [](tomoto::IMGLDAModel& self, std::vector<std::string> words, std::string delimiter) {
27
27
  auto doc = buildDoc(words);
28
28
  doc.misc["delimiter"] = delimiter;
29
29
  return self.addDoc(doc);
30
30
  })
31
31
  .define_method(
32
32
  "alpha_g",
33
- *[](tomoto::IMGLDAModel& self) {
33
+ [](tomoto::IMGLDAModel& self) {
34
34
  return self.getAlpha();
35
35
  })
36
36
  .define_method(
37
37
  "alpha_l",
38
- *[](tomoto::IMGLDAModel& self) {
38
+ [](tomoto::IMGLDAModel& self) {
39
39
  return self.getAlphaL();
40
40
  })
41
41
  .define_method(
42
42
  "alpha_mg",
43
- *[](tomoto::IMGLDAModel& self) {
43
+ [](tomoto::IMGLDAModel& self) {
44
44
  return self.getAlphaM();
45
45
  })
46
46
  .define_method(
47
47
  "alpha_ml",
48
- *[](tomoto::IMGLDAModel& self) {
48
+ [](tomoto::IMGLDAModel& self) {
49
49
  return self.getAlphaML();
50
50
  })
51
51
  .define_method(
52
52
  "eta_g",
53
- *[](tomoto::IMGLDAModel& self) {
53
+ [](tomoto::IMGLDAModel& self) {
54
54
  return self.getEta();
55
55
  })
56
56
  .define_method(
57
57
  "eta_l",
58
- *[](tomoto::IMGLDAModel& self) {
58
+ [](tomoto::IMGLDAModel& self) {
59
59
  return self.getEtaL();
60
60
  })
61
61
  .define_method(
62
62
  "gamma",
63
- *[](tomoto::IMGLDAModel& self) {
63
+ [](tomoto::IMGLDAModel& self) {
64
64
  return self.getGamma();
65
65
  })
66
66
  .define_method(
67
67
  "k_g",
68
- *[](tomoto::IMGLDAModel& self) {
68
+ [](tomoto::IMGLDAModel& self) {
69
69
  return self.getK();
70
70
  })
71
71
  .define_method(
72
72
  "k_l",
73
- *[](tomoto::IMGLDAModel& self) {
73
+ [](tomoto::IMGLDAModel& self) {
74
74
  return self.getKL();
75
75
  })
76
76
  .define_method(
77
77
  "t",
78
- *[](tomoto::IMGLDAModel& self) {
78
+ [](tomoto::IMGLDAModel& self) {
79
79
  return self.getT();
80
80
  });
81
81
  }
data/ext/tomoto/pa.cpp CHANGED
@@ -1,14 +1,14 @@
1
1
  #include <PA.h>
2
2
 
3
- #include <rice/Module.hpp>
3
+ #include <rice/rice.hpp>
4
4
 
5
5
  #include "utils.h"
6
6
 
7
7
  void init_pa(Rice::Module& m) {
8
8
  Rice::define_class_under<tomoto::IPAModel, tomoto::ILDAModel>(m, "PA")
9
- .define_singleton_method(
9
+ .define_singleton_function(
10
10
  "_new",
11
- *[](size_t tw, size_t k1, size_t k2, tomoto::Float alpha, tomoto::Float eta, size_t seed) {
11
+ [](size_t tw, size_t k1, size_t k2, tomoto::Float alpha, tomoto::Float eta, size_t seed) {
12
12
  tomoto::PAArgs args;
13
13
  args.k = k1;
14
14
  args.k2 = k2;
@@ -18,15 +18,15 @@ void init_pa(Rice::Module& m) {
18
18
  args.seed = seed;
19
19
  }
20
20
  return tomoto::IPAModel::create((tomoto::TermWeight)tw, args);
21
- })
21
+ }, Rice::Return().takeOwnership())
22
22
  .define_method(
23
23
  "k1",
24
- *[](tomoto::IPAModel& self) {
24
+ [](tomoto::IPAModel& self) {
25
25
  return self.getK();
26
26
  })
27
27
  .define_method(
28
28
  "k2",
29
- *[](tomoto::IPAModel& self) {
29
+ [](tomoto::IPAModel& self) {
30
30
  return self.getK2();
31
31
  });
32
32
  }
data/ext/tomoto/plda.cpp CHANGED
@@ -1,14 +1,14 @@
1
1
  #include <PLDA.h>
2
2
 
3
- #include <rice/Module.hpp>
3
+ #include <rice/rice.hpp>
4
4
 
5
5
  #include "utils.h"
6
6
 
7
7
  void init_plda(Rice::Module& m) {
8
8
  Rice::define_class_under<tomoto::IPLDAModel, tomoto::ILLDAModel>(m, "PLDA")
9
- .define_singleton_method(
9
+ .define_singleton_function(
10
10
  "_new",
11
- *[](size_t tw, size_t latent_topics, tomoto::Float alpha, tomoto::Float eta, size_t seed) {
11
+ [](size_t tw, size_t latent_topics, tomoto::Float alpha, tomoto::Float eta, size_t seed) {
12
12
  tomoto::PLDAArgs args;
13
13
  args.numLatentTopics = latent_topics;
14
14
  args.alpha = {alpha};
@@ -17,17 +17,17 @@ void init_plda(Rice::Module& m) {
17
17
  args.seed = seed;
18
18
  }
19
19
  return tomoto::IPLDAModel::create((tomoto::TermWeight)tw, args);
20
- })
20
+ }, Rice::Return().takeOwnership())
21
21
  .define_method(
22
22
  "_add_doc",
23
- *[](tomoto::IPLDAModel& self, std::vector<std::string> words, std::vector<std::string> labels) {
23
+ [](tomoto::IPLDAModel& self, std::vector<std::string> words, std::vector<std::string> labels) {
24
24
  auto doc = buildDoc(words);
25
25
  doc.misc["labels"] = labels;
26
26
  return self.addDoc(doc);
27
27
  })
28
28
  .define_method(
29
29
  "latent_topics",
30
- *[](tomoto::IPLDAModel& self) {
30
+ [](tomoto::IPLDAModel& self) {
31
31
  return self.getNumLatentTopics();
32
32
  });
33
33
  }
data/ext/tomoto/slda.cpp CHANGED
@@ -1,18 +1,18 @@
1
1
  #include <SLDA.h>
2
2
 
3
- #include <rice/Module.hpp>
3
+ #include <rice/rice.hpp>
4
4
 
5
5
  #include "utils.h"
6
6
 
7
7
  void init_slda(Rice::Module& m) {
8
8
  Rice::define_class_under<tomoto::ISLDAModel, tomoto::ILDAModel>(m, "SLDA")
9
- .define_singleton_method(
9
+ .define_singleton_function(
10
10
  "_new",
11
- *[](size_t tw, size_t k, Array rb_vars, tomoto::Float alpha, tomoto::Float eta, std::vector<tomoto::Float> mu, std::vector<tomoto::Float> nu_sq, std::vector<tomoto::Float> glm_param, size_t seed) {
11
+ [](size_t tw, size_t k, Array rb_vars, tomoto::Float alpha, tomoto::Float eta, std::vector<tomoto::Float> mu, std::vector<tomoto::Float> nu_sq, std::vector<tomoto::Float> glm_param, size_t seed) {
12
12
  std::vector<tomoto::ISLDAModel::GLM> vars;
13
13
  vars.reserve(rb_vars.size());
14
14
  for (auto const& v : rb_vars) {
15
- vars.push_back((tomoto::ISLDAModel::GLM) from_ruby<int>(v));
15
+ vars.push_back((tomoto::ISLDAModel::GLM) Rice::detail::From_Ruby<int>().convert(v.value()));
16
16
  }
17
17
  tomoto::SLDAArgs args;
18
18
  args.k = k;
@@ -26,22 +26,22 @@ void init_slda(Rice::Module& m) {
26
26
  args.seed = seed;
27
27
  }
28
28
  return tomoto::ISLDAModel::create((tomoto::TermWeight)tw, args);
29
- })
29
+ }, Rice::Return().takeOwnership())
30
30
  .define_method(
31
31
  "_add_doc",
32
- *[](tomoto::ISLDAModel& self, std::vector<std::string> words, std::vector<tomoto::Float> y) {
32
+ [](tomoto::ISLDAModel& self, std::vector<std::string> words, std::vector<tomoto::Float> y) {
33
33
  auto doc = buildDoc(words);
34
34
  doc.misc["y"] = y;
35
35
  return self.addDoc(doc);
36
36
  })
37
37
  .define_method(
38
38
  "f",
39
- *[](tomoto::ISLDAModel& self) {
39
+ [](tomoto::ISLDAModel& self) {
40
40
  return self.getF();
41
41
  })
42
42
  .define_method(
43
43
  "_var_type",
44
- *[](tomoto::ISLDAModel& self, size_t var_id) {
44
+ [](tomoto::ISLDAModel& self, size_t var_id) {
45
45
  if (var_id >= self.getF()) throw std::runtime_error{ "'var_id' must be < 'f'" };
46
46
  return self.getTypeOfVar(var_id) == tomoto::ISLDAModel::GLM::linear ? "l" : "b";
47
47
  });
@@ -1,4 +1,4 @@
1
- #include <rice/Module.hpp>
1
+ #include <rice/rice.hpp>
2
2
 
3
3
  void init_lda(Rice::Module& m);
4
4
  void init_ct(Rice::Module& m);
@@ -15,20 +15,20 @@ void init_plda(Rice::Module& m);
15
15
  void init_slda(Rice::Module& m);
16
16
 
17
17
  extern "C"
18
- void Init_ext()
18
+ void Init_tomoto()
19
19
  {
20
20
  auto m = Rice::define_module("Tomoto")
21
- .define_singleton_method(
21
+ .define_singleton_function(
22
22
  "isa",
23
- *[]() {
23
+ []() {
24
24
  #ifdef __AVX2__
25
- return "avx2";
25
+ return Rice::String("avx2");
26
26
  #elif defined(__AVX__)
27
- return "avx";
27
+ return Rice::String("avx");
28
28
  #elif defined(__SSE2__) || defined(__x86_64__) || defined(_WIN64)
29
- return "sse2";
29
+ return Rice::String("sse2");
30
30
  #else
31
- return "none";
31
+ return Rice::String("none");
32
32
  #endif
33
33
  });
34
34
 
data/ext/tomoto/utils.h CHANGED
@@ -1,80 +1,26 @@
1
1
  #pragma once
2
2
 
3
- #include <rice/Array.hpp>
3
+ #include <rice/rice.hpp>
4
+ #include <rice/stl.hpp>
4
5
 
5
6
  using Rice::Array;
6
7
  using Rice::Object;
7
8
 
8
- template<>
9
- inline
10
- Object to_ruby<std::vector<tomoto::Float>>(std::vector<tomoto::Float> const & x)
9
+ namespace Rice::detail
11
10
  {
12
- Array res;
13
- for (auto const& v : x) {
14
- res.push(v);
15
- }
16
- return res;
17
- }
18
-
19
- template<>
20
- inline
21
- Object to_ruby<std::vector<uint32_t>>(std::vector<uint32_t> const & x)
22
- {
23
- Array res;
24
- for (auto const& v : x) {
25
- res.push(v);
26
- }
27
- return res;
28
- }
29
-
30
- template<>
31
- inline
32
- Object to_ruby<std::vector<uint64_t>>(std::vector<uint64_t> const & x)
33
- {
34
- Array res;
35
- for (auto const& v : x) {
36
- res.push(v);
37
- }
38
- return res;
39
- }
40
-
41
- template<>
42
- inline
43
- std::vector<std::string> from_ruby<std::vector<std::string>>(Object x)
44
- {
45
- Array a = Array(x);
46
- std::vector<std::string> res;
47
- res.reserve(a.size());
48
- for (auto const& v : a) {
49
- res.push_back(from_ruby<std::string>(v));
50
- }
51
- return res;
52
- }
53
-
54
- template<>
55
- inline
56
- std::vector<tomoto::Float> from_ruby<std::vector<tomoto::Float>>(Object x)
57
- {
58
- Array a = Array(x);
59
- std::vector<tomoto::Float> res;
60
- res.reserve(a.size());
61
- for (auto const& v : a) {
62
- res.push_back(from_ruby<tomoto::Float>(v));
63
- }
64
- return res;
65
- }
66
-
67
- template<>
68
- inline
69
- std::vector<uint64_t> from_ruby<std::vector<uint64_t>>(Object x)
70
- {
71
- Array a = Array(x);
72
- std::vector<uint64_t> res;
73
- res.reserve(a.size());
74
- for (auto const& v : a) {
75
- res.push_back(from_ruby<uint64_t>(v));
76
- }
77
- return res;
11
+ template<typename T>
12
+ class To_Ruby<std::vector<T>>
13
+ {
14
+ public:
15
+ VALUE convert(std::vector<T> const & x)
16
+ {
17
+ auto a = rb_ary_new2(x.size());
18
+ for (const auto& v : x) {
19
+ detail::protect(rb_ary_push, a, To_Ruby<T>().convert(v));
20
+ }
21
+ return a;
22
+ }
23
+ };
78
24
  }
79
25
 
80
26
  inline tomoto::RawDoc buildDoc(std::vector<std::string>& words) {
@@ -1,3 +1,3 @@
1
1
  module Tomoto
2
- VERSION = "0.2.2"
2
+ VERSION = "0.3.0"
3
3
  end
data/lib/tomoto.rb CHANGED
@@ -1,5 +1,9 @@
1
1
  # ext
2
- require "tomoto/ext"
2
+ begin
3
+ require "tomoto/#{RUBY_VERSION.to_f}/tomoto"
4
+ rescue LoadError
5
+ require "tomoto/tomoto"
6
+ end
3
7
 
4
8
  # modules
5
9
  require "tomoto/ct"