aboutsummaryrefslogtreecommitdiffhomepage
path: root/Eigen/src/Cholesky
Commit message (Collapse)AuthorAge
...
* Propagate all five matrix template parameters to members and temporaries of ↵Gravatar Adolfo Rodriguez Tsouroukdissian2010-03-08
| | | | | | | decomposition classes. One particular advantage of this is that decomposing matrices with max sizes known at compile time will not allocate. NOTE: The ComplexEigenSolver class currently _does_ allocate (line 135 of Eigenvalues/ComplexEigenSolver.h), but the reason appears to be in the implementation of matrix-matrix products, and not in the decomposition itself. The nomalloc unit test has been extended to verify that decompositions do not allocate when max sizes are specified. There are currently two workarounds to prevent the test from failing (see comments in test/nomalloc.cpp), both of which are related to matrix products that allocate on the stack.
* add reconstructedMatrix() to LLT, and LUsGravatar Gael Guennebaud2010-02-24
| | | | | => they show that some improvements have still to be done for permutations, tr*tr, trapezoidal matrices
* * fix LDLT's default ctor useGravatar Gael Guennebaud2010-02-24
| | | | * add a reconstructedMatrix() function to LDLT for debug purpose
* actually, this is not even meant to be a termination criterion. so the ↵Gravatar Benoit Jacob2010-02-23
| | | | proper fix is this.
* LDLT:Gravatar Benoit Jacob2010-02-23
| | | | | * fix bug thanks to Ben Goodrich: we were terminating at the wrong place, leaving some matrix coefficients with wrong values. * don't use Higham's formula here: we're not trying to be rank-revealing.
* compilation fix in ldlt() for non matrix typesGravatar Gael Guennebaud2010-02-21
|
* Renamed PlainMatrixType to PlainObject (Array != Matrix).Gravatar Hauke Heibel2010-02-20
| | | | Renamed ReturnByValue::ReturnMatrixType ReturnByValue::ReturnType (again, Array != Matrix).
* * move dummy_precision and epsilon to NumTraitsGravatar Gael Guennebaud2010-02-10
| | | | * make NumTraits inherits std::numeric_limits
* Introduced NestParentByRefBit and NestByRefBit - this should fix temporaries ↵Gravatar Hauke Heibel2010-02-06
| | | | | | related to nested products. Fixed a few typos and a few warnings.
* remove the Triangular suffix to Upper, Lower, UnitLower, etc,Gravatar Gael Guennebaud2010-01-07
| | | | and remove the respective bit flags
* merge and add start/end to Eigen2SupportGravatar Gael Guennebaud2010-01-05
|\
| * Big renaming:Gravatar Benoit Jacob2010-01-04
| | | | | | | | | | | | start ---> head end ---> tail Much frustration with sed syntax. Need to learn perl some day.
* | merge with default branchGravatar Gael Guennebaud2009-12-22
|\|
| * Warning 4512 (assignment operators could not be generated) is now simply ↵Gravatar Hauke Heibel2009-12-14
| | | | | | | | | | | | disabled. All unimplemented assignment operators have been removed.
| * Correct type of ei_solve_retval<LLT,...>::operator=Gravatar Jitse Niesen2009-12-12
| |
| * Added to possibility to compile unit tests at maximum warning level.Gravatar Hauke Heibel2009-12-12
| | | | | | | | Silenced (amongst others) many conversion related warnings.
| * Removed NestByValue dependency from Cholesky, Eigenvalues, LU and QR.Gravatar Hauke Heibel2009-12-01
| |
* | Another big refactoring change:Gravatar Gael Guennebaud2009-11-18
|/ | | | | | * add a new Eigen2Support module including Cwise, Flagged, and some other deprecated stuff * add a few cwiseXxx functions * adapt a few modules to use cwiseXxx instead of the .cwise() prefix
* last round of changes, mainly to return derived types instead of base types, ↵Gravatar Benoit Jacob2009-11-09
| | | | and fix various compilation issues
* simplifications in the ei_solve_impl system, factor out some boilerplate codeGravatar Benoit Jacob2009-11-08
|
* move cholesky to ei_xxx_return_valueGravatar Benoit Jacob2009-11-03
|
* *port the Cholesky module to the new solve() APIGravatar Benoit Jacob2009-10-29
| | | | *improve documentation
* big huge changes, so i dont remember everything.Gravatar Benoit Jacob2009-10-28
| | | | | | | | | | * renaming, e.g. LU ---> FullPivLU * split tests framework: more robust, e.g. dont generate empty tests if a number is skipped * make all remaining tests use that splitting, as needed. * Fix 4x4 inversion (see stable branch) * Transform::inverse() and geo_transform test : adapt to new inverse() API, it was also trying to instantiate inverse() for 3x4 matrices. * CMakeLists: more robust regexp to parse the version number * misc fixes in unit tests
* * Introduce make targets btest (build tests), blas (build blas lib), demos ↵Gravatar Benoit Jacob2009-09-27
| | | | | | | | | (build demos). * remove EIGEN_BUILD_TESTS and siblings * add summary at the end of cmake run, hopefully not too verbose * fix build of quaternion demo * kill remnants of old binary library option
* Adaptions from .lazy() towards .noalias().Gravatar Hauke Heibel2009-08-31
| | | | Added missing casts.
* in all decs, make the compute() methods return *thisGravatar Benoit Jacob2009-08-15
| | | | (implements feature request #18)
* machine_epsilon -> epsilon as wrapper around numeric_traitsGravatar Benoit Jacob2009-08-14
|
* fix a couple of compilations issuesGravatar Gael Guennebaud2009-08-06
|
* bugfix in trsmGravatar Gael Guennebaud2009-08-02
|
* add blocked LLT, and bugfix in trsm assertsGravatar Gael Guennebaud2009-08-01
|
* synch with main branchGravatar Gael Guennebaud2009-07-28
|\
| * minor compilation fixes for Sun CC and ICCGravatar Gael Guennebaud2009-07-20
| |
* | * change the nesting order of adjoint_return_type toGravatar Gael Guennebaud2009-07-07
| | | | | | | | | | | | 1 - make it easier to catch conjugate expressions 2 - make sure there is no unecessary copy (we had NestByValue<Derived> which seems to be very bad) * update eigensolver wrt recent changes
* | * take advantage of new possibilies in LLT (mat -= product)Gravatar Gael Guennebaud2009-07-07
| | | | | | | | | | * fix Block::operator+= product which was not optimized * fix some compilation issues
* | add a generic mechanism to copy a special matrix to a dense matrix so thatGravatar Gael Guennebaud2009-07-07
| | | | | | | | | | we don't need to add other specialization of MatrixBase::operator=, Matrix::=, and Matrix::Matrix(...)
* | * bybye Part, welcome TriangularView and SelfAdjointView.Gravatar Gael Guennebaud2009-07-06
|/ | | | | | | | * move solveTriangular*() to TriangularView::solve*() * move .llt() to SelfAdjointView * add a high level wrapper to the efficient selfadjoint * vector product * improve LLT so that we can specify which triangular part is meaningless => there are still many things to do (doc, cleaning, improve the matrix products, etc.)
* fix #14: make llt::solve() and also ldlt::solve() work with uninitialized resultGravatar Benoit Jacob2009-06-19
|
* mergeGravatar Benoit Jacob2009-05-22
|\
| * remove sentence "Eigen itself is part of the KDE project."Gravatar Benoit Jacob2009-05-22
| | | | | | | | it never made very precise sense. but now does it still make any?
* | Oops, here the actual LLT and LDLT patch.Gravatar Hauke Heibel2009-05-22
|/
* constant Diagonal ---> DiagonalBitsGravatar Benoit Jacob2009-05-10
| | | | | | introduce ei_is_diagonal to check for it DiagonalCoeffs ---> Diagonal and allow Index to by Dynamic -> add MatrixBase::diagonal(int) with unittest and doc
* More Cholesky fixes.Gravatar Benoit Jacob2009-04-01
| | | | | | | | | * Cholesky decs are NOT rank revealing so remove all the rank/isPositiveDefinite etc stuff. * fix bug in LLT: s/return/continue/ * introduce machine_epsilon constants, they are actually needed for Higman's formula determining the cutoff in Cholesky. Btw fix the page reference to his book (chat with Keir). * solve methods always return true, since this isn't a rank revealing dec. Actually... they already did always return true!! Now it's explicit. * updated dox and unit-test
* add adjointInPlace() and add documentation warnings on adjoint() and ↵Gravatar Benoit Jacob2009-03-31
| | | | transpose() about aliasing effects.
* Many improvements in LLT and LDLT:Gravatar Benoit Jacob2009-03-30
| | | | | | | | | * in LDLT, support the negative semidefinite case * fix bad floating-point comparisons, improves greatly the accuracy of methods like isPositiveDefinite() and rank() * simplifications * identify (but not resolve) bug: claim that only triangular part is used, is inaccurate * expanded unit-tests
* * allows fixed size matrix with size==0 (via a specialization ofGravatar Gael Guennebaud2009-03-23
| | | | | | MatrixStorage returning a null pointer). For instance this is very useful to make Tridiagonalization compile for 1x1 matrices * fix LLT and eigensolver for 1x1 matrix
* Add COMPONENT DevelGravatar Laurent Montel2009-02-23
|
* s/cholesky/llt in precompiled lib and BTLGravatar Gael Guennebaud2009-02-06
|
* bugfix in LDLt for size==1Gravatar Gael Guennebaud2009-02-04
|
* Add full pivoting to LDLT decomposition.Gravatar Keir Mierle2009-02-03
|
* LLT: makes the non positive definite test less strict, but we still needGravatar Gael Guennebaud2009-01-27
| | | | something better.
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-04 03:56:59 +0000