tomoto 0.1.0

Sign up to get free protection for your applications and to get access to all the features.
Files changed (420) hide show
  1. checksums.yaml +7 -0
  2. data/CHANGELOG.md +3 -0
  3. data/LICENSE.txt +22 -0
  4. data/README.md +123 -0
  5. data/ext/tomoto/ext.cpp +245 -0
  6. data/ext/tomoto/extconf.rb +28 -0
  7. data/lib/tomoto.rb +12 -0
  8. data/lib/tomoto/ct.rb +11 -0
  9. data/lib/tomoto/hdp.rb +11 -0
  10. data/lib/tomoto/lda.rb +67 -0
  11. data/lib/tomoto/version.rb +3 -0
  12. data/vendor/EigenRand/EigenRand/Core.h +1139 -0
  13. data/vendor/EigenRand/EigenRand/Dists/Basic.h +111 -0
  14. data/vendor/EigenRand/EigenRand/Dists/Discrete.h +877 -0
  15. data/vendor/EigenRand/EigenRand/Dists/GammaPoisson.h +108 -0
  16. data/vendor/EigenRand/EigenRand/Dists/NormalExp.h +626 -0
  17. data/vendor/EigenRand/EigenRand/EigenRand +19 -0
  18. data/vendor/EigenRand/EigenRand/Macro.h +24 -0
  19. data/vendor/EigenRand/EigenRand/MorePacketMath.h +978 -0
  20. data/vendor/EigenRand/EigenRand/PacketFilter.h +286 -0
  21. data/vendor/EigenRand/EigenRand/PacketRandomEngine.h +624 -0
  22. data/vendor/EigenRand/EigenRand/RandUtils.h +413 -0
  23. data/vendor/EigenRand/EigenRand/doc.h +220 -0
  24. data/vendor/EigenRand/LICENSE +21 -0
  25. data/vendor/EigenRand/README.md +288 -0
  26. data/vendor/eigen/COPYING.BSD +26 -0
  27. data/vendor/eigen/COPYING.GPL +674 -0
  28. data/vendor/eigen/COPYING.LGPL +502 -0
  29. data/vendor/eigen/COPYING.MINPACK +52 -0
  30. data/vendor/eigen/COPYING.MPL2 +373 -0
  31. data/vendor/eigen/COPYING.README +18 -0
  32. data/vendor/eigen/Eigen/CMakeLists.txt +19 -0
  33. data/vendor/eigen/Eigen/Cholesky +46 -0
  34. data/vendor/eigen/Eigen/CholmodSupport +48 -0
  35. data/vendor/eigen/Eigen/Core +537 -0
  36. data/vendor/eigen/Eigen/Dense +7 -0
  37. data/vendor/eigen/Eigen/Eigen +2 -0
  38. data/vendor/eigen/Eigen/Eigenvalues +61 -0
  39. data/vendor/eigen/Eigen/Geometry +62 -0
  40. data/vendor/eigen/Eigen/Householder +30 -0
  41. data/vendor/eigen/Eigen/IterativeLinearSolvers +48 -0
  42. data/vendor/eigen/Eigen/Jacobi +33 -0
  43. data/vendor/eigen/Eigen/LU +50 -0
  44. data/vendor/eigen/Eigen/MetisSupport +35 -0
  45. data/vendor/eigen/Eigen/OrderingMethods +73 -0
  46. data/vendor/eigen/Eigen/PaStiXSupport +48 -0
  47. data/vendor/eigen/Eigen/PardisoSupport +35 -0
  48. data/vendor/eigen/Eigen/QR +51 -0
  49. data/vendor/eigen/Eigen/QtAlignedMalloc +40 -0
  50. data/vendor/eigen/Eigen/SPQRSupport +34 -0
  51. data/vendor/eigen/Eigen/SVD +51 -0
  52. data/vendor/eigen/Eigen/Sparse +36 -0
  53. data/vendor/eigen/Eigen/SparseCholesky +45 -0
  54. data/vendor/eigen/Eigen/SparseCore +69 -0
  55. data/vendor/eigen/Eigen/SparseLU +46 -0
  56. data/vendor/eigen/Eigen/SparseQR +37 -0
  57. data/vendor/eigen/Eigen/StdDeque +27 -0
  58. data/vendor/eigen/Eigen/StdList +26 -0
  59. data/vendor/eigen/Eigen/StdVector +27 -0
  60. data/vendor/eigen/Eigen/SuperLUSupport +64 -0
  61. data/vendor/eigen/Eigen/UmfPackSupport +40 -0
  62. data/vendor/eigen/Eigen/src/Cholesky/LDLT.h +673 -0
  63. data/vendor/eigen/Eigen/src/Cholesky/LLT.h +542 -0
  64. data/vendor/eigen/Eigen/src/Cholesky/LLT_LAPACKE.h +99 -0
  65. data/vendor/eigen/Eigen/src/CholmodSupport/CholmodSupport.h +639 -0
  66. data/vendor/eigen/Eigen/src/Core/Array.h +329 -0
  67. data/vendor/eigen/Eigen/src/Core/ArrayBase.h +226 -0
  68. data/vendor/eigen/Eigen/src/Core/ArrayWrapper.h +209 -0
  69. data/vendor/eigen/Eigen/src/Core/Assign.h +90 -0
  70. data/vendor/eigen/Eigen/src/Core/AssignEvaluator.h +935 -0
  71. data/vendor/eigen/Eigen/src/Core/Assign_MKL.h +178 -0
  72. data/vendor/eigen/Eigen/src/Core/BandMatrix.h +353 -0
  73. data/vendor/eigen/Eigen/src/Core/Block.h +452 -0
  74. data/vendor/eigen/Eigen/src/Core/BooleanRedux.h +164 -0
  75. data/vendor/eigen/Eigen/src/Core/CommaInitializer.h +160 -0
  76. data/vendor/eigen/Eigen/src/Core/ConditionEstimator.h +175 -0
  77. data/vendor/eigen/Eigen/src/Core/CoreEvaluators.h +1688 -0
  78. data/vendor/eigen/Eigen/src/Core/CoreIterators.h +127 -0
  79. data/vendor/eigen/Eigen/src/Core/CwiseBinaryOp.h +184 -0
  80. data/vendor/eigen/Eigen/src/Core/CwiseNullaryOp.h +866 -0
  81. data/vendor/eigen/Eigen/src/Core/CwiseTernaryOp.h +197 -0
  82. data/vendor/eigen/Eigen/src/Core/CwiseUnaryOp.h +103 -0
  83. data/vendor/eigen/Eigen/src/Core/CwiseUnaryView.h +128 -0
  84. data/vendor/eigen/Eigen/src/Core/DenseBase.h +611 -0
  85. data/vendor/eigen/Eigen/src/Core/DenseCoeffsBase.h +681 -0
  86. data/vendor/eigen/Eigen/src/Core/DenseStorage.h +570 -0
  87. data/vendor/eigen/Eigen/src/Core/Diagonal.h +260 -0
  88. data/vendor/eigen/Eigen/src/Core/DiagonalMatrix.h +343 -0
  89. data/vendor/eigen/Eigen/src/Core/DiagonalProduct.h +28 -0
  90. data/vendor/eigen/Eigen/src/Core/Dot.h +318 -0
  91. data/vendor/eigen/Eigen/src/Core/EigenBase.h +159 -0
  92. data/vendor/eigen/Eigen/src/Core/ForceAlignedAccess.h +146 -0
  93. data/vendor/eigen/Eigen/src/Core/Fuzzy.h +155 -0
  94. data/vendor/eigen/Eigen/src/Core/GeneralProduct.h +455 -0
  95. data/vendor/eigen/Eigen/src/Core/GenericPacketMath.h +593 -0
  96. data/vendor/eigen/Eigen/src/Core/GlobalFunctions.h +187 -0
  97. data/vendor/eigen/Eigen/src/Core/IO.h +225 -0
  98. data/vendor/eigen/Eigen/src/Core/Inverse.h +118 -0
  99. data/vendor/eigen/Eigen/src/Core/Map.h +171 -0
  100. data/vendor/eigen/Eigen/src/Core/MapBase.h +303 -0
  101. data/vendor/eigen/Eigen/src/Core/MathFunctions.h +1415 -0
  102. data/vendor/eigen/Eigen/src/Core/MathFunctionsImpl.h +101 -0
  103. data/vendor/eigen/Eigen/src/Core/Matrix.h +459 -0
  104. data/vendor/eigen/Eigen/src/Core/MatrixBase.h +529 -0
  105. data/vendor/eigen/Eigen/src/Core/NestByValue.h +110 -0
  106. data/vendor/eigen/Eigen/src/Core/NoAlias.h +108 -0
  107. data/vendor/eigen/Eigen/src/Core/NumTraits.h +248 -0
  108. data/vendor/eigen/Eigen/src/Core/PermutationMatrix.h +633 -0
  109. data/vendor/eigen/Eigen/src/Core/PlainObjectBase.h +1035 -0
  110. data/vendor/eigen/Eigen/src/Core/Product.h +186 -0
  111. data/vendor/eigen/Eigen/src/Core/ProductEvaluators.h +1112 -0
  112. data/vendor/eigen/Eigen/src/Core/Random.h +182 -0
  113. data/vendor/eigen/Eigen/src/Core/Redux.h +505 -0
  114. data/vendor/eigen/Eigen/src/Core/Ref.h +283 -0
  115. data/vendor/eigen/Eigen/src/Core/Replicate.h +142 -0
  116. data/vendor/eigen/Eigen/src/Core/ReturnByValue.h +117 -0
  117. data/vendor/eigen/Eigen/src/Core/Reverse.h +211 -0
  118. data/vendor/eigen/Eigen/src/Core/Select.h +162 -0
  119. data/vendor/eigen/Eigen/src/Core/SelfAdjointView.h +352 -0
  120. data/vendor/eigen/Eigen/src/Core/SelfCwiseBinaryOp.h +47 -0
  121. data/vendor/eigen/Eigen/src/Core/Solve.h +188 -0
  122. data/vendor/eigen/Eigen/src/Core/SolveTriangular.h +235 -0
  123. data/vendor/eigen/Eigen/src/Core/SolverBase.h +130 -0
  124. data/vendor/eigen/Eigen/src/Core/StableNorm.h +221 -0
  125. data/vendor/eigen/Eigen/src/Core/Stride.h +111 -0
  126. data/vendor/eigen/Eigen/src/Core/Swap.h +67 -0
  127. data/vendor/eigen/Eigen/src/Core/Transpose.h +403 -0
  128. data/vendor/eigen/Eigen/src/Core/Transpositions.h +407 -0
  129. data/vendor/eigen/Eigen/src/Core/TriangularMatrix.h +983 -0
  130. data/vendor/eigen/Eigen/src/Core/VectorBlock.h +96 -0
  131. data/vendor/eigen/Eigen/src/Core/VectorwiseOp.h +695 -0
  132. data/vendor/eigen/Eigen/src/Core/Visitor.h +273 -0
  133. data/vendor/eigen/Eigen/src/Core/arch/AVX/Complex.h +451 -0
  134. data/vendor/eigen/Eigen/src/Core/arch/AVX/MathFunctions.h +439 -0
  135. data/vendor/eigen/Eigen/src/Core/arch/AVX/PacketMath.h +637 -0
  136. data/vendor/eigen/Eigen/src/Core/arch/AVX/TypeCasting.h +51 -0
  137. data/vendor/eigen/Eigen/src/Core/arch/AVX512/MathFunctions.h +391 -0
  138. data/vendor/eigen/Eigen/src/Core/arch/AVX512/PacketMath.h +1316 -0
  139. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/Complex.h +430 -0
  140. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/MathFunctions.h +322 -0
  141. data/vendor/eigen/Eigen/src/Core/arch/AltiVec/PacketMath.h +1061 -0
  142. data/vendor/eigen/Eigen/src/Core/arch/CUDA/Complex.h +103 -0
  143. data/vendor/eigen/Eigen/src/Core/arch/CUDA/Half.h +674 -0
  144. data/vendor/eigen/Eigen/src/Core/arch/CUDA/MathFunctions.h +91 -0
  145. data/vendor/eigen/Eigen/src/Core/arch/CUDA/PacketMath.h +333 -0
  146. data/vendor/eigen/Eigen/src/Core/arch/CUDA/PacketMathHalf.h +1124 -0
  147. data/vendor/eigen/Eigen/src/Core/arch/CUDA/TypeCasting.h +212 -0
  148. data/vendor/eigen/Eigen/src/Core/arch/Default/ConjHelper.h +29 -0
  149. data/vendor/eigen/Eigen/src/Core/arch/Default/Settings.h +49 -0
  150. data/vendor/eigen/Eigen/src/Core/arch/NEON/Complex.h +490 -0
  151. data/vendor/eigen/Eigen/src/Core/arch/NEON/MathFunctions.h +91 -0
  152. data/vendor/eigen/Eigen/src/Core/arch/NEON/PacketMath.h +760 -0
  153. data/vendor/eigen/Eigen/src/Core/arch/SSE/Complex.h +471 -0
  154. data/vendor/eigen/Eigen/src/Core/arch/SSE/MathFunctions.h +562 -0
  155. data/vendor/eigen/Eigen/src/Core/arch/SSE/PacketMath.h +895 -0
  156. data/vendor/eigen/Eigen/src/Core/arch/SSE/TypeCasting.h +77 -0
  157. data/vendor/eigen/Eigen/src/Core/arch/ZVector/Complex.h +397 -0
  158. data/vendor/eigen/Eigen/src/Core/arch/ZVector/MathFunctions.h +137 -0
  159. data/vendor/eigen/Eigen/src/Core/arch/ZVector/PacketMath.h +945 -0
  160. data/vendor/eigen/Eigen/src/Core/functors/AssignmentFunctors.h +168 -0
  161. data/vendor/eigen/Eigen/src/Core/functors/BinaryFunctors.h +475 -0
  162. data/vendor/eigen/Eigen/src/Core/functors/NullaryFunctors.h +188 -0
  163. data/vendor/eigen/Eigen/src/Core/functors/StlFunctors.h +136 -0
  164. data/vendor/eigen/Eigen/src/Core/functors/TernaryFunctors.h +25 -0
  165. data/vendor/eigen/Eigen/src/Core/functors/UnaryFunctors.h +792 -0
  166. data/vendor/eigen/Eigen/src/Core/products/GeneralBlockPanelKernel.h +2156 -0
  167. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrix.h +492 -0
  168. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrixTriangular.h +311 -0
  169. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h +145 -0
  170. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrix_BLAS.h +122 -0
  171. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixVector.h +619 -0
  172. data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixVector_BLAS.h +136 -0
  173. data/vendor/eigen/Eigen/src/Core/products/Parallelizer.h +163 -0
  174. data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixMatrix.h +521 -0
  175. data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixMatrix_BLAS.h +287 -0
  176. data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixVector.h +260 -0
  177. data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixVector_BLAS.h +118 -0
  178. data/vendor/eigen/Eigen/src/Core/products/SelfadjointProduct.h +133 -0
  179. data/vendor/eigen/Eigen/src/Core/products/SelfadjointRank2Update.h +93 -0
  180. data/vendor/eigen/Eigen/src/Core/products/TriangularMatrixMatrix.h +466 -0
  181. data/vendor/eigen/Eigen/src/Core/products/TriangularMatrixMatrix_BLAS.h +315 -0
  182. data/vendor/eigen/Eigen/src/Core/products/TriangularMatrixVector.h +350 -0
  183. data/vendor/eigen/Eigen/src/Core/products/TriangularMatrixVector_BLAS.h +255 -0
  184. data/vendor/eigen/Eigen/src/Core/products/TriangularSolverMatrix.h +335 -0
  185. data/vendor/eigen/Eigen/src/Core/products/TriangularSolverMatrix_BLAS.h +163 -0
  186. data/vendor/eigen/Eigen/src/Core/products/TriangularSolverVector.h +145 -0
  187. data/vendor/eigen/Eigen/src/Core/util/BlasUtil.h +398 -0
  188. data/vendor/eigen/Eigen/src/Core/util/Constants.h +547 -0
  189. data/vendor/eigen/Eigen/src/Core/util/DisableStupidWarnings.h +83 -0
  190. data/vendor/eigen/Eigen/src/Core/util/ForwardDeclarations.h +302 -0
  191. data/vendor/eigen/Eigen/src/Core/util/MKL_support.h +130 -0
  192. data/vendor/eigen/Eigen/src/Core/util/Macros.h +1001 -0
  193. data/vendor/eigen/Eigen/src/Core/util/Memory.h +993 -0
  194. data/vendor/eigen/Eigen/src/Core/util/Meta.h +534 -0
  195. data/vendor/eigen/Eigen/src/Core/util/NonMPL2.h +3 -0
  196. data/vendor/eigen/Eigen/src/Core/util/ReenableStupidWarnings.h +27 -0
  197. data/vendor/eigen/Eigen/src/Core/util/StaticAssert.h +218 -0
  198. data/vendor/eigen/Eigen/src/Core/util/XprHelper.h +821 -0
  199. data/vendor/eigen/Eigen/src/Eigenvalues/ComplexEigenSolver.h +346 -0
  200. data/vendor/eigen/Eigen/src/Eigenvalues/ComplexSchur.h +459 -0
  201. data/vendor/eigen/Eigen/src/Eigenvalues/ComplexSchur_LAPACKE.h +91 -0
  202. data/vendor/eigen/Eigen/src/Eigenvalues/EigenSolver.h +622 -0
  203. data/vendor/eigen/Eigen/src/Eigenvalues/GeneralizedEigenSolver.h +418 -0
  204. data/vendor/eigen/Eigen/src/Eigenvalues/GeneralizedSelfAdjointEigenSolver.h +226 -0
  205. data/vendor/eigen/Eigen/src/Eigenvalues/HessenbergDecomposition.h +374 -0
  206. data/vendor/eigen/Eigen/src/Eigenvalues/MatrixBaseEigenvalues.h +158 -0
  207. data/vendor/eigen/Eigen/src/Eigenvalues/RealQZ.h +654 -0
  208. data/vendor/eigen/Eigen/src/Eigenvalues/RealSchur.h +546 -0
  209. data/vendor/eigen/Eigen/src/Eigenvalues/RealSchur_LAPACKE.h +77 -0
  210. data/vendor/eigen/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h +870 -0
  211. data/vendor/eigen/Eigen/src/Eigenvalues/SelfAdjointEigenSolver_LAPACKE.h +87 -0
  212. data/vendor/eigen/Eigen/src/Eigenvalues/Tridiagonalization.h +556 -0
  213. data/vendor/eigen/Eigen/src/Geometry/AlignedBox.h +392 -0
  214. data/vendor/eigen/Eigen/src/Geometry/AngleAxis.h +247 -0
  215. data/vendor/eigen/Eigen/src/Geometry/EulerAngles.h +114 -0
  216. data/vendor/eigen/Eigen/src/Geometry/Homogeneous.h +497 -0
  217. data/vendor/eigen/Eigen/src/Geometry/Hyperplane.h +282 -0
  218. data/vendor/eigen/Eigen/src/Geometry/OrthoMethods.h +234 -0
  219. data/vendor/eigen/Eigen/src/Geometry/ParametrizedLine.h +195 -0
  220. data/vendor/eigen/Eigen/src/Geometry/Quaternion.h +814 -0
  221. data/vendor/eigen/Eigen/src/Geometry/Rotation2D.h +199 -0
  222. data/vendor/eigen/Eigen/src/Geometry/RotationBase.h +206 -0
  223. data/vendor/eigen/Eigen/src/Geometry/Scaling.h +170 -0
  224. data/vendor/eigen/Eigen/src/Geometry/Transform.h +1542 -0
  225. data/vendor/eigen/Eigen/src/Geometry/Translation.h +208 -0
  226. data/vendor/eigen/Eigen/src/Geometry/Umeyama.h +166 -0
  227. data/vendor/eigen/Eigen/src/Geometry/arch/Geometry_SSE.h +161 -0
  228. data/vendor/eigen/Eigen/src/Householder/BlockHouseholder.h +103 -0
  229. data/vendor/eigen/Eigen/src/Householder/Householder.h +172 -0
  230. data/vendor/eigen/Eigen/src/Householder/HouseholderSequence.h +470 -0
  231. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/BasicPreconditioners.h +226 -0
  232. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/BiCGSTAB.h +228 -0
  233. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h +246 -0
  234. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/IncompleteCholesky.h +400 -0
  235. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/IncompleteLUT.h +462 -0
  236. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/IterativeSolverBase.h +394 -0
  237. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/LeastSquareConjugateGradient.h +216 -0
  238. data/vendor/eigen/Eigen/src/IterativeLinearSolvers/SolveWithGuess.h +115 -0
  239. data/vendor/eigen/Eigen/src/Jacobi/Jacobi.h +462 -0
  240. data/vendor/eigen/Eigen/src/LU/Determinant.h +101 -0
  241. data/vendor/eigen/Eigen/src/LU/FullPivLU.h +891 -0
  242. data/vendor/eigen/Eigen/src/LU/InverseImpl.h +415 -0
  243. data/vendor/eigen/Eigen/src/LU/PartialPivLU.h +611 -0
  244. data/vendor/eigen/Eigen/src/LU/PartialPivLU_LAPACKE.h +83 -0
  245. data/vendor/eigen/Eigen/src/LU/arch/Inverse_SSE.h +338 -0
  246. data/vendor/eigen/Eigen/src/MetisSupport/MetisSupport.h +137 -0
  247. data/vendor/eigen/Eigen/src/OrderingMethods/Amd.h +445 -0
  248. data/vendor/eigen/Eigen/src/OrderingMethods/Eigen_Colamd.h +1843 -0
  249. data/vendor/eigen/Eigen/src/OrderingMethods/Ordering.h +157 -0
  250. data/vendor/eigen/Eigen/src/PaStiXSupport/PaStiXSupport.h +678 -0
  251. data/vendor/eigen/Eigen/src/PardisoSupport/PardisoSupport.h +543 -0
  252. data/vendor/eigen/Eigen/src/QR/ColPivHouseholderQR.h +653 -0
  253. data/vendor/eigen/Eigen/src/QR/ColPivHouseholderQR_LAPACKE.h +97 -0
  254. data/vendor/eigen/Eigen/src/QR/CompleteOrthogonalDecomposition.h +562 -0
  255. data/vendor/eigen/Eigen/src/QR/FullPivHouseholderQR.h +676 -0
  256. data/vendor/eigen/Eigen/src/QR/HouseholderQR.h +409 -0
  257. data/vendor/eigen/Eigen/src/QR/HouseholderQR_LAPACKE.h +68 -0
  258. data/vendor/eigen/Eigen/src/SPQRSupport/SuiteSparseQRSupport.h +313 -0
  259. data/vendor/eigen/Eigen/src/SVD/BDCSVD.h +1246 -0
  260. data/vendor/eigen/Eigen/src/SVD/JacobiSVD.h +804 -0
  261. data/vendor/eigen/Eigen/src/SVD/JacobiSVD_LAPACKE.h +91 -0
  262. data/vendor/eigen/Eigen/src/SVD/SVDBase.h +315 -0
  263. data/vendor/eigen/Eigen/src/SVD/UpperBidiagonalization.h +414 -0
  264. data/vendor/eigen/Eigen/src/SparseCholesky/SimplicialCholesky.h +689 -0
  265. data/vendor/eigen/Eigen/src/SparseCholesky/SimplicialCholesky_impl.h +199 -0
  266. data/vendor/eigen/Eigen/src/SparseCore/AmbiVector.h +377 -0
  267. data/vendor/eigen/Eigen/src/SparseCore/CompressedStorage.h +258 -0
  268. data/vendor/eigen/Eigen/src/SparseCore/ConservativeSparseSparseProduct.h +352 -0
  269. data/vendor/eigen/Eigen/src/SparseCore/MappedSparseMatrix.h +67 -0
  270. data/vendor/eigen/Eigen/src/SparseCore/SparseAssign.h +216 -0
  271. data/vendor/eigen/Eigen/src/SparseCore/SparseBlock.h +603 -0
  272. data/vendor/eigen/Eigen/src/SparseCore/SparseColEtree.h +206 -0
  273. data/vendor/eigen/Eigen/src/SparseCore/SparseCompressedBase.h +341 -0
  274. data/vendor/eigen/Eigen/src/SparseCore/SparseCwiseBinaryOp.h +726 -0
  275. data/vendor/eigen/Eigen/src/SparseCore/SparseCwiseUnaryOp.h +148 -0
  276. data/vendor/eigen/Eigen/src/SparseCore/SparseDenseProduct.h +320 -0
  277. data/vendor/eigen/Eigen/src/SparseCore/SparseDiagonalProduct.h +138 -0
  278. data/vendor/eigen/Eigen/src/SparseCore/SparseDot.h +98 -0
  279. data/vendor/eigen/Eigen/src/SparseCore/SparseFuzzy.h +29 -0
  280. data/vendor/eigen/Eigen/src/SparseCore/SparseMap.h +305 -0
  281. data/vendor/eigen/Eigen/src/SparseCore/SparseMatrix.h +1403 -0
  282. data/vendor/eigen/Eigen/src/SparseCore/SparseMatrixBase.h +405 -0
  283. data/vendor/eigen/Eigen/src/SparseCore/SparsePermutation.h +178 -0
  284. data/vendor/eigen/Eigen/src/SparseCore/SparseProduct.h +169 -0
  285. data/vendor/eigen/Eigen/src/SparseCore/SparseRedux.h +49 -0
  286. data/vendor/eigen/Eigen/src/SparseCore/SparseRef.h +397 -0
  287. data/vendor/eigen/Eigen/src/SparseCore/SparseSelfAdjointView.h +656 -0
  288. data/vendor/eigen/Eigen/src/SparseCore/SparseSolverBase.h +124 -0
  289. data/vendor/eigen/Eigen/src/SparseCore/SparseSparseProductWithPruning.h +198 -0
  290. data/vendor/eigen/Eigen/src/SparseCore/SparseTranspose.h +92 -0
  291. data/vendor/eigen/Eigen/src/SparseCore/SparseTriangularView.h +189 -0
  292. data/vendor/eigen/Eigen/src/SparseCore/SparseUtil.h +178 -0
  293. data/vendor/eigen/Eigen/src/SparseCore/SparseVector.h +478 -0
  294. data/vendor/eigen/Eigen/src/SparseCore/SparseView.h +253 -0
  295. data/vendor/eigen/Eigen/src/SparseCore/TriangularSolver.h +315 -0
  296. data/vendor/eigen/Eigen/src/SparseLU/SparseLU.h +773 -0
  297. data/vendor/eigen/Eigen/src/SparseLU/SparseLUImpl.h +66 -0
  298. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_Memory.h +226 -0
  299. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_Structs.h +110 -0
  300. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_SupernodalMatrix.h +301 -0
  301. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_Utils.h +80 -0
  302. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_column_bmod.h +181 -0
  303. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_column_dfs.h +179 -0
  304. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_copy_to_ucol.h +107 -0
  305. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_gemm_kernel.h +280 -0
  306. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_heap_relax_snode.h +126 -0
  307. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_kernel_bmod.h +130 -0
  308. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_panel_bmod.h +223 -0
  309. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_panel_dfs.h +258 -0
  310. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_pivotL.h +137 -0
  311. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_pruneL.h +136 -0
  312. data/vendor/eigen/Eigen/src/SparseLU/SparseLU_relax_snode.h +83 -0
  313. data/vendor/eigen/Eigen/src/SparseQR/SparseQR.h +745 -0
  314. data/vendor/eigen/Eigen/src/StlSupport/StdDeque.h +126 -0
  315. data/vendor/eigen/Eigen/src/StlSupport/StdList.h +106 -0
  316. data/vendor/eigen/Eigen/src/StlSupport/StdVector.h +131 -0
  317. data/vendor/eigen/Eigen/src/StlSupport/details.h +84 -0
  318. data/vendor/eigen/Eigen/src/SuperLUSupport/SuperLUSupport.h +1027 -0
  319. data/vendor/eigen/Eigen/src/UmfPackSupport/UmfPackSupport.h +506 -0
  320. data/vendor/eigen/Eigen/src/misc/Image.h +82 -0
  321. data/vendor/eigen/Eigen/src/misc/Kernel.h +79 -0
  322. data/vendor/eigen/Eigen/src/misc/RealSvd2x2.h +55 -0
  323. data/vendor/eigen/Eigen/src/misc/blas.h +440 -0
  324. data/vendor/eigen/Eigen/src/misc/lapack.h +152 -0
  325. data/vendor/eigen/Eigen/src/misc/lapacke.h +16291 -0
  326. data/vendor/eigen/Eigen/src/misc/lapacke_mangling.h +17 -0
  327. data/vendor/eigen/Eigen/src/plugins/ArrayCwiseBinaryOps.h +332 -0
  328. data/vendor/eigen/Eigen/src/plugins/ArrayCwiseUnaryOps.h +552 -0
  329. data/vendor/eigen/Eigen/src/plugins/BlockMethods.h +1058 -0
  330. data/vendor/eigen/Eigen/src/plugins/CommonCwiseBinaryOps.h +115 -0
  331. data/vendor/eigen/Eigen/src/plugins/CommonCwiseUnaryOps.h +163 -0
  332. data/vendor/eigen/Eigen/src/plugins/MatrixCwiseBinaryOps.h +152 -0
  333. data/vendor/eigen/Eigen/src/plugins/MatrixCwiseUnaryOps.h +85 -0
  334. data/vendor/eigen/README.md +3 -0
  335. data/vendor/eigen/bench/README.txt +55 -0
  336. data/vendor/eigen/bench/btl/COPYING +340 -0
  337. data/vendor/eigen/bench/btl/README +154 -0
  338. data/vendor/eigen/bench/tensors/README +21 -0
  339. data/vendor/eigen/blas/README.txt +6 -0
  340. data/vendor/eigen/demos/mandelbrot/README +10 -0
  341. data/vendor/eigen/demos/mix_eigen_and_c/README +9 -0
  342. data/vendor/eigen/demos/opengl/README +13 -0
  343. data/vendor/eigen/unsupported/Eigen/CXX11/src/Tensor/README.md +1760 -0
  344. data/vendor/eigen/unsupported/README.txt +50 -0
  345. data/vendor/tomotopy/LICENSE +21 -0
  346. data/vendor/tomotopy/README.kr.rst +375 -0
  347. data/vendor/tomotopy/README.rst +382 -0
  348. data/vendor/tomotopy/src/Labeling/FoRelevance.cpp +362 -0
  349. data/vendor/tomotopy/src/Labeling/FoRelevance.h +88 -0
  350. data/vendor/tomotopy/src/Labeling/Labeler.h +50 -0
  351. data/vendor/tomotopy/src/TopicModel/CT.h +37 -0
  352. data/vendor/tomotopy/src/TopicModel/CTModel.cpp +13 -0
  353. data/vendor/tomotopy/src/TopicModel/CTModel.hpp +293 -0
  354. data/vendor/tomotopy/src/TopicModel/DMR.h +51 -0
  355. data/vendor/tomotopy/src/TopicModel/DMRModel.cpp +13 -0
  356. data/vendor/tomotopy/src/TopicModel/DMRModel.hpp +374 -0
  357. data/vendor/tomotopy/src/TopicModel/DT.h +65 -0
  358. data/vendor/tomotopy/src/TopicModel/DTM.h +22 -0
  359. data/vendor/tomotopy/src/TopicModel/DTModel.cpp +15 -0
  360. data/vendor/tomotopy/src/TopicModel/DTModel.hpp +572 -0
  361. data/vendor/tomotopy/src/TopicModel/GDMR.h +37 -0
  362. data/vendor/tomotopy/src/TopicModel/GDMRModel.cpp +14 -0
  363. data/vendor/tomotopy/src/TopicModel/GDMRModel.hpp +485 -0
  364. data/vendor/tomotopy/src/TopicModel/HDP.h +74 -0
  365. data/vendor/tomotopy/src/TopicModel/HDPModel.cpp +13 -0
  366. data/vendor/tomotopy/src/TopicModel/HDPModel.hpp +592 -0
  367. data/vendor/tomotopy/src/TopicModel/HLDA.h +40 -0
  368. data/vendor/tomotopy/src/TopicModel/HLDAModel.cpp +13 -0
  369. data/vendor/tomotopy/src/TopicModel/HLDAModel.hpp +681 -0
  370. data/vendor/tomotopy/src/TopicModel/HPA.h +27 -0
  371. data/vendor/tomotopy/src/TopicModel/HPAModel.cpp +21 -0
  372. data/vendor/tomotopy/src/TopicModel/HPAModel.hpp +588 -0
  373. data/vendor/tomotopy/src/TopicModel/LDA.h +144 -0
  374. data/vendor/tomotopy/src/TopicModel/LDACVB0Model.hpp +442 -0
  375. data/vendor/tomotopy/src/TopicModel/LDAModel.cpp +13 -0
  376. data/vendor/tomotopy/src/TopicModel/LDAModel.hpp +1058 -0
  377. data/vendor/tomotopy/src/TopicModel/LLDA.h +45 -0
  378. data/vendor/tomotopy/src/TopicModel/LLDAModel.cpp +13 -0
  379. data/vendor/tomotopy/src/TopicModel/LLDAModel.hpp +203 -0
  380. data/vendor/tomotopy/src/TopicModel/MGLDA.h +63 -0
  381. data/vendor/tomotopy/src/TopicModel/MGLDAModel.cpp +17 -0
  382. data/vendor/tomotopy/src/TopicModel/MGLDAModel.hpp +558 -0
  383. data/vendor/tomotopy/src/TopicModel/PA.h +43 -0
  384. data/vendor/tomotopy/src/TopicModel/PAModel.cpp +13 -0
  385. data/vendor/tomotopy/src/TopicModel/PAModel.hpp +467 -0
  386. data/vendor/tomotopy/src/TopicModel/PLDA.h +17 -0
  387. data/vendor/tomotopy/src/TopicModel/PLDAModel.cpp +13 -0
  388. data/vendor/tomotopy/src/TopicModel/PLDAModel.hpp +214 -0
  389. data/vendor/tomotopy/src/TopicModel/SLDA.h +54 -0
  390. data/vendor/tomotopy/src/TopicModel/SLDAModel.cpp +17 -0
  391. data/vendor/tomotopy/src/TopicModel/SLDAModel.hpp +456 -0
  392. data/vendor/tomotopy/src/TopicModel/TopicModel.hpp +692 -0
  393. data/vendor/tomotopy/src/Utils/AliasMethod.hpp +169 -0
  394. data/vendor/tomotopy/src/Utils/Dictionary.h +80 -0
  395. data/vendor/tomotopy/src/Utils/EigenAddonOps.hpp +181 -0
  396. data/vendor/tomotopy/src/Utils/LBFGS.h +202 -0
  397. data/vendor/tomotopy/src/Utils/LBFGS/LineSearchBacktracking.h +120 -0
  398. data/vendor/tomotopy/src/Utils/LBFGS/LineSearchBracketing.h +122 -0
  399. data/vendor/tomotopy/src/Utils/LBFGS/Param.h +213 -0
  400. data/vendor/tomotopy/src/Utils/LUT.hpp +82 -0
  401. data/vendor/tomotopy/src/Utils/MultiNormalDistribution.hpp +69 -0
  402. data/vendor/tomotopy/src/Utils/PolyaGamma.hpp +200 -0
  403. data/vendor/tomotopy/src/Utils/PolyaGammaHybrid.hpp +672 -0
  404. data/vendor/tomotopy/src/Utils/ThreadPool.hpp +150 -0
  405. data/vendor/tomotopy/src/Utils/Trie.hpp +220 -0
  406. data/vendor/tomotopy/src/Utils/TruncMultiNormal.hpp +94 -0
  407. data/vendor/tomotopy/src/Utils/Utils.hpp +337 -0
  408. data/vendor/tomotopy/src/Utils/avx_gamma.h +46 -0
  409. data/vendor/tomotopy/src/Utils/avx_mathfun.h +736 -0
  410. data/vendor/tomotopy/src/Utils/exception.h +28 -0
  411. data/vendor/tomotopy/src/Utils/math.h +281 -0
  412. data/vendor/tomotopy/src/Utils/rtnorm.hpp +2690 -0
  413. data/vendor/tomotopy/src/Utils/sample.hpp +192 -0
  414. data/vendor/tomotopy/src/Utils/serializer.hpp +695 -0
  415. data/vendor/tomotopy/src/Utils/slp.hpp +131 -0
  416. data/vendor/tomotopy/src/Utils/sse_gamma.h +48 -0
  417. data/vendor/tomotopy/src/Utils/sse_mathfun.h +710 -0
  418. data/vendor/tomotopy/src/Utils/text.hpp +49 -0
  419. data/vendor/tomotopy/src/Utils/tvector.hpp +543 -0
  420. metadata +531 -0
