tomoto 0.1.0
Sign up to get free protection for your applications and to get access to all the features.
- checksums.yaml +7 -0
- data/CHANGELOG.md +3 -0
- data/LICENSE.txt +22 -0
- data/README.md +123 -0
- data/ext/tomoto/ext.cpp +245 -0
- data/ext/tomoto/extconf.rb +28 -0
- data/lib/tomoto.rb +12 -0
- data/lib/tomoto/ct.rb +11 -0
- data/lib/tomoto/hdp.rb +11 -0
- data/lib/tomoto/lda.rb +67 -0
- data/lib/tomoto/version.rb +3 -0
- data/vendor/EigenRand/EigenRand/Core.h +1139 -0
- data/vendor/EigenRand/EigenRand/Dists/Basic.h +111 -0
- data/vendor/EigenRand/EigenRand/Dists/Discrete.h +877 -0
- data/vendor/EigenRand/EigenRand/Dists/GammaPoisson.h +108 -0
- data/vendor/EigenRand/EigenRand/Dists/NormalExp.h +626 -0
- data/vendor/EigenRand/EigenRand/EigenRand +19 -0
- data/vendor/EigenRand/EigenRand/Macro.h +24 -0
- data/vendor/EigenRand/EigenRand/MorePacketMath.h +978 -0
- data/vendor/EigenRand/EigenRand/PacketFilter.h +286 -0
- data/vendor/EigenRand/EigenRand/PacketRandomEngine.h +624 -0
- data/vendor/EigenRand/EigenRand/RandUtils.h +413 -0
- data/vendor/EigenRand/EigenRand/doc.h +220 -0
- data/vendor/EigenRand/LICENSE +21 -0
- data/vendor/EigenRand/README.md +288 -0
- data/vendor/eigen/COPYING.BSD +26 -0
- data/vendor/eigen/COPYING.GPL +674 -0
- data/vendor/eigen/COPYING.LGPL +502 -0
- data/vendor/eigen/COPYING.MINPACK +52 -0
- data/vendor/eigen/COPYING.MPL2 +373 -0
- data/vendor/eigen/COPYING.README +18 -0
- data/vendor/eigen/Eigen/CMakeLists.txt +19 -0
- data/vendor/eigen/Eigen/Cholesky +46 -0
- data/vendor/eigen/Eigen/CholmodSupport +48 -0
- data/vendor/eigen/Eigen/Core +537 -0
- data/vendor/eigen/Eigen/Dense +7 -0
- data/vendor/eigen/Eigen/Eigen +2 -0
- data/vendor/eigen/Eigen/Eigenvalues +61 -0
- data/vendor/eigen/Eigen/Geometry +62 -0
- data/vendor/eigen/Eigen/Householder +30 -0
- data/vendor/eigen/Eigen/IterativeLinearSolvers +48 -0
- data/vendor/eigen/Eigen/Jacobi +33 -0
- data/vendor/eigen/Eigen/LU +50 -0
- data/vendor/eigen/Eigen/MetisSupport +35 -0
- data/vendor/eigen/Eigen/OrderingMethods +73 -0
- data/vendor/eigen/Eigen/PaStiXSupport +48 -0
- data/vendor/eigen/Eigen/PardisoSupport +35 -0
- data/vendor/eigen/Eigen/QR +51 -0
- data/vendor/eigen/Eigen/QtAlignedMalloc +40 -0
- data/vendor/eigen/Eigen/SPQRSupport +34 -0
- data/vendor/eigen/Eigen/SVD +51 -0
- data/vendor/eigen/Eigen/Sparse +36 -0
- data/vendor/eigen/Eigen/SparseCholesky +45 -0
- data/vendor/eigen/Eigen/SparseCore +69 -0
- data/vendor/eigen/Eigen/SparseLU +46 -0
- data/vendor/eigen/Eigen/SparseQR +37 -0
- data/vendor/eigen/Eigen/StdDeque +27 -0
- data/vendor/eigen/Eigen/StdList +26 -0
- data/vendor/eigen/Eigen/StdVector +27 -0
- data/vendor/eigen/Eigen/SuperLUSupport +64 -0
- data/vendor/eigen/Eigen/UmfPackSupport +40 -0
- data/vendor/eigen/Eigen/src/Cholesky/LDLT.h +673 -0
- data/vendor/eigen/Eigen/src/Cholesky/LLT.h +542 -0
- data/vendor/eigen/Eigen/src/Cholesky/LLT_LAPACKE.h +99 -0
- data/vendor/eigen/Eigen/src/CholmodSupport/CholmodSupport.h +639 -0
- data/vendor/eigen/Eigen/src/Core/Array.h +329 -0
- data/vendor/eigen/Eigen/src/Core/ArrayBase.h +226 -0
- data/vendor/eigen/Eigen/src/Core/ArrayWrapper.h +209 -0
- data/vendor/eigen/Eigen/src/Core/Assign.h +90 -0
- data/vendor/eigen/Eigen/src/Core/AssignEvaluator.h +935 -0
- data/vendor/eigen/Eigen/src/Core/Assign_MKL.h +178 -0
- data/vendor/eigen/Eigen/src/Core/BandMatrix.h +353 -0
- data/vendor/eigen/Eigen/src/Core/Block.h +452 -0
- data/vendor/eigen/Eigen/src/Core/BooleanRedux.h +164 -0
- data/vendor/eigen/Eigen/src/Core/CommaInitializer.h +160 -0
- data/vendor/eigen/Eigen/src/Core/ConditionEstimator.h +175 -0
- data/vendor/eigen/Eigen/src/Core/CoreEvaluators.h +1688 -0
- data/vendor/eigen/Eigen/src/Core/CoreIterators.h +127 -0
- data/vendor/eigen/Eigen/src/Core/CwiseBinaryOp.h +184 -0
- data/vendor/eigen/Eigen/src/Core/CwiseNullaryOp.h +866 -0
- data/vendor/eigen/Eigen/src/Core/CwiseTernaryOp.h +197 -0
- data/vendor/eigen/Eigen/src/Core/CwiseUnaryOp.h +103 -0
- data/vendor/eigen/Eigen/src/Core/CwiseUnaryView.h +128 -0
- data/vendor/eigen/Eigen/src/Core/DenseBase.h +611 -0
- data/vendor/eigen/Eigen/src/Core/DenseCoeffsBase.h +681 -0
- data/vendor/eigen/Eigen/src/Core/DenseStorage.h +570 -0
- data/vendor/eigen/Eigen/src/Core/Diagonal.h +260 -0
- data/vendor/eigen/Eigen/src/Core/DiagonalMatrix.h +343 -0
- data/vendor/eigen/Eigen/src/Core/DiagonalProduct.h +28 -0
- data/vendor/eigen/Eigen/src/Core/Dot.h +318 -0
- data/vendor/eigen/Eigen/src/Core/EigenBase.h +159 -0
- data/vendor/eigen/Eigen/src/Core/ForceAlignedAccess.h +146 -0
- data/vendor/eigen/Eigen/src/Core/Fuzzy.h +155 -0
- data/vendor/eigen/Eigen/src/Core/GeneralProduct.h +455 -0
- data/vendor/eigen/Eigen/src/Core/GenericPacketMath.h +593 -0
- data/vendor/eigen/Eigen/src/Core/GlobalFunctions.h +187 -0
- data/vendor/eigen/Eigen/src/Core/IO.h +225 -0
- data/vendor/eigen/Eigen/src/Core/Inverse.h +118 -0
- data/vendor/eigen/Eigen/src/Core/Map.h +171 -0
- data/vendor/eigen/Eigen/src/Core/MapBase.h +303 -0
- data/vendor/eigen/Eigen/src/Core/MathFunctions.h +1415 -0
- data/vendor/eigen/Eigen/src/Core/MathFunctionsImpl.h +101 -0
- data/vendor/eigen/Eigen/src/Core/Matrix.h +459 -0
- data/vendor/eigen/Eigen/src/Core/MatrixBase.h +529 -0
- data/vendor/eigen/Eigen/src/Core/NestByValue.h +110 -0
- data/vendor/eigen/Eigen/src/Core/NoAlias.h +108 -0
- data/vendor/eigen/Eigen/src/Core/NumTraits.h +248 -0
- data/vendor/eigen/Eigen/src/Core/PermutationMatrix.h +633 -0
- data/vendor/eigen/Eigen/src/Core/PlainObjectBase.h +1035 -0
- data/vendor/eigen/Eigen/src/Core/Product.h +186 -0
- data/vendor/eigen/Eigen/src/Core/ProductEvaluators.h +1112 -0
- data/vendor/eigen/Eigen/src/Core/Random.h +182 -0
- data/vendor/eigen/Eigen/src/Core/Redux.h +505 -0
- data/vendor/eigen/Eigen/src/Core/Ref.h +283 -0
- data/vendor/eigen/Eigen/src/Core/Replicate.h +142 -0
- data/vendor/eigen/Eigen/src/Core/ReturnByValue.h +117 -0
- data/vendor/eigen/Eigen/src/Core/Reverse.h +211 -0
- data/vendor/eigen/Eigen/src/Core/Select.h +162 -0
- data/vendor/eigen/Eigen/src/Core/SelfAdjointView.h +352 -0
- data/vendor/eigen/Eigen/src/Core/SelfCwiseBinaryOp.h +47 -0
- data/vendor/eigen/Eigen/src/Core/Solve.h +188 -0
- data/vendor/eigen/Eigen/src/Core/SolveTriangular.h +235 -0
- data/vendor/eigen/Eigen/src/Core/SolverBase.h +130 -0
- data/vendor/eigen/Eigen/src/Core/StableNorm.h +221 -0
- data/vendor/eigen/Eigen/src/Core/Stride.h +111 -0
- data/vendor/eigen/Eigen/src/Core/Swap.h +67 -0
- data/vendor/eigen/Eigen/src/Core/Transpose.h +403 -0
- data/vendor/eigen/Eigen/src/Core/Transpositions.h +407 -0
- data/vendor/eigen/Eigen/src/Core/TriangularMatrix.h +983 -0
- data/vendor/eigen/Eigen/src/Core/VectorBlock.h +96 -0
- data/vendor/eigen/Eigen/src/Core/VectorwiseOp.h +695 -0
- data/vendor/eigen/Eigen/src/Core/Visitor.h +273 -0
- data/vendor/eigen/Eigen/src/Core/arch/AVX/Complex.h +451 -0
- data/vendor/eigen/Eigen/src/Core/arch/AVX/MathFunctions.h +439 -0
- data/vendor/eigen/Eigen/src/Core/arch/AVX/PacketMath.h +637 -0
- data/vendor/eigen/Eigen/src/Core/arch/AVX/TypeCasting.h +51 -0
- data/vendor/eigen/Eigen/src/Core/arch/AVX512/MathFunctions.h +391 -0
- data/vendor/eigen/Eigen/src/Core/arch/AVX512/PacketMath.h +1316 -0
- data/vendor/eigen/Eigen/src/Core/arch/AltiVec/Complex.h +430 -0
- data/vendor/eigen/Eigen/src/Core/arch/AltiVec/MathFunctions.h +322 -0
- data/vendor/eigen/Eigen/src/Core/arch/AltiVec/PacketMath.h +1061 -0
- data/vendor/eigen/Eigen/src/Core/arch/CUDA/Complex.h +103 -0
- data/vendor/eigen/Eigen/src/Core/arch/CUDA/Half.h +674 -0
- data/vendor/eigen/Eigen/src/Core/arch/CUDA/MathFunctions.h +91 -0
- data/vendor/eigen/Eigen/src/Core/arch/CUDA/PacketMath.h +333 -0
- data/vendor/eigen/Eigen/src/Core/arch/CUDA/PacketMathHalf.h +1124 -0
- data/vendor/eigen/Eigen/src/Core/arch/CUDA/TypeCasting.h +212 -0
- data/vendor/eigen/Eigen/src/Core/arch/Default/ConjHelper.h +29 -0
- data/vendor/eigen/Eigen/src/Core/arch/Default/Settings.h +49 -0
- data/vendor/eigen/Eigen/src/Core/arch/NEON/Complex.h +490 -0
- data/vendor/eigen/Eigen/src/Core/arch/NEON/MathFunctions.h +91 -0
- data/vendor/eigen/Eigen/src/Core/arch/NEON/PacketMath.h +760 -0
- data/vendor/eigen/Eigen/src/Core/arch/SSE/Complex.h +471 -0
- data/vendor/eigen/Eigen/src/Core/arch/SSE/MathFunctions.h +562 -0
- data/vendor/eigen/Eigen/src/Core/arch/SSE/PacketMath.h +895 -0
- data/vendor/eigen/Eigen/src/Core/arch/SSE/TypeCasting.h +77 -0
- data/vendor/eigen/Eigen/src/Core/arch/ZVector/Complex.h +397 -0
- data/vendor/eigen/Eigen/src/Core/arch/ZVector/MathFunctions.h +137 -0
- data/vendor/eigen/Eigen/src/Core/arch/ZVector/PacketMath.h +945 -0
- data/vendor/eigen/Eigen/src/Core/functors/AssignmentFunctors.h +168 -0
- data/vendor/eigen/Eigen/src/Core/functors/BinaryFunctors.h +475 -0
- data/vendor/eigen/Eigen/src/Core/functors/NullaryFunctors.h +188 -0
- data/vendor/eigen/Eigen/src/Core/functors/StlFunctors.h +136 -0
- data/vendor/eigen/Eigen/src/Core/functors/TernaryFunctors.h +25 -0
- data/vendor/eigen/Eigen/src/Core/functors/UnaryFunctors.h +792 -0
- data/vendor/eigen/Eigen/src/Core/products/GeneralBlockPanelKernel.h +2156 -0
- data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrix.h +492 -0
- data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrixTriangular.h +311 -0
- data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h +145 -0
- data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixMatrix_BLAS.h +122 -0
- data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixVector.h +619 -0
- data/vendor/eigen/Eigen/src/Core/products/GeneralMatrixVector_BLAS.h +136 -0
- data/vendor/eigen/Eigen/src/Core/products/Parallelizer.h +163 -0
- data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixMatrix.h +521 -0
- data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixMatrix_BLAS.h +287 -0
- data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixVector.h +260 -0
- data/vendor/eigen/Eigen/src/Core/products/SelfadjointMatrixVector_BLAS.h +118 -0
- data/vendor/eigen/Eigen/src/Core/products/SelfadjointProduct.h +133 -0
- data/vendor/eigen/Eigen/src/Core/products/SelfadjointRank2Update.h +93 -0
- data/vendor/eigen/Eigen/src/Core/products/TriangularMatrixMatrix.h +466 -0
- data/vendor/eigen/Eigen/src/Core/products/TriangularMatrixMatrix_BLAS.h +315 -0
- data/vendor/eigen/Eigen/src/Core/products/TriangularMatrixVector.h +350 -0
- data/vendor/eigen/Eigen/src/Core/products/TriangularMatrixVector_BLAS.h +255 -0
- data/vendor/eigen/Eigen/src/Core/products/TriangularSolverMatrix.h +335 -0
- data/vendor/eigen/Eigen/src/Core/products/TriangularSolverMatrix_BLAS.h +163 -0
- data/vendor/eigen/Eigen/src/Core/products/TriangularSolverVector.h +145 -0
- data/vendor/eigen/Eigen/src/Core/util/BlasUtil.h +398 -0
- data/vendor/eigen/Eigen/src/Core/util/Constants.h +547 -0
- data/vendor/eigen/Eigen/src/Core/util/DisableStupidWarnings.h +83 -0
- data/vendor/eigen/Eigen/src/Core/util/ForwardDeclarations.h +302 -0
- data/vendor/eigen/Eigen/src/Core/util/MKL_support.h +130 -0
- data/vendor/eigen/Eigen/src/Core/util/Macros.h +1001 -0
- data/vendor/eigen/Eigen/src/Core/util/Memory.h +993 -0
- data/vendor/eigen/Eigen/src/Core/util/Meta.h +534 -0
- data/vendor/eigen/Eigen/src/Core/util/NonMPL2.h +3 -0
- data/vendor/eigen/Eigen/src/Core/util/ReenableStupidWarnings.h +27 -0
- data/vendor/eigen/Eigen/src/Core/util/StaticAssert.h +218 -0
- data/vendor/eigen/Eigen/src/Core/util/XprHelper.h +821 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/ComplexEigenSolver.h +346 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/ComplexSchur.h +459 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/ComplexSchur_LAPACKE.h +91 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/EigenSolver.h +622 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/GeneralizedEigenSolver.h +418 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/GeneralizedSelfAdjointEigenSolver.h +226 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/HessenbergDecomposition.h +374 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/MatrixBaseEigenvalues.h +158 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/RealQZ.h +654 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/RealSchur.h +546 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/RealSchur_LAPACKE.h +77 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h +870 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/SelfAdjointEigenSolver_LAPACKE.h +87 -0
- data/vendor/eigen/Eigen/src/Eigenvalues/Tridiagonalization.h +556 -0
- data/vendor/eigen/Eigen/src/Geometry/AlignedBox.h +392 -0
- data/vendor/eigen/Eigen/src/Geometry/AngleAxis.h +247 -0
- data/vendor/eigen/Eigen/src/Geometry/EulerAngles.h +114 -0
- data/vendor/eigen/Eigen/src/Geometry/Homogeneous.h +497 -0
- data/vendor/eigen/Eigen/src/Geometry/Hyperplane.h +282 -0
- data/vendor/eigen/Eigen/src/Geometry/OrthoMethods.h +234 -0
- data/vendor/eigen/Eigen/src/Geometry/ParametrizedLine.h +195 -0
- data/vendor/eigen/Eigen/src/Geometry/Quaternion.h +814 -0
- data/vendor/eigen/Eigen/src/Geometry/Rotation2D.h +199 -0
- data/vendor/eigen/Eigen/src/Geometry/RotationBase.h +206 -0
- data/vendor/eigen/Eigen/src/Geometry/Scaling.h +170 -0
- data/vendor/eigen/Eigen/src/Geometry/Transform.h +1542 -0
- data/vendor/eigen/Eigen/src/Geometry/Translation.h +208 -0
- data/vendor/eigen/Eigen/src/Geometry/Umeyama.h +166 -0
- data/vendor/eigen/Eigen/src/Geometry/arch/Geometry_SSE.h +161 -0
- data/vendor/eigen/Eigen/src/Householder/BlockHouseholder.h +103 -0
- data/vendor/eigen/Eigen/src/Householder/Householder.h +172 -0
- data/vendor/eigen/Eigen/src/Householder/HouseholderSequence.h +470 -0
- data/vendor/eigen/Eigen/src/IterativeLinearSolvers/BasicPreconditioners.h +226 -0
- data/vendor/eigen/Eigen/src/IterativeLinearSolvers/BiCGSTAB.h +228 -0
- data/vendor/eigen/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h +246 -0
- data/vendor/eigen/Eigen/src/IterativeLinearSolvers/IncompleteCholesky.h +400 -0
- data/vendor/eigen/Eigen/src/IterativeLinearSolvers/IncompleteLUT.h +462 -0
- data/vendor/eigen/Eigen/src/IterativeLinearSolvers/IterativeSolverBase.h +394 -0
- data/vendor/eigen/Eigen/src/IterativeLinearSolvers/LeastSquareConjugateGradient.h +216 -0
- data/vendor/eigen/Eigen/src/IterativeLinearSolvers/SolveWithGuess.h +115 -0
- data/vendor/eigen/Eigen/src/Jacobi/Jacobi.h +462 -0
- data/vendor/eigen/Eigen/src/LU/Determinant.h +101 -0
- data/vendor/eigen/Eigen/src/LU/FullPivLU.h +891 -0
- data/vendor/eigen/Eigen/src/LU/InverseImpl.h +415 -0
- data/vendor/eigen/Eigen/src/LU/PartialPivLU.h +611 -0
- data/vendor/eigen/Eigen/src/LU/PartialPivLU_LAPACKE.h +83 -0
- data/vendor/eigen/Eigen/src/LU/arch/Inverse_SSE.h +338 -0
- data/vendor/eigen/Eigen/src/MetisSupport/MetisSupport.h +137 -0
- data/vendor/eigen/Eigen/src/OrderingMethods/Amd.h +445 -0
- data/vendor/eigen/Eigen/src/OrderingMethods/Eigen_Colamd.h +1843 -0
- data/vendor/eigen/Eigen/src/OrderingMethods/Ordering.h +157 -0
- data/vendor/eigen/Eigen/src/PaStiXSupport/PaStiXSupport.h +678 -0
- data/vendor/eigen/Eigen/src/PardisoSupport/PardisoSupport.h +543 -0
- data/vendor/eigen/Eigen/src/QR/ColPivHouseholderQR.h +653 -0
- data/vendor/eigen/Eigen/src/QR/ColPivHouseholderQR_LAPACKE.h +97 -0
- data/vendor/eigen/Eigen/src/QR/CompleteOrthogonalDecomposition.h +562 -0
- data/vendor/eigen/Eigen/src/QR/FullPivHouseholderQR.h +676 -0
- data/vendor/eigen/Eigen/src/QR/HouseholderQR.h +409 -0
- data/vendor/eigen/Eigen/src/QR/HouseholderQR_LAPACKE.h +68 -0
- data/vendor/eigen/Eigen/src/SPQRSupport/SuiteSparseQRSupport.h +313 -0
- data/vendor/eigen/Eigen/src/SVD/BDCSVD.h +1246 -0
- data/vendor/eigen/Eigen/src/SVD/JacobiSVD.h +804 -0
- data/vendor/eigen/Eigen/src/SVD/JacobiSVD_LAPACKE.h +91 -0
- data/vendor/eigen/Eigen/src/SVD/SVDBase.h +315 -0
- data/vendor/eigen/Eigen/src/SVD/UpperBidiagonalization.h +414 -0
- data/vendor/eigen/Eigen/src/SparseCholesky/SimplicialCholesky.h +689 -0
- data/vendor/eigen/Eigen/src/SparseCholesky/SimplicialCholesky_impl.h +199 -0
- data/vendor/eigen/Eigen/src/SparseCore/AmbiVector.h +377 -0
- data/vendor/eigen/Eigen/src/SparseCore/CompressedStorage.h +258 -0
- data/vendor/eigen/Eigen/src/SparseCore/ConservativeSparseSparseProduct.h +352 -0
- data/vendor/eigen/Eigen/src/SparseCore/MappedSparseMatrix.h +67 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseAssign.h +216 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseBlock.h +603 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseColEtree.h +206 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseCompressedBase.h +341 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseCwiseBinaryOp.h +726 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseCwiseUnaryOp.h +148 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseDenseProduct.h +320 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseDiagonalProduct.h +138 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseDot.h +98 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseFuzzy.h +29 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseMap.h +305 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseMatrix.h +1403 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseMatrixBase.h +405 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparsePermutation.h +178 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseProduct.h +169 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseRedux.h +49 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseRef.h +397 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseSelfAdjointView.h +656 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseSolverBase.h +124 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseSparseProductWithPruning.h +198 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseTranspose.h +92 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseTriangularView.h +189 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseUtil.h +178 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseVector.h +478 -0
- data/vendor/eigen/Eigen/src/SparseCore/SparseView.h +253 -0
- data/vendor/eigen/Eigen/src/SparseCore/TriangularSolver.h +315 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU.h +773 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLUImpl.h +66 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_Memory.h +226 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_Structs.h +110 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_SupernodalMatrix.h +301 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_Utils.h +80 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_column_bmod.h +181 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_column_dfs.h +179 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_copy_to_ucol.h +107 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_gemm_kernel.h +280 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_heap_relax_snode.h +126 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_kernel_bmod.h +130 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_panel_bmod.h +223 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_panel_dfs.h +258 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_pivotL.h +137 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_pruneL.h +136 -0
- data/vendor/eigen/Eigen/src/SparseLU/SparseLU_relax_snode.h +83 -0
- data/vendor/eigen/Eigen/src/SparseQR/SparseQR.h +745 -0
- data/vendor/eigen/Eigen/src/StlSupport/StdDeque.h +126 -0
- data/vendor/eigen/Eigen/src/StlSupport/StdList.h +106 -0
- data/vendor/eigen/Eigen/src/StlSupport/StdVector.h +131 -0
- data/vendor/eigen/Eigen/src/StlSupport/details.h +84 -0
- data/vendor/eigen/Eigen/src/SuperLUSupport/SuperLUSupport.h +1027 -0
- data/vendor/eigen/Eigen/src/UmfPackSupport/UmfPackSupport.h +506 -0
- data/vendor/eigen/Eigen/src/misc/Image.h +82 -0
- data/vendor/eigen/Eigen/src/misc/Kernel.h +79 -0
- data/vendor/eigen/Eigen/src/misc/RealSvd2x2.h +55 -0
- data/vendor/eigen/Eigen/src/misc/blas.h +440 -0
- data/vendor/eigen/Eigen/src/misc/lapack.h +152 -0
- data/vendor/eigen/Eigen/src/misc/lapacke.h +16291 -0
- data/vendor/eigen/Eigen/src/misc/lapacke_mangling.h +17 -0
- data/vendor/eigen/Eigen/src/plugins/ArrayCwiseBinaryOps.h +332 -0
- data/vendor/eigen/Eigen/src/plugins/ArrayCwiseUnaryOps.h +552 -0
- data/vendor/eigen/Eigen/src/plugins/BlockMethods.h +1058 -0
- data/vendor/eigen/Eigen/src/plugins/CommonCwiseBinaryOps.h +115 -0
- data/vendor/eigen/Eigen/src/plugins/CommonCwiseUnaryOps.h +163 -0
- data/vendor/eigen/Eigen/src/plugins/MatrixCwiseBinaryOps.h +152 -0
- data/vendor/eigen/Eigen/src/plugins/MatrixCwiseUnaryOps.h +85 -0
- data/vendor/eigen/README.md +3 -0
- data/vendor/eigen/bench/README.txt +55 -0
- data/vendor/eigen/bench/btl/COPYING +340 -0
- data/vendor/eigen/bench/btl/README +154 -0
- data/vendor/eigen/bench/tensors/README +21 -0
- data/vendor/eigen/blas/README.txt +6 -0
- data/vendor/eigen/demos/mandelbrot/README +10 -0
- data/vendor/eigen/demos/mix_eigen_and_c/README +9 -0
- data/vendor/eigen/demos/opengl/README +13 -0
- data/vendor/eigen/unsupported/Eigen/CXX11/src/Tensor/README.md +1760 -0
- data/vendor/eigen/unsupported/README.txt +50 -0
- data/vendor/tomotopy/LICENSE +21 -0
- data/vendor/tomotopy/README.kr.rst +375 -0
- data/vendor/tomotopy/README.rst +382 -0
- data/vendor/tomotopy/src/Labeling/FoRelevance.cpp +362 -0
- data/vendor/tomotopy/src/Labeling/FoRelevance.h +88 -0
- data/vendor/tomotopy/src/Labeling/Labeler.h +50 -0
- data/vendor/tomotopy/src/TopicModel/CT.h +37 -0
- data/vendor/tomotopy/src/TopicModel/CTModel.cpp +13 -0
- data/vendor/tomotopy/src/TopicModel/CTModel.hpp +293 -0
- data/vendor/tomotopy/src/TopicModel/DMR.h +51 -0
- data/vendor/tomotopy/src/TopicModel/DMRModel.cpp +13 -0
- data/vendor/tomotopy/src/TopicModel/DMRModel.hpp +374 -0
- data/vendor/tomotopy/src/TopicModel/DT.h +65 -0
- data/vendor/tomotopy/src/TopicModel/DTM.h +22 -0
- data/vendor/tomotopy/src/TopicModel/DTModel.cpp +15 -0
- data/vendor/tomotopy/src/TopicModel/DTModel.hpp +572 -0
- data/vendor/tomotopy/src/TopicModel/GDMR.h +37 -0
- data/vendor/tomotopy/src/TopicModel/GDMRModel.cpp +14 -0
- data/vendor/tomotopy/src/TopicModel/GDMRModel.hpp +485 -0
- data/vendor/tomotopy/src/TopicModel/HDP.h +74 -0
- data/vendor/tomotopy/src/TopicModel/HDPModel.cpp +13 -0
- data/vendor/tomotopy/src/TopicModel/HDPModel.hpp +592 -0
- data/vendor/tomotopy/src/TopicModel/HLDA.h +40 -0
- data/vendor/tomotopy/src/TopicModel/HLDAModel.cpp +13 -0
- data/vendor/tomotopy/src/TopicModel/HLDAModel.hpp +681 -0
- data/vendor/tomotopy/src/TopicModel/HPA.h +27 -0
- data/vendor/tomotopy/src/TopicModel/HPAModel.cpp +21 -0
- data/vendor/tomotopy/src/TopicModel/HPAModel.hpp +588 -0
- data/vendor/tomotopy/src/TopicModel/LDA.h +144 -0
- data/vendor/tomotopy/src/TopicModel/LDACVB0Model.hpp +442 -0
- data/vendor/tomotopy/src/TopicModel/LDAModel.cpp +13 -0
- data/vendor/tomotopy/src/TopicModel/LDAModel.hpp +1058 -0
- data/vendor/tomotopy/src/TopicModel/LLDA.h +45 -0
- data/vendor/tomotopy/src/TopicModel/LLDAModel.cpp +13 -0
- data/vendor/tomotopy/src/TopicModel/LLDAModel.hpp +203 -0
- data/vendor/tomotopy/src/TopicModel/MGLDA.h +63 -0
- data/vendor/tomotopy/src/TopicModel/MGLDAModel.cpp +17 -0
- data/vendor/tomotopy/src/TopicModel/MGLDAModel.hpp +558 -0
- data/vendor/tomotopy/src/TopicModel/PA.h +43 -0
- data/vendor/tomotopy/src/TopicModel/PAModel.cpp +13 -0
- data/vendor/tomotopy/src/TopicModel/PAModel.hpp +467 -0
- data/vendor/tomotopy/src/TopicModel/PLDA.h +17 -0
- data/vendor/tomotopy/src/TopicModel/PLDAModel.cpp +13 -0
- data/vendor/tomotopy/src/TopicModel/PLDAModel.hpp +214 -0
- data/vendor/tomotopy/src/TopicModel/SLDA.h +54 -0
- data/vendor/tomotopy/src/TopicModel/SLDAModel.cpp +17 -0
- data/vendor/tomotopy/src/TopicModel/SLDAModel.hpp +456 -0
- data/vendor/tomotopy/src/TopicModel/TopicModel.hpp +692 -0
- data/vendor/tomotopy/src/Utils/AliasMethod.hpp +169 -0
- data/vendor/tomotopy/src/Utils/Dictionary.h +80 -0
- data/vendor/tomotopy/src/Utils/EigenAddonOps.hpp +181 -0
- data/vendor/tomotopy/src/Utils/LBFGS.h +202 -0
- data/vendor/tomotopy/src/Utils/LBFGS/LineSearchBacktracking.h +120 -0
- data/vendor/tomotopy/src/Utils/LBFGS/LineSearchBracketing.h +122 -0
- data/vendor/tomotopy/src/Utils/LBFGS/Param.h +213 -0
- data/vendor/tomotopy/src/Utils/LUT.hpp +82 -0
- data/vendor/tomotopy/src/Utils/MultiNormalDistribution.hpp +69 -0
- data/vendor/tomotopy/src/Utils/PolyaGamma.hpp +200 -0
- data/vendor/tomotopy/src/Utils/PolyaGammaHybrid.hpp +672 -0
- data/vendor/tomotopy/src/Utils/ThreadPool.hpp +150 -0
- data/vendor/tomotopy/src/Utils/Trie.hpp +220 -0
- data/vendor/tomotopy/src/Utils/TruncMultiNormal.hpp +94 -0
- data/vendor/tomotopy/src/Utils/Utils.hpp +337 -0
- data/vendor/tomotopy/src/Utils/avx_gamma.h +46 -0
- data/vendor/tomotopy/src/Utils/avx_mathfun.h +736 -0
- data/vendor/tomotopy/src/Utils/exception.h +28 -0
- data/vendor/tomotopy/src/Utils/math.h +281 -0
- data/vendor/tomotopy/src/Utils/rtnorm.hpp +2690 -0
- data/vendor/tomotopy/src/Utils/sample.hpp +192 -0
- data/vendor/tomotopy/src/Utils/serializer.hpp +695 -0
- data/vendor/tomotopy/src/Utils/slp.hpp +131 -0
- data/vendor/tomotopy/src/Utils/sse_gamma.h +48 -0
- data/vendor/tomotopy/src/Utils/sse_mathfun.h +710 -0
- data/vendor/tomotopy/src/Utils/text.hpp +49 -0
- data/vendor/tomotopy/src/Utils/tvector.hpp +543 -0
- 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
|
+
}
|