aboutsummaryrefslogtreecommitdiffhomepage
Commit message (Collapse)AuthorAge
* doc fixes, and extended Basic Linear Algebra and Reductions sectionsGravatar Gael Guennebaud2008-08-20
|
* some updated of the quick start guideGravatar Gael Guennebaud2008-08-20
|
* * it's \returns not \ReturnsGravatar Benoit Jacob2008-08-20
| | | | | | | | | * add some explanations in the typedefs page * expand a bit the new QuickStartGuide. Some placeholders (not a pb since it's not even yet linked to from other pages). The point I want to make is that it's super important to have fully compilable short programs (even with compile instructions for the first one) not just small snippets, at least at the beginning. Let's start with examples of compilable programs.
* fix a few typosGravatar Gael Guennebaud2008-08-20
|
* * Added .all() and .any() members to PartialReduxGravatar Gael Guennebaud2008-08-20
| | | | | * Bug fixes in euler angle snippet, Assign and MapBase * Started a "quick start guide" (draft state)
* missing SVD header fileGravatar Gael Guennebaud2008-08-19
|
* Added a SVD module:Gravatar Gael Guennebaud2008-08-19
| | | | | | | - the decompostion code has been adfapted from JAMA - handles non square matrices of size MxN with M>=N - does not work for complex matrices - includes a solver where the parts corresponding to zero singular values are set to zero
* * fix bug found by Boudewijn Rempt: no CMakeLists in arch/ subdirGravatar Benoit Jacob2008-08-19
| | | | * fix warning in SolveTriangular
* Various compilation fixes for MSVC 9. All tests compile but someGravatar Gael Guennebaud2008-08-19
| | | | still fail at runtime in ei_aligned_free() (even without vectorization).
* * doc improvements in Cwise and PartialRedux:Gravatar Benoit Jacob2008-08-19
| | | | | | | | | - 33 new snippets - unfuck doxygen output in Cwise (issues with function macros) - more see-also links from outside, making Cwise more discoverable * rename matrixNorm() to operatorNorm(). There are many matrix norms (the L2 is another one) but only one is called the operator norm. Risk of confusion with keyword operator is not too scary after all.
* * revert the previous interface change in solveTriangular (pointer vs reference)Gravatar Gael Guennebaud2008-08-18
| | | | | | | | | * remove the cast operators in the Geometry module: they are replaced by constructors and new operator= in Matrix * extended the operations supported by Rotation2D * rewrite in solveTriangular: - merge the Upper and Lower specializations - big optimization of the path for row-major triangular matrices
* only use alloca on linux. Use malloc on other platforms. Needed for mingwGravatar Benoit Jacob2008-08-18
| | | | compatibility.
* *add operator()(int) for vectors, synonymous to operator[](int).Gravatar Benoit Jacob2008-08-18
| | | | | | I don't see any reason not to allow it, it doesn't add much code, and it makes porting from eigen1 easier. *expand tests/basicstuff to first test coefficient access methods
* fix compilation of unit-tests when cmake is run on kdesupport/ parentGravatar Benoit Jacob2008-08-17
| | | | directory
* *make Eigen2 install to prefix/include/eigen2/EigenGravatar Benoit Jacob2008-08-16
| | | | *Add FindEigen2.cmake module in kdelibs
* bugfix in CholeskyWithoutSquareRoot::solve found by Timothy HunterGravatar Gael Guennebaud2008-08-16
|
* updated tutorialGravatar Benoit Jacob2008-08-15
|
* Copy eigen2 into kdesupport, alongside eigen1 for now, co-installable,Gravatar Benoit Jacob2008-08-14
| | | | will remove eigen1 once no code uses it anymore (soon).
* change solveTriangularInPlace() to take a pointer as input (as discussed on ↵Gravatar Gael Guennebaud2008-08-12
| | | | | | IRC). extended the documentation of the triangular solver.
* last small fixes, this is alpha6, eigen2 is now ready for eigen1 apps toGravatar Benoit Jacob2008-08-12
| | | | port to.
* Complete LU documentationGravatar Benoit Jacob2008-08-11
|
* * add Regression module, from eigen1, improved, with doc and unit-test.Gravatar Benoit Jacob2008-08-11
| | | | | | | | | * fix .normalized() so that Random().normalized() works; since the return type became complicated to write down i just let it return an actual vector, perhaps not optimal. * add Sparse/CMakeLists.txt. I suppose that it was intentional that it didn't have CMakeLists, but in <=2.0 releases I'll just manually remove Sparse.
* Renamed allowAligned() => forceAligned() and added the constants ForceAlignedGravatar Gael Guennebaud2008-08-09
| | | | | and AsRequested for the PacketAccess parameter of MapBase. Updated respective documentation.
* renamed inverseProduct => solveTriangularGravatar Gael Guennebaud2008-08-09
|
* * add LU unit-test. Seems like we have very good numerical stability!Gravatar Benoit Jacob2008-08-09
| | | | * some cleanup, and grant me a copyright line on the determinant test.
* * Big change in Block and Map:Gravatar Gael Guennebaud2008-08-09
| | | | | | | | | | | | | | - added a MapBase base xpr on top of which Map and the specialization of Block are implemented - MapBase forces both aligned loads (and aligned stores, see below) in expressions such as "x.block(...) += other_expr" * Significant vectorization improvement: - added a AlignedBit flag meaning the first coeff/packet is aligned, this allows to not generate extra code to deal with the first unaligned part - removed all unaligned stores when no unrolling - removed unaligned loads in Sum when the input as the DirectAccessBit flag * Some code simplification in CacheFriendly product * Some minor documentation improvements
* * reimplement the general case of inverse() on top of LU. Advantages:Gravatar Benoit Jacob2008-08-09
| | | | | | | | | | | | - removes much code - 2.5x faster (even though LU uses complete pivoting contrary to what inverse used to do!) - there _were_ numeric stability problems with the partial pivoting approach of inverse(), with 200x200 matrices inversion failed almost half of the times (overflow). Now these problems are solved thanks to complete pivoting. * remove some useless stuff in LU
* *implement LU solver (solves any rectangular system)Gravatar Benoit Jacob2008-08-09
| | | | | | *in test/CMakeLists : modify EI_ADD_TEST so that 2nd argument is additional compiler flags. used to add -O2 to test_product_large so it doesn't take forever.
* * fix bug in SwapWrapper : store the wrapped expression by referenceGravatar Benoit Jacob2008-08-09
| | | | | | * optimize setIdentity: when the matrix is large enough it is better to setZero() and overwrite the diagonal * start of LU solver, disabled for now
* forgot to remove that old codeGravatar Benoit Jacob2008-08-07
|
* - add kernel computation using the triangular solverGravatar Benoit Jacob2008-08-07
| | | | | | - take advantage of the fact that our LU dec sorts the eigenvalues of U in decreasing order - add meta selector in determinant
* LU: remove partial-pivoting path (moderately useful since it's doesGravatar Benoit Jacob2008-08-07
| | | | | | | | not allow to easily get the rank), fix a bug (which could have been triggered by matrices having coefficients of very different magnitudes). Part: add an assert to prevent hard to find bugs Swap: update comments
* New implementation of Swap as discussed, reusing Assign. Makes LU runGravatar Benoit Jacob2008-08-05
| | | | 10% faster overall.
* introduce copyCoeff and copyPacket methods in MatrixBase, used byGravatar Benoit Jacob2008-08-05
| | | | | Assign, in preparation for new Swap impl reusing Assign code. remove last remnant of old Inverse class in Transform.
* Add partial pivoting runtime option to LU.Gravatar Benoit Jacob2008-08-05
| | | | | | | Note: in fact, inverse() always uses partial pivoting because the algo currently used doesn't make sense with complete pivoting. No num stability issue so far even with size 200x200. If there is any problem we can of course reimplement inverse on top of LU.
* further big perf improvement in InverseGravatar Benoit Jacob2008-08-04
|
* big performance improvement in inverse and LUGravatar Benoit Jacob2008-08-04
|
* Add a LU decomposition action in BTL and various cleaning in BTL. For instanceGravatar Gael Guennebaud2008-08-04
| | | | | | all per plot settings have been moved to a single file, go_mean now takes an optional second argument "tiny" to generate plots for tiny matrices, and output of comparison information wrt to previous benchs (if any).
* * LU decomposition, supporting all rectangular matrices, with fullGravatar Benoit Jacob2008-08-04
| | | | | | | | | | pivoting for better numerical stability. For now the only application is determinant. * New determinant unit-test. * Disable most of Swap.h for now as it makes LU fail (mysterious). Anyway Swap needs a big overhaul as proposed on IRC. * Remnants of old class Inverse removed. * Some warnings fixed.
* fix two perf issues in product.Gravatar Gael Guennebaud2008-08-03
| | | | | fix positive definite test in Cholesky. remove #include <cstring> in CoreDeclaration.
* fix compile errors with gcc 4.3: unresolved func call toGravatar Benoit Jacob2008-08-03
| | | | ei_cache_friendly_product, and undeclared memcpy
* Added a ei_palign function align a packet from two others.Gravatar Gael Guennebaud2008-08-03
| | | | | This allows much faster code dealing with unligned as in the updated matrix-vector product functions.
* Optimizations:Gravatar Gael Guennebaud2008-08-01
| | | | | | | | * faster matrix-matrix and matrix-vector products (especially for not aligned cases) * faster tridiagonalization (make it using our matrix-vector impl.) Others: * fix Flags of Map * split the test_product to two smaller ones
* removed the packet specializations of some functorsGravatar Gael Guennebaud2008-07-31
| | | | (GCC generates better code without those "optimizations")
* Several compilation fixes for MSVC and NVCC, basically:Gravatar Gael Guennebaud2008-07-29
| | | | | | | | - added explicit enum to int conversion where needed - if a function is not defined as declared and the return type is "tricky" then the type must be typedefined somewhere. A "tricky return type" can be: * a template class with a default parameter which depends on another template parameter * a nested template class, or type of a nested template class
* BTL: - added tridiagonalization and hessenberg decomposition benchGravatar Gael Guennebaud2008-07-28
| | | | - added GOTO library
* fix some internal asserts in CacheFrinedlyProductGravatar Gael Guennebaud2008-07-27
|
* forgot to include this file in previous commitGravatar Gael Guennebaud2008-07-27
|
* * updated benchmark files according to recent renamingsGravatar Gael Guennebaud2008-07-27
| | | | * various improvements in BTL including trisolver and cholesky bench
* Fix a couple issues introduced in the previous commit:Gravatar Gael Guennebaud2008-07-26
| | | | | * removed DirectAccessBit from Part * use a template specialization in inverseProduct() to transform a Part xpr to a Flagged xpr