@@ -0,0 +1,49 @@
1
+ #pragma once
2
+
3
+ #include <vector>
4
+ #include <string>
5
+ #include <iostream>
6
+ #include <iterator>
7
+ #include <algorithm>
8
+ #include <cstdio>
9
+
10
+ namespace tomoto
11
+ {
12
+ namespace text
13
+ {
14
+ template<typename ... _Args>
15
+ std::string format(const std::string& format, _Args ... args)
16
+ {
17
+ size_t size = snprintf(nullptr, 0, format.c_str(), args ...) + 1;
18
+ std::vector<char> buf(size);
19
+ snprintf(buf.data(), size, format.c_str(), args ...);
20
+ return std::string{ buf.data(), buf.data() + size - 1 };
21
+ }
22
+
23
+ template<class _Iter, class _Target = decltype(*_Iter{}) >
24
+ std::string join(_Iter first, _Iter last, const std::string& delimiter = ",")
25
+ {
26
+ if (first == last) return "";
27
+ std::ostringstream stream;
28
+ std::copy(first, last, std::ostream_iterator<_Target>(stream, delimiter.c_str()));
29
+ std::string s = stream.str();
30
+ s.erase(s.end() - delimiter.size(), s.end());
31
+ return s;
32
+ }
33
+
34
+ inline std::vector<std::string> split(const std::string& str, const std::string& delim)
35
+ {
36
+ std::vector<std::string> tokens;
37
+ size_t prev = 0, pos = 0;
38
+ do
39
+ {
40
+ pos = str.find(delim, prev);
41
+ if (pos == std::string::npos) pos = str.length();
42
+ std::string token = str.substr(prev, pos - prev);
43
+ if (!token.empty()) tokens.push_back(token);
44
+ prev = pos + delim.length();
45
+ } while (pos < str.length() && prev < str.length());
46
+ return tokens;
47
+ }
48
+ }
49
+ }
@@ -0,0 +1,543 @@
1
+ #pragma once
2
+ #include <algorithm>
3
+ #include <cstring>
4
+ namespace tomoto
5
+ {
6
+ /**
7
+ tvector : std::vector-like container which can handle non-owning std::array as well as traditional std::vector function
8
+ */
9
+ template <typename T, class _Alloc = std::allocator<T>>
10
+ class tvector
11
+ {
12
+ public:
13
+ typedef T value_type;
14
+ typedef T & reference;
15
+ typedef const T & const_reference;
16
+ typedef T * pointer;
17
+ typedef const T * const_pointer;
18
+ typedef T * iterator;
19
+ typedef const T * const_iterator;
20
+ typedef std::reverse_iterator<iterator> reverse_iterator;
21
+ typedef std::reverse_iterator<const_iterator> const_reverse_iterator;
22
+ typedef std::ptrdiff_t difference_type;
23
+ typedef std::size_t size_type;
24
+
25
+ tvector() noexcept
26
+ {
27
+ _first = _Alloc{}.allocate(4);
28
+ _last = _first;
29
+ _rsvEnd = _first + 4;
30
+ }
31
+
32
+ tvector(std::nullptr_t) noexcept
33
+ {
34
+
35
+ }
36
+
37
+ // non-owning, just pointing constructor
38
+ tvector(pointer ptr, size_type size) noexcept
39
+ {
40
+ _first = ptr;
41
+ _last = ptr + size;
42
+ _rsvEnd = nullptr;
43
+ }
44
+
45
+ explicit tvector(size_type n, const T &val = T())
46
+ {
47
+ _first = _Alloc{}.allocate(n);
48
+ for (size_type i = 0; i < n; ++i) _first[i] = val;
49
+ _last = _first + n;
50
+ _rsvEnd = _last;
51
+ }
52
+
53
+ tvector(iterator first, iterator last)
54
+ {
55
+ size_type count = last - first;
56
+ _first = _Alloc{}.allocate(count);
57
+ for (size_type i = 0; i < count; ++i, ++first)
58
+ _first[i] = *first;
59
+ _last = _first + count;
60
+ _rsvEnd = _last;
61
+ }
62
+
63
+ tvector(std::initializer_list<T> lst)
64
+ {
65
+ _first = _Alloc{}.allocate(lst.size());
66
+ _last = _first;
67
+ for (auto &item : lst) *_last++ = item;
68
+ _rsvEnd = _last;
69
+ }
70
+
71
+ tvector(const tvector<T> & other)
72
+ {
73
+ if (other._rsvEnd)
74
+ {
75
+ _first = _Alloc{}.allocate(other.capacity());
76
+ _last = _first;
77
+ for (size_type i = 0; i < other.size(); ++i) *_last++ = other._first[i];
78
+ _rsvEnd = _first + other.capacity();
79
+ }
80
+ else
81
+ {
82
+ if (_rsvEnd) _Alloc{}.deallocate(_first, capacity());
83
+ _first = other._first;
84
+ _last = other._last;
85
+ _rsvEnd = nullptr;
86
+ }
87
+ }
88
+
89
+ tvector(tvector<T> && other) noexcept
90
+ {
91
+ swap(other);
92
+ }
93
+
94
+ ~tvector()
95
+ {
96
+ if (_first && _rsvEnd) _Alloc{}.deallocate(_first, capacity());
97
+ }
98
+
99
+ tvector<T> & operator = (const tvector<T> &other)
100
+ {
101
+ if (other._rsvEnd)
102
+ {
103
+ buy(other.size());
104
+ for (size_type i = 0; i < other.size(); ++i)
105
+ _first[i] = other._first[i];
106
+ _last = _first + other.size();
107
+ }
108
+ else
109
+ {
110
+ if (_rsvEnd) _Alloc{}.deallocate(_first, capacity());
111
+ _first = other._first;
112
+ _last = other._last;
113
+ _rsvEnd = nullptr;
114
+ }
115
+ return *this;
116
+ }
117
+
118
+ tvector<T> & operator = (tvector<T> &&other)
119
+ {
120
+ swap(other);
121
+ return *this;
122
+ }
123
+
124
+ tvector<T> & operator = (std::initializer_list<T> lst)
125
+ {
126
+ buy(lst.size());
127
+ _last = _first;
128
+ for (auto &item : lst)
129
+ *_last++ = item;
130
+ return *this;
131
+ }
132
+
133
+ void assign(size_type count, const T &value)
134
+ {
135
+ size_type i;
136
+ buy(count);
137
+ for (i = 0; i < count; ++i)
138
+ _first[i] = value;
139
+ _last = _first + count;
140
+ }
141
+
142
+ void assign(iterator first, iterator last)
143
+ {
144
+ size_type count = last - first;
145
+ buy(count);
146
+ for (size_type i = 0; i < count; ++i, ++first)
147
+ _first[i] = *first;
148
+ _last = _first + count;
149
+ }
150
+
151
+ void assign(std::initializer_list<T> lst)
152
+ {
153
+ size_type count = lst.size();
154
+ buy(count);
155
+ _last = _first;
156
+ for (auto &item : lst)
157
+ *_last++ = item;
158
+ }
159
+
160
+ // iterators:
161
+ iterator begin() noexcept { return _first; }
162
+
163
+ const_iterator begin() const noexcept { return _first; }
164
+
165
+ const_iterator cbegin() const noexcept { return _first; }
166
+
167
+ iterator end() noexcept { return _last; }
168
+
169
+ const_iterator end() const noexcept { return _last; }
170
+
171
+ const_iterator cend() const noexcept { return _last; }
172
+
173
+ reverse_iterator rbegin() noexcept { return reverse_iterator(_last); }
174
+
175
+ const_reverse_iterator crbegin() const noexcept { return reverse_iterator(_last); }
176
+
177
+ reverse_iterator rend() noexcept { return reverse_iterator(_first); }
178
+
179
+ const_reverse_iterator crend() const noexcept { return reverse_iterator(_first); }
180
+
181
+ // 23.3.11.3, capacity:
182
+ bool empty() const noexcept
183
+ {
184
+ return _last == _first;
185
+ }
186
+
187
+ size_type size() const noexcept
188
+ {
189
+ return _last - _first;
190
+ }
191
+
192
+
193
+ size_type max_size() const noexcept
194
+ {
195
+ return 0xFFFFFFFF;
196
+ }
197
+
198
+
199
+ size_type capacity() const noexcept
200
+ {
201
+ return _rsvEnd - _first;
202
+ }
203
+
204
+ void resize(size_type sz, const T &c = T())
205
+ {
206
+ if (sz > size())
207
+ {
208
+ if (!isOwner()) throw std::out_of_range("cannot increase size of non-owning mode");
209
+ reallocate(sz);
210
+ for (size_type i = size(); i < sz; ++i) _first[i] = c;
211
+ }
212
+ else
213
+ {
214
+ for (size_type i = size(); i < sz; ++i) _first[i].~T();
215
+ }
216
+ _last = _first + sz;
217
+ }
218
+
219
+ void reserve(size_type _sz)
220
+ {
221
+ if (_sz > capacity()) {
222
+ reallocate(_sz);
223
+ }
224
+ }
225
+
226
+ void shrink_to_fit()
227
+ {
228
+ reallocate(size());
229
+ }
230
+
231
+ // element access
232
+ reference operator [](size_type idx)
233
+ {
234
+ return _first[idx];
235
+ }
236
+
237
+ const_reference operator [](size_type idx) const
238
+ {
239
+ return _first[idx];
240
+ }
241
+
242
+ reference at(size_type pos)
243
+ {
244
+ if (pos < size())
245
+ return _first[pos];
246
+ else
247
+ throw std::out_of_range("accessed position is out of range");
248
+ }
249
+
250
+ const_reference at(size_type pos) const
251
+ {
252
+ if (pos < size())
253
+ return _first[pos];
254
+ else
255
+ throw std::out_of_range("accessed position is out of range");
256
+ }
257
+
258
+ reference front()
259
+ {
260
+ return _first[0];
261
+ }
262
+
263
+ const_reference front() const
264
+ {
265
+ return _first[0];
266
+ }
267
+
268
+ reference back()
269
+ {
270
+ return _last[-1];
271
+ }
272
+
273
+ const_reference back() const
274
+ {
275
+ return _last[-1];
276
+ }
277
+
278
+ // 23.3.11.4, data access:
279
+ T * data() noexcept
280
+ {
281
+ return _first;
282
+ }
283
+
284
+ const T * data() const noexcept
285
+ {
286
+ return _first;
287
+ }
288
+
289
+ bool isOwner() const noexcept
290
+ {
291
+ return _rsvEnd;
292
+ }
293
+
294
+ // 23.3.11.5, modifiers:
295
+ template <class ... Args>
296
+ void emplace_back(Args && ... args)
297
+ {
298
+ buy(size() + 1);
299
+ *_last++ = std::move(T(std::forward<Args>(args) ...));
300
+ }
301
+
302
+ void push_back(const T &val)
303
+ {
304
+ buy(size() + 1);
305
+ *_last++ = val;
306
+ }
307
+
308
+ void push_back(T &&val)
309
+ {
310
+ buy(size() + 1);
311
+ *_last++ = std::move(val);
312
+ }
313
+
314
+ void pop_back()
315
+ {
316
+ (_last--)[-1].~T();
317
+ }
318
+
319
+ template <class ... Args>
320
+ iterator emplace(const_iterator it, Args && ... args)
321
+ {
322
+ iterator iit = &_first[it - _first];
323
+ buy(size() + 1);
324
+ memmove(iit + 1, iit, (size() - (it - _first)) * sizeof(T));
325
+ (*iit) = std::move(T(std::forward<Args>(args) ...));
326
+ _last++;
327
+ return iit;
328
+ }
329
+
330
+ iterator insert(const_iterator it, const T &val)
331
+ {
332
+ iterator iit = &_first[it - _first];
333
+ buy(size() + 1);
334
+ memmove(iit + 1, iit, (size() - (it - _first)) * sizeof(T));
335
+ (*iit) = val;
336
+ _last++;
337
+ return iit;
338
+ }
339
+
340
+
341
+ iterator insert(const_iterator it, T &&val)
342
+ {
343
+ iterator iit = &_first[it - _first];
344
+ buy(size() + 1);
345
+ memmove(iit + 1, iit, (size() - (it - _first)) * sizeof(T));
346
+ (*iit) = std::move(val);
347
+ _last++;
348
+ return iit;
349
+ }
350
+
351
+
352
+ iterator insert(const_iterator it, size_type cnt, const T &val)
353
+ {
354
+ iterator f = &_first[it - _first];
355
+ if (!cnt) return f;
356
+ buy(size() + cnt);
357
+ memmove(f + cnt, f, (size() - (it - _first)) * sizeof(T));
358
+ _last += cnt;
359
+ for (iterator it = f; cnt--; ++it)
360
+ (*it) = val;
361
+ return f;
362
+ }
363
+
364
+
365
+ template <class InputIt>
366
+ iterator insert(const_iterator it, InputIt first, InputIt last)
367
+ {
368
+ iterator f = &_first[it - _first];
369
+ size_type cnt = last - first;
370
+ if (!cnt) return f;
371
+ buy(size() + cnt);
372
+ memmove(f + cnt, f, (size() - (it - _first)) * sizeof(T));
373
+ for (iterator it = f; first != last; ++it, ++first)
374
+ (*it) = *first;
375
+ _last += cnt;
376
+ return f;
377
+ }
378
+
379
+ iterator insert(const_iterator it, std::initializer_list<T> lst)
380
+ {
381
+ size_type cnt = lst.size();
382
+ iterator f = &_first[it - _first];
383
+ if (!cnt) return f;
384
+ buy(size() + cnt);
385
+ memmove(f + cnt, f, (size() - (it - _first)) * sizeof(T));
386
+ iterator iit = f;
387
+ for (auto &item : lst) {
388
+ (*iit) = item;
389
+ ++iit;
390
+ }
391
+ _last += cnt;
392
+ return f;
393
+ }
394
+
395
+ iterator erase(const_iterator it)
396
+ {
397
+ iterator iit = &_first[it - _first];
398
+ (*iit).~T();
399
+ memmove(iit, iit + 1, (size() - (it - _first) - 1) * sizeof(T));
400
+ --_last;
401
+ return iit;
402
+ }
403
+
404
+ iterator erase(const_iterator first, const_iterator last)
405
+ {
406
+ iterator f = &_first[first - _first];
407
+ if (first == last) return f;
408
+ for (; first != last; ++first)
409
+ (*first).~T();
410
+ memmove(f, last, (size() - (last - _first)) * sizeof(T));
411
+ _last -= last - first;
412
+ return f;
413
+ }
414
+
415
+ void swap(tvector<T> &rhs)
416
+ {
417
+ std::swap(_first, rhs._first);
418
+ std::swap(_last, rhs._last);
419
+ std::swap(_rsvEnd, rhs._rsvEnd);
420
+ }
421
+
422
+ void clear() noexcept
423
+ {
424
+ for (size_type i = 0; i < size(); ++i) _first[i].~T();
425
+ _last = _first;
426
+ }
427
+
428
+ bool operator == (const tvector<T> &rhs) const
429
+ {
430
+ if (size() != rhs.size()) return false;
431
+ for (size_type i = 0; i < size(); ++i)
432
+ if (_first[i] != rhs._first[i])
433
+ return false;
434
+ return true;
435
+ }
436
+
437
+
438
+ bool operator != (const tvector<T> &rhs) const
439
+ {
440
+ return !operator==(rhs);
441
+ }
442
+
443
+
444
+ bool operator < (const tvector<T> &rhs) const
445
+ {
446
+ size_type i, ub = std::min(size(), rhs.size());
447
+ for (i = 0; i < ub; ++i)
448
+ if (_first[i] != rhs._first[i])
449
+ return _first[i] < rhs._first[i];
450
+ return size() < rhs.size();
451
+ }
452
+
453
+
454
+ bool operator <= (const tvector<T> &rhs) const
455
+ {
456
+ size_type i, ub = std::min(size(), rhs.size());
457
+ for (i = 0; i < ub; ++i)
458
+ if (_first[i] != rhs._first[i])
459
+ return _first[i] < rhs._first[i];
460
+ return size() < rhs.size();
461
+ }
462
+
463
+
464
+ bool operator > (const tvector<T> &rhs) const
465
+ {
466
+ size_type i, ub = std::min(size(), rhs.size());
467
+ for (i = 0; i < ub; ++i)
468
+ if (_first[i] != rhs._first[i])
469
+ return _first[i] > rhs._first[i];
470
+ return size() > rhs.size();
471
+ }
472
+
473
+
474
+ bool operator >= (const tvector<T> &rhs) const
475
+ {
476
+ size_type i, ub = std::min(size(), rhs.size());
477
+ for (i = 0; i < ub; ++i)
478
+ if (_first[i] != rhs._first[i])
479
+ return _first[i] > rhs._first[i];
480
+ return size() >= rhs.size();
481
+ }
482
+
483
+ template <class _Cont, class _Iter>
484
+ static void trade(_Cont& dest, _Iter srcBegin, _Iter srcEnd)
485
+ {
486
+ /*static_assert(std::is_same<typename std::remove_reference<std::iterator_traits<_Iter>::value_type>::type, tvector<T>*>::value,
487
+ "value_type of InputIt must be tvector<T>* type");*/
488
+ size_type totalLen = 0;
489
+ for (auto it = srcBegin; it != srcEnd; ++it) totalLen += (*it)->size();
490
+ auto dend = dest.size();
491
+ dest.resize(dest.size() + totalLen);
492
+ T* dp = dest.data() + dend;
493
+ for (auto it = srcBegin; it != srcEnd; ++it)
494
+ {
495
+ auto& tv = **it;
496
+ std::copy(tv.begin(), tv.end(), dp);
497
+ tv = tvector<T>{ dp, tv.size() };
498
+ dp += tv.size();
499
+ }
500
+ }
501
+
502
+ private:
503
+ pointer _first = nullptr; // ptr for first data
504
+ pointer _last = nullptr; // ptr last data
505
+ pointer _rsvEnd = nullptr; // ptr reserved last data
506
+
507
+ inline void buy(size_type needSize)
508
+ {
509
+ if (!isOwner()) throw std::out_of_range("cannot increase size of non-owning mode");
510
+ if (needSize > capacity())
511
+ {
512
+ reallocate(calcGrowth(needSize));
513
+ }
514
+ }
515
+
516
+ inline void reallocate(size_type newSize)
517
+ {
518
+ size_type s = size();
519
+ T *tarr = _Alloc{}.allocate(newSize);
520
+ memcpy(tarr, _first, s * sizeof(T));
521
+ _Alloc{}.deallocate(_first, capacity());
522
+ _first = tarr;
523
+ _last = _first + s;
524
+ _rsvEnd = _first + newSize;
525
+ }
526
+
527
+ inline size_type calcGrowth(size_type newSize)
528
+ {
529
+ const size_type oldcapacity = capacity();
530
+ if (oldcapacity > max_size() - oldcapacity / 2)
531
+ {
532
+ return newSize;
533
+ }
534
+
535
+ const size_type geometric = oldcapacity + oldcapacity / 2;
536
+ if (geometric < newSize)
537
+ {
538
+ return newSize;
539
+ }
540
+ return geometric;
541
+ }
542
+ };
543
+ }