Commit message (Collapse) | Author | Age | |
---|---|---|---|
* | Handle some Index to int conversions in BLAS/LAPACK support. | Gael Guennebaud | 2016-05-26 |
| | |||
* | Prevent crash in CompleteOrthogonalDecomposition if object was default ↵ | Rasmus Munk Larsen | 2016-04-21 |
| | | | | constructed. | ||
* | Removed executable bit from header files | Benoit Steiner | 2016-03-23 |
| | |||
* | Rename inverse -> pseudoInverse. | Rasmus Munk Larsen | 2016-02-10 |
| | |||
* | Enable inverse() method for computing pseudo-inverse. | Rasmus Munk Larsen | 2016-02-09 |
| | |||
* | Make applyZAdjointOnTheLeftInPlace protected. | Rasmus Munk Larsen | 2016-02-08 |
| | |||
* | Add missing calls to tests of COD. | Rasmus Munk Larsen | 2016-02-08 |
| | | | | Fix a few mistakes in 3.2 -> 3.3 port. | ||
* | Nevermind. | Rasmus Munk Larsen | 2016-02-06 |
| | |||
* | Add my name to copyright notice in ColPivHouseholder.h, mostly for previous ↵ | Rasmus Munk Larsen | 2016-02-06 |
| | | | | work on stable norm downdate formula. | ||
* | Implement complete orthogonal decomposition in Eigen. | Rasmus Munk Larsen | 2016-02-06 |
| | |||
* | Fix bad line break. Don't repeat Kahan matrix test since it is deterministic. | Rasmus Munk Larsen | 2016-02-03 |
| | |||
* | Make the array of directly compute column norms a member to avoid allocation ↵ | Rasmus Munk Larsen | 2016-02-03 |
| | | | | in computeInPlace. | ||
* | Change Eigen's ColPivHouseholderQR to use the numerically stable norm ↵ | Rasmus Munk Larsen | 2016-01-28 |
| | | | | | | | | | | | | | | | | | | | | | | | | downdate formula from http://www.netlib.org/lapack/lawnspdf/lawn176.pdf, which has been used in LAPACK's xGEQPF and xGEQP3 since 2006. With the old formula, the code chooses the wrong pivots and fails to correctly determine rank on graded matrices. This change also adds additional checks for non-increasing diagonal in R11 to existing unit tests, and adds a new unit test with the Kahan matrix, which consistently fails for the original code. Benchmark timings on Intel(R) Xeon(R) CPU E5-1650 v3 @ 3.50GHz. Code compiled with AVX & FMA. I just ran on square matrices of 3 difference sizes. Benchmark Time(ns) CPU(ns) Iterations ------------------------------------------------------- Before: BM_EigencolPivQR/64 53677 53627 12890 BM_EigencolPivQR/512 15265408 15250784 46 BM_EigencolPivQR/4k 15403556228 15388788368 2 After (non-vectorized version): Benchmark Time(ns) CPU(ns) Iterations Degradation -------------------------------------------------------------------- BM_EigencolPivQR/64 63736 63669 10844 18.5% BM_EigencolPivQR/512 16052546 16037381 43 5.1% BM_EigencolPivQR/4k 15149263620 15132025316 2 -2.0% Performance-wise there seems to be a ~18.5% degradation for small (64x64) matrices, probably due to the cost of more O(min(m,n)^2) sqrt operations that are not needed for the unstable formula. | ||
* | Fix numerous doxygen shortcomings, and workaround some clang -Wdocumentation ↵ | Gael Guennebaud | 2016-01-01 |
| | | | | warnings | ||
* | Fix MKL compilation issue | Gael Guennebaud | 2015-12-11 |
| | |||
* | Fix compilation of MKL support. | Gael Guennebaud | 2015-12-11 |
| | |||
* | Generalize matrix ctor and compute() method of dense decomposition to 1) ↵ | Gael Guennebaud | 2015-09-07 |
| | | | | limit temporaries, 2) forward expressions to nested decompositions, 3) fix ambiguous ctor instanciation for square decomposition | ||
* | bug #949: add static assertion for incompatible scalar types in dense ↵ | Gael Guennebaud | 2015-03-13 |
| | | | | end-user decompositions. | ||
* | New scoring functor to select the pivot. | Marc Glisse | 2015-03-03 |
| | | | | This is can be useful for non-floating point scalars, where choosing the biggest element is generally not the best choice. | ||
* | Merged in chtz/eigen-indexconversion (pull request PR-92) | Gael Guennebaud | 2015-02-16 |
|\ | | | | | | | | | | | | | | | | | | | | | | | bug #877, bug #572: Get rid of Index conversion warnings, summary of changes: - Introduce a global typedef Eigen::Index making Eigen::DenseIndex and AnyExpr<>::Index deprecated (default is std::ptrdiff_t). - Eigen::Index is used throughout the API to represent indices, offsets, and sizes. - Classes storing an array of indices uses the type StorageIndex to store them. This is a template parameter of the class. Default is int. - Methods that *explicitly* set or return an element of such an array take or return a StorageIndex type. In all other cases, the Index type is used. | ||
| * | Remove deprecated usage of expr::Index. | Gael Guennebaud | 2015-02-16 |
| | | |||
* | | bug #952: Missing \endcode made doxygen fail to build ColPivHouseholderQR | Christoph Hertzberg | 2015-02-15 |
| | | |||
| * | Index refactoring: StorageIndex must be used for storage only (and locally ↵ | Gael Guennebaud | 2015-02-13 |
| | | | | | | | | when it make sense). In all other cases use the global Index type. | ||
| * | Merge Index-refactoring branch with default, fix PastixSupport, remove some ↵ | Gael Guennebaud | 2015-02-13 |
| |\ | |/ |/| | | | useless typedefs | ||
* | | bug #941: fix accuracy issue in ColPivHouseholderQR, do not stop ↵ | Gael Guennebaud | 2015-01-30 |
| | | | | | | | | decomposition on a small pivot | ||
| * | bug #877, bug #572: Introduce a global Index typedef. Rename Sparse*::Index ↵ | Christoph Hertzberg | 2014-12-04 |
|/ | | | | to StorageIndex, make Dense*::StorageIndex an alias to DenseIndex. Overall this commit gets rid of all Index conversion warnings. | ||
* | quieted more g++ warnings of the form: warning: typedef XXX locally defined ↵ | Mark Borgerding | 2014-10-16 |
| | | | | but not used [-Wunused-local-typedefs] | ||
* | Make constructors explicit if they could lead to unintended implicit conversion | Christoph Hertzberg | 2014-09-23 |
| | |||
* | Remove deprecated code not used by evaluators | Gael Guennebaud | 2014-09-18 |
| | |||
* | Optimization: enable cache-efficient application of HouseholderSequence. | Gael Guennebaud | 2014-09-04 |
| | |||
* | Refactoring of sparse solvers through a SparseSolverBase class and usage of ↵ | Gael Guennebaud | 2014-09-01 |
| | | | | the Solve<> expression. Introduce a SolveWithGuess expression on top of Solve. | ||
* | merge default and evaluator branches | Gael Guennebaud | 2014-03-12 |
|\ | |||
* | | Port QR module to Solve/Inverse | Gael Guennebaud | 2014-03-11 |
| | | |||
| * | Make pivoting HouseholderQR compatible with custom scalar types | Christoph Hertzberg | 2014-02-25 |
|/ | |||
* | Minor corrections in QR docs. | Jitse Niesen | 2013-12-31 |
| | |||
* | Applied patch from Richard JW Roberts, resolving bug #704 | Christoph Hertzberg | 2013-12-21 |
| | |||
* | Fix FullPivHouseholderQR ctors for non squared fixed size matrix types | Gael Guennebaud | 2013-11-19 |
| | |||
* | merge with default branch | Gael Guennebaud | 2013-11-05 |
|\ | |||
* | | Add a rank method with threshold control to JacobiSVD, and make solve uses ↵ | Gael Guennebaud | 2013-11-01 |
| | | | | | | | | it to return the minimal norm solution for rank-deficient problems | ||
* | | Fix changeset 2702788da71303afa71e53dfa6542a560115425a | Gael Guennebaud | 2013-11-01 |
| | | | | | | | | for fixed size matrices | ||
* | | Fix bug #678: vectors of row and columns transpositions were not properly ↵ | Gael Guennebaud | 2013-10-29 |
| | | | | | | | | resized in FullPivQR | ||
* | | Make FullPivHouseholderQR::solve returns the least-square solution instead ↵ | Gael Guennebaud | 2013-08-20 |
| | | | | | | | | of aborting if no exact solution exist | ||
| * | Disable some shortcuts with nvcc | Gael Guennebaud | 2013-07-31 |
|/ | |||
* | Fix various scalar type conversion warnings. | Gael Guennebaud | 2013-07-12 |
| | |||
* | Fix HouseholderSequence::conjugate() and ::adjoint() and add respective unit ↵ | Gael Guennebaud | 2013-06-17 |
| | | | | tests. | ||
* | Fix bug #314: move remaining math functions from internal to numext namespace | Gael Guennebaud | 2013-06-10 |
| | |||
* | Fix bug #599: add missing documentation of some members in QR module. | Gael Guennebaud | 2013-06-10 |
| | |||
* | Clean source code and unit tests with respect to -Wunused-local-typedefs | Gael Guennebaud | 2013-04-10 |
| | |||
* | Add assertion on the input matrix size in factorizations relying on ↵ | Gael Guennebaud | 2013-02-26 |
| | | | | permutations of 32bits int | ||
* | Remove superfluous cast. | Gael Guennebaud | 2013-02-24 |
| |