├── .gitignore ├── COPYING ├── README.mkd ├── eigen-patched ├── CMakeLists.txt ├── COPYING.BSD ├── COPYING.GPL ├── COPYING.LGPL ├── COPYING.MINPACK ├── COPYING.MPL2 ├── COPYING.README ├── CTestConfig.cmake ├── CTestCustom.cmake.in ├── Eigen │ ├── Array │ ├── CMakeLists.txt │ ├── Cholesky │ ├── CholmodSupport │ ├── Core │ ├── Dense │ ├── Eigen │ ├── Eigen2Support │ ├── Eigenvalues │ ├── Geometry │ ├── Householder │ ├── IterativeLinearSolvers │ ├── Jacobi │ ├── LU │ ├── LeastSquares │ ├── OrderingMethods │ ├── PaStiXSupport │ ├── PardisoSupport │ ├── QR │ ├── QtAlignedMalloc │ ├── SVD │ ├── Sparse │ ├── SparseCholesky │ ├── SparseCore │ ├── StdDeque │ ├── StdList │ ├── StdVector │ ├── SuperLUSupport │ ├── UmfPackSupport │ └── src │ │ ├── CMakeLists.txt │ │ ├── Cholesky │ │ ├── CMakeLists.txt │ │ ├── LDLT.h │ │ ├── LLT.h │ │ └── LLT_MKL.h │ │ ├── CholmodSupport │ │ ├── CMakeLists.txt │ │ └── CholmodSupport.h │ │ ├── Core │ │ ├── Array.h │ │ ├── ArrayBase.h │ │ ├── ArrayWrapper.h │ │ ├── Assign.h │ │ ├── Assign_MKL.h │ │ ├── BandMatrix.h │ │ ├── Block.h │ │ ├── BooleanRedux.h │ │ ├── CMakeLists.txt │ │ ├── CommaInitializer.h │ │ ├── CwiseBinaryOp.h │ │ ├── CwiseNullaryOp.h │ │ ├── CwiseUnaryOp.h │ │ ├── CwiseUnaryView.h │ │ ├── DenseBase.h │ │ ├── DenseCoeffsBase.h │ │ ├── DenseStorage.h │ │ ├── Diagonal.h │ │ ├── DiagonalMatrix.h │ │ ├── DiagonalProduct.h │ │ ├── Dot.h │ │ ├── EigenBase.h │ │ ├── Flagged.h │ │ ├── ForceAlignedAccess.h │ │ ├── Functors.h │ │ ├── Fuzzy.h │ │ ├── GeneralProduct.h │ │ ├── GenericPacketMath.h │ │ ├── GlobalFunctions.h │ │ ├── IO.h │ │ ├── Map.h │ │ ├── MapBase.h │ │ ├── MathFunctions.h │ │ ├── Matrix.h │ │ ├── MatrixBase.h │ │ ├── NestByValue.h │ │ ├── NoAlias.h │ │ ├── NumTraits.h │ │ ├── PermutationMatrix.h │ │ ├── PlainObjectBase.h │ │ ├── Product.h │ │ ├── ProductBase.h │ │ ├── Random.h │ │ ├── Redux.h │ │ ├── Replicate.h │ │ ├── ReturnByValue.h │ │ ├── Reverse.h │ │ ├── Select.h │ │ ├── SelfAdjointView.h │ │ ├── SelfCwiseBinaryOp.h │ │ ├── SolveTriangular.h │ │ ├── StableNorm.h │ │ ├── Stride.h │ │ ├── Swap.h │ │ ├── Transpose.h │ │ ├── Transpositions.h │ │ ├── TriangularMatrix.h │ │ ├── VectorBlock.h │ │ ├── VectorwiseOp.h │ │ ├── Visitor.h │ │ ├── arch │ │ │ ├── AltiVec │ │ │ │ ├── CMakeLists.txt │ │ │ │ ├── Complex.h │ │ │ │ └── PacketMath.h │ │ │ ├── CMakeLists.txt │ │ │ ├── Default │ │ │ │ ├── CMakeLists.txt │ │ │ │ └── Settings.h │ │ │ ├── NEON │ │ │ │ ├── CMakeLists.txt │ │ │ │ ├── Complex.h │ │ │ │ └── PacketMath.h │ │ │ └── SSE │ │ │ │ ├── CMakeLists.txt │ │ │ │ ├── Complex.h │ │ │ │ ├── MathFunctions.h │ │ │ │ └── PacketMath.h │ │ ├── products │ │ │ ├── CMakeLists.txt │ │ │ ├── CoeffBasedProduct.h │ │ │ ├── GeneralBlockPanelKernel.h │ │ │ ├── GeneralMatrixMatrix.h │ │ │ ├── GeneralMatrixMatrixTriangular.h │ │ │ ├── GeneralMatrixMatrixTriangular_MKL.h │ │ │ ├── GeneralMatrixMatrix_MKL.h │ │ │ ├── GeneralMatrixVector.h │ │ │ ├── GeneralMatrixVector_MKL.h │ │ │ ├── Parallelizer.h │ │ │ ├── SelfadjointMatrixMatrix.h │ │ │ ├── SelfadjointMatrixMatrix_MKL.h │ │ │ ├── SelfadjointMatrixVector.h │ │ │ ├── SelfadjointMatrixVector_MKL.h │ │ │ ├── SelfadjointProduct.h │ │ │ ├── SelfadjointRank2Update.h │ │ │ ├── TriangularMatrixMatrix.h │ │ │ ├── TriangularMatrixMatrix_MKL.h │ │ │ ├── TriangularMatrixVector.h │ │ │ ├── TriangularMatrixVector_MKL.h │ │ │ ├── TriangularSolverMatrix.h │ │ │ ├── TriangularSolverMatrix_MKL.h │ │ │ └── TriangularSolverVector.h │ │ └── util │ │ │ ├── BlasUtil.h │ │ │ ├── CMakeLists.txt │ │ │ ├── Constants.h │ │ │ ├── DisableStupidWarnings.h │ │ │ ├── ForwardDeclarations.h │ │ │ ├── MKL_support.h │ │ │ ├── Macros.h │ │ │ ├── Memory.h │ │ │ ├── Meta.h │ │ │ ├── NonMPL2.h │ │ │ ├── ReenableStupidWarnings.h │ │ │ ├── StaticAssert.h │ │ │ └── XprHelper.h │ │ ├── Eigen2Support │ │ ├── Block.h │ │ ├── CMakeLists.txt │ │ ├── Cwise.h │ │ ├── CwiseOperators.h │ │ ├── Geometry │ │ │ ├── AlignedBox.h │ │ │ ├── All.h │ │ │ ├── AngleAxis.h │ │ │ ├── CMakeLists.txt │ │ │ ├── Hyperplane.h │ │ │ ├── ParametrizedLine.h │ │ │ ├── Quaternion.h │ │ │ ├── Rotation2D.h │ │ │ ├── RotationBase.h │ │ │ ├── Scaling.h │ │ │ ├── Transform.h │ │ │ └── Translation.h │ │ ├── LU.h │ │ ├── Lazy.h │ │ ├── LeastSquares.h │ │ ├── Macros.h │ │ ├── MathFunctions.h │ │ ├── Memory.h │ │ ├── Meta.h │ │ ├── Minor.h │ │ ├── QR.h │ │ ├── SVD.h │ │ ├── TriangularSolver.h │ │ └── VectorBlock.h │ │ ├── Eigenvalues │ │ ├── CMakeLists.txt │ │ ├── ComplexEigenSolver.h │ │ ├── ComplexSchur.h │ │ ├── ComplexSchur_MKL.h │ │ ├── EigenSolver.h │ │ ├── GeneralizedSelfAdjointEigenSolver.h │ │ ├── HessenbergDecomposition.h │ │ ├── MatrixBaseEigenvalues.h │ │ ├── RealSchur.h │ │ ├── RealSchur_MKL.h │ │ ├── SelfAdjointEigenSolver.h │ │ ├── SelfAdjointEigenSolver_MKL.h │ │ └── Tridiagonalization.h │ │ ├── Geometry │ │ ├── AlignedBox.h │ │ ├── AngleAxis.h │ │ ├── CMakeLists.txt │ │ ├── EulerAngles.h │ │ ├── Homogeneous.h │ │ ├── Hyperplane.h │ │ ├── OrthoMethods.h │ │ ├── ParametrizedLine.h │ │ ├── Quaternion.h │ │ ├── Rotation2D.h │ │ ├── RotationBase.h │ │ ├── Scaling.h │ │ ├── Transform.h │ │ ├── Translation.h │ │ ├── Umeyama.h │ │ └── arch │ │ │ ├── CMakeLists.txt │ │ │ └── Geometry_SSE.h │ │ ├── Householder │ │ ├── BlockHouseholder.h │ │ ├── CMakeLists.txt │ │ ├── Householder.h │ │ └── HouseholderSequence.h │ │ ├── IterativeLinearSolvers │ │ ├── BasicPreconditioners.h │ │ ├── BiCGSTAB.h │ │ ├── CMakeLists.txt │ │ ├── ConjugateGradient.h │ │ ├── IncompleteLUT.h │ │ └── IterativeSolverBase.h │ │ ├── Jacobi │ │ ├── CMakeLists.txt │ │ └── Jacobi.h │ │ ├── LU │ │ ├── CMakeLists.txt │ │ ├── Determinant.h │ │ ├── FullPivLU.h │ │ ├── Inverse.h │ │ ├── PartialPivLU.h │ │ ├── PartialPivLU_MKL.h │ │ └── arch │ │ │ ├── CMakeLists.txt │ │ │ └── Inverse_SSE.h │ │ ├── OrderingMethods │ │ ├── Amd.h │ │ └── CMakeLists.txt │ │ ├── PaStiXSupport │ │ ├── CMakeLists.txt │ │ └── PaStiXSupport.h │ │ ├── PardisoSupport │ │ ├── CMakeLists.txt │ │ └── PardisoSupport.h │ │ ├── QR │ │ ├── CMakeLists.txt │ │ ├── ColPivHouseholderQR.h │ │ ├── ColPivHouseholderQR_MKL.h │ │ ├── FullPivHouseholderQR.h │ │ ├── HouseholderQR.h │ │ └── HouseholderQR_MKL.h │ │ ├── SVD │ │ ├── CMakeLists.txt │ │ ├── JacobiSVD.h │ │ ├── JacobiSVD_MKL.h │ │ └── UpperBidiagonalization.h │ │ ├── SparseCholesky │ │ ├── CMakeLists.txt │ │ └── SimplicialCholesky.h │ │ ├── SparseCore │ │ ├── AmbiVector.h │ │ ├── CMakeLists.txt │ │ ├── CompressedStorage.h │ │ ├── ConservativeSparseSparseProduct.h │ │ ├── CoreIterators.h │ │ ├── MappedSparseMatrix.h │ │ ├── SparseAssign.h │ │ ├── SparseBlock.h │ │ ├── SparseCwiseBinaryOp.h │ │ ├── SparseCwiseUnaryOp.h │ │ ├── SparseDenseProduct.h │ │ ├── SparseDiagonalProduct.h │ │ ├── SparseDot.h │ │ ├── SparseFuzzy.h │ │ ├── SparseMatrix.h │ │ ├── SparseMatrixBase.h │ │ ├── SparsePermutation.h │ │ ├── SparseProduct.h │ │ ├── SparseRedux.h │ │ ├── SparseSelfAdjointView.h │ │ ├── SparseSparseProductWithPruning.h │ │ ├── SparseTranspose.h │ │ ├── SparseTriangularView.h │ │ ├── SparseUtil.h │ │ ├── SparseVector.h │ │ ├── SparseView.h │ │ └── TriangularSolver.h │ │ ├── StlSupport │ │ ├── CMakeLists.txt │ │ ├── StdDeque.h │ │ ├── StdList.h │ │ ├── StdVector.h │ │ └── details.h │ │ ├── SuperLUSupport │ │ ├── CMakeLists.txt │ │ └── SuperLUSupport.h │ │ ├── UmfPackSupport │ │ ├── CMakeLists.txt │ │ └── UmfPackSupport.h │ │ ├── misc │ │ ├── CMakeLists.txt │ │ ├── Image.h │ │ ├── Kernel.h │ │ ├── Solve.h │ │ ├── SparseSolve.h │ │ └── blas.h │ │ └── plugins │ │ ├── ArrayCwiseBinaryOps.h │ │ ├── ArrayCwiseUnaryOps.h │ │ ├── BlockMethods.h │ │ ├── CMakeLists.txt │ │ ├── CommonCwiseBinaryOps.h │ │ ├── CommonCwiseUnaryOps.h │ │ ├── MatrixCwiseBinaryOps.h │ │ └── MatrixCwiseUnaryOps.h ├── INSTALL ├── bench │ ├── BenchSparseUtil.h │ ├── BenchTimer.h │ ├── BenchUtil.h │ ├── README.txt │ ├── basicbench.cxxlist │ ├── basicbenchmark.cpp │ ├── basicbenchmark.h │ ├── benchBlasGemm.cpp │ ├── benchCholesky.cpp │ ├── benchEigenSolver.cpp │ ├── benchFFT.cpp │ ├── benchVecAdd.cpp │ ├── bench_gemm.cpp │ ├── bench_multi_compilers.sh │ ├── bench_norm.cpp │ ├── bench_reverse.cpp │ ├── bench_sum.cpp │ ├── bench_unrolling │ ├── benchmark.cpp │ ├── benchmarkSlice.cpp │ ├── benchmarkX.cpp │ ├── benchmarkXcwise.cpp │ ├── benchmark_suite │ ├── btl │ │ ├── CMakeLists.txt │ │ ├── COPYING │ │ ├── README │ │ ├── actions │ │ │ ├── action_aat_product.hh │ │ │ ├── action_ata_product.hh │ │ │ ├── action_atv_product.hh │ │ │ ├── action_axpby.hh │ │ │ ├── action_axpy.hh │ │ │ ├── action_cholesky.hh │ │ │ ├── action_ger.hh │ │ │ ├── action_hessenberg.hh │ │ │ ├── action_lu_decomp.hh │ │ │ ├── action_lu_solve.hh │ │ │ ├── action_matrix_matrix_product.hh │ │ │ ├── action_matrix_matrix_product_bis.hh │ │ │ ├── action_matrix_vector_product.hh │ │ │ ├── action_partial_lu.hh │ │ │ ├── action_rot.hh │ │ │ ├── action_symv.hh │ │ │ ├── action_syr2.hh │ │ │ ├── action_trisolve.hh │ │ │ ├── action_trisolve_matrix.hh │ │ │ ├── action_trmm.hh │ │ │ └── basic_actions.hh │ │ ├── cmake │ │ │ ├── FindACML.cmake │ │ │ ├── FindATLAS.cmake │ │ │ ├── FindBlitz.cmake │ │ │ ├── FindCBLAS.cmake │ │ │ ├── FindGMM.cmake │ │ │ ├── FindGOTO.cmake │ │ │ ├── FindGOTO2.cmake │ │ │ ├── FindMKL.cmake │ │ │ ├── FindMTL4.cmake │ │ │ ├── FindPackageHandleStandardArgs.cmake │ │ │ ├── FindTvmet.cmake │ │ │ └── MacroOptionalAddSubdirectory.cmake │ │ ├── data │ │ │ ├── CMakeLists.txt │ │ │ ├── action_settings.txt │ │ │ ├── gnuplot_common_settings.hh │ │ │ ├── go_mean │ │ │ ├── mean.cxx │ │ │ ├── mk_gnuplot_script.sh │ │ │ ├── mk_mean_script.sh │ │ │ ├── mk_new_gnuplot.sh │ │ │ ├── perlib_plot_settings.txt │ │ │ ├── regularize.cxx │ │ │ ├── smooth.cxx │ │ │ └── smooth_all.sh │ │ ├── generic_bench │ │ │ ├── bench.hh │ │ │ ├── bench_parameter.hh │ │ │ ├── btl.hh │ │ │ ├── init │ │ │ │ ├── init_function.hh │ │ │ │ ├── init_matrix.hh │ │ │ │ └── init_vector.hh │ │ │ ├── static │ │ │ │ ├── bench_static.hh │ │ │ │ ├── intel_bench_fixed_size.hh │ │ │ │ └── static_size_generator.hh │ │ │ ├── timers │ │ │ │ ├── STL_perf_analyzer.hh │ │ │ │ ├── STL_timer.hh │ │ │ │ ├── mixed_perf_analyzer.hh │ │ │ │ ├── portable_perf_analyzer.hh │ │ │ │ ├── portable_perf_analyzer_old.hh │ │ │ │ ├── portable_timer.hh │ │ │ │ ├── x86_perf_analyzer.hh │ │ │ │ └── x86_timer.hh │ │ │ └── utils │ │ │ │ ├── size_lin_log.hh │ │ │ │ ├── size_log.hh │ │ │ │ ├── utilities.h │ │ │ │ └── xy_file.hh │ │ └── libs │ │ │ ├── BLAS │ │ │ ├── CMakeLists.txt │ │ │ ├── blas.h │ │ │ ├── blas_interface.hh │ │ │ ├── blas_interface_impl.hh │ │ │ ├── c_interface_base.h │ │ │ └── main.cpp │ │ │ ├── STL │ │ │ ├── CMakeLists.txt │ │ │ ├── STL_interface.hh │ │ │ └── main.cpp │ │ │ ├── blitz │ │ │ ├── CMakeLists.txt │ │ │ ├── blitz_LU_solve_interface.hh │ │ │ ├── blitz_interface.hh │ │ │ ├── btl_blitz.cpp │ │ │ ├── btl_tiny_blitz.cpp │ │ │ └── tiny_blitz_interface.hh │ │ │ ├── eigen2 │ │ │ ├── CMakeLists.txt │ │ │ ├── btl_tiny_eigen2.cpp │ │ │ ├── eigen2_interface.hh │ │ │ ├── main_adv.cpp │ │ │ ├── main_linear.cpp │ │ │ ├── main_matmat.cpp │ │ │ └── main_vecmat.cpp │ │ │ ├── eigen3 │ │ │ ├── CMakeLists.txt │ │ │ ├── btl_tiny_eigen3.cpp │ │ │ ├── eigen3_interface.hh │ │ │ ├── main_adv.cpp │ │ │ ├── main_linear.cpp │ │ │ ├── main_matmat.cpp │ │ │ └── main_vecmat.cpp │ │ │ ├── gmm │ │ │ ├── CMakeLists.txt │ │ │ ├── gmm_LU_solve_interface.hh │ │ │ ├── gmm_interface.hh │ │ │ └── main.cpp │ │ │ ├── mtl4 │ │ │ ├── .kdbgrc.main │ │ │ ├── CMakeLists.txt │ │ │ ├── main.cpp │ │ │ ├── mtl4_LU_solve_interface.hh │ │ │ └── mtl4_interface.hh │ │ │ ├── tvmet │ │ │ ├── CMakeLists.txt │ │ │ ├── main.cpp │ │ │ └── tvmet_interface.hh │ │ │ └── ublas │ │ │ ├── CMakeLists.txt │ │ │ ├── main.cpp │ │ │ └── ublas_interface.hh │ ├── check_cache_queries.cpp │ ├── eig33.cpp │ ├── geometry.cpp │ ├── product_threshold.cpp │ ├── quat_slerp.cpp │ ├── quatmul.cpp │ ├── sparse_cholesky.cpp │ ├── sparse_dense_product.cpp │ ├── sparse_lu.cpp │ ├── sparse_product.cpp │ ├── sparse_randomsetter.cpp │ ├── sparse_setter.cpp │ ├── sparse_transpose.cpp │ ├── sparse_trisolver.cpp │ ├── spbench │ │ ├── CMakeLists.txt │ │ ├── spbenchsolver.cpp │ │ └── spbenchsolver.h │ ├── spmv.cpp │ └── vdw_new.cpp ├── blas │ ├── BandTriangularSolver.h │ ├── CMakeLists.txt │ ├── README.txt │ ├── chbmv.f │ ├── chpmv.f │ ├── chpr.f │ ├── chpr2.f │ ├── common.h │ ├── complex_double.cpp │ ├── complex_single.cpp │ ├── complexdots.f │ ├── ctbmv.f │ ├── ctpmv.f │ ├── ctpsv.f │ ├── double.cpp │ ├── drotm.f │ ├── drotmg.f │ ├── dsbmv.f │ ├── dspmv.f │ ├── dspr.f │ ├── dspr2.f │ ├── dtbmv.f │ ├── dtpmv.f │ ├── dtpsv.f │ ├── level1_cplx_impl.h │ ├── level1_impl.h │ ├── level1_real_impl.h │ ├── level2_cplx_impl.h │ ├── level2_impl.h │ ├── level2_real_impl.h │ ├── level3_impl.h │ ├── lsame.f │ ├── single.cpp │ ├── srotm.f │ ├── srotmg.f │ ├── ssbmv.f │ ├── sspmv.f │ ├── sspr.f │ ├── sspr2.f │ ├── stbmv.f │ ├── stpmv.f │ ├── stpsv.f │ ├── testing │ │ ├── CMakeLists.txt │ │ ├── cblat1.f │ │ ├── cblat2.dat │ │ ├── cblat2.f │ │ ├── cblat3.dat │ │ ├── cblat3.f │ │ ├── dblat1.f │ │ ├── dblat2.dat │ │ ├── dblat2.f │ │ ├── dblat3.dat │ │ ├── dblat3.f │ │ ├── runblastest.sh │ │ ├── sblat1.f │ │ ├── sblat2.dat │ │ ├── sblat2.f │ │ ├── sblat3.dat │ │ ├── sblat3.f │ │ ├── zblat1.f │ │ ├── zblat2.dat │ │ ├── zblat2.f │ │ ├── zblat3.dat │ │ └── zblat3.f │ ├── xerbla.cpp │ ├── zhbmv.f │ ├── zhpmv.f │ ├── zhpr.f │ ├── zhpr2.f │ ├── ztbmv.f │ ├── ztpmv.f │ └── ztpsv.f ├── cmake │ ├── CMakeDetermineVSServicePack.cmake │ ├── EigenConfigureTesting.cmake │ ├── EigenDetermineOSVersion.cmake │ ├── EigenTesting.cmake │ ├── FindAdolc.cmake │ ├── FindBLAS.cmake │ ├── FindCholmod.cmake │ ├── FindEigen2.cmake │ ├── FindEigen3.cmake │ ├── FindFFTW.cmake │ ├── FindGLEW.cmake │ ├── FindGMP.cmake │ ├── FindGSL.cmake │ ├── FindGoogleHash.cmake │ ├── FindLAPACK.cmake │ ├── FindMPFR.cmake │ ├── FindMetis.cmake │ ├── FindPastix.cmake │ ├── FindScotch.cmake │ ├── FindStandardMathLibrary.cmake │ ├── FindSuperLU.cmake │ ├── FindUmfpack.cmake │ ├── RegexUtils.cmake │ └── language_support.cmake ├── debug │ ├── gdb │ │ ├── __init__.py │ │ └── printers.py │ └── msvc │ │ └── eigen_autoexp_part.dat ├── demos │ ├── CMakeLists.txt │ ├── mandelbrot │ │ ├── CMakeLists.txt │ │ ├── README │ │ ├── mandelbrot.cpp │ │ └── mandelbrot.h │ ├── mix_eigen_and_c │ │ ├── README │ │ ├── binary_library.cpp │ │ ├── binary_library.h │ │ └── example.c │ └── opengl │ │ ├── CMakeLists.txt │ │ ├── README │ │ ├── camera.cpp │ │ ├── camera.h │ │ ├── gpuhelper.cpp │ │ ├── gpuhelper.h │ │ ├── icosphere.cpp │ │ ├── icosphere.h │ │ ├── quaternion_demo.cpp │ │ ├── quaternion_demo.h │ │ ├── trackball.cpp │ │ └── trackball.h ├── doc │ ├── A05_PortingFrom2To3.dox │ ├── A10_Eigen2SupportModes.dox │ ├── AsciiQuickReference.txt │ ├── B01_Experimental.dox │ ├── C00_QuickStartGuide.dox │ ├── C01_TutorialMatrixClass.dox │ ├── C02_TutorialMatrixArithmetic.dox │ ├── C03_TutorialArrayClass.dox │ ├── C04_TutorialBlockOperations.dox │ ├── C05_TutorialAdvancedInitialization.dox │ ├── C06_TutorialLinearAlgebra.dox │ ├── C07_TutorialReductionsVisitorsBroadcasting.dox │ ├── C08_TutorialGeometry.dox │ ├── C09_TutorialSparse.dox │ ├── C10_TutorialMapClass.dox │ ├── CMakeLists.txt │ ├── D01_StlContainers.dox │ ├── D03_WrongStackAlignment.dox │ ├── D07_PassingByValue.dox │ ├── D09_StructHavingEigenMembers.dox │ ├── D11_UnalignedArrayAssert.dox │ ├── Doxyfile.in │ ├── Eigen_Silly_Professor_64x64.png │ ├── I00_CustomizingEigen.dox │ ├── I01_TopicLazyEvaluation.dox │ ├── I02_HiPerformance.dox │ ├── I03_InsideEigenExample.dox │ ├── I05_FixedSizeVectorizable.dox │ ├── I06_TopicEigenExpressionTemplates.dox │ ├── I07_TopicScalarTypes.dox │ ├── I08_Resizing.dox │ ├── I09_Vectorization.dox │ ├── I10_Assertions.dox │ ├── I11_Aliasing.dox │ ├── I12_ClassHierarchy.dox │ ├── I13_FunctionsTakingEigenTypes.dox │ ├── I14_PreprocessorDirectives.dox │ ├── I15_StorageOrders.dox │ ├── I16_TemplateKeyword.dox │ ├── Overview.dox │ ├── QuickReference.dox │ ├── SparseQuickReference.dox │ ├── TopicLinearAlgebraDecompositions.dox │ ├── TopicMultithreading.dox │ ├── TutorialSparse_example_details.dox │ ├── UsingIntelMKL.dox │ ├── eigendoxy.css │ ├── eigendoxy_footer.html.in │ ├── eigendoxy_header.html.in │ ├── eigendoxy_tabs.css │ ├── examples │ │ ├── .krazy │ │ ├── CMakeLists.txt │ │ ├── DenseBase_middleCols_int.cpp │ │ ├── DenseBase_middleRows_int.cpp │ │ ├── DenseBase_template_int_middleCols.cpp │ │ ├── DenseBase_template_int_middleRows.cpp │ │ ├── MatrixBase_cwise_const.cpp │ │ ├── QuickStart_example.cpp │ │ ├── QuickStart_example2_dynamic.cpp │ │ ├── QuickStart_example2_fixed.cpp │ │ ├── TemplateKeyword_flexible.cpp │ │ ├── TemplateKeyword_simple.cpp │ │ ├── TutorialLinAlgComputeTwice.cpp │ │ ├── TutorialLinAlgExComputeSolveError.cpp │ │ ├── TutorialLinAlgExSolveColPivHouseholderQR.cpp │ │ ├── TutorialLinAlgExSolveLDLT.cpp │ │ ├── TutorialLinAlgInverseDeterminant.cpp │ │ ├── TutorialLinAlgRankRevealing.cpp │ │ ├── TutorialLinAlgSVDSolve.cpp │ │ ├── TutorialLinAlgSelfAdjointEigenSolver.cpp │ │ ├── TutorialLinAlgSetThreshold.cpp │ │ ├── Tutorial_ArrayClass_accessors.cpp │ │ ├── Tutorial_ArrayClass_addition.cpp │ │ ├── Tutorial_ArrayClass_cwise_other.cpp │ │ ├── Tutorial_ArrayClass_interop.cpp │ │ ├── Tutorial_ArrayClass_interop_matrix.cpp │ │ ├── Tutorial_ArrayClass_mult.cpp │ │ ├── Tutorial_BlockOperations_block_assignment.cpp │ │ ├── Tutorial_BlockOperations_colrow.cpp │ │ ├── Tutorial_BlockOperations_corner.cpp │ │ ├── Tutorial_BlockOperations_print_block.cpp │ │ ├── Tutorial_BlockOperations_vector.cpp │ │ ├── Tutorial_PartialLU_solve.cpp │ │ ├── Tutorial_ReductionsVisitorsBroadcasting_broadcast_1nn.cpp │ │ ├── Tutorial_ReductionsVisitorsBroadcasting_broadcast_simple.cpp │ │ ├── Tutorial_ReductionsVisitorsBroadcasting_broadcast_simple_rowwise.cpp │ │ ├── Tutorial_ReductionsVisitorsBroadcasting_colwise.cpp │ │ ├── Tutorial_ReductionsVisitorsBroadcasting_maxnorm.cpp │ │ ├── Tutorial_ReductionsVisitorsBroadcasting_reductions_bool.cpp │ │ ├── Tutorial_ReductionsVisitorsBroadcasting_reductions_norm.cpp │ │ ├── Tutorial_ReductionsVisitorsBroadcasting_rowwise.cpp │ │ ├── Tutorial_ReductionsVisitorsBroadcasting_visitors.cpp │ │ ├── Tutorial_simple_example_dynamic_size.cpp │ │ ├── Tutorial_simple_example_fixed_size.cpp │ │ ├── class_Block.cpp │ │ ├── class_CwiseBinaryOp.cpp │ │ ├── class_CwiseUnaryOp.cpp │ │ ├── class_CwiseUnaryOp_ptrfun.cpp │ │ ├── class_FixedBlock.cpp │ │ ├── class_FixedVectorBlock.cpp │ │ ├── class_VectorBlock.cpp │ │ ├── function_taking_eigenbase.cpp │ │ ├── tut_arithmetic_add_sub.cpp │ │ ├── tut_arithmetic_dot_cross.cpp │ │ ├── tut_arithmetic_matrix_mul.cpp │ │ ├── tut_arithmetic_redux_basic.cpp │ │ ├── tut_arithmetic_scalar_mul_div.cpp │ │ ├── tut_matrix_coefficient_accessors.cpp │ │ ├── tut_matrix_resize.cpp │ │ └── tut_matrix_resize_fixed_size.cpp │ ├── snippets │ │ ├── .krazy │ │ ├── AngleAxis_mimic_euler.cpp │ │ ├── CMakeLists.txt │ │ ├── ColPivHouseholderQR_solve.cpp │ │ ├── ComplexEigenSolver_compute.cpp │ │ ├── ComplexEigenSolver_eigenvalues.cpp │ │ ├── ComplexEigenSolver_eigenvectors.cpp │ │ ├── ComplexSchur_compute.cpp │ │ ├── ComplexSchur_matrixT.cpp │ │ ├── ComplexSchur_matrixU.cpp │ │ ├── Cwise_abs.cpp │ │ ├── Cwise_abs2.cpp │ │ ├── Cwise_acos.cpp │ │ ├── Cwise_boolean_and.cpp │ │ ├── Cwise_boolean_or.cpp │ │ ├── Cwise_cos.cpp │ │ ├── Cwise_cube.cpp │ │ ├── Cwise_equal_equal.cpp │ │ ├── Cwise_exp.cpp │ │ ├── Cwise_greater.cpp │ │ ├── Cwise_greater_equal.cpp │ │ ├── Cwise_inverse.cpp │ │ ├── Cwise_less.cpp │ │ ├── Cwise_less_equal.cpp │ │ ├── Cwise_log.cpp │ │ ├── Cwise_max.cpp │ │ ├── Cwise_min.cpp │ │ ├── Cwise_minus.cpp │ │ ├── Cwise_minus_equal.cpp │ │ ├── Cwise_not_equal.cpp │ │ ├── Cwise_plus.cpp │ │ ├── Cwise_plus_equal.cpp │ │ ├── Cwise_pow.cpp │ │ ├── Cwise_product.cpp │ │ ├── Cwise_quotient.cpp │ │ ├── Cwise_sin.cpp │ │ ├── Cwise_slash_equal.cpp │ │ ├── Cwise_sqrt.cpp │ │ ├── Cwise_square.cpp │ │ ├── Cwise_tan.cpp │ │ ├── Cwise_times_equal.cpp │ │ ├── DenseBase_LinSpaced.cpp │ │ ├── DenseBase_LinSpaced_seq.cpp │ │ ├── DenseBase_setLinSpaced.cpp │ │ ├── DirectionWise_replicate.cpp │ │ ├── DirectionWise_replicate_int.cpp │ │ ├── EigenSolver_EigenSolver_MatrixType.cpp │ │ ├── EigenSolver_compute.cpp │ │ ├── EigenSolver_eigenvalues.cpp │ │ ├── EigenSolver_eigenvectors.cpp │ │ ├── EigenSolver_pseudoEigenvectors.cpp │ │ ├── FullPivHouseholderQR_solve.cpp │ │ ├── FullPivLU_image.cpp │ │ ├── FullPivLU_kernel.cpp │ │ ├── FullPivLU_solve.cpp │ │ ├── HessenbergDecomposition_compute.cpp │ │ ├── HessenbergDecomposition_matrixH.cpp │ │ ├── HessenbergDecomposition_packedMatrix.cpp │ │ ├── HouseholderQR_solve.cpp │ │ ├── HouseholderSequence_HouseholderSequence.cpp │ │ ├── IOFormat.cpp │ │ ├── JacobiSVD_basic.cpp │ │ ├── Jacobi_makeGivens.cpp │ │ ├── Jacobi_makeJacobi.cpp │ │ ├── LLT_example.cpp │ │ ├── LLT_solve.cpp │ │ ├── Map_general_stride.cpp │ │ ├── Map_inner_stride.cpp │ │ ├── Map_outer_stride.cpp │ │ ├── Map_placement_new.cpp │ │ ├── Map_simple.cpp │ │ ├── MatrixBase_adjoint.cpp │ │ ├── MatrixBase_all.cpp │ │ ├── MatrixBase_array.cpp │ │ ├── MatrixBase_array_const.cpp │ │ ├── MatrixBase_asDiagonal.cpp │ │ ├── MatrixBase_block_int_int.cpp │ │ ├── MatrixBase_block_int_int_int_int.cpp │ │ ├── MatrixBase_bottomLeftCorner_int_int.cpp │ │ ├── MatrixBase_bottomRightCorner_int_int.cpp │ │ ├── MatrixBase_bottomRows_int.cpp │ │ ├── MatrixBase_cast.cpp │ │ ├── MatrixBase_col.cpp │ │ ├── MatrixBase_colwise.cpp │ │ ├── MatrixBase_computeInverseAndDetWithCheck.cpp │ │ ├── MatrixBase_computeInverseWithCheck.cpp │ │ ├── MatrixBase_cwiseAbs.cpp │ │ ├── MatrixBase_cwiseAbs2.cpp │ │ ├── MatrixBase_cwiseEqual.cpp │ │ ├── MatrixBase_cwiseInverse.cpp │ │ ├── MatrixBase_cwiseMax.cpp │ │ ├── MatrixBase_cwiseMin.cpp │ │ ├── MatrixBase_cwiseNotEqual.cpp │ │ ├── MatrixBase_cwiseProduct.cpp │ │ ├── MatrixBase_cwiseQuotient.cpp │ │ ├── MatrixBase_cwiseSqrt.cpp │ │ ├── MatrixBase_diagonal.cpp │ │ ├── MatrixBase_diagonal_int.cpp │ │ ├── MatrixBase_diagonal_template_int.cpp │ │ ├── MatrixBase_eigenvalues.cpp │ │ ├── MatrixBase_end_int.cpp │ │ ├── MatrixBase_eval.cpp │ │ ├── MatrixBase_extract.cpp │ │ ├── MatrixBase_fixedBlock_int_int.cpp │ │ ├── MatrixBase_identity.cpp │ │ ├── MatrixBase_identity_int_int.cpp │ │ ├── MatrixBase_inverse.cpp │ │ ├── MatrixBase_isDiagonal.cpp │ │ ├── MatrixBase_isIdentity.cpp │ │ ├── MatrixBase_isOnes.cpp │ │ ├── MatrixBase_isOrthogonal.cpp │ │ ├── MatrixBase_isUnitary.cpp │ │ ├── MatrixBase_isZero.cpp │ │ ├── MatrixBase_leftCols_int.cpp │ │ ├── MatrixBase_marked.cpp │ │ ├── MatrixBase_noalias.cpp │ │ ├── MatrixBase_ones.cpp │ │ ├── MatrixBase_ones_int.cpp │ │ ├── MatrixBase_ones_int_int.cpp │ │ ├── MatrixBase_operatorNorm.cpp │ │ ├── MatrixBase_part.cpp │ │ ├── MatrixBase_prod.cpp │ │ ├── MatrixBase_random.cpp │ │ ├── MatrixBase_random_int.cpp │ │ ├── MatrixBase_random_int_int.cpp │ │ ├── MatrixBase_replicate.cpp │ │ ├── MatrixBase_replicate_int_int.cpp │ │ ├── MatrixBase_reverse.cpp │ │ ├── MatrixBase_rightCols_int.cpp │ │ ├── MatrixBase_row.cpp │ │ ├── MatrixBase_rowwise.cpp │ │ ├── MatrixBase_segment_int_int.cpp │ │ ├── MatrixBase_select.cpp │ │ ├── MatrixBase_set.cpp │ │ ├── MatrixBase_setIdentity.cpp │ │ ├── MatrixBase_setOnes.cpp │ │ ├── MatrixBase_setRandom.cpp │ │ ├── MatrixBase_setZero.cpp │ │ ├── MatrixBase_start_int.cpp │ │ ├── MatrixBase_template_int_bottomRows.cpp │ │ ├── MatrixBase_template_int_end.cpp │ │ ├── MatrixBase_template_int_int_bottomLeftCorner.cpp │ │ ├── MatrixBase_template_int_int_bottomRightCorner.cpp │ │ ├── MatrixBase_template_int_int_topLeftCorner.cpp │ │ ├── MatrixBase_template_int_int_topRightCorner.cpp │ │ ├── MatrixBase_template_int_leftCols.cpp │ │ ├── MatrixBase_template_int_rightCols.cpp │ │ ├── MatrixBase_template_int_segment.cpp │ │ ├── MatrixBase_template_int_start.cpp │ │ ├── MatrixBase_template_int_topRows.cpp │ │ ├── MatrixBase_topLeftCorner_int_int.cpp │ │ ├── MatrixBase_topRightCorner_int_int.cpp │ │ ├── MatrixBase_topRows_int.cpp │ │ ├── MatrixBase_transpose.cpp │ │ ├── MatrixBase_zero.cpp │ │ ├── MatrixBase_zero_int.cpp │ │ ├── MatrixBase_zero_int_int.cpp │ │ ├── Matrix_resize_NoChange_int.cpp │ │ ├── Matrix_resize_int.cpp │ │ ├── Matrix_resize_int_NoChange.cpp │ │ ├── Matrix_resize_int_int.cpp │ │ ├── Matrix_setConstant_int.cpp │ │ ├── Matrix_setConstant_int_int.cpp │ │ ├── Matrix_setIdentity_int_int.cpp │ │ ├── Matrix_setOnes_int.cpp │ │ ├── Matrix_setOnes_int_int.cpp │ │ ├── Matrix_setRandom_int.cpp │ │ ├── Matrix_setRandom_int_int.cpp │ │ ├── Matrix_setZero_int.cpp │ │ ├── Matrix_setZero_int_int.cpp │ │ ├── PartialPivLU_solve.cpp │ │ ├── PartialRedux_count.cpp │ │ ├── PartialRedux_maxCoeff.cpp │ │ ├── PartialRedux_minCoeff.cpp │ │ ├── PartialRedux_norm.cpp │ │ ├── PartialRedux_prod.cpp │ │ ├── PartialRedux_squaredNorm.cpp │ │ ├── PartialRedux_sum.cpp │ │ ├── RealSchur_RealSchur_MatrixType.cpp │ │ ├── RealSchur_compute.cpp │ │ ├── SelfAdjointEigenSolver_SelfAdjointEigenSolver.cpp │ │ ├── SelfAdjointEigenSolver_SelfAdjointEigenSolver_MatrixType.cpp │ │ ├── SelfAdjointEigenSolver_SelfAdjointEigenSolver_MatrixType2.cpp │ │ ├── SelfAdjointEigenSolver_compute_MatrixType.cpp │ │ ├── SelfAdjointEigenSolver_compute_MatrixType2.cpp │ │ ├── SelfAdjointEigenSolver_eigenvalues.cpp │ │ ├── SelfAdjointEigenSolver_eigenvectors.cpp │ │ ├── SelfAdjointEigenSolver_operatorInverseSqrt.cpp │ │ ├── SelfAdjointEigenSolver_operatorSqrt.cpp │ │ ├── SelfAdjointView_eigenvalues.cpp │ │ ├── SelfAdjointView_operatorNorm.cpp │ │ ├── TopicAliasing_block.cpp │ │ ├── TopicAliasing_block_correct.cpp │ │ ├── TopicAliasing_cwise.cpp │ │ ├── TopicAliasing_mult1.cpp │ │ ├── TopicAliasing_mult2.cpp │ │ ├── TopicAliasing_mult3.cpp │ │ ├── TopicStorageOrders_example.cpp │ │ ├── Tridiagonalization_Tridiagonalization_MatrixType.cpp │ │ ├── Tridiagonalization_compute.cpp │ │ ├── Tridiagonalization_decomposeInPlace.cpp │ │ ├── Tridiagonalization_diagonal.cpp │ │ ├── Tridiagonalization_householderCoefficients.cpp │ │ ├── Tridiagonalization_packedMatrix.cpp │ │ ├── Tutorial_AdvancedInitialization_Block.cpp │ │ ├── Tutorial_AdvancedInitialization_CommaTemporary.cpp │ │ ├── Tutorial_AdvancedInitialization_Join.cpp │ │ ├── Tutorial_AdvancedInitialization_LinSpaced.cpp │ │ ├── Tutorial_AdvancedInitialization_ThreeWays.cpp │ │ ├── Tutorial_AdvancedInitialization_Zero.cpp │ │ ├── Tutorial_Map_rowmajor.cpp │ │ ├── Tutorial_Map_using.cpp │ │ ├── Tutorial_commainit_01.cpp │ │ ├── Tutorial_commainit_01b.cpp │ │ ├── Tutorial_commainit_02.cpp │ │ ├── Tutorial_solve_matrix_inverse.cpp │ │ ├── Tutorial_solve_multiple_rhs.cpp │ │ ├── Tutorial_solve_reuse_decomposition.cpp │ │ ├── Tutorial_solve_singular.cpp │ │ ├── Tutorial_solve_triangular.cpp │ │ ├── Tutorial_solve_triangular_inplace.cpp │ │ ├── Vectorwise_reverse.cpp │ │ ├── class_FullPivLU.cpp │ │ ├── compile_snippet.cpp.in │ │ ├── tut_arithmetic_redux_minmax.cpp │ │ ├── tut_arithmetic_transpose_aliasing.cpp │ │ ├── tut_arithmetic_transpose_conjugate.cpp │ │ ├── tut_arithmetic_transpose_inplace.cpp │ │ └── tut_matrix_assignment_resizing.cpp │ ├── special_examples │ │ ├── CMakeLists.txt │ │ ├── Tutorial_sparse_example.cpp │ │ └── Tutorial_sparse_example_details.cpp │ └── tutorial.cpp ├── eigen3.pc.in ├── failtest │ ├── CMakeLists.txt │ ├── block_nonconst_ctor_on_const_xpr_0.cpp │ ├── block_nonconst_ctor_on_const_xpr_1.cpp │ ├── block_nonconst_ctor_on_const_xpr_2.cpp │ ├── block_on_const_type_actually_const_0.cpp │ ├── block_on_const_type_actually_const_1.cpp │ ├── const_qualified_block_method_retval_0.cpp │ ├── const_qualified_block_method_retval_1.cpp │ ├── const_qualified_diagonal_method_retval.cpp │ ├── const_qualified_transpose_method_retval.cpp │ ├── diagonal_nonconst_ctor_on_const_xpr.cpp │ ├── diagonal_on_const_type_actually_const.cpp │ ├── failtest_sanity_check.cpp │ ├── map_nonconst_ctor_on_const_ptr_0.cpp │ ├── map_nonconst_ctor_on_const_ptr_1.cpp │ ├── map_nonconst_ctor_on_const_ptr_2.cpp │ ├── map_nonconst_ctor_on_const_ptr_3.cpp │ ├── map_nonconst_ctor_on_const_ptr_4.cpp │ ├── map_on_const_type_actually_const_0.cpp │ ├── map_on_const_type_actually_const_1.cpp │ ├── transpose_nonconst_ctor_on_const_xpr.cpp │ └── transpose_on_const_type_actually_const.cpp ├── lapack │ ├── CMakeLists.txt │ ├── cholesky.cpp │ ├── complex_double.cpp │ ├── complex_single.cpp │ ├── double.cpp │ ├── eigenvalues.cpp │ ├── lapack_common.h │ ├── lu.cpp │ └── single.cpp ├── scripts │ ├── CMakeLists.txt │ ├── buildtests.in │ ├── check.in │ ├── debug.in │ ├── eigen_gen_credits.cpp │ ├── eigen_gen_docs │ ├── release.in │ └── relicense.py ├── signature_of_eigen3_matrix_library ├── test-patched ├── test-patched.cpp ├── test │ ├── CMakeLists.txt │ ├── adjoint.cpp │ ├── array.cpp │ ├── array_for_matrix.cpp │ ├── array_replicate.cpp │ ├── array_reverse.cpp │ ├── bandmatrix.cpp │ ├── basicstuff.cpp │ ├── bicgstab.cpp │ ├── block.cpp │ ├── cholesky.cpp │ ├── cholmod_support.cpp │ ├── commainitializer.cpp │ ├── conjugate_gradient.cpp │ ├── conservative_resize.cpp │ ├── corners.cpp │ ├── cwiseop.cpp │ ├── determinant.cpp │ ├── diagonal.cpp │ ├── diagonalmatrices.cpp │ ├── dontalign.cpp │ ├── dynalloc.cpp │ ├── eigen2 │ │ ├── CMakeLists.txt │ │ ├── eigen2_adjoint.cpp │ │ ├── eigen2_alignedbox.cpp │ │ ├── eigen2_array.cpp │ │ ├── eigen2_basicstuff.cpp │ │ ├── eigen2_bug_132.cpp │ │ ├── eigen2_cholesky.cpp │ │ ├── eigen2_commainitializer.cpp │ │ ├── eigen2_cwiseop.cpp │ │ ├── eigen2_determinant.cpp │ │ ├── eigen2_dynalloc.cpp │ │ ├── eigen2_eigensolver.cpp │ │ ├── eigen2_first_aligned.cpp │ │ ├── eigen2_geometry.cpp │ │ ├── eigen2_geometry_with_eigen2_prefix.cpp │ │ ├── eigen2_hyperplane.cpp │ │ ├── eigen2_inverse.cpp │ │ ├── eigen2_linearstructure.cpp │ │ ├── eigen2_lu.cpp │ │ ├── eigen2_map.cpp │ │ ├── eigen2_meta.cpp │ │ ├── eigen2_miscmatrices.cpp │ │ ├── eigen2_mixingtypes.cpp │ │ ├── eigen2_newstdvector.cpp │ │ ├── eigen2_nomalloc.cpp │ │ ├── eigen2_packetmath.cpp │ │ ├── eigen2_parametrizedline.cpp │ │ ├── eigen2_prec_inverse_4x4.cpp │ │ ├── eigen2_product_large.cpp │ │ ├── eigen2_product_small.cpp │ │ ├── eigen2_qr.cpp │ │ ├── eigen2_qtvector.cpp │ │ ├── eigen2_regression.cpp │ │ ├── eigen2_sizeof.cpp │ │ ├── eigen2_smallvectors.cpp │ │ ├── eigen2_sparse_basic.cpp │ │ ├── eigen2_sparse_product.cpp │ │ ├── eigen2_sparse_solvers.cpp │ │ ├── eigen2_sparse_vector.cpp │ │ ├── eigen2_stdvector.cpp │ │ ├── eigen2_submatrices.cpp │ │ ├── eigen2_sum.cpp │ │ ├── eigen2_svd.cpp │ │ ├── eigen2_swap.cpp │ │ ├── eigen2_triangular.cpp │ │ ├── eigen2_unalignedassert.cpp │ │ ├── eigen2_visitor.cpp │ │ ├── gsl_helper.h │ │ ├── main.h │ │ ├── product.h │ │ ├── runtest.sh │ │ ├── sparse.h │ │ └── testsuite.cmake │ ├── eigen2support.cpp │ ├── eigensolver_complex.cpp │ ├── eigensolver_generic.cpp │ ├── eigensolver_selfadjoint.cpp │ ├── exceptions.cpp │ ├── first_aligned.cpp │ ├── geo_alignedbox.cpp │ ├── geo_eulerangles.cpp │ ├── geo_homogeneous.cpp │ ├── geo_hyperplane.cpp │ ├── geo_orthomethods.cpp │ ├── geo_parametrizedline.cpp │ ├── geo_quaternion.cpp │ ├── geo_transformations.cpp │ ├── hessenberg.cpp │ ├── householder.cpp │ ├── integer_types.cpp │ ├── inverse.cpp │ ├── jacobi.cpp │ ├── jacobisvd.cpp │ ├── linearstructure.cpp │ ├── lu.cpp │ ├── main.h │ ├── map.cpp │ ├── mapstaticmethods.cpp │ ├── mapstride.cpp │ ├── meta.cpp │ ├── miscmatrices.cpp │ ├── mixingtypes.cpp │ ├── nesting_ops.cpp │ ├── nomalloc.cpp │ ├── nullary.cpp │ ├── packetmath.cpp │ ├── pardiso_support.cpp │ ├── pastix_support.cpp │ ├── permutationmatrices.cpp │ ├── prec_inverse_4x4.cpp │ ├── product.h │ ├── product_extra.cpp │ ├── product_large.cpp │ ├── product_mmtr.cpp │ ├── product_notemporary.cpp │ ├── product_selfadjoint.cpp │ ├── product_small.cpp │ ├── product_symm.cpp │ ├── product_syrk.cpp │ ├── product_trmm.cpp │ ├── product_trmv.cpp │ ├── product_trsolve.cpp │ ├── qr.cpp │ ├── qr_colpivoting.cpp │ ├── qr_fullpivoting.cpp │ ├── qtvector.cpp │ ├── redux.cpp │ ├── resize.cpp │ ├── runtest.sh │ ├── schur_complex.cpp │ ├── schur_real.cpp │ ├── selfadjoint.cpp │ ├── simplicial_cholesky.cpp │ ├── sizeof.cpp │ ├── sizeoverflow.cpp │ ├── smallvectors.cpp │ ├── sparse.h │ ├── sparse_basic.cpp │ ├── sparse_permutations.cpp │ ├── sparse_product.cpp │ ├── sparse_solver.h │ ├── sparse_solvers.cpp │ ├── sparse_vector.cpp │ ├── stable_norm.cpp │ ├── stddeque.cpp │ ├── stdlist.cpp │ ├── stdvector.cpp │ ├── stdvector_overload.cpp │ ├── superlu_support.cpp │ ├── swap.cpp │ ├── testsuite.cmake │ ├── triangular.cpp │ ├── umeyama.cpp │ ├── umfpack_support.cpp │ ├── unalignedassert.cpp │ ├── unalignedcount.cpp │ ├── upperbidiagonalization.cpp │ ├── vectorization_logic.cpp │ ├── vectorwiseop.cpp │ ├── visitor.cpp │ └── zerosized.cpp └── unsupported │ ├── CMakeLists.txt │ ├── Eigen │ ├── AdolcForward │ ├── AlignedVector3 │ ├── AutoDiff │ ├── BVH │ ├── CMakeLists.txt │ ├── FFT │ ├── IterativeSolvers │ ├── KroneckerProduct │ ├── MPRealSupport │ ├── MatrixFunctions │ ├── MoreVectorization │ ├── NonLinearOptimization │ ├── NumericalDiff │ ├── OpenGLSupport │ ├── Polynomials │ ├── Skyline │ ├── SparseExtra │ ├── Splines │ └── src │ │ ├── AutoDiff │ │ ├── AutoDiffJacobian.h │ │ ├── AutoDiffScalar.h │ │ ├── AutoDiffVector.h │ │ └── CMakeLists.txt │ │ ├── BVH │ │ ├── BVAlgorithms.h │ │ ├── CMakeLists.txt │ │ └── KdBVH.h │ │ ├── CMakeLists.txt │ │ ├── FFT │ │ ├── CMakeLists.txt │ │ ├── ei_fftw_impl.h │ │ └── ei_kissfft_impl.h │ │ ├── IterativeSolvers │ │ ├── CMakeLists.txt │ │ ├── ConstrainedConjGrad.h │ │ ├── GMRES.h │ │ ├── IncompleteLU.h │ │ ├── IterationController.h │ │ └── Scaling.h │ │ ├── KroneckerProduct │ │ ├── CMakeLists.txt │ │ └── KroneckerTensorProduct.h │ │ ├── MatrixFunctions │ │ ├── CMakeLists.txt │ │ ├── MatrixExponential.h │ │ ├── MatrixFunction.h │ │ ├── MatrixFunctionAtomic.h │ │ ├── MatrixLogarithm.h │ │ ├── MatrixSquareRoot.h │ │ └── StemFunction.h │ │ ├── MoreVectorization │ │ ├── CMakeLists.txt │ │ └── MathFunctions.h │ │ ├── NonLinearOptimization │ │ ├── CMakeLists.txt │ │ ├── HybridNonLinearSolver.h │ │ ├── LevenbergMarquardt.h │ │ ├── chkder.h │ │ ├── covar.h │ │ ├── dogleg.h │ │ ├── fdjac1.h │ │ ├── lmpar.h │ │ ├── qrsolv.h │ │ ├── r1mpyq.h │ │ ├── r1updt.h │ │ └── rwupdt.h │ │ ├── NumericalDiff │ │ ├── CMakeLists.txt │ │ └── NumericalDiff.h │ │ ├── Polynomials │ │ ├── CMakeLists.txt │ │ ├── Companion.h │ │ ├── PolynomialSolver.h │ │ └── PolynomialUtils.h │ │ ├── Skyline │ │ ├── CMakeLists.txt │ │ ├── SkylineInplaceLU.h │ │ ├── SkylineMatrix.h │ │ ├── SkylineMatrixBase.h │ │ ├── SkylineProduct.h │ │ ├── SkylineStorage.h │ │ └── SkylineUtil.h │ │ ├── SparseExtra │ │ ├── BlockOfDynamicSparseMatrix.h │ │ ├── CMakeLists.txt │ │ ├── DynamicSparseMatrix.h │ │ ├── MarketIO.h │ │ ├── MatrixMarketIterator.h │ │ └── RandomSetter.h │ │ └── Splines │ │ ├── CMakeLists.txt │ │ ├── Spline.h │ │ ├── SplineFitting.h │ │ └── SplineFwd.h │ ├── README.txt │ ├── doc │ ├── CMakeLists.txt │ ├── Doxyfile.in │ ├── Overview.dox │ ├── examples │ │ ├── BVH_Example.cpp │ │ ├── CMakeLists.txt │ │ ├── FFT.cpp │ │ ├── MatrixExponential.cpp │ │ ├── MatrixFunction.cpp │ │ ├── MatrixLogarithm.cpp │ │ ├── MatrixSine.cpp │ │ ├── MatrixSinh.cpp │ │ ├── MatrixSquareRoot.cpp │ │ ├── PolynomialSolver1.cpp │ │ └── PolynomialUtils1.cpp │ └── snippets │ │ └── CMakeLists.txt │ └── test │ ├── BVH.cpp │ ├── CMakeLists.txt │ ├── FFT.cpp │ ├── FFTW.cpp │ ├── NonLinearOptimization.cpp │ ├── NumericalDiff.cpp │ ├── alignedvector3.cpp │ ├── autodiff.cpp │ ├── forward_adolc.cpp │ ├── gmres.cpp │ ├── kronecker_product.cpp │ ├── matrix_exponential.cpp │ ├── matrix_function.cpp │ ├── matrix_square_root.cpp │ ├── mpreal │ ├── dlmalloc.c │ ├── dlmalloc.h │ ├── mpreal.cpp │ └── mpreal.h │ ├── mpreal_support.cpp │ ├── openglsupport.cpp │ ├── polynomialsolver.cpp │ ├── polynomialutils.cpp │ ├── sparse_extra.cpp │ └── splines.cpp ├── gui ├── codeinput.cpp ├── codeinput.h ├── codewidget.cpp ├── codewidget.h ├── consoleoutput.cpp ├── consoleoutput.h ├── widget.cpp └── widget.h ├── main.cpp ├── qling.pro ├── qling ├── jiteventlistener.cpp ├── jiteventlistener.h ├── qling.cpp ├── qling.h ├── tests.cpp └── tests.h ├── qt-hack ├── QtCore │ ├── QAtomicInt │ ├── QAtomicPointer │ ├── qatomic.h │ ├── qatomic_bootstrap.h │ ├── qatomic_x86_64.h │ └── qbasicatomic.h ├── qatomic.h ├── qatomic_arch.h ├── qatomic_sun.s ├── qatomic_x86_64.h └── qbasicatomic.h ├── setupenv.sh ├── streamlistener.cpp └── streamlistener.h /.gitignore: -------------------------------------------------------------------------------- 1 | *.o 2 | moc_* 3 | ui_* 4 | *.pro.user* 5 | log 6 | Makefile* 7 | bin/ 8 | core 9 | core.* 10 | tst_* 11 | .cling_history 12 | -------------------------------------------------------------------------------- /eigen-patched/COPYING.README: -------------------------------------------------------------------------------- 1 | Eigen is primarily MPL2 licensed. See COPYING.MPL2 and these links: 2 | http://www.mozilla.org/MPL/2.0/ 3 | http://www.mozilla.org/MPL/2.0/FAQ.html 4 | 5 | Some files contain third-party code under BSD or LGPL licenses, whence the other 6 | COPYING.* files here. 7 | 8 | All the LGPL code is either LGPL 2.1-only, or LGPL 2.1-or-later. 9 | For this reason, the COPYING.LGPL file contains the LGPL 2.1 text. 10 | 11 | If you want to guarantee that the Eigen code that you are #including is licensed 12 | under the MPL2 and possibly more permissive licenses (like BSD), #define this 13 | preprocessor symbol: 14 | EIGEN_MPL2_ONLY 15 | For example, with most compilers, you could add this to your project CXXFLAGS: 16 | -DEIGEN_MPL2_ONLY 17 | This will cause a compilation error to be generated if you #include any code that is 18 | LGPL licensed. 19 | -------------------------------------------------------------------------------- /eigen-patched/CTestConfig.cmake: -------------------------------------------------------------------------------- 1 | ## This file should be placed in the root directory of your project. 2 | ## Then modify the CMakeLists.txt file in the root directory of your 3 | ## project to incorporate the testing dashboard. 4 | ## # The following are required to uses Dart and the Cdash dashboard 5 | ## ENABLE_TESTING() 6 | ## INCLUDE(CTest) 7 | set(CTEST_PROJECT_NAME "Eigen") 8 | set(CTEST_NIGHTLY_START_TIME "00:00:00 UTC") 9 | 10 | set(CTEST_DROP_METHOD "http") 11 | set(CTEST_DROP_SITE "manao.inria.fr") 12 | set(CTEST_DROP_LOCATION "/CDash/submit.php?project=Eigen") 13 | set(CTEST_DROP_SITE_CDASH TRUE) 14 | -------------------------------------------------------------------------------- /eigen-patched/CTestCustom.cmake.in: -------------------------------------------------------------------------------- 1 | 2 | ## A tribute to Dynamic! 3 | set(CTEST_CUSTOM_MAXIMUM_NUMBER_OF_WARNINGS "33331") 4 | set(CTEST_CUSTOM_MAXIMUM_NUMBER_OF_ERRORS "33331") 5 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/Array: -------------------------------------------------------------------------------- 1 | #ifndef EIGEN_ARRAY_MODULE_H 2 | #define EIGEN_ARRAY_MODULE_H 3 | 4 | // include Core first to handle Eigen2 support macros 5 | #include "Core" 6 | 7 | #ifndef EIGEN2_SUPPORT 8 | #error The Eigen/Array header does no longer exist in Eigen3. All that functionality has moved to Eigen/Core. 9 | #endif 10 | 11 | #endif // EIGEN_ARRAY_MODULE_H 12 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | include(RegexUtils) 2 | test_escape_string_as_regex() 3 | 4 | file(GLOB Eigen_directory_files "*") 5 | 6 | escape_string_as_regex(ESCAPED_CMAKE_CURRENT_SOURCE_DIR "${CMAKE_CURRENT_SOURCE_DIR}") 7 | 8 | foreach(f ${Eigen_directory_files}) 9 | if(NOT f MATCHES "\\.txt" AND NOT f MATCHES "${ESCAPED_CMAKE_CURRENT_SOURCE_DIR}/[.].+" AND NOT f MATCHES "${ESCAPED_CMAKE_CURRENT_SOURCE_DIR}/src") 10 | list(APPEND Eigen_directory_files_to_install ${f}) 11 | endif() 12 | endforeach(f ${Eigen_directory_files}) 13 | 14 | install(FILES 15 | ${Eigen_directory_files_to_install} 16 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen COMPONENT Devel 17 | ) 18 | 19 | add_subdirectory(src) 20 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/Dense: -------------------------------------------------------------------------------- 1 | #include "Core" 2 | #include "LU" 3 | #include "Cholesky" 4 | #include "QR" 5 | #include "SVD" 6 | #include "Geometry" 7 | #include "Eigenvalues" 8 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/Eigen: -------------------------------------------------------------------------------- 1 | #include "Dense" 2 | //#include "Sparse" 3 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/Householder: -------------------------------------------------------------------------------- 1 | #ifndef EIGEN_HOUSEHOLDER_MODULE_H 2 | #define EIGEN_HOUSEHOLDER_MODULE_H 3 | 4 | #include "Core" 5 | 6 | #include "src/Core/util/DisableStupidWarnings.h" 7 | 8 | /** \defgroup Householder_Module Householder module 9 | * This module provides Householder transformations. 10 | * 11 | * \code 12 | * #include 13 | * \endcode 14 | */ 15 | 16 | #include "src/Householder/Householder.h" 17 | #include "src/Householder/HouseholderSequence.h" 18 | #include "src/Householder/BlockHouseholder.h" 19 | 20 | #include "src/Core/util/ReenableStupidWarnings.h" 21 | 22 | #endif // EIGEN_HOUSEHOLDER_MODULE_H 23 | /* vim: set filetype=cpp et sw=2 ts=2 ai: */ 24 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/Jacobi: -------------------------------------------------------------------------------- 1 | #ifndef EIGEN_JACOBI_MODULE_H 2 | #define EIGEN_JACOBI_MODULE_H 3 | 4 | #include "Core" 5 | 6 | #include "src/Core/util/DisableStupidWarnings.h" 7 | 8 | /** \defgroup Jacobi_Module Jacobi module 9 | * This module provides Jacobi and Givens rotations. 10 | * 11 | * \code 12 | * #include 13 | * \endcode 14 | * 15 | * In addition to listed classes, it defines the two following MatrixBase methods to apply a Jacobi or Givens rotation: 16 | * - MatrixBase::applyOnTheLeft() 17 | * - MatrixBase::applyOnTheRight(). 18 | */ 19 | 20 | #include "src/Jacobi/Jacobi.h" 21 | 22 | #include "src/Core/util/ReenableStupidWarnings.h" 23 | 24 | #endif // EIGEN_JACOBI_MODULE_H 25 | /* vim: set filetype=cpp et sw=2 ts=2 ai: */ 26 | 27 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/LeastSquares: -------------------------------------------------------------------------------- 1 | #ifndef EIGEN_REGRESSION_MODULE_H 2 | #define EIGEN_REGRESSION_MODULE_H 3 | 4 | #ifndef EIGEN2_SUPPORT 5 | #error LeastSquares is only available in Eigen2 support mode (define EIGEN2_SUPPORT) 6 | #endif 7 | 8 | // exclude from normal eigen3-only documentation 9 | #ifdef EIGEN2_SUPPORT 10 | 11 | #include "Core" 12 | 13 | #include "src/Core/util/DisableStupidWarnings.h" 14 | 15 | #include "Eigenvalues" 16 | #include "Geometry" 17 | 18 | /** \defgroup LeastSquares_Module LeastSquares module 19 | * This module provides linear regression and related features. 20 | * 21 | * \code 22 | * #include 23 | * \endcode 24 | */ 25 | 26 | #include "src/Eigen2Support/LeastSquares.h" 27 | 28 | #include "src/Core/util/ReenableStupidWarnings.h" 29 | 30 | #endif // EIGEN2_SUPPORT 31 | 32 | #endif // EIGEN_REGRESSION_MODULE_H 33 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/OrderingMethods: -------------------------------------------------------------------------------- 1 | #ifndef EIGEN_ORDERINGMETHODS_MODULE_H 2 | #define EIGEN_ORDERINGMETHODS_MODULE_H 3 | 4 | #include "SparseCore" 5 | 6 | #include "src/Core/util/DisableStupidWarnings.h" 7 | 8 | /** \ingroup Sparse_modules 9 | * \defgroup OrderingMethods_Module OrderingMethods module 10 | * 11 | * This module is currently for internal use only. 12 | * 13 | * 14 | * \code 15 | * #include 16 | * \endcode 17 | */ 18 | 19 | #include "src/OrderingMethods/Amd.h" 20 | 21 | #include "src/Core/util/ReenableStupidWarnings.h" 22 | 23 | #endif // EIGEN_ORDERINGMETHODS_MODULE_H 24 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/QtAlignedMalloc: -------------------------------------------------------------------------------- 1 | 2 | #ifndef EIGEN_QTMALLOC_MODULE_H 3 | #define EIGEN_QTMALLOC_MODULE_H 4 | 5 | #include "Core" 6 | 7 | #if (!EIGEN_MALLOC_ALREADY_ALIGNED) 8 | 9 | #include "src/Core/util/DisableStupidWarnings.h" 10 | 11 | void *qMalloc(size_t size) 12 | { 13 | return Eigen::internal::aligned_malloc(size); 14 | } 15 | 16 | void qFree(void *ptr) 17 | { 18 | Eigen::internal::aligned_free(ptr); 19 | } 20 | 21 | void *qRealloc(void *ptr, size_t size) 22 | { 23 | void* newPtr = Eigen::internal::aligned_malloc(size); 24 | memcpy(newPtr, ptr, size); 25 | Eigen::internal::aligned_free(ptr); 26 | return newPtr; 27 | } 28 | 29 | #include "src/Core/util/ReenableStupidWarnings.h" 30 | 31 | #endif 32 | 33 | #endif // EIGEN_QTMALLOC_MODULE_H 34 | /* vim: set filetype=cpp et sw=2 ts=2 ai: */ 35 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/Sparse: -------------------------------------------------------------------------------- 1 | #ifndef EIGEN_SPARSE_MODULE_H 2 | #define EIGEN_SPARSE_MODULE_H 3 | 4 | /** \defgroup Sparse_modules Sparse modules 5 | * 6 | * Meta-module including all related modules: 7 | * - SparseCore 8 | * - OrderingMethods 9 | * - SparseCholesky 10 | * - IterativeLinearSolvers 11 | * 12 | * \code 13 | * #include 14 | * \endcode 15 | */ 16 | 17 | #include "SparseCore" 18 | #include "OrderingMethods" 19 | #include "SparseCholesky" 20 | #include "IterativeLinearSolvers" 21 | 22 | #endif // EIGEN_SPARSE_MODULE_H 23 | 24 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/StdList: -------------------------------------------------------------------------------- 1 | // This file is part of Eigen, a lightweight C++ template library 2 | // for linear algebra. 3 | // 4 | // Copyright (C) 2009 Hauke Heibel 5 | // 6 | // This Source Code Form is subject to the terms of the Mozilla 7 | // Public License v. 2.0. If a copy of the MPL was not distributed 8 | // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9 | 10 | #ifndef EIGEN_STDLIST_MODULE_H 11 | #define EIGEN_STDLIST_MODULE_H 12 | 13 | #include "Core" 14 | #include 15 | 16 | #if (defined(_MSC_VER) && defined(_WIN64)) /* MSVC auto aligns in 64 bit builds */ 17 | 18 | #define EIGEN_DEFINE_STL_LIST_SPECIALIZATION(...) 19 | 20 | #else 21 | 22 | #include "src/StlSupport/StdList.h" 23 | 24 | #endif 25 | 26 | #endif // EIGEN_STDLIST_MODULE_H 27 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | file(GLOB Eigen_src_subdirectories "*") 2 | escape_string_as_regex(ESCAPED_CMAKE_CURRENT_SOURCE_DIR "${CMAKE_CURRENT_SOURCE_DIR}") 3 | foreach(f ${Eigen_src_subdirectories}) 4 | if(NOT f MATCHES "\\.txt" AND NOT f MATCHES "${ESCAPED_CMAKE_CURRENT_SOURCE_DIR}/[.].+" ) 5 | add_subdirectory(${f}) 6 | endif() 7 | endforeach() 8 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Cholesky/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Cholesky_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Cholesky_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Cholesky COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/CholmodSupport/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_CholmodSupport_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_CholmodSupport_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/CholmodSupport COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Core/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Core_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Core_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Core COMPONENT Devel 6 | ) 7 | 8 | ADD_SUBDIRECTORY(products) 9 | ADD_SUBDIRECTORY(util) 10 | ADD_SUBDIRECTORY(arch) 11 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Core/arch/AltiVec/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Core_arch_AltiVec_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Core_arch_AltiVec_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Core/arch/AltiVec COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Core/arch/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | ADD_SUBDIRECTORY(SSE) 2 | ADD_SUBDIRECTORY(AltiVec) 3 | ADD_SUBDIRECTORY(NEON) 4 | ADD_SUBDIRECTORY(Default) 5 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Core/arch/Default/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Core_arch_Default_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Core_arch_Default_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Core/arch/Default COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Core/arch/NEON/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Core_arch_NEON_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Core_arch_NEON_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Core/arch/NEON COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Core/arch/SSE/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Core_arch_SSE_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Core_arch_SSE_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Core/arch/SSE COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Core/products/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Core_Product_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Core_Product_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Core/products COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Core/util/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Core_util_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Core_util_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Core/util COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Core/util/NonMPL2.h: -------------------------------------------------------------------------------- 1 | #ifdef EIGEN_MPL2_ONLY 2 | #error Including non-MPL2 code in EIGEN_MPL2_ONLY mode 3 | #endif 4 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Core/util/ReenableStupidWarnings.h: -------------------------------------------------------------------------------- 1 | #ifdef EIGEN_WARNINGS_DISABLED 2 | #undef EIGEN_WARNINGS_DISABLED 3 | 4 | #ifndef EIGEN_PERMANENTLY_DISABLE_STUPID_WARNINGS 5 | #ifdef _MSC_VER 6 | #pragma warning( pop ) 7 | #elif defined __INTEL_COMPILER 8 | #pragma warning pop 9 | #elif defined __clang__ 10 | #pragma clang diagnostic pop 11 | #endif 12 | #endif 13 | 14 | #endif // EIGEN_WARNINGS_DISABLED 15 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Eigen2Support/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Eigen2Support_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Eigen2Support_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Eigen2Support COMPONENT Devel 6 | ) 7 | 8 | ADD_SUBDIRECTORY(Geometry) -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Eigen2Support/Geometry/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Eigen2Support_Geometry_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Eigen2Support_Geometry_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Eigen2Support/Geometry 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Eigen2Support/Macros.h: -------------------------------------------------------------------------------- 1 | // This file is part of Eigen, a lightweight C++ template library 2 | // for linear algebra. 3 | // 4 | // Copyright (C) 2011 Benoit Jacob 5 | // 6 | // This Source Code Form is subject to the terms of the Mozilla 7 | // Public License v. 2.0. If a copy of the MPL was not distributed 8 | // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9 | 10 | #ifndef EIGEN2_MACROS_H 11 | #define EIGEN2_MACROS_H 12 | 13 | #define ei_assert eigen_assert 14 | #define ei_internal_assert eigen_internal_assert 15 | 16 | #define EIGEN_ALIGN_128 EIGEN_ALIGN16 17 | 18 | #define EIGEN_ARCH_WANTS_ALIGNMENT EIGEN_ALIGN_STATICALLY 19 | 20 | #endif // EIGEN2_MACROS_H 21 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Eigenvalues/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_EIGENVALUES_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_EIGENVALUES_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Eigenvalues COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Geometry/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Geometry_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Geometry_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Geometry COMPONENT Devel 6 | ) 7 | 8 | ADD_SUBDIRECTORY(arch) 9 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Geometry/arch/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Geometry_arch_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Geometry_arch_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Geometry/arch COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Householder/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Householder_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Householder_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Householder COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/IterativeLinearSolvers/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_IterativeLinearSolvers_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_IterativeLinearSolvers_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/IterativeLinearSolvers COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/Jacobi/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Jacobi_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Jacobi_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/Jacobi COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/LU/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_LU_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_LU_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/LU COMPONENT Devel 6 | ) 7 | 8 | ADD_SUBDIRECTORY(arch) 9 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/LU/arch/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_LU_arch_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_LU_arch_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/LU/arch COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/OrderingMethods/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_OrderingMethods_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_OrderingMethods_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/OrderingMethods COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/PaStiXSupport/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_PastixSupport_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_PastixSupport_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/PaStiXSupport COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/PardisoSupport/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_PardisoSupport_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_PardisoSupport_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/PardisoSupport COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/QR/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_QR_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_QR_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/QR COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/SVD/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_SVD_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_SVD_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/SVD COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/SparseCholesky/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_SparseCholesky_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_SparseCholesky_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/SparseCholesky COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/SparseCore/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_SparseCore_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_SparseCore_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/SparseCore COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/SparseCore/SparseAssign.h: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/herrgahr/qling/22e56c4be0bbccb1d0437f610bfa37374b44b87f/eigen-patched/Eigen/src/SparseCore/SparseAssign.h -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/StlSupport/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_StlSupport_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_StlSupport_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/StlSupport COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/SuperLUSupport/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_SuperLUSupport_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_SuperLUSupport_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/SuperLUSupport COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/UmfPackSupport/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_UmfPackSupport_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_UmfPackSupport_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/UmfPackSupport COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/misc/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_misc_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_misc_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/misc COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/Eigen/src/plugins/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_plugins_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_plugins_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/Eigen/src/plugins COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/bench/bench_multi_compilers.sh: -------------------------------------------------------------------------------- 1 | #!/bin/bash 2 | 3 | if (($# < 2)); then 4 | echo "Usage: $0 compilerlist.txt benchfile.cpp" 5 | else 6 | 7 | compilerlist=$1 8 | benchfile=$2 9 | 10 | g=0 11 | source $compilerlist 12 | 13 | # for each compiler, compile benchfile and run the benchmark 14 | for (( i=0 ; i /dev/null 22 | echo "" 23 | else 24 | echo "compiler not found: $compiler" 25 | fi 26 | done 27 | 28 | fi 29 | -------------------------------------------------------------------------------- /eigen-patched/bench/bench_sum.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | using namespace Eigen; 4 | using namespace std; 5 | 6 | int main() 7 | { 8 | typedef Matrix Vec; 9 | Vec v(SIZE); 10 | v.setZero(); 11 | v[0] = 1; 12 | v[1] = 2; 13 | for(int i = 0; i < 1000000; i++) 14 | { 15 | v.coeffRef(0) += v.sum() * SCALAR(1e-20); 16 | } 17 | cout << v.sum() << endl; 18 | } 19 | -------------------------------------------------------------------------------- /eigen-patched/bench/bench_unrolling: -------------------------------------------------------------------------------- 1 | #!/bin/bash 2 | 3 | # gcc : CXX="g++ -finline-limit=10000 -ftemplate-depth-2000 --param max-inline-recursive-depth=2000" 4 | # icc : CXX="icpc -fast -no-inline-max-size -fno-exceptions" 5 | CXX=${CXX-g++ -finline-limit=10000 -ftemplate-depth-2000 --param max-inline-recursive-depth=2000} # default value 6 | 7 | for ((i=1; i<16; ++i)); do 8 | echo "Matrix size: $i x $i :" 9 | $CXX -O3 -I.. -DNDEBUG benchmark.cpp -DMATSIZE=$i -DEIGEN_UNROLLING_LIMIT=400 -o benchmark && time ./benchmark >/dev/null 10 | $CXX -O3 -I.. -DNDEBUG -finline-limit=10000 benchmark.cpp -DMATSIZE=$i -DEIGEN_DONT_USE_UNROLLED_LOOPS=1 -o benchmark && time ./benchmark >/dev/null 11 | echo " " 12 | done 13 | -------------------------------------------------------------------------------- /eigen-patched/bench/benchmarkX.cpp: -------------------------------------------------------------------------------- 1 | // g++ -fopenmp -I .. -O3 -DNDEBUG -finline-limit=1000 benchmarkX.cpp -o b && time ./b 2 | 3 | #include 4 | 5 | #include 6 | 7 | using namespace std; 8 | using namespace Eigen; 9 | 10 | #ifndef MATTYPE 11 | #define MATTYPE MatrixXLd 12 | #endif 13 | 14 | #ifndef MATSIZE 15 | #define MATSIZE 400 16 | #endif 17 | 18 | #ifndef REPEAT 19 | #define REPEAT 100 20 | #endif 21 | 22 | int main(int argc, char *argv[]) 23 | { 24 | MATTYPE I = MATTYPE::Ones(MATSIZE,MATSIZE); 25 | MATTYPE m(MATSIZE,MATSIZE); 26 | for(int i = 0; i < MATSIZE; i++) for(int j = 0; j < MATSIZE; j++) 27 | { 28 | m(i,j) = (i+j+1)/(MATSIZE*MATSIZE); 29 | } 30 | for(int a = 0; a < REPEAT; a++) 31 | { 32 | m = I + 0.0001 * (m + m*m); 33 | } 34 | cout << m(0,0) << endl; 35 | return 0; 36 | } 37 | -------------------------------------------------------------------------------- /eigen-patched/bench/benchmarkXcwise.cpp: -------------------------------------------------------------------------------- 1 | // g++ -O3 -DNDEBUG benchmarkX.cpp -o benchmarkX && time ./benchmarkX 2 | 3 | #include 4 | #include 5 | 6 | using namespace std; 7 | using namespace Eigen; 8 | 9 | #ifndef VECTYPE 10 | #define VECTYPE VectorXLd 11 | #endif 12 | 13 | #ifndef VECSIZE 14 | #define VECSIZE 1000000 15 | #endif 16 | 17 | #ifndef REPEAT 18 | #define REPEAT 1000 19 | #endif 20 | 21 | int main(int argc, char *argv[]) 22 | { 23 | VECTYPE I = VECTYPE::Ones(VECSIZE); 24 | VECTYPE m(VECSIZE,1); 25 | for(int i = 0; i < VECSIZE; i++) 26 | { 27 | m[i] = 0.1 * i/VECSIZE; 28 | } 29 | for(int a = 0; a < REPEAT; a++) 30 | { 31 | m = VECTYPE::Ones(VECSIZE) + 0.00005 * (m.cwise().square() + m/4); 32 | } 33 | cout << m[0] << endl; 34 | return 0; 35 | } 36 | -------------------------------------------------------------------------------- /eigen-patched/bench/btl/actions/basic_actions.hh: -------------------------------------------------------------------------------- 1 | 2 | #include "action_axpy.hh" 3 | #include "action_axpby.hh" 4 | 5 | #include "action_matrix_vector_product.hh" 6 | #include "action_atv_product.hh" 7 | 8 | #include "action_matrix_matrix_product.hh" 9 | // #include "action_ata_product.hh" 10 | #include "action_aat_product.hh" 11 | 12 | #include "action_trisolve.hh" 13 | #include "action_trmm.hh" 14 | #include "action_symv.hh" 15 | // #include "action_symm.hh" 16 | #include "action_syr2.hh" 17 | #include "action_ger.hh" 18 | #include "action_rot.hh" 19 | 20 | // #include "action_lu_solve.hh" 21 | 22 | -------------------------------------------------------------------------------- /eigen-patched/bench/btl/cmake/FindCBLAS.cmake: -------------------------------------------------------------------------------- 1 | # include(FindLibraryWithDebug) 2 | 3 | if (CBLAS_INCLUDES AND CBLAS_LIBRARIES) 4 | set(CBLAS_FIND_QUIETLY TRUE) 5 | endif (CBLAS_INCLUDES AND CBLAS_LIBRARIES) 6 | 7 | find_path(CBLAS_INCLUDES 8 | NAMES 9 | cblas.h 10 | PATHS 11 | $ENV{CBLASDIR}/include 12 | ${INCLUDE_INSTALL_DIR} 13 | ) 14 | 15 | find_library(CBLAS_LIBRARIES 16 | cblas 17 | PATHS 18 | $ENV{CBLASDIR}/lib 19 | ${LIB_INSTALL_DIR} 20 | ) 21 | 22 | find_file(CBLAS_LIBRARIES 23 | libcblas.so.3 24 | PATHS 25 | /usr/lib 26 | $ENV{CBLASDIR}/lib 27 | ${LIB_INSTALL_DIR} 28 | ) 29 | 30 | include(FindPackageHandleStandardArgs) 31 | find_package_handle_standard_args(CBLAS DEFAULT_MSG 32 | CBLAS_INCLUDES CBLAS_LIBRARIES) 33 | 34 | mark_as_advanced(CBLAS_INCLUDES CBLAS_LIBRARIES) 35 | -------------------------------------------------------------------------------- /eigen-patched/bench/btl/cmake/FindGMM.cmake: -------------------------------------------------------------------------------- 1 | if (GMM_INCLUDE_DIR) 2 | # in cache already 3 | set(GMM_FOUND TRUE) 4 | else (GMM_INCLUDE_DIR) 5 | 6 | find_path(GMM_INCLUDE_DIR NAMES gmm/gmm.h 7 | PATHS 8 | ${INCLUDE_INSTALL_DIR} 9 | ${GMM_INCLUDE_PATH} 10 | ) 11 | 12 | include(FindPackageHandleStandardArgs) 13 | FIND_PACKAGE_HANDLE_STANDARD_ARGS(GMM DEFAULT_MSG GMM_INCLUDE_DIR ) 14 | 15 | mark_as_advanced(GMM_INCLUDE_DIR) 16 | 17 | endif(GMM_INCLUDE_DIR) 18 | -------------------------------------------------------------------------------- /eigen-patched/bench/btl/cmake/FindGOTO.cmake: -------------------------------------------------------------------------------- 1 | 2 | if (GOTO_LIBRARIES) 3 | set(GOTO_FIND_QUIETLY TRUE) 4 | endif (GOTO_LIBRARIES) 5 | 6 | find_library(GOTO_LIBRARIES goto PATHS $ENV{GOTODIR} ${LIB_INSTALL_DIR}) 7 | 8 | if(GOTO_LIBRARIES AND CMAKE_COMPILER_IS_GNUCXX) 9 | set(GOTO_LIBRARIES ${GOTO_LIBRARIES} "-lpthread -lgfortran") 10 | endif() 11 | 12 | include(FindPackageHandleStandardArgs) 13 | find_package_handle_standard_args(GOTO DEFAULT_MSG GOTO_LIBRARIES) 14 | 15 | mark_as_advanced(GOTO_LIBRARIES) 16 | -------------------------------------------------------------------------------- /eigen-patched/bench/btl/cmake/FindGOTO2.cmake: -------------------------------------------------------------------------------- 1 | 2 | if (GOTO2_LIBRARIES) 3 | set(GOTO2_FIND_QUIETLY TRUE) 4 | endif (GOTO2_LIBRARIES) 5 | # 6 | # find_path(GOTO_INCLUDES 7 | # NAMES 8 | # cblas.h 9 | # PATHS 10 | # $ENV{GOTODIR}/include 11 | # ${INCLUDE_INSTALL_DIR} 12 | # ) 13 | 14 | find_file(GOTO2_LIBRARIES libgoto2.so PATHS /usr/lib $ENV{GOTO2DIR} ${LIB_INSTALL_DIR}) 15 | find_library(GOTO2_LIBRARIES goto2 PATHS $ENV{GOTO2DIR} ${LIB_INSTALL_DIR}) 16 | 17 | if(GOTO2_LIBRARIES AND CMAKE_COMPILER_IS_GNUCXX) 18 | set(GOTO2_LIBRARIES ${GOTO2_LIBRARIES} "-lpthread -lgfortran") 19 | endif() 20 | 21 | include(FindPackageHandleStandardArgs) 22 | find_package_handle_standard_args(GOTO2 DEFAULT_MSG 23 | GOTO2_LIBRARIES) 24 | 25 | mark_as_advanced(GOTO2_LIBRARIES) 26 | -------------------------------------------------------------------------------- /eigen-patched/bench/btl/generic_bench/static/intel_bench_fixed_size.hh: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/herrgahr/qling/22e56c4be0bbccb1d0437f610bfa37374b44b87f/eigen-patched/bench/btl/generic_bench/static/intel_bench_fixed_size.hh -------------------------------------------------------------------------------- /eigen-patched/bench/btl/generic_bench/static/static_size_generator.hh: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/herrgahr/qling/22e56c4be0bbccb1d0437f610bfa37374b44b87f/eigen-patched/bench/btl/generic_bench/static/static_size_generator.hh -------------------------------------------------------------------------------- /eigen-patched/bench/btl/generic_bench/timers/STL_perf_analyzer.hh: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/herrgahr/qling/22e56c4be0bbccb1d0437f610bfa37374b44b87f/eigen-patched/bench/btl/generic_bench/timers/STL_perf_analyzer.hh -------------------------------------------------------------------------------- /eigen-patched/bench/btl/generic_bench/timers/STL_timer.hh: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/herrgahr/qling/22e56c4be0bbccb1d0437f610bfa37374b44b87f/eigen-patched/bench/btl/generic_bench/timers/STL_timer.hh -------------------------------------------------------------------------------- /eigen-patched/bench/btl/generic_bench/timers/mixed_perf_analyzer.hh: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/herrgahr/qling/22e56c4be0bbccb1d0437f610bfa37374b44b87f/eigen-patched/bench/btl/generic_bench/timers/mixed_perf_analyzer.hh -------------------------------------------------------------------------------- /eigen-patched/bench/btl/generic_bench/utils/size_lin_log.hh: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/herrgahr/qling/22e56c4be0bbccb1d0437f610bfa37374b44b87f/eigen-patched/bench/btl/generic_bench/utils/size_lin_log.hh -------------------------------------------------------------------------------- /eigen-patched/bench/btl/libs/STL/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | 2 | btl_add_bench(btl_STL main.cpp OFF) 3 | -------------------------------------------------------------------------------- /eigen-patched/bench/btl/libs/blitz/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | 2 | find_package(Blitz) 3 | 4 | if (BLITZ_FOUND) 5 | include_directories(${BLITZ_INCLUDES}) 6 | 7 | btl_add_bench(btl_blitz btl_blitz.cpp) 8 | if (BUILD_btl_blitz) 9 | target_link_libraries(btl_blitz ${BLITZ_LIBRARIES}) 10 | endif (BUILD_btl_blitz) 11 | 12 | btl_add_bench(btl_tiny_blitz btl_tiny_blitz.cpp OFF) 13 | if (BUILD_btl_tiny_blitz) 14 | target_link_libraries(btl_tiny_blitz ${BLITZ_LIBRARIES}) 15 | endif (BUILD_btl_tiny_blitz) 16 | 17 | endif (BLITZ_FOUND) 18 | -------------------------------------------------------------------------------- /eigen-patched/bench/btl/libs/gmm/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | 2 | find_package(GMM) 3 | if (GMM_FOUND) 4 | include_directories(${GMM_INCLUDES}) 5 | btl_add_bench(btl_gmm main.cpp) 6 | endif (GMM_FOUND) 7 | -------------------------------------------------------------------------------- /eigen-patched/bench/btl/libs/mtl4/.kdbgrc.main: -------------------------------------------------------------------------------- 1 | [General] 2 | DebuggerCmdStr= 3 | DriverName=GDB 4 | FileVersion=1 5 | OptionsSelected= 6 | ProgramArgs= 7 | TTYLevel=7 8 | WorkingDirectory= 9 | 10 | [Memory] 11 | ColumnWidths=80,0 12 | NumExprs=0 13 | -------------------------------------------------------------------------------- /eigen-patched/bench/btl/libs/mtl4/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | 2 | find_package(MTL4) 3 | if (MTL4_FOUND) 4 | include_directories(${MTL4_INCLUDE_DIR}) 5 | btl_add_bench(btl_mtl4 main.cpp) 6 | endif (MTL4_FOUND) 7 | -------------------------------------------------------------------------------- /eigen-patched/bench/btl/libs/tvmet/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | 2 | find_package(Tvmet) 3 | if (TVMET_FOUND) 4 | include_directories(${TVMET_INCLUDE_DIR}) 5 | btl_add_bench(btl_tvmet main.cpp OFF) 6 | endif (TVMET_FOUND) 7 | -------------------------------------------------------------------------------- /eigen-patched/bench/btl/libs/ublas/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | 2 | find_package(Boost) 3 | if (Boost_FOUND) 4 | include_directories(${Boost_INCLUDE_DIRS}) 5 | include_directories(${Boost_INCLUDES}) 6 | btl_add_bench(btl_ublas main.cpp) 7 | endif (Boost_FOUND) 8 | -------------------------------------------------------------------------------- /eigen-patched/blas/README.txt: -------------------------------------------------------------------------------- 1 | 2 | This directory contains a BLAS library built on top of Eigen. 3 | 4 | This is currently a work in progress which is far to be ready for use, 5 | but feel free to contribute to it if you wish. 6 | 7 | This module is not built by default. In order to compile it, you need to 8 | type 'make blas' from within your build dir. 9 | 10 | -------------------------------------------------------------------------------- /eigen-patched/blas/complex_double.cpp: -------------------------------------------------------------------------------- 1 | // This file is part of Eigen, a lightweight C++ template library 2 | // for linear algebra. 3 | // 4 | // Copyright (C) 2009 Gael Guennebaud 5 | // 6 | // This Source Code Form is subject to the terms of the Mozilla 7 | // Public License v. 2.0. If a copy of the MPL was not distributed 8 | // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9 | 10 | #define SCALAR std::complex 11 | #define SCALAR_SUFFIX z 12 | #define SCALAR_SUFFIX_UP "Z" 13 | #define REAL_SCALAR_SUFFIX d 14 | #define ISCOMPLEX 1 15 | 16 | #include "level1_impl.h" 17 | #include "level1_cplx_impl.h" 18 | #include "level2_impl.h" 19 | #include "level2_cplx_impl.h" 20 | #include "level3_impl.h" 21 | -------------------------------------------------------------------------------- /eigen-patched/blas/complex_single.cpp: -------------------------------------------------------------------------------- 1 | // This file is part of Eigen, a lightweight C++ template library 2 | // for linear algebra. 3 | // 4 | // Copyright (C) 2009 Gael Guennebaud 5 | // 6 | // This Source Code Form is subject to the terms of the Mozilla 7 | // Public License v. 2.0. If a copy of the MPL was not distributed 8 | // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9 | 10 | #define SCALAR std::complex 11 | #define SCALAR_SUFFIX c 12 | #define SCALAR_SUFFIX_UP "C" 13 | #define REAL_SCALAR_SUFFIX s 14 | #define ISCOMPLEX 1 15 | 16 | #include "level1_impl.h" 17 | #include "level1_cplx_impl.h" 18 | #include "level2_impl.h" 19 | #include "level2_cplx_impl.h" 20 | #include "level3_impl.h" 21 | -------------------------------------------------------------------------------- /eigen-patched/blas/double.cpp: -------------------------------------------------------------------------------- 1 | // This file is part of Eigen, a lightweight C++ template library 2 | // for linear algebra. 3 | // 4 | // Copyright (C) 2009 Gael Guennebaud 5 | // 6 | // This Source Code Form is subject to the terms of the Mozilla 7 | // Public License v. 2.0. If a copy of the MPL was not distributed 8 | // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9 | 10 | #define SCALAR double 11 | #define SCALAR_SUFFIX d 12 | #define SCALAR_SUFFIX_UP "D" 13 | #define ISCOMPLEX 0 14 | 15 | #include "level1_impl.h" 16 | #include "level1_real_impl.h" 17 | #include "level2_impl.h" 18 | #include "level2_real_impl.h" 19 | #include "level3_impl.h" 20 | -------------------------------------------------------------------------------- /eigen-patched/blas/single.cpp: -------------------------------------------------------------------------------- 1 | // This file is part of Eigen, a lightweight C++ template library 2 | // for linear algebra. 3 | // 4 | // Copyright (C) 2009 Gael Guennebaud 5 | // 6 | // This Source Code Form is subject to the terms of the Mozilla 7 | // Public License v. 2.0. If a copy of the MPL was not distributed 8 | // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9 | 10 | #define SCALAR float 11 | #define SCALAR_SUFFIX s 12 | #define SCALAR_SUFFIX_UP "S" 13 | #define ISCOMPLEX 0 14 | 15 | #include "level1_impl.h" 16 | #include "level1_real_impl.h" 17 | #include "level2_impl.h" 18 | #include "level2_real_impl.h" 19 | #include "level3_impl.h" 20 | -------------------------------------------------------------------------------- /eigen-patched/blas/xerbla.cpp: -------------------------------------------------------------------------------- 1 | 2 | #include 3 | 4 | #if (defined __GNUC__) 5 | #define EIGEN_WEAK_LINKING __attribute__ ((weak)) 6 | #else 7 | #define EIGEN_WEAK_LINKING 8 | #endif 9 | 10 | #ifdef __cplusplus 11 | extern "C" 12 | { 13 | #endif 14 | 15 | EIGEN_WEAK_LINKING int xerbla_(const char * msg, int *info, int) 16 | { 17 | std::cerr << "Eigen BLAS ERROR #" << *info << ": " << msg << "\n"; 18 | return 0; 19 | } 20 | 21 | #ifdef __cplusplus 22 | } 23 | #endif 24 | -------------------------------------------------------------------------------- /eigen-patched/cmake/FindAdolc.cmake: -------------------------------------------------------------------------------- 1 | 2 | if (ADOLC_INCLUDES AND ADOLC_LIBRARIES) 3 | set(ADOLC_FIND_QUIETLY TRUE) 4 | endif (ADOLC_INCLUDES AND ADOLC_LIBRARIES) 5 | 6 | find_path(ADOLC_INCLUDES 7 | NAMES 8 | adolc/adouble.h 9 | PATHS 10 | $ENV{ADOLCDIR} 11 | ${INCLUDE_INSTALL_DIR} 12 | ) 13 | 14 | find_library(ADOLC_LIBRARIES adolc PATHS $ENV{ADOLCDIR} ${LIB_INSTALL_DIR}) 15 | 16 | include(FindPackageHandleStandardArgs) 17 | find_package_handle_standard_args(ADOLC DEFAULT_MSG 18 | ADOLC_INCLUDES ADOLC_LIBRARIES) 19 | 20 | mark_as_advanced(ADOLC_INCLUDES ADOLC_LIBRARIES) 21 | -------------------------------------------------------------------------------- /eigen-patched/cmake/FindGMP.cmake: -------------------------------------------------------------------------------- 1 | # Try to find the GNU Multiple Precision Arithmetic Library (GMP) 2 | # See http://gmplib.org/ 3 | 4 | if (GMP_INCLUDES AND GMP_LIBRARIES) 5 | set(GMP_FIND_QUIETLY TRUE) 6 | endif (GMP_INCLUDES AND GMP_LIBRARIES) 7 | 8 | find_path(GMP_INCLUDES 9 | NAMES 10 | gmp.h 11 | PATHS 12 | $ENV{GMPDIR} 13 | ${INCLUDE_INSTALL_DIR} 14 | ) 15 | 16 | find_library(GMP_LIBRARIES gmp PATHS $ENV{GMPDIR} ${LIB_INSTALL_DIR}) 17 | 18 | include(FindPackageHandleStandardArgs) 19 | find_package_handle_standard_args(GMP DEFAULT_MSG 20 | GMP_INCLUDES GMP_LIBRARIES) 21 | mark_as_advanced(GMP_INCLUDES GMP_LIBRARIES) 22 | -------------------------------------------------------------------------------- /eigen-patched/cmake/FindMetis.cmake: -------------------------------------------------------------------------------- 1 | # Pastix requires METIS or METIS (partitioning and reordering tools) 2 | 3 | if (METIS_INCLUDES AND METIS_LIBRARIES) 4 | set(METIS_FIND_QUIETLY TRUE) 5 | endif (METIS_INCLUDES AND METIS_LIBRARIES) 6 | 7 | find_path(METIS_INCLUDES 8 | NAMES 9 | metis.h 10 | PATHS 11 | $ENV{METISDIR} 12 | ${INCLUDE_INSTALL_DIR} 13 | PATH_SUFFIXES 14 | metis 15 | ) 16 | 17 | 18 | find_library(METIS_LIBRARIES metis PATHS $ENV{METISDIR} ${LIB_INSTALL_DIR}) 19 | 20 | include(FindPackageHandleStandardArgs) 21 | find_package_handle_standard_args(METIS DEFAULT_MSG 22 | METIS_INCLUDES METIS_LIBRARIES) 23 | 24 | mark_as_advanced(METIS_INCLUDES METIS_LIBRARIES) 25 | -------------------------------------------------------------------------------- /eigen-patched/cmake/FindPastix.cmake: -------------------------------------------------------------------------------- 1 | # Pastix lib requires linking to a blas library. 2 | # It is up to the user of this module to find a BLAS and link to it. 3 | # Pastix requires SCOTCH or METIS (partitioning and reordering tools) as well 4 | 5 | if (PASTIX_INCLUDES AND PASTIX_LIBRARIES) 6 | set(PASTIX_FIND_QUIETLY TRUE) 7 | endif (PASTIX_INCLUDES AND PASTIX_LIBRARIES) 8 | 9 | find_path(PASTIX_INCLUDES 10 | NAMES 11 | pastix_nompi.h 12 | PATHS 13 | $ENV{PASTIXDIR} 14 | ${INCLUDE_INSTALL_DIR} 15 | ) 16 | 17 | find_library(PASTIX_LIBRARIES pastix PATHS $ENV{PASTIXDIR} ${LIB_INSTALL_DIR}) 18 | 19 | 20 | 21 | include(FindPackageHandleStandardArgs) 22 | find_package_handle_standard_args(PASTIX DEFAULT_MSG 23 | PASTIX_INCLUDES PASTIX_LIBRARIES) 24 | 25 | mark_as_advanced(PASTIX_INCLUDES PASTIX_LIBRARIES) 26 | -------------------------------------------------------------------------------- /eigen-patched/cmake/FindScotch.cmake: -------------------------------------------------------------------------------- 1 | # Pastix requires SCOTCH or METIS (partitioning and reordering tools) 2 | 3 | if (SCOTCH_INCLUDES AND SCOTCH_LIBRARIES) 4 | set(SCOTCH_FIND_QUIETLY TRUE) 5 | endif (SCOTCH_INCLUDES AND SCOTCH_LIBRARIES) 6 | 7 | find_path(SCOTCH_INCLUDES 8 | NAMES 9 | scotch.h 10 | PATHS 11 | $ENV{SCOTCHDIR} 12 | ${INCLUDE_INSTALL_DIR} 13 | PATH_SUFFIXES 14 | scotch 15 | ) 16 | 17 | 18 | find_library(SCOTCH_LIBRARIES scotch PATHS $ENV{SCOTCHDIR} ${LIB_INSTALL_DIR}) 19 | 20 | include(FindPackageHandleStandardArgs) 21 | find_package_handle_standard_args(SCOTCH DEFAULT_MSG 22 | SCOTCH_INCLUDES SCOTCH_LIBRARIES) 23 | 24 | mark_as_advanced(SCOTCH_INCLUDES SCOTCH_LIBRARIES) 25 | -------------------------------------------------------------------------------- /eigen-patched/cmake/FindSuperLU.cmake: -------------------------------------------------------------------------------- 1 | 2 | # Umfpack lib usually requires linking to a blas library. 3 | # It is up to the user of this module to find a BLAS and link to it. 4 | 5 | if (SUPERLU_INCLUDES AND SUPERLU_LIBRARIES) 6 | set(SUPERLU_FIND_QUIETLY TRUE) 7 | endif (SUPERLU_INCLUDES AND SUPERLU_LIBRARIES) 8 | 9 | find_path(SUPERLU_INCLUDES 10 | NAMES 11 | supermatrix.h 12 | PATHS 13 | $ENV{SUPERLUDIR} 14 | ${INCLUDE_INSTALL_DIR} 15 | PATH_SUFFIXES 16 | superlu 17 | ) 18 | 19 | find_library(SUPERLU_LIBRARIES superlu PATHS $ENV{SUPERLUDIR} ${LIB_INSTALL_DIR}) 20 | 21 | include(FindPackageHandleStandardArgs) 22 | find_package_handle_standard_args(SUPERLU DEFAULT_MSG 23 | SUPERLU_INCLUDES SUPERLU_LIBRARIES) 24 | 25 | mark_as_advanced(SUPERLU_INCLUDES SUPERLU_LIBRARIES) 26 | -------------------------------------------------------------------------------- /eigen-patched/debug/gdb/__init__.py: -------------------------------------------------------------------------------- 1 | # Intentionally empty 2 | -------------------------------------------------------------------------------- /eigen-patched/demos/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | project(EigenDemos) 2 | 3 | add_custom_target(demos) 4 | 5 | if(NOT EIGEN_TEST_NOQT) 6 | find_package(Qt4) 7 | if(QT4_FOUND) 8 | add_subdirectory(mandelbrot) 9 | add_subdirectory(opengl) 10 | else(QT4_FOUND) 11 | message(STATUS "Qt4 not found, so disabling the mandelbrot and opengl demos") 12 | endif(QT4_FOUND) 13 | endif() 14 | -------------------------------------------------------------------------------- /eigen-patched/demos/mandelbrot/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | find_package(Qt4 REQUIRED) 2 | 3 | set(CMAKE_INCLUDE_CURRENT_DIR ON) 4 | 5 | if (CMAKE_COMPILER_IS_GNUCXX) 6 | set ( CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -O2") 7 | add_definitions ( "-DNDEBUG" ) 8 | endif (CMAKE_COMPILER_IS_GNUCXX) 9 | 10 | include_directories( ${QT_INCLUDE_DIR} ) 11 | 12 | set(mandelbrot_SRCS 13 | mandelbrot.cpp 14 | ) 15 | 16 | qt4_automoc(${mandelbrot_SRCS}) 17 | 18 | add_executable(mandelbrot ${mandelbrot_SRCS}) 19 | add_dependencies(demos mandelbrot) 20 | 21 | target_link_libraries(mandelbrot ${QT_QTCORE_LIBRARY} ${QT_QTGUI_LIBRARY}) 22 | -------------------------------------------------------------------------------- /eigen-patched/demos/mandelbrot/README: -------------------------------------------------------------------------------- 1 | *** Mandelbrot demo *** 2 | 3 | Controls: 4 | * Left mouse button to center view at a point. 5 | * Drag vertically with left mouse button to zoom in and out. 6 | 7 | Be sure to enable SSE2 or AltiVec to improve performance. 8 | 9 | The number of iterations, and the choice between single and double precision, are 10 | determined at runtime depending on the zoom level. 11 | -------------------------------------------------------------------------------- /eigen-patched/demos/mix_eigen_and_c/README: -------------------------------------------------------------------------------- 1 | This is an example of how one can wrap some of Eigen into a C library. 2 | 3 | To try this with GCC, do: 4 | 5 | g++ -c binary_library.cpp -O2 -msse2 -I ../.. 6 | gcc example.c binary_library.o -o example -lstdc++ 7 | ./example 8 | 9 | TODO: add CMakeLists, add more explanations here -------------------------------------------------------------------------------- /eigen-patched/demos/opengl/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | find_package(Qt4) 2 | find_package(OpenGL) 3 | 4 | if(QT4_FOUND AND OPENGL_FOUND) 5 | 6 | set(QT_USE_QTOPENGL TRUE) 7 | include(${QT_USE_FILE}) 8 | 9 | set(CMAKE_INCLUDE_CURRENT_DIR ON) 10 | 11 | include_directories( ${QT_INCLUDE_DIR} ) 12 | 13 | set(quaternion_demo_SRCS gpuhelper.cpp icosphere.cpp camera.cpp trackball.cpp quaternion_demo.cpp) 14 | 15 | qt4_automoc(${quaternion_demo_SRCS}) 16 | 17 | add_executable(quaternion_demo ${quaternion_demo_SRCS}) 18 | add_dependencies(demos quaternion_demo) 19 | 20 | target_link_libraries(quaternion_demo 21 | ${QT_QTCORE_LIBRARY} ${QT_QTGUI_LIBRARY} 22 | ${QT_QTOPENGL_LIBRARY} ${OPENGL_LIBRARIES} ) 23 | 24 | else() 25 | 26 | message(STATUS "OpenGL demo disabled because Qt4 and/or OpenGL have not been found.") 27 | 28 | endif() -------------------------------------------------------------------------------- /eigen-patched/demos/opengl/README: -------------------------------------------------------------------------------- 1 | 2 | Navigation: 3 | left button: rotate around the target 4 | middle button: zoom 5 | left button + ctrl quake rotate (rotate around camera position) 6 | middle button + ctrl walk (progress along camera's z direction) 7 | left button: pan (translate in the XY camera's plane) 8 | 9 | R : move the camera to initial position 10 | A : start/stop animation 11 | C : clear the animation 12 | G : add a key frame 13 | 14 | -------------------------------------------------------------------------------- /eigen-patched/doc/Eigen_Silly_Professor_64x64.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/herrgahr/qling/22e56c4be0bbccb1d0437f610bfa37374b44b87f/eigen-patched/doc/Eigen_Silly_Professor_64x64.png -------------------------------------------------------------------------------- /eigen-patched/doc/I06_TopicEigenExpressionTemplates.dox: -------------------------------------------------------------------------------- 1 | namespace Eigen { 2 | 3 | /** \page TopicEigenExpressionTemplates Expression templates in Eigen 4 | 5 | 6 | TODO: write this dox page! 7 | 8 | Is linked from the tutorial on arithmetic ops. 9 | 10 | */ 11 | 12 | } 13 | -------------------------------------------------------------------------------- /eigen-patched/doc/I07_TopicScalarTypes.dox: -------------------------------------------------------------------------------- 1 | namespace Eigen { 2 | 3 | /** \page TopicScalarTypes Scalar types 4 | 5 | 6 | TODO: write this dox page! 7 | 8 | Is linked from the tutorial on the Matrix class. 9 | 10 | */ 11 | 12 | } 13 | -------------------------------------------------------------------------------- /eigen-patched/doc/I08_Resizing.dox: -------------------------------------------------------------------------------- 1 | namespace Eigen { 2 | 3 | /** \page TopicResizing Resizing 4 | 5 | 6 | TODO: write this dox page! 7 | 8 | Is linked from the tutorial on the Matrix class. 9 | 10 | */ 11 | } 12 | -------------------------------------------------------------------------------- /eigen-patched/doc/I09_Vectorization.dox: -------------------------------------------------------------------------------- 1 | namespace Eigen { 2 | 3 | /** \page TopicVectorization Vectorization 4 | 5 | 6 | TODO: write this dox page! 7 | 8 | */ 9 | } 10 | -------------------------------------------------------------------------------- /eigen-patched/doc/I10_Assertions.dox: -------------------------------------------------------------------------------- 1 | namespace Eigen { 2 | 3 | /** \page TopicAssertions Assertions 4 | 5 | 6 | TODO: write this dox page! 7 | 8 | Is linked from the tutorial on matrix arithmetic. 9 | 10 | \sa Section \ref TopicPreprocessorDirectivesAssertions on page \ref TopicPreprocessorDirectives. 11 | 12 | */ 13 | } 14 | -------------------------------------------------------------------------------- /eigen-patched/doc/TutorialSparse_example_details.dox: -------------------------------------------------------------------------------- 1 | /** 2 | \page TutorialSparse_example_details 3 | \include Tutorial_sparse_example_details.cpp 4 | */ 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/eigendoxy_footer.html.in: -------------------------------------------------------------------------------- 1 |
2 | 3 | 4 | 14 | 15 | 16 | 17 | 18 | -------------------------------------------------------------------------------- /eigen-patched/doc/eigendoxy_header.html.in: -------------------------------------------------------------------------------- 1 | 2 | 3 | 4 | 5 | $title 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 15 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/.krazy: -------------------------------------------------------------------------------- 1 | EXCLUDE copyright 2 | EXCLUDE license 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | file(GLOB examples_SRCS "*.cpp") 2 | 3 | foreach(example_src ${examples_SRCS}) 4 | get_filename_component(example ${example_src} NAME_WE) 5 | add_executable(${example} ${example_src}) 6 | if(EIGEN_STANDARD_LIBRARIES_TO_LINK_TO) 7 | target_link_libraries(${example} ${EIGEN_STANDARD_LIBRARIES_TO_LINK_TO}) 8 | endif() 9 | get_target_property(example_executable 10 | ${example} LOCATION) 11 | add_custom_command( 12 | TARGET ${example} 13 | POST_BUILD 14 | COMMAND ${example_executable} 15 | ARGS >${CMAKE_CURRENT_BINARY_DIR}/${example}.out 16 | ) 17 | add_dependencies(all_examples ${example}) 18 | endforeach(example_src) 19 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/DenseBase_middleCols_int.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main(void) 8 | { 9 | int const N = 5; 10 | MatrixXi A(N,N); 11 | A.setRandom(); 12 | cout << "A =\n" << A << '\n' << endl; 13 | cout << "A(1..3,:) =\n" << A.middleCols(1,3) << endl; 14 | return 0; 15 | } 16 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/DenseBase_middleRows_int.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main(void) 8 | { 9 | int const N = 5; 10 | MatrixXi A(N,N); 11 | A.setRandom(); 12 | cout << "A =\n" << A << '\n' << endl; 13 | cout << "A(2..3,:) =\n" << A.middleRows(2,2) << endl; 14 | return 0; 15 | } 16 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/DenseBase_template_int_middleCols.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main(void) 8 | { 9 | int const N = 5; 10 | MatrixXi A(N,N); 11 | A.setRandom(); 12 | cout << "A =\n" << A << '\n' << endl; 13 | cout << "A(:,1..3) =\n" << A.middleCols<3>(1) << endl; 14 | return 0; 15 | } 16 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/DenseBase_template_int_middleRows.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main(void) 8 | { 9 | int const N = 5; 10 | MatrixXi A(N,N); 11 | A.setRandom(); 12 | cout << "A =\n" << A << '\n' << endl; 13 | cout << "A(1..3,:) =\n" << A.middleRows<3>(1) << endl; 14 | return 0; 15 | } 16 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/MatrixBase_cwise_const.cpp: -------------------------------------------------------------------------------- 1 | #define EIGEN2_SUPPORT 2 | #include 3 | #include 4 | 5 | using namespace Eigen; 6 | using namespace std; 7 | 8 | int main() 9 | { 10 | Matrix3i m = Matrix3i::Random(); 11 | cout << "Here is the matrix m:" << endl << m << endl; 12 | Matrix3i n = Matrix3i::Random(); 13 | cout << "And here is the matrix n:" << endl << n << endl; 14 | cout << "The coefficient-wise product of m and n is:" << endl; 15 | cout << m.cwise() * n << endl; 16 | cout << "Taking the cube of the coefficients of m yields:" << endl; 17 | cout << m.cwise().pow(3) << endl; 18 | } 19 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/QuickStart_example.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using Eigen::MatrixXd; 5 | 6 | int main() 7 | { 8 | MatrixXd m(2,2); 9 | m(0,0) = 3; 10 | m(1,0) = 2.5; 11 | m(0,1) = -1; 12 | m(1,1) = m(1,0) + m(0,1); 13 | std::cout << m << std::endl; 14 | } 15 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/QuickStart_example2_dynamic.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main() 8 | { 9 | MatrixXd m = MatrixXd::Random(3,3); 10 | m = (m + MatrixXd::Constant(3,3,1.2)) * 50; 11 | cout << "m =" << endl << m << endl; 12 | VectorXd v(3); 13 | v << 1, 2, 3; 14 | cout << "m * v =" << endl << m * v << endl; 15 | } 16 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/QuickStart_example2_fixed.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main() 8 | { 9 | Matrix3d m = Matrix3d::Random(); 10 | m = (m + Matrix3d::Constant(1.2)) * 50; 11 | cout << "m =" << endl << m << endl; 12 | Vector3d v(1,2,3); 13 | 14 | cout << "m * v =" << endl << m * v << endl; 15 | } 16 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/TemplateKeyword_flexible.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | template 7 | void copyUpperTriangularPart(MatrixBase& dst, const MatrixBase& src) 8 | { 9 | /* Note the 'template' keywords in the following line! */ 10 | dst.template triangularView() = src.template triangularView(); 11 | } 12 | 13 | int main() 14 | { 15 | MatrixXi m1 = MatrixXi::Ones(5,5); 16 | MatrixXi m2 = MatrixXi::Random(4,4); 17 | std::cout << "m2 before copy:" << std::endl; 18 | std::cout << m2 << std::endl << std::endl; 19 | copyUpperTriangularPart(m2, m1.topLeftCorner(4,4)); 20 | std::cout << "m2 after copy:" << std::endl; 21 | std::cout << m2 << std::endl << std::endl; 22 | } 23 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/TemplateKeyword_simple.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | void copyUpperTriangularPart(MatrixXf& dst, const MatrixXf& src) 7 | { 8 | dst.triangularView() = src.triangularView(); 9 | } 10 | 11 | int main() 12 | { 13 | MatrixXf m1 = MatrixXf::Ones(4,4); 14 | MatrixXf m2 = MatrixXf::Random(4,4); 15 | std::cout << "m2 before copy:" << std::endl; 16 | std::cout << m2 << std::endl << std::endl; 17 | copyUpperTriangularPart(m2, m1); 18 | std::cout << "m2 after copy:" << std::endl; 19 | std::cout << m2 << std::endl << std::endl; 20 | } 21 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/TutorialLinAlgComputeTwice.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | Matrix2f A, b; 10 | LLT llt; 11 | A << 2, -1, -1, 3; 12 | b << 1, 2, 3, 1; 13 | cout << "Here is the matrix A:\n" << A << endl; 14 | cout << "Here is the right hand side b:\n" << b << endl; 15 | cout << "Computing LLT decomposition..." << endl; 16 | llt.compute(A); 17 | cout << "The solution is:\n" << llt.solve(b) << endl; 18 | A(1,1)++; 19 | cout << "The matrix A is now:\n" << A << endl; 20 | cout << "Computing LLT decomposition..." << endl; 21 | llt.compute(A); 22 | cout << "The solution is now:\n" << llt.solve(b) << endl; 23 | } 24 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/TutorialLinAlgExComputeSolveError.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | MatrixXd A = MatrixXd::Random(100,100); 10 | MatrixXd b = MatrixXd::Random(100,50); 11 | MatrixXd x = A.fullPivLu().solve(b); 12 | double relative_error = (A*x - b).norm() / b.norm(); // norm() is L2 norm 13 | cout << "The relative error is:\n" << relative_error << endl; 14 | } 15 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/TutorialLinAlgExSolveColPivHouseholderQR.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | Matrix3f A; 10 | Vector3f b; 11 | A << 1,2,3, 4,5,6, 7,8,10; 12 | b << 3, 3, 4; 13 | cout << "Here is the matrix A:\n" << A << endl; 14 | cout << "Here is the vector b:\n" << b << endl; 15 | Vector3f x = A.colPivHouseholderQr().solve(b); 16 | cout << "The solution is:\n" << x << endl; 17 | } 18 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/TutorialLinAlgExSolveLDLT.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | Matrix2f A, b; 10 | A << 2, -1, -1, 3; 11 | b << 1, 2, 3, 1; 12 | cout << "Here is the matrix A:\n" << A << endl; 13 | cout << "Here is the right hand side b:\n" << b << endl; 14 | Matrix2f x = A.ldlt().solve(b); 15 | cout << "The solution is:\n" << x << endl; 16 | } 17 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/TutorialLinAlgInverseDeterminant.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | Matrix3f A; 10 | A << 1, 2, 1, 11 | 2, 1, 0, 12 | -1, 1, 2; 13 | cout << "Here is the matrix A:\n" << A << endl; 14 | cout << "The determinant of A is " << A.determinant() << endl; 15 | cout << "The inverse of A is:\n" << A.inverse() << endl; 16 | } -------------------------------------------------------------------------------- /eigen-patched/doc/examples/TutorialLinAlgRankRevealing.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | Matrix3f A; 10 | A << 1, 2, 5, 11 | 2, 1, 4, 12 | 3, 0, 3; 13 | cout << "Here is the matrix A:\n" << A << endl; 14 | FullPivLU lu_decomp(A); 15 | cout << "The rank of A is " << lu_decomp.rank() << endl; 16 | cout << "Here is a matrix whose columns form a basis of the null-space of A:\n" 17 | << lu_decomp.kernel() << endl; 18 | cout << "Here is a matrix whose columns form a basis of the column-space of A:\n" 19 | << lu_decomp.image(A) << endl; // yes, have to pass the original A 20 | } 21 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/TutorialLinAlgSVDSolve.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | MatrixXf A = MatrixXf::Random(3, 2); 10 | cout << "Here is the matrix A:\n" << A << endl; 11 | VectorXf b = VectorXf::Random(3); 12 | cout << "Here is the right hand side b:\n" << b << endl; 13 | cout << "The least-squares solution is:\n" 14 | << A.jacobiSvd(ComputeThinU | ComputeThinV).solve(b) << endl; 15 | } 16 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/TutorialLinAlgSelfAdjointEigenSolver.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | Matrix2f A; 10 | A << 1, 2, 2, 3; 11 | cout << "Here is the matrix A:\n" << A << endl; 12 | SelfAdjointEigenSolver eigensolver(A); 13 | if (eigensolver.info() != Success) abort(); 14 | cout << "The eigenvalues of A are:\n" << eigensolver.eigenvalues() << endl; 15 | cout << "Here's a matrix whose columns are eigenvectors of A \n" 16 | << "corresponding to these eigenvalues:\n" 17 | << eigensolver.eigenvectors() << endl; 18 | } 19 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/TutorialLinAlgSetThreshold.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | Matrix2d A; 10 | A << 2, 1, 11 | 2, 0.9999999999; 12 | FullPivLU lu(A); 13 | cout << "By default, the rank of A is found to be " << lu.rank() << endl; 14 | lu.setThreshold(1e-5); 15 | cout << "With threshold 1e-5, the rank of A is found to be " << lu.rank() << endl; 16 | } 17 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ArrayClass_accessors.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main() 8 | { 9 | ArrayXXf m(2,2); 10 | 11 | // assign some values coefficient by coefficient 12 | m(0,0) = 1.0; m(0,1) = 2.0; 13 | m(1,0) = 3.0; m(1,1) = m(0,1) + m(1,0); 14 | 15 | // print values to standard output 16 | cout << m << endl << endl; 17 | 18 | // using the comma-initializer is also allowed 19 | m << 1.0,2.0, 20 | 3.0,4.0; 21 | 22 | // print values to standard output 23 | cout << m << endl; 24 | } 25 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ArrayClass_addition.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main() 8 | { 9 | ArrayXXf a(3,3); 10 | ArrayXXf b(3,3); 11 | a << 1,2,3, 12 | 4,5,6, 13 | 7,8,9; 14 | b << 1,2,3, 15 | 1,2,3, 16 | 1,2,3; 17 | 18 | // Adding two arrays 19 | cout << "a + b = " << endl << a + b << endl << endl; 20 | 21 | // Subtracting a scalar from an array 22 | cout << "a - 2 = " << endl << a - 2 << endl; 23 | } 24 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ArrayClass_cwise_other.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main() 8 | { 9 | ArrayXf a = ArrayXf::Random(5); 10 | a *= 2; 11 | cout << "a =" << endl 12 | << a << endl; 13 | cout << "a.abs() =" << endl 14 | << a.abs() << endl; 15 | cout << "a.abs().sqrt() =" << endl 16 | << a.abs().sqrt() << endl; 17 | cout << "a.min(a.abs().sqrt()) =" << endl 18 | << a.min(a.abs().sqrt()) << endl; 19 | } 20 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ArrayClass_interop.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main() 8 | { 9 | MatrixXf m(2,2); 10 | MatrixXf n(2,2); 11 | MatrixXf result(2,2); 12 | 13 | m << 1,2, 14 | 3,4; 15 | n << 5,6, 16 | 7,8; 17 | 18 | result = (m.array() + 4).matrix() * m; 19 | cout << "-- Combination 1: --" << endl << result << endl << endl; 20 | result = (m.array() * n.array()).matrix() * m; 21 | cout << "-- Combination 2: --" << endl << result << endl << endl; 22 | } 23 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ArrayClass_interop_matrix.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main() 8 | { 9 | MatrixXf m(2,2); 10 | MatrixXf n(2,2); 11 | MatrixXf result(2,2); 12 | 13 | m << 1,2, 14 | 3,4; 15 | n << 5,6, 16 | 7,8; 17 | 18 | result = m * n; 19 | cout << "-- Matrix m*n: --" << endl << result << endl << endl; 20 | result = m.array() * n.array(); 21 | cout << "-- Array m*n: --" << endl << result << endl << endl; 22 | result = m.cwiseProduct(n); 23 | cout << "-- With cwiseProduct: --" << endl << result << endl << endl; 24 | result = m.array() + 4; 25 | cout << "-- Array m + 4: --" << endl << result << endl << endl; 26 | } 27 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ArrayClass_mult.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main() 8 | { 9 | ArrayXXf a(2,2); 10 | ArrayXXf b(2,2); 11 | a << 1,2, 12 | 3,4; 13 | b << 5,6, 14 | 7,8; 15 | cout << "a * b = " << endl << a * b << endl; 16 | } 17 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_BlockOperations_block_assignment.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | Array22f m; 10 | m << 1,2, 11 | 3,4; 12 | Array44f a = Array44f::Constant(0.6); 13 | cout << "Here is the array a:" << endl << a << endl << endl; 14 | a.block<2,2>(1,1) = m; 15 | cout << "Here is now a with m copied into its central 2x2 block:" << endl << a << endl << endl; 16 | a.block(0,0,2,3) = a.block(2,1,2,3); 17 | cout << "Here is now a with bottom-right 2x3 block copied into top-left 2x2 block:" << endl << a << endl << endl; 18 | } 19 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_BlockOperations_colrow.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | 6 | int main() 7 | { 8 | Eigen::MatrixXf m(3,3); 9 | m << 1,2,3, 10 | 4,5,6, 11 | 7,8,9; 12 | cout << "Here is the matrix m:" << endl << m << endl; 13 | cout << "2nd Row: " << m.row(1) << endl; 14 | m.col(2) += 3 * m.col(0); 15 | cout << "After adding 3 times the first column into the third column, the matrix m is:\n"; 16 | cout << m << endl; 17 | } 18 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_BlockOperations_corner.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | 6 | int main() 7 | { 8 | Eigen::Matrix4f m; 9 | m << 1, 2, 3, 4, 10 | 5, 6, 7, 8, 11 | 9, 10,11,12, 12 | 13,14,15,16; 13 | cout << "m.leftCols(2) =" << endl << m.leftCols(2) << endl << endl; 14 | cout << "m.bottomRows<2>() =" << endl << m.bottomRows<2>() << endl << endl; 15 | m.topLeftCorner(1,3) = m.bottomRightCorner(3,1).transpose(); 16 | cout << "After assignment, m = " << endl << m << endl; 17 | } 18 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_BlockOperations_print_block.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | 6 | int main() 7 | { 8 | Eigen::MatrixXf m(4,4); 9 | m << 1, 2, 3, 4, 10 | 5, 6, 7, 8, 11 | 9,10,11,12, 12 | 13,14,15,16; 13 | cout << "Block in the middle" << endl; 14 | cout << m.block<2,2>(1,1) << endl << endl; 15 | for (int i = 1; i <= 3; ++i) 16 | { 17 | cout << "Block of size " << i << "x" << i << endl; 18 | cout << m.block(0,0,i,i) << endl << endl; 19 | } 20 | } 21 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_BlockOperations_vector.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | 6 | int main() 7 | { 8 | Eigen::ArrayXf v(6); 9 | v << 1, 2, 3, 4, 5, 6; 10 | cout << "v.head(3) =" << endl << v.head(3) << endl << endl; 11 | cout << "v.tail<3>() = " << endl << v.tail<3>() << endl << endl; 12 | v.segment(1,4) *= 2; 13 | cout << "after 'v.segment(1,4) *= 2', v =" << endl << v << endl; 14 | } 15 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_PartialLU_solve.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | #include 4 | 5 | using namespace std; 6 | using namespace Eigen; 7 | 8 | int main() 9 | { 10 | Matrix3f A; 11 | Vector3f b; 12 | A << 1,2,3, 4,5,6, 7,8,10; 13 | b << 3, 3, 4; 14 | cout << "Here is the matrix A:" << endl << A << endl; 15 | cout << "Here is the vector b:" << endl << b << endl; 16 | Vector3f x = A.lu().solve(b); 17 | cout << "The solution is:" << endl << x << endl; 18 | } 19 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ReductionsVisitorsBroadcasting_broadcast_1nn.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | Eigen::MatrixXf m(2,4); 10 | Eigen::VectorXf v(2); 11 | 12 | m << 1, 23, 6, 9, 13 | 3, 11, 7, 2; 14 | 15 | v << 2, 16 | 3; 17 | 18 | MatrixXf::Index index; 19 | // find nearest neighbour 20 | (m.colwise() - v).colwise().squaredNorm().minCoeff(&index); 21 | 22 | cout << "Nearest neighbour is column " << index << ":" << endl; 23 | cout << m.col(index) << endl; 24 | } 25 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ReductionsVisitorsBroadcasting_broadcast_simple.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | int main() 6 | { 7 | Eigen::MatrixXf mat(2,4); 8 | Eigen::VectorXf v(2); 9 | 10 | mat << 1, 2, 6, 9, 11 | 3, 1, 7, 2; 12 | 13 | v << 0, 14 | 1; 15 | 16 | //add v to each column of m 17 | mat.colwise() += v; 18 | 19 | std::cout << "Broadcasting result: " << std::endl; 20 | std::cout << mat << std::endl; 21 | } 22 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ReductionsVisitorsBroadcasting_broadcast_simple_rowwise.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | int main() 6 | { 7 | Eigen::MatrixXf mat(2,4); 8 | Eigen::VectorXf v(4); 9 | 10 | mat << 1, 2, 6, 9, 11 | 3, 1, 7, 2; 12 | 13 | v << 0,1,2,3; 14 | 15 | //add v to each row of m 16 | mat.rowwise() += v.transpose(); 17 | 18 | std::cout << "Broadcasting result: " << std::endl; 19 | std::cout << mat << std::endl; 20 | } 21 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ReductionsVisitorsBroadcasting_colwise.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | int main() 6 | { 7 | Eigen::MatrixXf mat(2,4); 8 | mat << 1, 2, 6, 9, 9 | 3, 1, 7, 2; 10 | 11 | std::cout << "Column's maximum: " << std::endl 12 | << mat.colwise().maxCoeff() << std::endl; 13 | } 14 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ReductionsVisitorsBroadcasting_maxnorm.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | int main() 7 | { 8 | MatrixXf mat(2,4); 9 | mat << 1, 2, 6, 9, 10 | 3, 1, 7, 2; 11 | 12 | MatrixXf::Index maxIndex; 13 | float maxNorm = mat.colwise().sum().maxCoeff(&maxIndex); 14 | 15 | std::cout << "Maximum sum at position " << maxIndex << std::endl; 16 | 17 | std::cout << "The corresponding vector is: " << std::endl; 18 | std::cout << mat.col( maxIndex ) << std::endl; 19 | std::cout << "And its sum is is: " << maxNorm << std::endl; 20 | } 21 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ReductionsVisitorsBroadcasting_reductions_bool.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | ArrayXXf a(2,2); 10 | 11 | a << 1,2, 12 | 3,4; 13 | 14 | cout << "(a > 0).all() = " << (a > 0).all() << endl; 15 | cout << "(a > 0).any() = " << (a > 0).any() << endl; 16 | cout << "(a > 0).count() = " << (a > 0).count() << endl; 17 | cout << endl; 18 | cout << "(a > 2).all() = " << (a > 2).all() << endl; 19 | cout << "(a > 2).any() = " << (a > 2).any() << endl; 20 | cout << "(a > 2).count() = " << (a > 2).count() << endl; 21 | } 22 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ReductionsVisitorsBroadcasting_reductions_norm.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | VectorXf v(2); 10 | MatrixXf m(2,2), n(2,2); 11 | 12 | v << -1, 13 | 2; 14 | 15 | m << 1,-2, 16 | -3,4; 17 | 18 | cout << "v.squaredNorm() = " << v.squaredNorm() << endl; 19 | cout << "v.norm() = " << v.norm() << endl; 20 | cout << "v.lpNorm<1>() = " << v.lpNorm<1>() << endl; 21 | cout << "v.lpNorm() = " << v.lpNorm() << endl; 22 | 23 | cout << endl; 24 | cout << "m.squaredNorm() = " << m.squaredNorm() << endl; 25 | cout << "m.norm() = " << m.norm() << endl; 26 | cout << "m.lpNorm<1>() = " << m.lpNorm<1>() << endl; 27 | cout << "m.lpNorm() = " << m.lpNorm() << endl; 28 | } 29 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ReductionsVisitorsBroadcasting_rowwise.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | int main() 6 | { 7 | Eigen::MatrixXf mat(2,4); 8 | mat << 1, 2, 6, 9, 9 | 3, 1, 7, 2; 10 | 11 | std::cout << "Row's maximum: " << std::endl 12 | << mat.rowwise().maxCoeff() << std::endl; 13 | } 14 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_ReductionsVisitorsBroadcasting_visitors.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | using namespace Eigen; 6 | 7 | int main() 8 | { 9 | Eigen::MatrixXf m(2,2); 10 | 11 | m << 1, 2, 12 | 3, 4; 13 | 14 | //get location of maximum 15 | MatrixXf::Index maxRow, maxCol; 16 | float max = m.maxCoeff(&maxRow, &maxCol); 17 | 18 | //get location of minimum 19 | MatrixXf::Index minRow, minCol; 20 | float min = m.minCoeff(&minRow, &minCol); 21 | 22 | cout << "Max: " << max << ", at: " << 23 | maxRow << "," << maxCol << endl; 24 | cout << "Min: " << min << ", at: " << 25 | minRow << "," << minCol << endl; 26 | } 27 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/Tutorial_simple_example_dynamic_size.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | int main() 7 | { 8 | for (int size=1; size<=4; ++size) 9 | { 10 | MatrixXi m(size,size+1); // a (size)x(size+1)-matrix of int's 11 | for (int j=0; j 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | int main() 7 | { 8 | Matrix3f m3; 9 | m3 << 1, 2, 3, 4, 5, 6, 7, 8, 9; 10 | Matrix4f m4 = Matrix4f::Identity(); 11 | Vector4i v4(1, 2, 3, 4); 12 | 13 | std::cout << "m3\n" << m3 << "\nm4:\n" 14 | << m4 << "\nv4:\n" << v4 << std::endl; 15 | } 16 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/class_CwiseBinaryOp.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | using namespace Eigen; 4 | using namespace std; 5 | 6 | // define a custom template binary functor 7 | template struct MakeComplexOp { 8 | EIGEN_EMPTY_STRUCT_CTOR(MakeComplexOp) 9 | typedef complex result_type; 10 | complex operator()(const Scalar& a, const Scalar& b) const { return complex(a,b); } 11 | }; 12 | 13 | int main(int, char**) 14 | { 15 | Matrix4d m1 = Matrix4d::Random(), m2 = Matrix4d::Random(); 16 | cout << m1.binaryExpr(m2, MakeComplexOp()) << endl; 17 | return 0; 18 | } 19 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/class_CwiseUnaryOp.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | using namespace Eigen; 4 | using namespace std; 5 | 6 | // define a custom template unary functor 7 | template 8 | struct CwiseClampOp { 9 | CwiseClampOp(const Scalar& inf, const Scalar& sup) : m_inf(inf), m_sup(sup) {} 10 | const Scalar operator()(const Scalar& x) const { return xm_sup ? m_sup : x); } 11 | Scalar m_inf, m_sup; 12 | }; 13 | 14 | int main(int, char**) 15 | { 16 | Matrix4d m1 = Matrix4d::Random(); 17 | cout << m1 << endl << "becomes: " << endl << m1.unaryExpr(CwiseClampOp(-0.5,0.5)) << endl; 18 | return 0; 19 | } 20 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/class_CwiseUnaryOp_ptrfun.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | using namespace Eigen; 4 | using namespace std; 5 | 6 | // define function to be applied coefficient-wise 7 | double ramp(double x) 8 | { 9 | if (x > 0) 10 | return x; 11 | else 12 | return 0; 13 | } 14 | 15 | int main(int, char**) 16 | { 17 | Matrix4d m1 = Matrix4d::Random(); 18 | cout << m1 << endl << "becomes: " << endl << m1.unaryExpr(ptr_fun(ramp)) << endl; 19 | return 0; 20 | } 21 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/class_FixedBlock.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | using namespace Eigen; 4 | using namespace std; 5 | 6 | template 7 | Eigen::Block 8 | topLeft2x2Corner(MatrixBase& m) 9 | { 10 | return Eigen::Block(m.derived(), 0, 0); 11 | } 12 | 13 | template 14 | const Eigen::Block 15 | topLeft2x2Corner(const MatrixBase& m) 16 | { 17 | return Eigen::Block(m.derived(), 0, 0); 18 | } 19 | 20 | int main(int, char**) 21 | { 22 | Matrix3d m = Matrix3d::Identity(); 23 | cout << topLeft2x2Corner(4*m) << endl; // calls the const version 24 | topLeft2x2Corner(m) *= 2; // calls the non-const version 25 | cout << "Now the matrix m is:" << endl << m << endl; 26 | return 0; 27 | } 28 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/class_FixedVectorBlock.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | using namespace Eigen; 4 | using namespace std; 5 | 6 | template 7 | Eigen::VectorBlock 8 | firstTwo(MatrixBase& v) 9 | { 10 | return Eigen::VectorBlock(v.derived(), 0); 11 | } 12 | 13 | template 14 | const Eigen::VectorBlock 15 | firstTwo(const MatrixBase& v) 16 | { 17 | return Eigen::VectorBlock(v.derived(), 0); 18 | } 19 | 20 | int main(int, char**) 21 | { 22 | Matrix v; v << 1,2,3,4,5,6; 23 | cout << firstTwo(4*v) << endl; // calls the const version 24 | firstTwo(v) *= 2; // calls the non-const version 25 | cout << "Now the vector v is:" << endl << v << endl; 26 | return 0; 27 | } 28 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/function_taking_eigenbase.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | using namespace Eigen; 4 | 5 | template 6 | void print_size(const EigenBase& b) 7 | { 8 | std::cout << "size (rows, cols): " << b.size() << " (" << b.rows() 9 | << ", " << b.cols() << ")" << std::endl; 10 | } 11 | 12 | int main() 13 | { 14 | Vector3f v; 15 | print_size(v); 16 | // v.asDiagonal() returns a 3x3 diagonal matrix pseudo-expression 17 | print_size(v.asDiagonal()); 18 | } 19 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/tut_arithmetic_add_sub.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | int main() 7 | { 8 | Matrix2d a; 9 | a << 1, 2, 10 | 3, 4; 11 | MatrixXd b(2,2); 12 | b << 2, 3, 13 | 1, 4; 14 | std::cout << "a + b =\n" << a + b << std::endl; 15 | std::cout << "a - b =\n" << a - b << std::endl; 16 | std::cout << "Doing a += b;" << std::endl; 17 | a += b; 18 | std::cout << "Now a =\n" << a << std::endl; 19 | Vector3d v(1,2,3); 20 | Vector3d w(1,0,0); 21 | std::cout << "-v + w - v =\n" << -v + w - v << std::endl; 22 | } 23 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/tut_arithmetic_dot_cross.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | int main() 7 | { 8 | Vector3d v(1,2,3); 9 | Vector3d w(0,1,2); 10 | 11 | cout << "Dot product: " << v.dot(w) << endl; 12 | double dp = v.adjoint()*w; // automatic conversion of the inner product to a scalar 13 | cout << "Dot product via a matrix product: " << dp << endl; 14 | cout << "Cross product:\n" << v.cross(w) << endl; 15 | } 16 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/tut_arithmetic_matrix_mul.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | int main() 6 | { 7 | Matrix2d mat; 8 | mat << 1, 2, 9 | 3, 4; 10 | Vector2d u(-1,1), v(2,0); 11 | std::cout << "Here is mat*mat:\n" << mat*mat << std::endl; 12 | std::cout << "Here is mat*u:\n" << mat*u << std::endl; 13 | std::cout << "Here is u^T*mat:\n" << u.transpose()*mat << std::endl; 14 | std::cout << "Here is u^T*v:\n" << u.transpose()*v << std::endl; 15 | std::cout << "Here is u*v^T:\n" << u*v.transpose() << std::endl; 16 | std::cout << "Let's multiply mat by itself" << std::endl; 17 | mat = mat*mat; 18 | std::cout << "Now mat is mat:\n" << mat << std::endl; 19 | } 20 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/tut_arithmetic_redux_basic.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace std; 5 | int main() 6 | { 7 | Eigen::Matrix2d mat; 8 | mat << 1, 2, 9 | 3, 4; 10 | cout << "Here is mat.sum(): " << mat.sum() << endl; 11 | cout << "Here is mat.prod(): " << mat.prod() << endl; 12 | cout << "Here is mat.mean(): " << mat.mean() << endl; 13 | cout << "Here is mat.minCoeff(): " << mat.minCoeff() << endl; 14 | cout << "Here is mat.maxCoeff(): " << mat.maxCoeff() << endl; 15 | cout << "Here is mat.trace(): " << mat.trace() << endl; 16 | } 17 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/tut_arithmetic_scalar_mul_div.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | int main() 7 | { 8 | Matrix2d a; 9 | a << 1, 2, 10 | 3, 4; 11 | Vector3d v(1,2,3); 12 | std::cout << "a * 2.5 =\n" << a * 2.5 << std::endl; 13 | std::cout << "0.1 * v =\n" << 0.1 * v << std::endl; 14 | std::cout << "Doing v *= 2;" << std::endl; 15 | v *= 2; 16 | std::cout << "Now v =\n" << v << std::endl; 17 | } 18 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/tut_matrix_coefficient_accessors.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | int main() 7 | { 8 | MatrixXd m(2,2); 9 | m(0,0) = 3; 10 | m(1,0) = 2.5; 11 | m(0,1) = -1; 12 | m(1,1) = m(1,0) + m(0,1); 13 | std::cout << "Here is the matrix m:\n" << m << std::endl; 14 | VectorXd v(2); 15 | v(0) = 4; 16 | v(1) = v(0) - 1; 17 | std::cout << "Here is the vector v:\n" << v << std::endl; 18 | } 19 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/tut_matrix_resize.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | int main() 7 | { 8 | MatrixXd m(2,5); 9 | m.resize(4,3); 10 | std::cout << "The matrix m is of size " 11 | << m.rows() << "x" << m.cols() << std::endl; 12 | std::cout << "It has " << m.size() << " coefficients" << std::endl; 13 | VectorXd v(2); 14 | v.resize(5); 15 | std::cout << "The vector v is of size " << v.size() << std::endl; 16 | std::cout << "As a matrix, v is of size " 17 | << v.rows() << "x" << v.cols() << std::endl; 18 | } 19 | -------------------------------------------------------------------------------- /eigen-patched/doc/examples/tut_matrix_resize_fixed_size.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | int main() 7 | { 8 | Matrix4d m; 9 | m.resize(4,4); // no operation 10 | std::cout << "The matrix m is of size " 11 | << m.rows() << "x" << m.cols() << std::endl; 12 | } 13 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/.krazy: -------------------------------------------------------------------------------- 1 | EXCLUDE copyright 2 | EXCLUDE license 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/AngleAxis_mimic_euler.cpp: -------------------------------------------------------------------------------- 1 | Matrix3f m; 2 | m = AngleAxisf(0.25*M_PI, Vector3f::UnitX()) 3 | * AngleAxisf(0.5*M_PI, Vector3f::UnitY()) 4 | * AngleAxisf(0.33*M_PI, Vector3f::UnitZ()); 5 | cout << m << endl << "is unitary: " << m.isUnitary() << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/ColPivHouseholderQR_solve.cpp: -------------------------------------------------------------------------------- 1 | Matrix3f m = Matrix3f::Random(); 2 | Matrix3f y = Matrix3f::Random(); 3 | cout << "Here is the matrix m:" << endl << m << endl; 4 | cout << "Here is the matrix y:" << endl << y << endl; 5 | Matrix3f x; 6 | x = m.colPivHouseholderQr().solve(y); 7 | assert(y.isApprox(m*x)); 8 | cout << "Here is a solution x to the equation mx=y:" << endl << x << endl; 9 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/ComplexEigenSolver_eigenvalues.cpp: -------------------------------------------------------------------------------- 1 | MatrixXcf ones = MatrixXcf::Ones(3,3); 2 | ComplexEigenSolver ces(ones, /* computeEigenvectors = */ false); 3 | cout << "The eigenvalues of the 3x3 matrix of ones are:" 4 | << endl << ces.eigenvalues() << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/ComplexEigenSolver_eigenvectors.cpp: -------------------------------------------------------------------------------- 1 | MatrixXcf ones = MatrixXcf::Ones(3,3); 2 | ComplexEigenSolver ces(ones); 3 | cout << "The first eigenvector of the 3x3 matrix of ones is:" 4 | << endl << ces.eigenvectors().col(1) << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/ComplexSchur_compute.cpp: -------------------------------------------------------------------------------- 1 | MatrixXcf A = MatrixXcf::Random(4,4); 2 | ComplexSchur schur(4); 3 | schur.compute(A); 4 | cout << "The matrix T in the decomposition of A is:" << endl << schur.matrixT() << endl; 5 | schur.compute(A.inverse()); 6 | cout << "The matrix T in the decomposition of A^(-1) is:" << endl << schur.matrixT() << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/ComplexSchur_matrixT.cpp: -------------------------------------------------------------------------------- 1 | MatrixXcf A = MatrixXcf::Random(4,4); 2 | cout << "Here is a random 4x4 matrix, A:" << endl << A << endl << endl; 3 | ComplexSchur schurOfA(A, false); // false means do not compute U 4 | cout << "The triangular matrix T is:" << endl << schurOfA.matrixT() << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/ComplexSchur_matrixU.cpp: -------------------------------------------------------------------------------- 1 | MatrixXcf A = MatrixXcf::Random(4,4); 2 | cout << "Here is a random 4x4 matrix, A:" << endl << A << endl << endl; 3 | ComplexSchur schurOfA(A); 4 | cout << "The unitary matrix U is:" << endl << schurOfA.matrixU() << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Cwise_abs.cpp: -------------------------------------------------------------------------------- 1 | Array3d v(1,-2,-3); 2 | cout << v.abs() << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Cwise_abs2.cpp: -------------------------------------------------------------------------------- 1 | Array3d v(1,-2,-3); 2 | cout << v.abs2() << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Cwise_acos.cpp: -------------------------------------------------------------------------------- 1 | Array3d v(0, sqrt(2.)/2, 1); 2 | cout << v.acos() << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Cwise_boolean_and.cpp: -------------------------------------------------------------------------------- 1 | Array3d v(-1,2,1), w(-3,2,3); 2 | cout << ((vw) << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Cwise_greater_equal.cpp: -------------------------------------------------------------------------------- 1 | Array3d v(1,2,3), w(3,2,1); 2 | cout << (v>=w) << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Cwise_inverse.cpp: -------------------------------------------------------------------------------- 1 | Array3d v(2,3,4); 2 | cout << v.inverse() << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Cwise_less.cpp: -------------------------------------------------------------------------------- 1 | Array3d v(1,2,3), w(3,2,1); 2 | cout << (v() = ..." << endl; 4 | cout << m.colwise().replicate<3>() << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/DirectionWise_replicate_int.cpp: -------------------------------------------------------------------------------- 1 | Vector3i v = Vector3i::Random(); 2 | cout << "Here is the vector v:" << endl << v << endl; 3 | cout << "v.rowwise().replicate(5) = ..." << endl; 4 | cout << v.rowwise().replicate(5) << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/EigenSolver_compute.cpp: -------------------------------------------------------------------------------- 1 | EigenSolver es; 2 | MatrixXf A = MatrixXf::Random(4,4); 3 | es.compute(A, /* computeEigenvectors = */ false); 4 | cout << "The eigenvalues of A are: " << es.eigenvalues().transpose() << endl; 5 | es.compute(A + MatrixXf::Identity(4,4), false); // re-use es to compute eigenvalues of A+I 6 | cout << "The eigenvalues of A+I are: " << es.eigenvalues().transpose() << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/EigenSolver_eigenvalues.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd ones = MatrixXd::Ones(3,3); 2 | EigenSolver es(ones, false); 3 | cout << "The eigenvalues of the 3x3 matrix of ones are:" 4 | << endl << es.eigenvalues() << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/EigenSolver_eigenvectors.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd ones = MatrixXd::Ones(3,3); 2 | EigenSolver es(ones); 3 | cout << "The first eigenvector of the 3x3 matrix of ones is:" 4 | << endl << es.eigenvectors().col(1) << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/EigenSolver_pseudoEigenvectors.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd A = MatrixXd::Random(6,6); 2 | cout << "Here is a random 6x6 matrix, A:" << endl << A << endl << endl; 3 | 4 | EigenSolver es(A); 5 | MatrixXd D = es.pseudoEigenvalueMatrix(); 6 | MatrixXd V = es.pseudoEigenvectors(); 7 | cout << "The pseudo-eigenvalue matrix D is:" << endl << D << endl; 8 | cout << "The pseudo-eigenvector matrix V is:" << endl << V << endl; 9 | cout << "Finally, V * D * V^(-1) = " << endl << V * D * V.inverse() << endl; 10 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/FullPivHouseholderQR_solve.cpp: -------------------------------------------------------------------------------- 1 | Matrix3f m = Matrix3f::Random(); 2 | Matrix3f y = Matrix3f::Random(); 3 | cout << "Here is the matrix m:" << endl << m << endl; 4 | cout << "Here is the matrix y:" << endl << y << endl; 5 | Matrix3f x; 6 | x = m.fullPivHouseholderQr().solve(y); 7 | assert(y.isApprox(m*x)); 8 | cout << "Here is a solution x to the equation mx=y:" << endl << x << endl; 9 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/FullPivLU_image.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m; 2 | m << 1,1,0, 3 | 1,3,2, 4 | 0,1,1; 5 | cout << "Here is the matrix m:" << endl << m << endl; 6 | cout << "Notice that the middle column is the sum of the two others, so the " 7 | << "columns are linearly dependent." << endl; 8 | cout << "Here is a matrix whose columns have the same span but are linearly independent:" 9 | << endl << m.fullPivLu().image(m) << endl; 10 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/FullPivLU_kernel.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf m = MatrixXf::Random(3,5); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | MatrixXf ker = m.fullPivLu().kernel(); 4 | cout << "Here is a matrix whose columns form a basis of the kernel of m:" 5 | << endl << ker << endl; 6 | cout << "By definition of the kernel, m*ker is zero:" 7 | << endl << m*ker << endl; 8 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/FullPivLU_solve.cpp: -------------------------------------------------------------------------------- 1 | Matrix m = Matrix::Random(); 2 | Matrix2f y = Matrix2f::Random(); 3 | cout << "Here is the matrix m:" << endl << m << endl; 4 | cout << "Here is the matrix y:" << endl << y << endl; 5 | Matrix x = m.fullPivLu().solve(y); 6 | if((m*x).isApprox(y)) 7 | { 8 | cout << "Here is a solution x to the equation mx=y:" << endl << x << endl; 9 | } 10 | else 11 | cout << "The equation mx=y does not have any solution." << endl; 12 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/HessenbergDecomposition_compute.cpp: -------------------------------------------------------------------------------- 1 | MatrixXcf A = MatrixXcf::Random(4,4); 2 | HessenbergDecomposition hd(4); 3 | hd.compute(A); 4 | cout << "The matrix H in the decomposition of A is:" << endl << hd.matrixH() << endl; 5 | hd.compute(2*A); // re-use hd to compute and store decomposition of 2A 6 | cout << "The matrix H in the decomposition of 2A is:" << endl << hd.matrixH() << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/HessenbergDecomposition_matrixH.cpp: -------------------------------------------------------------------------------- 1 | Matrix4f A = MatrixXf::Random(4,4); 2 | cout << "Here is a random 4x4 matrix:" << endl << A << endl; 3 | HessenbergDecomposition hessOfA(A); 4 | MatrixXf H = hessOfA.matrixH(); 5 | cout << "The Hessenberg matrix H is:" << endl << H << endl; 6 | MatrixXf Q = hessOfA.matrixQ(); 7 | cout << "The orthogonal matrix Q is:" << endl << Q << endl; 8 | cout << "Q H Q^T is:" << endl << Q * H * Q.transpose() << endl; 9 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/HessenbergDecomposition_packedMatrix.cpp: -------------------------------------------------------------------------------- 1 | Matrix4d A = Matrix4d::Random(4,4); 2 | cout << "Here is a random 4x4 matrix:" << endl << A << endl; 3 | HessenbergDecomposition hessOfA(A); 4 | Matrix4d pm = hessOfA.packedMatrix(); 5 | cout << "The packed matrix M is:" << endl << pm << endl; 6 | cout << "The upper Hessenberg part corresponds to the matrix H, which is:" 7 | << endl << hessOfA.matrixH() << endl; 8 | Vector3d hc = hessOfA.householderCoefficients(); 9 | cout << "The vector of Householder coefficients is:" << endl << hc << endl; 10 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/HouseholderQR_solve.cpp: -------------------------------------------------------------------------------- 1 | typedef Matrix Matrix3x3; 2 | Matrix3x3 m = Matrix3x3::Random(); 3 | Matrix3f y = Matrix3f::Random(); 4 | cout << "Here is the matrix m:" << endl << m << endl; 5 | cout << "Here is the matrix y:" << endl << y << endl; 6 | Matrix3f x; 7 | x = m.householderQr().solve(y); 8 | assert(y.isApprox(m*x)); 9 | cout << "Here is a solution x to the equation mx=y:" << endl << x << endl; 10 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/IOFormat.cpp: -------------------------------------------------------------------------------- 1 | std::string sep = "\n----------------------------------------\n"; 2 | Matrix3d m1; 3 | m1 << 1.111111, 2, 3.33333, 4, 5, 6, 7, 8.888888, 9; 4 | 5 | IOFormat CommaInitFmt(StreamPrecision, DontAlignCols, ", ", ", ", "", "", " << ", ";"); 6 | IOFormat CleanFmt(4, 0, ", ", "\n", "[", "]"); 7 | IOFormat OctaveFmt(StreamPrecision, 0, ", ", ";\n", "", "", "[", "]"); 8 | IOFormat HeavyFmt(FullPrecision, 0, ", ", ";\n", "[", "]", "[", "]"); 9 | 10 | std::cout << m1 << sep; 11 | std::cout << m1.format(CommaInitFmt) << sep; 12 | std::cout << m1.format(CleanFmt) << sep; 13 | std::cout << m1.format(OctaveFmt) << sep; 14 | std::cout << m1.format(HeavyFmt) << sep; 15 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/JacobiSVD_basic.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf m = MatrixXf::Random(3,2); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | JacobiSVD svd(m, ComputeThinU | ComputeThinV); 4 | cout << "Its singular values are:" << endl << svd.singularValues() << endl; 5 | cout << "Its left singular vectors are the columns of the thin U matrix:" << endl << svd.matrixU() << endl; 6 | cout << "Its right singular vectors are the columns of the thin V matrix:" << endl << svd.matrixV() << endl; 7 | Vector3f rhs(1, 0, 0); 8 | cout << "Now consider this rhs vector:" << endl << rhs << endl; 9 | cout << "A least-squares solution of m*x = rhs is:" << endl << svd.solve(rhs) << endl; 10 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Jacobi_makeGivens.cpp: -------------------------------------------------------------------------------- 1 | Vector2f v = Vector2f::Random(); 2 | JacobiRotation G; 3 | G.makeGivens(v.x(), v.y()); 4 | cout << "Here is the vector v:" << endl << v << endl; 5 | v.applyOnTheLeft(0, 1, G.adjoint()); 6 | cout << "Here is the vector J' * v:" << endl << v << endl; -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Jacobi_makeJacobi.cpp: -------------------------------------------------------------------------------- 1 | Matrix2f m = Matrix2f::Random(); 2 | m = (m + m.adjoint()).eval(); 3 | JacobiRotation J; 4 | J.makeJacobi(m, 0, 1); 5 | cout << "Here is the matrix m:" << endl << m << endl; 6 | m.applyOnTheLeft(0, 1, J.adjoint()); 7 | m.applyOnTheRight(0, 1, J); 8 | cout << "Here is the matrix J' * m * J:" << endl << m << endl; -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/LLT_example.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd A(3,3); 2 | A << 4,-1,2, -1,6,0, 2,0,5; 3 | cout << "The matrix A is" << endl << A << endl; 4 | 5 | LLT lltOfA(A); // compute the Cholesky decomposition of A 6 | MatrixXd L = lltOfA.matrixL(); // retrieve factor L in the decomposition 7 | // The previous two lines can also be written as "L = A.llt().matrixL()" 8 | 9 | cout << "The Cholesky factor L is" << endl << L << endl; 10 | cout << "To check this, let us compute L * L.transpose()" << endl; 11 | cout << L * L.transpose() << endl; 12 | cout << "This should equal the matrix A" << endl; 13 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/LLT_solve.cpp: -------------------------------------------------------------------------------- 1 | typedef Matrix DataMatrix; 2 | // let's generate some samples on the 3D plane of equation z = 2x+3y (with some noise) 3 | DataMatrix samples = DataMatrix::Random(12,2); 4 | VectorXf elevations = 2*samples.col(0) + 3*samples.col(1) + VectorXf::Random(12)*0.1; 5 | // and let's solve samples * [x y]^T = elevations in least square sense: 6 | Matrix xy 7 | = (samples.adjoint() * samples).llt().solve((samples.adjoint()*elevations)); 8 | cout << xy << endl; 9 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Map_general_stride.cpp: -------------------------------------------------------------------------------- 1 | int array[24]; 2 | for(int i = 0; i < 24; ++i) array[i] = i; 3 | cout << Map > 4 | (array, 3, 3, Stride(8, 2)) 5 | << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Map_inner_stride.cpp: -------------------------------------------------------------------------------- 1 | int array[12]; 2 | for(int i = 0; i < 12; ++i) array[i] = i; 3 | cout << Map > 4 | (array, 6) // the inner stride has already been passed as template parameter 5 | << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Map_outer_stride.cpp: -------------------------------------------------------------------------------- 1 | int array[12]; 2 | for(int i = 0; i < 12; ++i) array[i] = i; 3 | cout << Map >(array, 3, 3, OuterStride<>(4)) << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Map_placement_new.cpp: -------------------------------------------------------------------------------- 1 | int data[] = {1,2,3,4,5,6,7,8,9}; 2 | Map v(data,4); 3 | cout << "The mapped vector v is: " << v << "\n"; 4 | new (&v) Map(data+4,5); 5 | cout << "Now v is: " << v << "\n"; -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Map_simple.cpp: -------------------------------------------------------------------------------- 1 | int array[9]; 2 | for(int i = 0; i < 9; ++i) array[i] = i; 3 | cout << Map(array) << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_adjoint.cpp: -------------------------------------------------------------------------------- 1 | Matrix2cf m = Matrix2cf::Random(); 2 | cout << "Here is the 2x2 complex matrix m:" << endl << m << endl; 3 | cout << "Here is the adjoint of m:" << endl << m.adjoint() << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_all.cpp: -------------------------------------------------------------------------------- 1 | Vector3f boxMin(Vector3f::Zero()), boxMax(Vector3f::Ones()); 2 | Vector3f p0 = Vector3f::Random(), p1 = Vector3f::Random().cwiseAbs(); 3 | // let's check if p0 and p1 are inside the axis aligned box defined by the corners boxMin,boxMax: 4 | cout << "Is (" << p0.transpose() << ") inside the box: " 5 | << ((boxMin.array()p0.array()).all()) << endl; 6 | cout << "Is (" << p1.transpose() << ") inside the box: " 7 | << ((boxMin.array()p1.array()).all()) << endl; 8 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_array.cpp: -------------------------------------------------------------------------------- 1 | Vector3d v(1,2,3); 2 | v.array() += 3; 3 | v.array() -= 2; 4 | cout << v << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_array_const.cpp: -------------------------------------------------------------------------------- 1 | Vector3d v(-1,2,-3); 2 | cout << "the absolute values:" << endl << v.array().abs() << endl; 3 | cout << "the absolute values plus one:" << endl << v.array().abs()+1 << endl; 4 | cout << "sum of the squares: " << v.array().square().sum() << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_asDiagonal.cpp: -------------------------------------------------------------------------------- 1 | cout << Matrix3i(Vector3i(2,5,6).asDiagonal()) << endl; 2 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_block_int_int.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is m.block<2,2>(1,1):" << endl << m.block<2,2>(1,1) << endl; 4 | m.block<2,2>(1,1).setZero(); 5 | cout << "Now the matrix m is:" << endl << m << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_block_int_int_int_int.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is m.block(1, 1, 2, 2):" << endl << m.block(1, 1, 2, 2) << endl; 4 | m.block(1, 1, 2, 2).setZero(); 5 | cout << "Now the matrix m is:" << endl << m << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_bottomLeftCorner_int_int.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is m.bottomLeftCorner(2, 2):" << endl; 4 | cout << m.bottomLeftCorner(2, 2) << endl; 5 | m.bottomLeftCorner(2, 2).setZero(); 6 | cout << "Now the matrix m is:" << endl << m << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_bottomRightCorner_int_int.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is m.bottomRightCorner(2, 2):" << endl; 4 | cout << m.bottomRightCorner(2, 2) << endl; 5 | m.bottomRightCorner(2, 2).setZero(); 6 | cout << "Now the matrix m is:" << endl << m << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_bottomRows_int.cpp: -------------------------------------------------------------------------------- 1 | Array44i a = Array44i::Random(); 2 | cout << "Here is the array a:" << endl << a << endl; 3 | cout << "Here is a.bottomRows(2):" << endl; 4 | cout << a.bottomRows(2) << endl; 5 | a.bottomRows(2).setZero(); 6 | cout << "Now the array a is:" << endl << a << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_cast.cpp: -------------------------------------------------------------------------------- 1 | Matrix2d md = Matrix2d::Identity() * 0.45; 2 | Matrix2f mf = Matrix2f::Identity(); 3 | cout << md + mf.cast() << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_col.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Identity(); 2 | m.col(1) = Vector3d(4,5,6); 3 | cout << m << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_colwise.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is the sum of each column:" << endl << m.colwise().sum() << endl; 4 | cout << "Here is the maximum absolute value of each column:" 5 | << endl << m.cwiseAbs().colwise().maxCoeff() << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_computeInverseAndDetWithCheck.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | Matrix3d inverse; 4 | bool invertible; 5 | double determinant; 6 | m.computeInverseAndDetWithCheck(inverse,determinant,invertible); 7 | cout << "Its determinant is " << determinant << endl; 8 | if(invertible) { 9 | cout << "It is invertible, and its inverse is:" << endl << inverse << endl; 10 | } 11 | else { 12 | cout << "It is not invertible." << endl; 13 | } 14 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_computeInverseWithCheck.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | Matrix3d inverse; 4 | bool invertible; 5 | m.computeInverseWithCheck(inverse,invertible); 6 | if(invertible) { 7 | cout << "It is invertible, and its inverse is:" << endl << inverse << endl; 8 | } 9 | else { 10 | cout << "It is not invertible." << endl; 11 | } 12 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_cwiseAbs.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd m(2,3); 2 | m << 2, -4, 6, 3 | -5, 1, 0; 4 | cout << m.cwiseAbs() << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_cwiseAbs2.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd m(2,3); 2 | m << 2, -4, 6, 3 | -5, 1, 0; 4 | cout << m.cwiseAbs2() << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_cwiseEqual.cpp: -------------------------------------------------------------------------------- 1 | MatrixXi m(2,2); 2 | m << 1, 0, 3 | 1, 1; 4 | cout << "Comparing m with identity matrix:" << endl; 5 | cout << m.cwiseEqual(MatrixXi::Identity(2,2)) << endl; 6 | int count = m.cwiseEqual(MatrixXi::Identity(2,2)).count(); 7 | cout << "Number of coefficients that are equal: " << count << endl; 8 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_cwiseInverse.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd m(2,3); 2 | m << 2, 0.5, 1, 3 | 3, 0.25, 1; 4 | cout << m.cwiseInverse() << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_cwiseMax.cpp: -------------------------------------------------------------------------------- 1 | Vector3d v(2,3,4), w(4,2,3); 2 | cout << v.cwiseMax(w) << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_cwiseMin.cpp: -------------------------------------------------------------------------------- 1 | Vector3d v(2,3,4), w(4,2,3); 2 | cout << v.cwiseMin(w) << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_cwiseNotEqual.cpp: -------------------------------------------------------------------------------- 1 | MatrixXi m(2,2); 2 | m << 1, 0, 3 | 1, 1; 4 | cout << "Comparing m with identity matrix:" << endl; 5 | cout << m.cwiseNotEqual(MatrixXi::Identity(2,2)) << endl; 6 | int count = m.cwiseNotEqual(MatrixXi::Identity(2,2)).count(); 7 | cout << "Number of coefficients that are not equal: " << count << endl; 8 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_cwiseProduct.cpp: -------------------------------------------------------------------------------- 1 | Matrix3i a = Matrix3i::Random(), b = Matrix3i::Random(); 2 | Matrix3i c = a.cwiseProduct(b); 3 | cout << "a:\n" << a << "\nb:\n" << b << "\nc:\n" << c << endl; 4 | 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_cwiseQuotient.cpp: -------------------------------------------------------------------------------- 1 | Vector3d v(2,3,4), w(4,2,3); 2 | cout << v.cwiseQuotient(w) << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_cwiseSqrt.cpp: -------------------------------------------------------------------------------- 1 | Vector3d v(1,2,4); 2 | cout << v.cwiseSqrt() << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_diagonal.cpp: -------------------------------------------------------------------------------- 1 | Matrix3i m = Matrix3i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here are the coefficients on the main diagonal of m:" << endl 4 | << m.diagonal() << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_diagonal_int.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here are the coefficients on the 1st super-diagonal and 2nd sub-diagonal of m:" << endl 4 | << m.diagonal(1).transpose() << endl 5 | << m.diagonal(-2).transpose() << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_diagonal_template_int.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here are the coefficients on the 1st super-diagonal and 2nd sub-diagonal of m:" << endl 4 | << m.diagonal<1>().transpose() << endl 5 | << m.diagonal<-2>().transpose() << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_eigenvalues.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd ones = MatrixXd::Ones(3,3); 2 | VectorXcd eivals = ones.eigenvalues(); 3 | cout << "The eigenvalues of the 3x3 matrix of ones are:" << endl << eivals << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_end_int.cpp: -------------------------------------------------------------------------------- 1 | RowVector4i v = RowVector4i::Random(); 2 | cout << "Here is the vector v:" << endl << v << endl; 3 | cout << "Here is v.tail(2):" << endl << v.tail(2) << endl; 4 | v.tail(2).setZero(); 5 | cout << "Now the vector v is:" << endl << v << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_eval.cpp: -------------------------------------------------------------------------------- 1 | Matrix2f M = Matrix2f::Random(); 2 | Matrix2f m; 3 | m = M; 4 | cout << "Here is the matrix m:" << endl << m << endl; 5 | cout << "Now we want to copy a column into a row." << endl; 6 | cout << "If we do m.col(1) = m.row(0), then m becomes:" << endl; 7 | m.col(1) = m.row(0); 8 | cout << m << endl << "which is wrong!" << endl; 9 | cout << "Now let us instead do m.col(1) = m.row(0).eval(). Then m becomes" << endl; 10 | m = M; 11 | m.col(1) = m.row(0).eval(); 12 | cout << m << endl << "which is right." << endl; 13 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_extract.cpp: -------------------------------------------------------------------------------- 1 | #ifndef _MSC_VER 2 | #warning deprecated 3 | #endif 4 | /* deprecated 5 | Matrix3i m = Matrix3i::Random(); 6 | cout << "Here is the matrix m:" << endl << m << endl; 7 | cout << "Here is the upper-triangular matrix extracted from m:" << endl 8 | << m.part() << endl; 9 | cout << "Here is the strictly-upper-triangular matrix extracted from m:" << endl 10 | << m.part() << endl; 11 | cout << "Here is the unit-lower-triangular matrix extracted from m:" << endl 12 | << m.part() << endl; 13 | */ -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_fixedBlock_int_int.cpp: -------------------------------------------------------------------------------- 1 | Matrix4d m = Vector4d(1,2,3,4).asDiagonal(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is m.fixed<2, 2>(2, 2):" << endl << m.block<2, 2>(2, 2) << endl; 4 | m.block<2, 2>(2, 0) = m.block<2, 2>(2, 2); 5 | cout << "Now the matrix m is:" << endl << m << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_identity.cpp: -------------------------------------------------------------------------------- 1 | cout << Matrix::Identity() << endl; 2 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_identity_int_int.cpp: -------------------------------------------------------------------------------- 1 | cout << MatrixXd::Identity(4, 3) << endl; 2 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_inverse.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Its inverse is:" << endl << m.inverse() << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_isDiagonal.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = 10000 * Matrix3d::Identity(); 2 | m(0,2) = 1; 3 | cout << "Here's the matrix m:" << endl << m << endl; 4 | cout << "m.isDiagonal() returns: " << m.isDiagonal() << endl; 5 | cout << "m.isDiagonal(1e-3) returns: " << m.isDiagonal(1e-3) << endl; 6 | 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_isIdentity.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Identity(); 2 | m(0,2) = 1e-4; 3 | cout << "Here's the matrix m:" << endl << m << endl; 4 | cout << "m.isIdentity() returns: " << m.isIdentity() << endl; 5 | cout << "m.isIdentity(1e-3) returns: " << m.isIdentity(1e-3) << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_isOnes.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Ones(); 2 | m(0,2) += 1e-4; 3 | cout << "Here's the matrix m:" << endl << m << endl; 4 | cout << "m.isOnes() returns: " << m.isOnes() << endl; 5 | cout << "m.isOnes(1e-3) returns: " << m.isOnes(1e-3) << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_isOrthogonal.cpp: -------------------------------------------------------------------------------- 1 | Vector3d v(1,0,0); 2 | Vector3d w(1e-4,0,1); 3 | cout << "Here's the vector v:" << endl << v << endl; 4 | cout << "Here's the vector w:" << endl << w << endl; 5 | cout << "v.isOrthogonal(w) returns: " << v.isOrthogonal(w) << endl; 6 | cout << "v.isOrthogonal(w,1e-3) returns: " << v.isOrthogonal(w,1e-3) << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_isUnitary.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Identity(); 2 | m(0,2) = 1e-4; 3 | cout << "Here's the matrix m:" << endl << m << endl; 4 | cout << "m.isUnitary() returns: " << m.isUnitary() << endl; 5 | cout << "m.isUnitary(1e-3) returns: " << m.isUnitary(1e-3) << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_isZero.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Zero(); 2 | m(0,2) = 1e-4; 3 | cout << "Here's the matrix m:" << endl << m << endl; 4 | cout << "m.isZero() returns: " << m.isZero() << endl; 5 | cout << "m.isZero(1e-3) returns: " << m.isZero(1e-3) << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_leftCols_int.cpp: -------------------------------------------------------------------------------- 1 | Array44i a = Array44i::Random(); 2 | cout << "Here is the array a:" << endl << a << endl; 3 | cout << "Here is a.leftCols(2):" << endl; 4 | cout << a.leftCols(2) << endl; 5 | a.leftCols(2).setZero(); 6 | cout << "Now the array a is:" << endl << a << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_marked.cpp: -------------------------------------------------------------------------------- 1 | #ifndef _MSC_VER 2 | #warning deprecated 3 | #endif 4 | /* 5 | Matrix3d m = Matrix3d::Zero(); 6 | m.part().setOnes(); 7 | cout << "Here is the matrix m:" << endl << m << endl; 8 | Matrix3d n = Matrix3d::Ones(); 9 | n.part() *= 2; 10 | cout << "Here is the matrix n:" << endl << n << endl; 11 | cout << "And now here is m.inverse()*n, taking advantage of the fact that" 12 | " m is upper-triangular:" << endl 13 | << m.marked().solveTriangular(n); 14 | */ -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_noalias.cpp: -------------------------------------------------------------------------------- 1 | Matrix2d a, b, c; a << 1,2,3,4; b << 5,6,7,8; 2 | c.noalias() = a * b; // this computes the product directly to c 3 | cout << c << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_ones.cpp: -------------------------------------------------------------------------------- 1 | cout << Matrix2d::Ones() << endl; 2 | cout << 6 * RowVector4i::Ones() << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_ones_int.cpp: -------------------------------------------------------------------------------- 1 | cout << 6 * RowVectorXi::Ones(4) << endl; 2 | cout << VectorXf::Ones(2) << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_ones_int_int.cpp: -------------------------------------------------------------------------------- 1 | cout << MatrixXi::Ones(2,3) << endl; 2 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_operatorNorm.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd ones = MatrixXd::Ones(3,3); 2 | cout << "The operator norm of the 3x3 matrix of ones is " 3 | << ones.operatorNorm() << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_part.cpp: -------------------------------------------------------------------------------- 1 | #ifndef _MSC_VER 2 | #warning deprecated 3 | #endif 4 | /* 5 | Matrix3d m = Matrix3d::Zero(); 6 | m.part().setOnes(); 7 | cout << "Here is the matrix m:" << endl << m << endl; 8 | cout << "And let us now compute m*m.adjoint() in a very optimized way" << endl 9 | << "taking advantage of the symmetry." << endl; 10 | Matrix3d n; 11 | n.part() = (m*m.adjoint()).lazy(); 12 | cout << "The result is:" << endl << n << endl; 13 | */ -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_prod.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is the product of all the coefficients:" << endl << m.prod() << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_random.cpp: -------------------------------------------------------------------------------- 1 | cout << 100 * Matrix2i::Random() << endl; 2 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_random_int.cpp: -------------------------------------------------------------------------------- 1 | cout << VectorXi::Random(2) << endl; 2 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_random_int_int.cpp: -------------------------------------------------------------------------------- 1 | cout << MatrixXi::Random(2,3) << endl; 2 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_replicate.cpp: -------------------------------------------------------------------------------- 1 | MatrixXi m = MatrixXi::Random(2,3); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "m.replicate<3,2>() = ..." << endl; 4 | cout << m.replicate<3,2>() << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_replicate_int_int.cpp: -------------------------------------------------------------------------------- 1 | Vector3i v = Vector3i::Random(); 2 | cout << "Here is the vector v:" << endl << v << endl; 3 | cout << "v.replicate(2,5) = ..." << endl; 4 | cout << v.replicate(2,5) << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_reverse.cpp: -------------------------------------------------------------------------------- 1 | MatrixXi m = MatrixXi::Random(3,4); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is the reverse of m:" << endl << m.reverse() << endl; 4 | cout << "Here is the coefficient (1,0) in the reverse of m:" << endl 5 | << m.reverse()(1,0) << endl; 6 | cout << "Let us overwrite this coefficient with the value 4." << endl; 7 | m.reverse()(1,0) = 4; 8 | cout << "Now the matrix m is:" << endl << m << endl; 9 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_rightCols_int.cpp: -------------------------------------------------------------------------------- 1 | Array44i a = Array44i::Random(); 2 | cout << "Here is the array a:" << endl << a << endl; 3 | cout << "Here is a.rightCols(2):" << endl; 4 | cout << a.rightCols(2) << endl; 5 | a.rightCols(2).setZero(); 6 | cout << "Now the array a is:" << endl << a << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_row.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Identity(); 2 | m.row(1) = Vector3d(4,5,6); 3 | cout << m << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_rowwise.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is the sum of each row:" << endl << m.rowwise().sum() << endl; 4 | cout << "Here is the maximum absolute value of each row:" 5 | << endl << m.cwiseAbs().rowwise().maxCoeff() << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_segment_int_int.cpp: -------------------------------------------------------------------------------- 1 | RowVector4i v = RowVector4i::Random(); 2 | cout << "Here is the vector v:" << endl << v << endl; 3 | cout << "Here is v.segment(1, 2):" << endl << v.segment(1, 2) << endl; 4 | v.segment(1, 2).setZero(); 5 | cout << "Now the vector v is:" << endl << v << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_select.cpp: -------------------------------------------------------------------------------- 1 | MatrixXi m(3, 3); 2 | m << 1, 2, 3, 3 | 4, 5, 6, 4 | 7, 8, 9; 5 | m = (m.array() >= 5).select(-m, m); 6 | cout << m << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_set.cpp: -------------------------------------------------------------------------------- 1 | Matrix3i m1; 2 | m1 << 1, 2, 3, 3 | 4, 5, 6, 4 | 7, 8, 9; 5 | cout << m1 << endl << endl; 6 | Matrix3i m2 = Matrix3i::Identity(); 7 | m2.block(0,0, 2,2) << 10, 11, 12, 13; 8 | cout << m2 << endl << endl; 9 | Vector2i v1; 10 | v1 << 14, 15; 11 | m2 << v1.transpose(), 16, 12 | v1, m1.block(1,1,2,2); 13 | cout << m2 << endl; 14 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_setIdentity.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Zero(); 2 | m.block<3,3>(1,0).setIdentity(); 3 | cout << m << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_setOnes.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | m.row(1).setOnes(); 3 | cout << m << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_setRandom.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Zero(); 2 | m.col(1).setRandom(); 3 | cout << m << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_setZero.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | m.row(1).setZero(); 3 | cout << m << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_start_int.cpp: -------------------------------------------------------------------------------- 1 | RowVector4i v = RowVector4i::Random(); 2 | cout << "Here is the vector v:" << endl << v << endl; 3 | cout << "Here is v.head(2):" << endl << v.head(2) << endl; 4 | v.head(2).setZero(); 5 | cout << "Now the vector v is:" << endl << v << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_template_int_bottomRows.cpp: -------------------------------------------------------------------------------- 1 | Array44i a = Array44i::Random(); 2 | cout << "Here is the array a:" << endl << a << endl; 3 | cout << "Here is a.bottomRows<2>():" << endl; 4 | cout << a.bottomRows<2>() << endl; 5 | a.bottomRows<2>().setZero(); 6 | cout << "Now the array a is:" << endl << a << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_template_int_end.cpp: -------------------------------------------------------------------------------- 1 | RowVector4i v = RowVector4i::Random(); 2 | cout << "Here is the vector v:" << endl << v << endl; 3 | cout << "Here is v.tail(2):" << endl << v.tail<2>() << endl; 4 | v.tail<2>().setZero(); 5 | cout << "Now the vector v is:" << endl << v << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_template_int_int_bottomLeftCorner.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is m.bottomLeftCorner<2,2>():" << endl; 4 | cout << m.bottomLeftCorner<2,2>() << endl; 5 | m.bottomLeftCorner<2,2>().setZero(); 6 | cout << "Now the matrix m is:" << endl << m << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_template_int_int_bottomRightCorner.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is m.bottomRightCorner<2,2>():" << endl; 4 | cout << m.bottomRightCorner<2,2>() << endl; 5 | m.bottomRightCorner<2,2>().setZero(); 6 | cout << "Now the matrix m is:" << endl << m << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_template_int_int_topLeftCorner.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is m.topLeftCorner<2,2>():" << endl; 4 | cout << m.topLeftCorner<2,2>() << endl; 5 | m.topLeftCorner<2,2>().setZero(); 6 | cout << "Now the matrix m is:" << endl << m << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_template_int_int_topRightCorner.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is m.topRightCorner<2,2>():" << endl; 4 | cout << m.topRightCorner<2,2>() << endl; 5 | m.topRightCorner<2,2>().setZero(); 6 | cout << "Now the matrix m is:" << endl << m << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_template_int_leftCols.cpp: -------------------------------------------------------------------------------- 1 | Array44i a = Array44i::Random(); 2 | cout << "Here is the array a:" << endl << a << endl; 3 | cout << "Here is a.leftCols<2>():" << endl; 4 | cout << a.leftCols<2>() << endl; 5 | a.leftCols<2>().setZero(); 6 | cout << "Now the array a is:" << endl << a << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_template_int_rightCols.cpp: -------------------------------------------------------------------------------- 1 | Array44i a = Array44i::Random(); 2 | cout << "Here is the array a:" << endl << a << endl; 3 | cout << "Here is a.rightCols<2>():" << endl; 4 | cout << a.rightCols<2>() << endl; 5 | a.rightCols<2>().setZero(); 6 | cout << "Now the array a is:" << endl << a << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_template_int_segment.cpp: -------------------------------------------------------------------------------- 1 | RowVector4i v = RowVector4i::Random(); 2 | cout << "Here is the vector v:" << endl << v << endl; 3 | cout << "Here is v.segment<2>(1):" << endl << v.segment<2>(1) << endl; 4 | v.segment<2>(2).setZero(); 5 | cout << "Now the vector v is:" << endl << v << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_template_int_start.cpp: -------------------------------------------------------------------------------- 1 | RowVector4i v = RowVector4i::Random(); 2 | cout << "Here is the vector v:" << endl << v << endl; 3 | cout << "Here is v.head(2):" << endl << v.head<2>() << endl; 4 | v.head<2>().setZero(); 5 | cout << "Now the vector v is:" << endl << v << endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_template_int_topRows.cpp: -------------------------------------------------------------------------------- 1 | Array44i a = Array44i::Random(); 2 | cout << "Here is the array a:" << endl << a << endl; 3 | cout << "Here is a.topRows<2>():" << endl; 4 | cout << a.topRows<2>() << endl; 5 | a.topRows<2>().setZero(); 6 | cout << "Now the array a is:" << endl << a << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_topLeftCorner_int_int.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is m.topLeftCorner(2, 2):" << endl; 4 | cout << m.topLeftCorner(2, 2) << endl; 5 | m.topLeftCorner(2, 2).setZero(); 6 | cout << "Now the matrix m is:" << endl << m << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_topRightCorner_int_int.cpp: -------------------------------------------------------------------------------- 1 | Matrix4i m = Matrix4i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is m.topRightCorner(2, 2):" << endl; 4 | cout << m.topRightCorner(2, 2) << endl; 5 | m.topRightCorner(2, 2).setZero(); 6 | cout << "Now the matrix m is:" << endl << m << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_topRows_int.cpp: -------------------------------------------------------------------------------- 1 | Array44i a = Array44i::Random(); 2 | cout << "Here is the array a:" << endl << a << endl; 3 | cout << "Here is a.topRows(2):" << endl; 4 | cout << a.topRows(2) << endl; 5 | a.topRows(2).setZero(); 6 | cout << "Now the array a is:" << endl << a << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_transpose.cpp: -------------------------------------------------------------------------------- 1 | Matrix2i m = Matrix2i::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is the transpose of m:" << endl << m.transpose() << endl; 4 | cout << "Here is the coefficient (1,0) in the transpose of m:" << endl 5 | << m.transpose()(1,0) << endl; 6 | cout << "Let us overwrite this coefficient with the value 0." << endl; 7 | m.transpose()(1,0) = 0; 8 | cout << "Now the matrix m is:" << endl << m << endl; 9 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_zero.cpp: -------------------------------------------------------------------------------- 1 | cout << Matrix2d::Zero() << endl; 2 | cout << RowVector4i::Zero() << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_zero_int.cpp: -------------------------------------------------------------------------------- 1 | cout << RowVectorXi::Zero(4) << endl; 2 | cout << VectorXf::Zero(2) << endl; 3 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/MatrixBase_zero_int_int.cpp: -------------------------------------------------------------------------------- 1 | cout << MatrixXi::Zero(2,3) << endl; 2 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Matrix_resize_NoChange_int.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd m(3,4); 2 | m.resize(NoChange, 5); 3 | cout << "m: " << m.rows() << " rows, " << m.cols() << " cols" << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Matrix_resize_int.cpp: -------------------------------------------------------------------------------- 1 | VectorXd v(10); 2 | v.resize(3); 3 | RowVector3d w; 4 | w.resize(3); // this is legal, but has no effect 5 | cout << "v: " << v.rows() << " rows, " << v.cols() << " cols" << endl; 6 | cout << "w: " << w.rows() << " rows, " << w.cols() << " cols" << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Matrix_resize_int_NoChange.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd m(3,4); 2 | m.resize(5, NoChange); 3 | cout << "m: " << m.rows() << " rows, " << m.cols() << " cols" << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Matrix_resize_int_int.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd m(2,3); 2 | m << 1,2,3,4,5,6; 3 | cout << "here's the 2x3 matrix m:" << endl << m << endl; 4 | cout << "let's resize m to 3x2. This is a conservative resizing because 2*3==3*2." << endl; 5 | m.resize(3,2); 6 | cout << "here's the 3x2 matrix m:" << endl << m << endl; 7 | cout << "now let's resize m to size 2x2. This is NOT a conservative resizing, so it becomes uninitialized:" << endl; 8 | m.resize(2,2); 9 | cout << m << endl; 10 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Matrix_setConstant_int.cpp: -------------------------------------------------------------------------------- 1 | VectorXf v; 2 | v.setConstant(3, 5); 3 | cout << v << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Matrix_setConstant_int_int.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf m; 2 | m.setConstant(3, 3, 5); 3 | cout << m << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Matrix_setIdentity_int_int.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf m; 2 | m.setIdentity(3, 3); 3 | cout << m << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Matrix_setOnes_int.cpp: -------------------------------------------------------------------------------- 1 | VectorXf v; 2 | v.setOnes(3); 3 | cout << v << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Matrix_setOnes_int_int.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf m; 2 | m.setOnes(3, 3); 3 | cout << m << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Matrix_setRandom_int.cpp: -------------------------------------------------------------------------------- 1 | VectorXf v; 2 | v.setRandom(3); 3 | cout << v << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Matrix_setRandom_int_int.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf m; 2 | m.setRandom(3, 3); 3 | cout << m << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Matrix_setZero_int.cpp: -------------------------------------------------------------------------------- 1 | VectorXf v; 2 | v.setZero(3); 3 | cout << v << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Matrix_setZero_int_int.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf m; 2 | m.setZero(3, 3); 3 | cout << m << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/PartialPivLU_solve.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd A = MatrixXd::Random(3,3); 2 | MatrixXd B = MatrixXd::Random(3,2); 3 | cout << "Here is the invertible matrix A:" << endl << A << endl; 4 | cout << "Here is the matrix B:" << endl << B << endl; 5 | MatrixXd X = A.lu().solve(B); 6 | cout << "Here is the (unique) solution X to the equation AX=B:" << endl << X << endl; 7 | cout << "Relative error: " << (A*X-B).norm() / B.norm() << endl; 8 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/PartialRedux_count.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is the count of elements larger or equal than 0.5 of each row:" << endl << (m.array() >= 0.5).rowwise().count() << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/PartialRedux_maxCoeff.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is the maximum of each column:" << endl << m.colwise().maxCoeff() << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/PartialRedux_minCoeff.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is the minimum of each column:" << endl << m.colwise().minCoeff() << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/PartialRedux_norm.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is the norm of each column:" << endl << m.colwise().norm() << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/PartialRedux_prod.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is the product of each row:" << endl << m.rowwise().prod() << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/PartialRedux_squaredNorm.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is the square norm of each row:" << endl << m.rowwise().squaredNorm() << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/PartialRedux_sum.cpp: -------------------------------------------------------------------------------- 1 | Matrix3d m = Matrix3d::Random(); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is the sum of each row:" << endl << m.rowwise().sum() << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/RealSchur_RealSchur_MatrixType.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd A = MatrixXd::Random(6,6); 2 | cout << "Here is a random 6x6 matrix, A:" << endl << A << endl << endl; 3 | 4 | RealSchur schur(A); 5 | cout << "The orthogonal matrix U is:" << endl << schur.matrixU() << endl; 6 | cout << "The quasi-triangular matrix T is:" << endl << schur.matrixT() << endl << endl; 7 | 8 | MatrixXd U = schur.matrixU(); 9 | MatrixXd T = schur.matrixT(); 10 | cout << "U * T * U^T = " << endl << U * T * U.transpose() << endl; 11 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/RealSchur_compute.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf A = MatrixXf::Random(4,4); 2 | RealSchur schur(4); 3 | schur.compute(A, /* computeU = */ false); 4 | cout << "The matrix T in the decomposition of A is:" << endl << schur.matrixT() << endl; 5 | schur.compute(A.inverse(), /* computeU = */ false); 6 | cout << "The matrix T in the decomposition of A^(-1) is:" << endl << schur.matrixT() << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/SelfAdjointEigenSolver_SelfAdjointEigenSolver.cpp: -------------------------------------------------------------------------------- 1 | SelfAdjointEigenSolver es; 2 | Matrix4f X = Matrix4f::Random(4,4); 3 | Matrix4f A = X + X.transpose(); 4 | es.compute(A); 5 | cout << "The eigenvalues of A are: " << es.eigenvalues().transpose() << endl; 6 | es.compute(A + Matrix4f::Identity(4,4)); // re-use es to compute eigenvalues of A+I 7 | cout << "The eigenvalues of A+I are: " << es.eigenvalues().transpose() << endl; 8 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/SelfAdjointEigenSolver_compute_MatrixType.cpp: -------------------------------------------------------------------------------- 1 | SelfAdjointEigenSolver es(4); 2 | MatrixXf X = MatrixXf::Random(4,4); 3 | MatrixXf A = X + X.transpose(); 4 | es.compute(A); 5 | cout << "The eigenvalues of A are: " << es.eigenvalues().transpose() << endl; 6 | es.compute(A + MatrixXf::Identity(4,4)); // re-use es to compute eigenvalues of A+I 7 | cout << "The eigenvalues of A+I are: " << es.eigenvalues().transpose() << endl; 8 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/SelfAdjointEigenSolver_compute_MatrixType2.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd X = MatrixXd::Random(5,5); 2 | MatrixXd A = X * X.transpose(); 3 | X = MatrixXd::Random(5,5); 4 | MatrixXd B = X * X.transpose(); 5 | 6 | GeneralizedSelfAdjointEigenSolver es(A,B,EigenvaluesOnly); 7 | cout << "The eigenvalues of the pencil (A,B) are:" << endl << es.eigenvalues() << endl; 8 | es.compute(B,A,false); 9 | cout << "The eigenvalues of the pencil (B,A) are:" << endl << es.eigenvalues() << endl; 10 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/SelfAdjointEigenSolver_eigenvalues.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd ones = MatrixXd::Ones(3,3); 2 | SelfAdjointEigenSolver es(ones); 3 | cout << "The eigenvalues of the 3x3 matrix of ones are:" 4 | << endl << es.eigenvalues() << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/SelfAdjointEigenSolver_eigenvectors.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd ones = MatrixXd::Ones(3,3); 2 | SelfAdjointEigenSolver es(ones); 3 | cout << "The first eigenvector of the 3x3 matrix of ones is:" 4 | << endl << es.eigenvectors().col(1) << endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/SelfAdjointEigenSolver_operatorInverseSqrt.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd X = MatrixXd::Random(4,4); 2 | MatrixXd A = X * X.transpose(); 3 | cout << "Here is a random positive-definite matrix, A:" << endl << A << endl << endl; 4 | 5 | SelfAdjointEigenSolver es(A); 6 | cout << "The inverse square root of A is: " << endl; 7 | cout << es.operatorInverseSqrt() << endl; 8 | cout << "We can also compute it with operatorSqrt() and inverse(). That yields: " << endl; 9 | cout << es.operatorSqrt().inverse() << endl; 10 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/SelfAdjointEigenSolver_operatorSqrt.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd X = MatrixXd::Random(4,4); 2 | MatrixXd A = X * X.transpose(); 3 | cout << "Here is a random positive-definite matrix, A:" << endl << A << endl << endl; 4 | 5 | SelfAdjointEigenSolver es(A); 6 | MatrixXd sqrtA = es.operatorSqrt(); 7 | cout << "The square root of A is: " << endl << sqrtA << endl; 8 | cout << "If we square this, we get: " << endl << sqrtA*sqrtA << endl; 9 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/SelfAdjointView_eigenvalues.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd ones = MatrixXd::Ones(3,3); 2 | VectorXd eivals = ones.selfadjointView().eigenvalues(); 3 | cout << "The eigenvalues of the 3x3 matrix of ones are:" << endl << eivals << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/SelfAdjointView_operatorNorm.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd ones = MatrixXd::Ones(3,3); 2 | cout << "The operator norm of the 3x3 matrix of ones is " 3 | << ones.selfadjointView().operatorNorm() << endl; 4 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/TopicAliasing_block.cpp: -------------------------------------------------------------------------------- 1 | MatrixXi mat(3,3); 2 | mat << 1, 2, 3, 4, 5, 6, 7, 8, 9; 3 | cout << "Here is the matrix mat:\n" << mat << endl; 4 | 5 | // This assignment shows the aliasing problem 6 | mat.bottomRightCorner(2,2) = mat.topLeftCorner(2,2); 7 | cout << "After the assignment, mat = \n" << mat << endl; 8 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/TopicAliasing_block_correct.cpp: -------------------------------------------------------------------------------- 1 | MatrixXi mat(3,3); 2 | mat << 1, 2, 3, 4, 5, 6, 7, 8, 9; 3 | cout << "Here is the matrix mat:\n" << mat << endl; 4 | 5 | // The eval() solves the aliasing problem 6 | mat.bottomRightCorner(2,2) = mat.topLeftCorner(2,2).eval(); 7 | cout << "After the assignment, mat = \n" << mat << endl; 8 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/TopicAliasing_cwise.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf mat(2,2); 2 | mat << 1, 2, 4, 7; 3 | cout << "Here is the matrix mat:\n" << mat << endl << endl; 4 | 5 | mat = 2 * mat; 6 | cout << "After 'mat = 2 * mat', mat = \n" << mat << endl << endl; 7 | 8 | 9 | mat = mat - MatrixXf::Identity(2,2); 10 | cout << "After the subtraction, it becomes\n" << mat << endl << endl; 11 | 12 | 13 | ArrayXXf arr = mat; 14 | arr = arr.square(); 15 | cout << "After squaring, it becomes\n" << arr << endl << endl; 16 | 17 | // Combining all operations in one statement: 18 | mat << 1, 2, 4, 7; 19 | mat = (2 * mat - MatrixXf::Identity(2,2)).array().square(); 20 | cout << "Doing everything at once yields\n" << mat << endl << endl; 21 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/TopicAliasing_mult1.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf matA(2,2); 2 | matA << 2, 0, 0, 2; 3 | matA = matA * matA; 4 | cout << matA; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/TopicAliasing_mult2.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf matA(2,2), matB(2,2); 2 | matA << 2, 0, 0, 2; 3 | 4 | // Simple but not quite as efficient 5 | matB = matA * matA; 6 | cout << matB << endl << endl; 7 | 8 | // More complicated but also more efficient 9 | matB.noalias() = matA * matA; 10 | cout << matB; 11 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/TopicAliasing_mult3.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf matA(2,2); 2 | matA << 2, 0, 0, 2; 3 | matA.noalias() = matA * matA; 4 | cout << matA; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/TopicStorageOrders_example.cpp: -------------------------------------------------------------------------------- 1 | Matrix Acolmajor; 2 | Acolmajor << 8, 2, 2, 9, 3 | 9, 1, 4, 4, 4 | 3, 5, 4, 5; 5 | cout << "The matrix A:" << endl; 6 | cout << Acolmajor << endl << endl; 7 | 8 | cout << "In memory (column-major):" << endl; 9 | for (int i = 0; i < Acolmajor.size(); i++) 10 | cout << *(Acolmajor.data() + i) << " "; 11 | cout << endl << endl; 12 | 13 | Matrix Arowmajor = Acolmajor; 14 | cout << "In memory (row-major):" << endl; 15 | for (int i = 0; i < Arowmajor.size(); i++) 16 | cout << *(Arowmajor.data() + i) << " "; 17 | cout << endl; 18 | 19 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tridiagonalization_Tridiagonalization_MatrixType.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd X = MatrixXd::Random(5,5); 2 | MatrixXd A = X + X.transpose(); 3 | cout << "Here is a random symmetric 5x5 matrix:" << endl << A << endl << endl; 4 | Tridiagonalization triOfA(A); 5 | MatrixXd Q = triOfA.matrixQ(); 6 | cout << "The orthogonal matrix Q is:" << endl << Q << endl; 7 | MatrixXd T = triOfA.matrixT(); 8 | cout << "The tridiagonal matrix T is:" << endl << T << endl << endl; 9 | cout << "Q * T * Q^T = " << endl << Q * T * Q.transpose() << endl; 10 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tridiagonalization_compute.cpp: -------------------------------------------------------------------------------- 1 | Tridiagonalization tri; 2 | MatrixXf X = MatrixXf::Random(4,4); 3 | MatrixXf A = X + X.transpose(); 4 | tri.compute(A); 5 | cout << "The matrix T in the tridiagonal decomposition of A is: " << endl; 6 | cout << tri.matrixT() << endl; 7 | tri.compute(2*A); // re-use tri to compute eigenvalues of 2A 8 | cout << "The matrix T in the tridiagonal decomposition of 2A is: " << endl; 9 | cout << tri.matrixT() << endl; 10 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tridiagonalization_decomposeInPlace.cpp: -------------------------------------------------------------------------------- 1 | MatrixXd X = MatrixXd::Random(5,5); 2 | MatrixXd A = X + X.transpose(); 3 | cout << "Here is a random symmetric 5x5 matrix:" << endl << A << endl << endl; 4 | 5 | VectorXd diag(5); 6 | VectorXd subdiag(4); 7 | internal::tridiagonalization_inplace(A, diag, subdiag, true); 8 | cout << "The orthogonal matrix Q is:" << endl << A << endl; 9 | cout << "The diagonal of the tridiagonal matrix T is:" << endl << diag << endl; 10 | cout << "The subdiagonal of the tridiagonal matrix T is:" << endl << subdiag << endl; 11 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tridiagonalization_diagonal.cpp: -------------------------------------------------------------------------------- 1 | MatrixXcd X = MatrixXcd::Random(4,4); 2 | MatrixXcd A = X + X.adjoint(); 3 | cout << "Here is a random self-adjoint 4x4 matrix:" << endl << A << endl << endl; 4 | 5 | Tridiagonalization triOfA(A); 6 | MatrixXd T = triOfA.matrixT(); 7 | cout << "The tridiagonal matrix T is:" << endl << T << endl << endl; 8 | 9 | cout << "We can also extract the diagonals of T directly ..." << endl; 10 | VectorXd diag = triOfA.diagonal(); 11 | cout << "The diagonal is:" << endl << diag << endl; 12 | VectorXd subdiag = triOfA.subDiagonal(); 13 | cout << "The subdiagonal is:" << endl << subdiag << endl; 14 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tridiagonalization_householderCoefficients.cpp: -------------------------------------------------------------------------------- 1 | Matrix4d X = Matrix4d::Random(4,4); 2 | Matrix4d A = X + X.transpose(); 3 | cout << "Here is a random symmetric 4x4 matrix:" << endl << A << endl; 4 | Tridiagonalization triOfA(A); 5 | Vector3d hc = triOfA.householderCoefficients(); 6 | cout << "The vector of Householder coefficients is:" << endl << hc << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tridiagonalization_packedMatrix.cpp: -------------------------------------------------------------------------------- 1 | Matrix4d X = Matrix4d::Random(4,4); 2 | Matrix4d A = X + X.transpose(); 3 | cout << "Here is a random symmetric 4x4 matrix:" << endl << A << endl; 4 | Tridiagonalization triOfA(A); 5 | Matrix4d pm = triOfA.packedMatrix(); 6 | cout << "The packed matrix M is:" << endl << pm << endl; 7 | cout << "The diagonal and subdiagonal corresponds to the matrix T, which is:" 8 | << endl << triOfA.matrixT() << endl; 9 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_AdvancedInitialization_Block.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf matA(2, 2); 2 | matA << 1, 2, 3, 4; 3 | MatrixXf matB(4, 4); 4 | matB << matA, matA/10, matA/10, matA; 5 | std::cout << matB << std::endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_AdvancedInitialization_CommaTemporary.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf mat = MatrixXf::Random(2, 3); 2 | std::cout << mat << std::endl << std::endl; 3 | mat = (MatrixXf(2,2) << 0, 1, 1, 0).finished() * mat; 4 | std::cout << mat << std::endl; 5 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_AdvancedInitialization_Join.cpp: -------------------------------------------------------------------------------- 1 | RowVectorXd vec1(3); 2 | vec1 << 1, 2, 3; 3 | std::cout << "vec1 = " << vec1 << std::endl; 4 | 5 | RowVectorXd vec2(4); 6 | vec2 << 1, 4, 9, 16;; 7 | std::cout << "vec2 = " << vec2 << std::endl; 8 | 9 | RowVectorXd joined(7); 10 | joined << vec1, vec2; 11 | std::cout << "joined = " << joined << std::endl; 12 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_AdvancedInitialization_LinSpaced.cpp: -------------------------------------------------------------------------------- 1 | ArrayXXf table(10, 4); 2 | table.col(0) = ArrayXf::LinSpaced(10, 0, 90); 3 | table.col(1) = M_PI / 180 * table.col(0); 4 | table.col(2) = table.col(1).sin(); 5 | table.col(3) = table.col(1).cos(); 6 | std::cout << " Degrees Radians Sine Cosine\n"; 7 | std::cout << table << std::endl; 8 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_AdvancedInitialization_Zero.cpp: -------------------------------------------------------------------------------- 1 | std::cout << "A fixed-size array:\n"; 2 | Array33f a1 = Array33f::Zero(); 3 | std::cout << a1 << "\n\n"; 4 | 5 | 6 | std::cout << "A one-dimensional dynamic-size array:\n"; 7 | ArrayXf a2 = ArrayXf::Zero(3); 8 | std::cout << a2 << "\n\n"; 9 | 10 | 11 | std::cout << "A two-dimensional dynamic-size array:\n"; 12 | ArrayXXf a3 = ArrayXXf::Zero(3, 4); 13 | std::cout << a3 << "\n"; 14 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_Map_rowmajor.cpp: -------------------------------------------------------------------------------- 1 | int array[8]; 2 | for(int i = 0; i < 8; ++i) array[i] = i; 3 | cout << "Column-major:\n" << Map >(array) << endl; 4 | cout << "Row-major:\n" << Map >(array) << endl; 5 | cout << "Row-major using stride:\n" << 6 | Map, Unaligned, Stride<1,4> >(array) << endl; 7 | 8 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_commainit_01.cpp: -------------------------------------------------------------------------------- 1 | Matrix3f m; 2 | m << 1, 2, 3, 3 | 4, 5, 6, 4 | 7, 8, 9; 5 | std::cout << m; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_commainit_01b.cpp: -------------------------------------------------------------------------------- 1 | Matrix3f m; 2 | m.row(0) << 1, 2, 3; 3 | m.block(1,0,2,2) << 4, 5, 7, 8; 4 | m.col(2).tail(2) << 6, 9; 5 | std::cout << m; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_commainit_02.cpp: -------------------------------------------------------------------------------- 1 | int rows=5, cols=5; 2 | MatrixXf m(rows,cols); 3 | m << (Matrix3f() << 1, 2, 3, 4, 5, 6, 7, 8, 9).finished(), 4 | MatrixXf::Zero(3,cols-3), 5 | MatrixXf::Zero(rows-3,3), 6 | MatrixXf::Identity(rows-3,cols-3); 7 | cout << m; 8 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_solve_matrix_inverse.cpp: -------------------------------------------------------------------------------- 1 | Matrix3f A; 2 | Vector3f b; 3 | A << 1,2,3, 4,5,6, 7,8,10; 4 | b << 3, 3, 4; 5 | Vector3f x = A.inverse() * b; 6 | cout << "The solution is:" << endl << x << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_solve_multiple_rhs.cpp: -------------------------------------------------------------------------------- 1 | Matrix3f A(3,3); 2 | A << 1,2,3, 4,5,6, 7,8,10; 3 | Matrix B; 4 | B << 3,1, 3,1, 4,1; 5 | Matrix X; 6 | X = A.fullPivLu().solve(B); 7 | cout << "The solution with right-hand side (3,3,4) is:" << endl; 8 | cout << X.col(0) << endl; 9 | cout << "The solution with right-hand side (1,1,1) is:" << endl; 10 | cout << X.col(1) << endl; 11 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_solve_reuse_decomposition.cpp: -------------------------------------------------------------------------------- 1 | Matrix3f A(3,3); 2 | A << 1,2,3, 4,5,6, 7,8,10; 3 | PartialPivLU luOfA(A); // compute LU decomposition of A 4 | Vector3f b; 5 | b << 3,3,4; 6 | Vector3f x; 7 | x = luOfA.solve(b); 8 | cout << "The solution with right-hand side (3,3,4) is:" << endl; 9 | cout << x << endl; 10 | b << 1,1,1; 11 | x = luOfA.solve(b); 12 | cout << "The solution with right-hand side (1,1,1) is:" << endl; 13 | cout << x << endl; 14 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_solve_singular.cpp: -------------------------------------------------------------------------------- 1 | Matrix3f A; 2 | Vector3f b; 3 | A << 1,2,3, 4,5,6, 7,8,9; 4 | b << 3, 3, 4; 5 | cout << "Here is the matrix A:" << endl << A << endl; 6 | cout << "Here is the vector b:" << endl << b << endl; 7 | Vector3f x; 8 | x = A.lu().solve(b); 9 | cout << "The solution is:" << endl << x << endl; 10 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_solve_triangular.cpp: -------------------------------------------------------------------------------- 1 | Matrix3f A; 2 | Vector3f b; 3 | A << 1,2,3, 0,5,6, 0,0,10; 4 | b << 3, 3, 4; 5 | cout << "Here is the matrix A:" << endl << A << endl; 6 | cout << "Here is the vector b:" << endl << b << endl; 7 | Vector3f x = A.triangularView().solve(b); 8 | cout << "The solution is:" << endl << x << endl; 9 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Tutorial_solve_triangular_inplace.cpp: -------------------------------------------------------------------------------- 1 | Matrix3f A; 2 | Vector3f b; 3 | A << 1,2,3, 0,5,6, 0,0,10; 4 | b << 3, 3, 4; 5 | A.triangularView().solveInPlace(b); 6 | cout << "The solution is:" << endl << b << endl; 7 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/Vectorwise_reverse.cpp: -------------------------------------------------------------------------------- 1 | MatrixXi m = MatrixXi::Random(3,4); 2 | cout << "Here is the matrix m:" << endl << m << endl; 3 | cout << "Here is the rowwise reverse of m:" << endl << m.rowwise().reverse() << endl; 4 | cout << "Here is the colwise reverse of m:" << endl << m.colwise().reverse() << endl; 5 | 6 | cout << "Here is the coefficient (1,0) in the rowise reverse of m:" << endl 7 | << m.rowwise().reverse()(1,0) << endl; 8 | cout << "Let us overwrite this coefficient with the value 4." << endl; 9 | //m.colwise().reverse()(1,0) = 4; 10 | cout << "Now the matrix m is:" << endl << m << endl; 11 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/class_FullPivLU.cpp: -------------------------------------------------------------------------------- 1 | typedef Matrix Matrix5x3; 2 | typedef Matrix Matrix5x5; 3 | Matrix5x3 m = Matrix5x3::Random(); 4 | cout << "Here is the matrix m:" << endl << m << endl; 5 | Eigen::FullPivLU lu(m); 6 | cout << "Here is, up to permutations, its LU decomposition matrix:" 7 | << endl << lu.matrixLU() << endl; 8 | cout << "Here is the L part:" << endl; 9 | Matrix5x5 l = Matrix5x5::Identity(); 10 | l.block<5,3>(0,0).triangularView() = lu.matrixLU(); 11 | cout << l << endl; 12 | cout << "Here is the U part:" << endl; 13 | Matrix5x3 u = lu.matrixLU().triangularView(); 14 | cout << u << endl; 15 | cout << "Let us now reconstruct the original matrix m:" << endl; 16 | cout << lu.permutationP().inverse() * l * u * lu.permutationQ().inverse() << endl; 17 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/compile_snippet.cpp.in: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main(int, char**) 8 | { 9 | cout.precision(3); 10 | ${snippet_source_code} 11 | return 0; 12 | } 13 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/tut_arithmetic_redux_minmax.cpp: -------------------------------------------------------------------------------- 1 | Matrix3f m = Matrix3f::Random(); 2 | std::ptrdiff_t i, j; 3 | float minOfM = m.minCoeff(&i,&j); 4 | cout << "Here is the matrix m:\n" << m << endl; 5 | cout << "Its minimum coefficient (" << minOfM 6 | << ") is at position (" << i << "," << j << ")\n\n"; 7 | 8 | RowVector4i v = RowVector4i::Random(); 9 | int maxOfV = v.maxCoeff(&i); 10 | cout << "Here is the vector v: " << v << endl; 11 | cout << "Its maximum coefficient (" << maxOfV 12 | << ") is at position " << i << endl; 13 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/tut_arithmetic_transpose_aliasing.cpp: -------------------------------------------------------------------------------- 1 | Matrix2i a; a << 1, 2, 3, 4; 2 | cout << "Here is the matrix a:\n" << a << endl; 3 | 4 | a = a.transpose(); // !!! do NOT do this !!! 5 | cout << "and the result of the aliasing effect:\n" << a << endl; -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/tut_arithmetic_transpose_conjugate.cpp: -------------------------------------------------------------------------------- 1 | MatrixXcf a = MatrixXcf::Random(2,2); 2 | cout << "Here is the matrix a\n" << a << endl; 3 | 4 | cout << "Here is the matrix a^T\n" << a.transpose() << endl; 5 | 6 | 7 | cout << "Here is the conjugate of a\n" << a.conjugate() << endl; 8 | 9 | 10 | cout << "Here is the matrix a^*\n" << a.adjoint() << endl; 11 | 12 | 13 | -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/tut_arithmetic_transpose_inplace.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf a(2,3); a << 1, 2, 3, 4, 5, 6; 2 | cout << "Here is the initial matrix a:\n" << a << endl; 3 | 4 | 5 | a.transposeInPlace(); 6 | cout << "and after being transposed:\n" << a << endl; -------------------------------------------------------------------------------- /eigen-patched/doc/snippets/tut_matrix_assignment_resizing.cpp: -------------------------------------------------------------------------------- 1 | MatrixXf a(2,2); 2 | std::cout << "a is of size " << a.rows() << "x" << a.cols() << std::endl; 3 | MatrixXf b(3,3); 4 | a = b; 5 | std::cout << "a is now of size " << a.rows() << "x" << a.cols() << std::endl; 6 | -------------------------------------------------------------------------------- /eigen-patched/doc/special_examples/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | 2 | if(NOT EIGEN_TEST_NOQT) 3 | find_package(Qt4) 4 | if(QT4_FOUND) 5 | include(${QT_USE_FILE}) 6 | endif() 7 | endif(NOT EIGEN_TEST_NOQT) 8 | 9 | 10 | if(QT4_FOUND) 11 | add_executable(Tutorial_sparse_example Tutorial_sparse_example.cpp Tutorial_sparse_example_details.cpp) 12 | target_link_libraries(Tutorial_sparse_example ${EIGEN_STANDARD_LIBRARIES_TO_LINK_TO} ${QT_QTCORE_LIBRARY} ${QT_QTGUI_LIBRARY}) 13 | 14 | add_custom_command( 15 | TARGET Tutorial_sparse_example 16 | POST_BUILD 17 | COMMAND Tutorial_sparse_example 18 | ARGS ${CMAKE_CURRENT_BINARY_DIR}/../html/Tutorial_sparse_example.jpeg 19 | ) 20 | add_dependencies(all_examples Tutorial_sparse_example) 21 | endif(QT4_FOUND) 22 | -------------------------------------------------------------------------------- /eigen-patched/eigen3.pc.in: -------------------------------------------------------------------------------- 1 | Name: Eigen3 2 | Description: A C++ template library for linear algebra: vectors, matrices, and related algorithms 3 | Requires: 4 | Version: ${EIGEN_VERSION_NUMBER} 5 | Libs: 6 | Cflags: -I${INCLUDE_INSTALL_DIR} 7 | -------------------------------------------------------------------------------- /eigen-patched/failtest/block_nonconst_ctor_on_const_xpr_0.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(CV_QUALIFIER Matrix3d &m){ 12 | Block b(m,0,0); 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/block_nonconst_ctor_on_const_xpr_1.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(CV_QUALIFIER Matrix3d &m){ 12 | Block b(m,0,0,3,3); 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/block_nonconst_ctor_on_const_xpr_2.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(CV_QUALIFIER Matrix3d &m){ 12 | // row/column constructor 13 | Block b(m,0); 14 | } 15 | 16 | int main() {} 17 | -------------------------------------------------------------------------------- /eigen-patched/failtest/block_on_const_type_actually_const_0.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(){ 12 | Matrix3f m; 13 | Block(m, 0, 0, 3, 3).coeffRef(0, 0) = 1.0f; 14 | } 15 | 16 | int main() {} 17 | -------------------------------------------------------------------------------- /eigen-patched/failtest/block_on_const_type_actually_const_1.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(){ 12 | MatrixXf m; 13 | Block(m, 0, 0).coeffRef(0, 0) = 1.0f; 14 | } 15 | 16 | int main() {} 17 | -------------------------------------------------------------------------------- /eigen-patched/failtest/const_qualified_block_method_retval_0.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(CV_QUALIFIER Matrix3d &m){ 12 | Block b(m.block<3,3>(0,0)); 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/const_qualified_block_method_retval_1.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(CV_QUALIFIER Matrix3d &m){ 12 | Block b(m.block(0,0,3,3)); 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/const_qualified_diagonal_method_retval.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(CV_QUALIFIER Matrix3d &m){ 12 | Diagonal b(m.diagonal()); 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/const_qualified_transpose_method_retval.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(CV_QUALIFIER Matrix3d &m){ 12 | Transpose b(m.transpose()); 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/diagonal_nonconst_ctor_on_const_xpr.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(CV_QUALIFIER Matrix3d &m){ 12 | Diagonal d(m); 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/diagonal_on_const_type_actually_const.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(){ 12 | MatrixXf m; 13 | Diagonal(m).coeffRef(0) = 1.0f; 14 | } 15 | 16 | int main() {} 17 | -------------------------------------------------------------------------------- /eigen-patched/failtest/failtest_sanity_check.cpp: -------------------------------------------------------------------------------- 1 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 2 | This is just some text that won't compile as a C++ file, as a basic sanity check for failtest. 3 | #else 4 | int main() {} 5 | #endif 6 | -------------------------------------------------------------------------------- /eigen-patched/failtest/map_nonconst_ctor_on_const_ptr_0.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(CV_QUALIFIER float *ptr){ 12 | Map m(ptr); 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/map_nonconst_ctor_on_const_ptr_1.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(CV_QUALIFIER float *ptr, DenseIndex size){ 12 | Map m(ptr, size); 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/map_nonconst_ctor_on_const_ptr_2.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(CV_QUALIFIER float *ptr, DenseIndex rows, DenseIndex cols){ 12 | Map m(ptr, rows, cols); 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/map_nonconst_ctor_on_const_ptr_3.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(CV_QUALIFIER float *ptr, DenseIndex rows, DenseIndex cols){ 12 | Map > m(ptr, rows, cols, InnerStride<2>()); 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/map_nonconst_ctor_on_const_ptr_4.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER 5 | #else 6 | #define CV_QUALIFIER const 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(const float *ptr, DenseIndex rows, DenseIndex cols){ 12 | Map > m(ptr, rows, cols, OuterStride<>(2)); 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/map_on_const_type_actually_const_0.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(float *ptr){ 12 | Map(ptr, 1, 1).coeffRef(0,0) = 1.0f; 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/map_on_const_type_actually_const_1.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(float *ptr){ 12 | Map(ptr).coeffRef(0) = 1.0f; 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/transpose_nonconst_ctor_on_const_xpr.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(CV_QUALIFIER Matrix3d &m){ 12 | Transpose t(m); 13 | } 14 | 15 | int main() {} 16 | -------------------------------------------------------------------------------- /eigen-patched/failtest/transpose_on_const_type_actually_const.cpp: -------------------------------------------------------------------------------- 1 | #include "../Eigen/Core" 2 | 3 | #ifdef EIGEN_SHOULD_FAIL_TO_BUILD 4 | #define CV_QUALIFIER const 5 | #else 6 | #define CV_QUALIFIER 7 | #endif 8 | 9 | using namespace Eigen; 10 | 11 | void foo(){ 12 | MatrixXf m; 13 | Transpose(m).coeffRef(0, 0) = 1.0f; 14 | } 15 | 16 | int main() {} 17 | -------------------------------------------------------------------------------- /eigen-patched/lapack/complex_double.cpp: -------------------------------------------------------------------------------- 1 | // This file is part of Eigen, a lightweight C++ template library 2 | // for linear algebra. 3 | // 4 | // Copyright (C) 2009-2011 Gael Guennebaud 5 | // 6 | // This Source Code Form is subject to the terms of the Mozilla 7 | // Public License v. 2.0. If a copy of the MPL was not distributed 8 | // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9 | 10 | #define SCALAR std::complex 11 | #define SCALAR_SUFFIX z 12 | #define SCALAR_SUFFIX_UP "Z" 13 | #define REAL_SCALAR_SUFFIX d 14 | #define ISCOMPLEX 1 15 | 16 | #include "cholesky.cpp" 17 | #include "lu.cpp" 18 | -------------------------------------------------------------------------------- /eigen-patched/lapack/complex_single.cpp: -------------------------------------------------------------------------------- 1 | // This file is part of Eigen, a lightweight C++ template library 2 | // for linear algebra. 3 | // 4 | // Copyright (C) 2009-2011 Gael Guennebaud 5 | // 6 | // This Source Code Form is subject to the terms of the Mozilla 7 | // Public License v. 2.0. If a copy of the MPL was not distributed 8 | // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9 | 10 | #define SCALAR std::complex 11 | #define SCALAR_SUFFIX c 12 | #define SCALAR_SUFFIX_UP "C" 13 | #define REAL_SCALAR_SUFFIX s 14 | #define ISCOMPLEX 1 15 | 16 | #include "cholesky.cpp" 17 | #include "lu.cpp" 18 | -------------------------------------------------------------------------------- /eigen-patched/lapack/double.cpp: -------------------------------------------------------------------------------- 1 | // This file is part of Eigen, a lightweight C++ template library 2 | // for linear algebra. 3 | // 4 | // Copyright (C) 2009-2011 Gael Guennebaud 5 | // 6 | // This Source Code Form is subject to the terms of the Mozilla 7 | // Public License v. 2.0. If a copy of the MPL was not distributed 8 | // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9 | 10 | #define SCALAR double 11 | #define SCALAR_SUFFIX d 12 | #define SCALAR_SUFFIX_UP "D" 13 | #define ISCOMPLEX 0 14 | 15 | #include "cholesky.cpp" 16 | #include "lu.cpp" 17 | #include "eigenvalues.cpp" 18 | -------------------------------------------------------------------------------- /eigen-patched/lapack/lapack_common.h: -------------------------------------------------------------------------------- 1 | // This file is part of Eigen, a lightweight C++ template library 2 | // for linear algebra. 3 | // 4 | // Copyright (C) 2010-2011 Gael Guennebaud 5 | // 6 | // This Source Code Form is subject to the terms of the Mozilla 7 | // Public License v. 2.0. If a copy of the MPL was not distributed 8 | // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9 | 10 | #ifndef EIGEN_LAPACK_COMMON_H 11 | #define EIGEN_LAPACK_COMMON_H 12 | 13 | #include "../blas/common.h" 14 | 15 | #define EIGEN_LAPACK_FUNC(FUNC,ARGLIST) \ 16 | extern "C" { int EIGEN_BLAS_FUNC(FUNC) ARGLIST; } \ 17 | int EIGEN_BLAS_FUNC(FUNC) ARGLIST 18 | 19 | typedef Eigen::Map > PivotsType; 20 | 21 | 22 | 23 | #endif // EIGEN_LAPACK_COMMON_H 24 | -------------------------------------------------------------------------------- /eigen-patched/lapack/single.cpp: -------------------------------------------------------------------------------- 1 | // This file is part of Eigen, a lightweight C++ template library 2 | // for linear algebra. 3 | // 4 | // Copyright (C) 2009-2011 Gael Guennebaud 5 | // 6 | // This Source Code Form is subject to the terms of the Mozilla 7 | // Public License v. 2.0. If a copy of the MPL was not distributed 8 | // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9 | 10 | #define SCALAR float 11 | #define SCALAR_SUFFIX s 12 | #define SCALAR_SUFFIX_UP "S" 13 | #define ISCOMPLEX 0 14 | 15 | #include "cholesky.cpp" 16 | #include "lu.cpp" 17 | #include "eigenvalues.cpp" 18 | -------------------------------------------------------------------------------- /eigen-patched/scripts/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | get_property(EIGEN_TESTS_LIST GLOBAL PROPERTY EIGEN_TESTS_LIST) 2 | configure_file(buildtests.in ${CMAKE_BINARY_DIR}/buildtests.sh @ONLY) 3 | 4 | configure_file(check.in ${CMAKE_BINARY_DIR}/check.sh COPYONLY) 5 | configure_file(debug.in ${CMAKE_BINARY_DIR}/debug.sh COPYONLY) 6 | configure_file(release.in ${CMAKE_BINARY_DIR}/release.sh COPYONLY) 7 | -------------------------------------------------------------------------------- /eigen-patched/scripts/buildtests.in: -------------------------------------------------------------------------------- 1 | #!/bin/bash 2 | 3 | if [[ $# != 1 || $1 == *help ]] 4 | then 5 | echo "usage: ./check regexp" 6 | echo " Builds tests matching the regexp." 7 | echo " The EIGEN_MAKE_ARGS environment variable allows to pass args to 'make'." 8 | echo " For example, to launch 5 concurrent builds, use EIGEN_MAKE_ARGS='-j5'" 9 | exit 0 10 | fi 11 | 12 | TESTSLIST="@EIGEN_TESTS_LIST@" 13 | targets_to_make=`echo "$TESTSLIST" | egrep "$1" | xargs echo` 14 | 15 | if [ -n "${EIGEN_MAKE_ARGS:+x}" ] 16 | then 17 | make $targets_to_make ${EIGEN_MAKE_ARGS} 18 | else 19 | make $targets_to_make 20 | fi 21 | exit $? 22 | 23 | -------------------------------------------------------------------------------- /eigen-patched/scripts/check.in: -------------------------------------------------------------------------------- 1 | #!/bin/bash 2 | # check : shorthand for make and ctest -R 3 | 4 | if [[ $# != 1 || $1 == *help ]] 5 | then 6 | echo "usage: ./check regexp" 7 | echo " Builds and runs tests matching the regexp." 8 | echo " The EIGEN_MAKE_ARGS environment variable allows to pass args to 'make'." 9 | echo " For example, to launch 5 concurrent builds, use EIGEN_MAKE_ARGS='-j5'" 10 | echo " The EIGEN_CTEST_ARGS environment variable allows to pass args to 'ctest'." 11 | echo " For example, with CTest 2.8, you can use EIGEN_CTEST_ARGS='-j5'." 12 | exit 0 13 | fi 14 | 15 | if [ -n "${EIGEN_CTEST_ARGS:+x}" ] 16 | then 17 | ./buildtests.sh "$1" && ctest -R "$1" ${EIGEN_CTEST_ARGS} 18 | else 19 | ./buildtests.sh "$1" && ctest -R "$1" 20 | fi 21 | exit $? 22 | -------------------------------------------------------------------------------- /eigen-patched/scripts/debug.in: -------------------------------------------------------------------------------- 1 | #!/bin/sh 2 | 3 | cmake -DCMAKE_BUILD_TYPE=Debug . 4 | -------------------------------------------------------------------------------- /eigen-patched/scripts/eigen_gen_docs: -------------------------------------------------------------------------------- 1 | #!/bin/sh 2 | 3 | # configuration 4 | # You should call this script with USER set as you want, else some default 5 | # will be used 6 | USER=${USER:-'orzel'} 7 | 8 | #ulimit -v 1024000 9 | 10 | # step 1 : build 11 | mkdir build -p 12 | (cd build && cmake .. && make doc) || { echo "make failed"; exit 1; } 13 | 14 | #step 2 : upload 15 | # (the '/' at the end of path is very important, see rsync documentation) 16 | rsync -az --no-p --delete build/doc/html/ $USER@ssh.tuxfamily.org:eigen/eigen.tuxfamily.org-web/htdocs/dox/ || { echo "upload failed"; exit 1; } 17 | 18 | #step 3 : fix the perm 19 | ssh $USER@ssh.tuxfamily.org 'chmod -R g+w /home/eigen/eigen.tuxfamily.org-web/htdocs/dox-devel' || { echo "perm failed"; exit 1; } 20 | 21 | echo "Uploaded successfully" 22 | 23 | -------------------------------------------------------------------------------- /eigen-patched/scripts/release.in: -------------------------------------------------------------------------------- 1 | #!/bin/sh 2 | 3 | cmake -DCMAKE_BUILD_TYPE=Release . 4 | -------------------------------------------------------------------------------- /eigen-patched/signature_of_eigen3_matrix_library: -------------------------------------------------------------------------------- 1 | This file is just there as a signature to help identify directories containing Eigen3. When writing a script looking for Eigen3, just look for this file. This is especially useful to help disambiguate with Eigen2... 2 | -------------------------------------------------------------------------------- /eigen-patched/test-patched: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/herrgahr/qling/22e56c4be0bbccb1d0437f610bfa37374b44b87f/eigen-patched/test-patched -------------------------------------------------------------------------------- /eigen-patched/test-patched.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | using Eigen::MatrixXd; 4 | int main() 5 | { 6 | MatrixXd m(2,2); 7 | m(0,0) = 3; 8 | m(1,0) = 2.5; 9 | m(0,1) = -1; 10 | m(1,1) = m(1,0) + m(0,1); 11 | std::cout << m << std::endl; 12 | } 13 | 14 | -------------------------------------------------------------------------------- /eigen-patched/test/eigen2/runtest.sh: -------------------------------------------------------------------------------- 1 | #!/bin/bash 2 | 3 | black='\E[30m' 4 | red='\E[31m' 5 | green='\E[32m' 6 | yellow='\E[33m' 7 | blue='\E[34m' 8 | magenta='\E[35m' 9 | cyan='\E[36m' 10 | white='\E[37m' 11 | 12 | if make test_$1 > /dev/null 2> .runtest.log ; then 13 | if ! ./test_$1 r20 > /dev/null 2> .runtest.log ; then 14 | echo -e $red Test $1 failed: $black 15 | echo -e $blue 16 | cat .runtest.log 17 | echo -e $black 18 | exit 1 19 | else 20 | echo -e $green Test $1 passed$black 21 | fi 22 | else 23 | echo -e $red Build of target $1 failed: $black 24 | echo -e $blue 25 | cat .runtest.log 26 | echo -e $black 27 | exit 1 28 | fi 29 | -------------------------------------------------------------------------------- /eigen-patched/test/runtest.sh: -------------------------------------------------------------------------------- 1 | #!/bin/bash 2 | 3 | black='\E[30m' 4 | red='\E[31m' 5 | green='\E[32m' 6 | yellow='\E[33m' 7 | blue='\E[34m' 8 | magenta='\E[35m' 9 | cyan='\E[36m' 10 | white='\E[37m' 11 | 12 | if ! ./$1 > /dev/null 2> .runtest.log ; then 13 | echo -e $red Test $1 failed: $black 14 | echo -e $blue 15 | cat .runtest.log 16 | echo -e $black 17 | exit 1 18 | else 19 | echo -e $green Test $1 passed$black 20 | fi 21 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | add_subdirectory(Eigen) 2 | add_subdirectory(doc EXCLUDE_FROM_ALL) 3 | if(EIGEN_LEAVE_TEST_IN_ALL_TARGET) 4 | add_subdirectory(test) # can't do EXCLUDE_FROM_ALL here, breaks CTest 5 | else() 6 | add_subdirectory(test EXCLUDE_FROM_ALL) 7 | endif() 8 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | set(Eigen_HEADERS AdolcForward BVH IterativeSolvers MatrixFunctions MoreVectorization AutoDiff AlignedVector3 Polynomials 2 | FFT NonLinearOptimization SparseExtra IterativeSolvers 3 | NumericalDiff Skyline MPRealSupport OpenGLSupport KroneckerProduct Splines 4 | ) 5 | 6 | install(FILES 7 | ${Eigen_HEADERS} 8 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen COMPONENT Devel 9 | ) 10 | 11 | add_subdirectory(src) 12 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/KroneckerProduct: -------------------------------------------------------------------------------- 1 | #ifndef EIGEN_KRONECKER_PRODUCT_MODULE_H 2 | #define EIGEN_KRONECKER_PRODUCT_MODULE_H 3 | 4 | #include "../../Eigen/Core" 5 | 6 | #include "../../Eigen/src/Core/util/DisableStupidWarnings.h" 7 | 8 | namespace Eigen { 9 | 10 | /** \ingroup Unsupported_modules 11 | * \defgroup KroneckerProduct_Module KroneckerProduct module 12 | * 13 | * This module contains an experimental Kronecker product implementation. 14 | * 15 | * \code 16 | * #include 17 | * \endcode 18 | */ 19 | 20 | } // namespace Eigen 21 | 22 | #include "src/KroneckerProduct/KroneckerTensorProduct.h" 23 | 24 | #include "../../Eigen/src/Core/util/ReenableStupidWarnings.h" 25 | 26 | #endif // EIGEN_KRONECKER_PRODUCT_MODULE_H 27 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/MoreVectorization: -------------------------------------------------------------------------------- 1 | #ifndef EIGEN_MOREVECTORIZATION_MODULE_H 2 | #define EIGEN_MOREVECTORIZATION_MODULE_H 3 | 4 | #include 5 | 6 | namespace Eigen { 7 | 8 | /** \ingroup Unsupported_modules 9 | * \defgroup MoreVectorization More vectorization module 10 | */ 11 | 12 | } 13 | 14 | #include "src/MoreVectorization/MathFunctions.h" 15 | 16 | #endif // EIGEN_MOREVECTORIZATION_MODULE_H 17 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/Skyline: -------------------------------------------------------------------------------- 1 | #ifndef EIGEN_SKYLINE_MODULE_H 2 | #define EIGEN_SKYLINE_MODULE_H 3 | 4 | 5 | #include "Eigen/Core" 6 | 7 | #include "Eigen/src/Core/util/DisableStupidWarnings.h" 8 | 9 | #include 10 | #include 11 | #include 12 | #include 13 | 14 | /** \ingroup Unsupported_modules 15 | * \defgroup Skyline_Module Skyline module 16 | * 17 | * 18 | * 19 | * 20 | */ 21 | 22 | #include "src/Skyline/SkylineUtil.h" 23 | #include "src/Skyline/SkylineMatrixBase.h" 24 | #include "src/Skyline/SkylineStorage.h" 25 | #include "src/Skyline/SkylineMatrix.h" 26 | #include "src/Skyline/SkylineInplaceLU.h" 27 | #include "src/Skyline/SkylineProduct.h" 28 | 29 | #include "Eigen/src/Core/util/ReenableStupidWarnings.h" 30 | 31 | #endif // EIGEN_SKYLINE_MODULE_H 32 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/AutoDiff/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_AutoDiff_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_AutoDiff_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen/src/AutoDiff COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/BVH/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_BVH_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_BVH_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen/src/BVH COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | ADD_SUBDIRECTORY(AutoDiff) 2 | ADD_SUBDIRECTORY(BVH) 3 | ADD_SUBDIRECTORY(FFT) 4 | ADD_SUBDIRECTORY(IterativeSolvers) 5 | ADD_SUBDIRECTORY(MatrixFunctions) 6 | ADD_SUBDIRECTORY(MoreVectorization) 7 | ADD_SUBDIRECTORY(NonLinearOptimization) 8 | ADD_SUBDIRECTORY(NumericalDiff) 9 | ADD_SUBDIRECTORY(Polynomials) 10 | ADD_SUBDIRECTORY(Skyline) 11 | ADD_SUBDIRECTORY(SparseExtra) 12 | ADD_SUBDIRECTORY(KroneckerProduct) 13 | ADD_SUBDIRECTORY(Splines) 14 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/FFT/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_FFT_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_FFT_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen/src/FFT COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/IterativeSolvers/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_IterativeSolvers_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_IterativeSolvers_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen/src/IterativeSolvers COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/KroneckerProduct/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_KroneckerProduct_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_KroneckerProduct_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen/src/KroneckerProduct COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/MatrixFunctions/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_MatrixFunctions_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_MatrixFunctions_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen/src/MatrixFunctions COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/MoreVectorization/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_MoreVectorization_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_MoreVectorization_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen/src/MoreVectorization COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/NonLinearOptimization/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_NonLinearOptimization_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_NonLinearOptimization_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen/src/NonLinearOptimization COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/NumericalDiff/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_NumericalDiff_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_NumericalDiff_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen/src/NumericalDiff COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/Polynomials/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Polynomials_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Polynomials_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen/src/Polynomials COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/Skyline/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Skyline_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Skyline_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen/src/Skyline COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/SparseExtra/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_SparseExtra_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_SparseExtra_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen/src/SparseExtra COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/Eigen/src/Splines/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | FILE(GLOB Eigen_Splines_SRCS "*.h") 2 | 3 | INSTALL(FILES 4 | ${Eigen_Splines_SRCS} 5 | DESTINATION ${INCLUDE_INSTALL_DIR}/unsupported/Eigen/src/Splines COMPONENT Devel 6 | ) 7 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/doc/CMakeLists.txt: -------------------------------------------------------------------------------- 1 | set_directory_properties(PROPERTIES EXCLUDE_FROM_ALL TRUE) 2 | 3 | add_subdirectory(examples) 4 | add_subdirectory(snippets) 5 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/doc/Overview.dox: -------------------------------------------------------------------------------- 1 | namespace Eigen { 2 | 3 | /** \mainpage Eigen's unsupported modules 4 | 5 | This is the API documentation for Eigen's unsupported modules. 6 | 7 | These modules are contributions from various users. They are provided "as is", without any support. 8 | 9 | Click on the \e Modules tab at the top of this page to get a list of all unsupported modules. 10 | 11 | Don't miss the official Eigen documentation. 12 | 13 | 14 | \defgroup Unsupported_modules Unsupported modules 15 | 16 | The unsupported modules are contributions from various users. They are 17 | provided "as is", without any support. Nevertheless, some of them are 18 | subject to be included in Eigen in the future. 19 | 20 | */ 21 | 22 | } 23 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/doc/examples/MatrixExponential.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | int main() 7 | { 8 | const double pi = std::acos(-1.0); 9 | 10 | MatrixXd A(3,3); 11 | A << 0, -pi/4, 0, 12 | pi/4, 0, 0, 13 | 0, 0, 0; 14 | std::cout << "The matrix A is:\n" << A << "\n\n"; 15 | std::cout << "The matrix exponential of A is:\n" << A.exp() << "\n\n"; 16 | } 17 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/doc/examples/MatrixFunction.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | std::complex expfn(std::complex x, int) 7 | { 8 | return std::exp(x); 9 | } 10 | 11 | int main() 12 | { 13 | const double pi = std::acos(-1.0); 14 | 15 | MatrixXd A(3,3); 16 | A << 0, -pi/4, 0, 17 | pi/4, 0, 0, 18 | 0, 0, 0; 19 | 20 | std::cout << "The matrix A is:\n" << A << "\n\n"; 21 | std::cout << "The matrix exponential of A is:\n" 22 | << A.matrixFunction(expfn) << "\n\n"; 23 | } 24 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/doc/examples/MatrixLogarithm.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | int main() 7 | { 8 | using std::sqrt; 9 | MatrixXd A(3,3); 10 | A << 0.5*sqrt(2), -0.5*sqrt(2), 0, 11 | 0.5*sqrt(2), 0.5*sqrt(2), 0, 12 | 0, 0, 1; 13 | std::cout << "The matrix A is:\n" << A << "\n\n"; 14 | std::cout << "The matrix logarithm of A is:\n" << A.log() << "\n"; 15 | } 16 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/doc/examples/MatrixSine.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | int main() 7 | { 8 | MatrixXd A = MatrixXd::Random(3,3); 9 | std::cout << "A = \n" << A << "\n\n"; 10 | 11 | MatrixXd sinA = A.sin(); 12 | std::cout << "sin(A) = \n" << sinA << "\n\n"; 13 | 14 | MatrixXd cosA = A.cos(); 15 | std::cout << "cos(A) = \n" << cosA << "\n\n"; 16 | 17 | // The matrix functions satisfy sin^2(A) + cos^2(A) = I, 18 | // like the scalar functions. 19 | std::cout << "sin^2(A) + cos^2(A) = \n" << sinA*sinA + cosA*cosA << "\n\n"; 20 | } 21 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/doc/examples/MatrixSinh.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | int main() 7 | { 8 | MatrixXf A = MatrixXf::Random(3,3); 9 | std::cout << "A = \n" << A << "\n\n"; 10 | 11 | MatrixXf sinhA = A.sinh(); 12 | std::cout << "sinh(A) = \n" << sinhA << "\n\n"; 13 | 14 | MatrixXf coshA = A.cosh(); 15 | std::cout << "cosh(A) = \n" << coshA << "\n\n"; 16 | 17 | // The matrix functions satisfy cosh^2(A) - sinh^2(A) = I, 18 | // like the scalar functions. 19 | std::cout << "cosh^2(A) - sinh^2(A) = \n" << coshA*coshA - sinhA*sinhA << "\n\n"; 20 | } 21 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/doc/examples/MatrixSquareRoot.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | 6 | int main() 7 | { 8 | const double pi = std::acos(-1.0); 9 | 10 | MatrixXd A(2,2); 11 | A << cos(pi/3), -sin(pi/3), 12 | sin(pi/3), cos(pi/3); 13 | std::cout << "The matrix A is:\n" << A << "\n\n"; 14 | std::cout << "The matrix square root of A is:\n" << A.sqrt() << "\n\n"; 15 | std::cout << "The square of the last matrix is:\n" << A.sqrt() * A.sqrt() << "\n"; 16 | } 17 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/doc/examples/PolynomialUtils1.cpp: -------------------------------------------------------------------------------- 1 | #include 2 | #include 3 | 4 | using namespace Eigen; 5 | using namespace std; 6 | 7 | int main() 8 | { 9 | Vector4d roots = Vector4d::Random(); 10 | cout << "Roots: " << roots.transpose() << endl; 11 | Eigen::Matrix polynomial; 12 | roots_to_monicPolynomial( roots, polynomial ); 13 | cout << "Polynomial: "; 14 | for( int i=0; i<4; ++i ){ cout << polynomial[i] << ".x^" << i << "+ "; } 15 | cout << polynomial[4] << ".x^4" << endl; 16 | Vector4d evaluation; 17 | for( int i=0; i<4; ++i ){ 18 | evaluation[i] = poly_eval( polynomial, roots[i] ); } 19 | cout << "Evaluation of the polynomial at the roots: " << evaluation.transpose(); 20 | } 21 | -------------------------------------------------------------------------------- /eigen-patched/unsupported/test/FFT.cpp: -------------------------------------------------------------------------------- 1 | #define test_FFTW test_FFT 2 | #include "FFTW.cpp" 3 | -------------------------------------------------------------------------------- /qling/jiteventlistener.cpp: -------------------------------------------------------------------------------- 1 | #include "jiteventlistener.h" 2 | #include 3 | #include "llvm/IR/Function.h" 4 | #include "llvm/Support/raw_ostream.h" 5 | 6 | void JitEventListener::NotifyFunctionEmitted(const llvm::Function &f, void *, size_t, const llvm::JITEventListener::EmittedFunctionDetails &) 7 | { 8 | if(f.getName().find("__cling_Un1Qu3")!=llvm::StringRef::npos) 9 | emit aboutToExecWrappedFunction(); 10 | } 11 | 12 | void JitEventListener::NotifyFreeingMachineCode(void *) 13 | { 14 | } 15 | -------------------------------------------------------------------------------- /qt-hack/QtCore/QAtomicInt: -------------------------------------------------------------------------------- 1 | #include "qatomic.h" 2 | -------------------------------------------------------------------------------- /qt-hack/QtCore/QAtomicPointer: -------------------------------------------------------------------------------- 1 | #include "qatomic.h" 2 | -------------------------------------------------------------------------------- /setupenv.sh: -------------------------------------------------------------------------------- 1 | #!/bin/sh 2 | #source this file to setup the environment for building and running qling 3 | 4 | #change to match your install 5 | export LLVM_INSTALL=/home/thomas/opt/llvm-debug 6 | 7 | 8 | export LD_LIBRARY_PATH=${LLVM_INSTALL}/lib 9 | export PATH=${LLVM_INSTALL}/bin:$PATH 10 | --------------------------------------------------------------------------------