devtools/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2026-04-10 11:34:33 +08:00
Code Issues Packages Projects Releases Wiki Activity

Commit Graph

  • fbbd8afe30 Started a Transform class in the Geometry module to represent homography. Fix indentation in Quaternion.h Gael Guennebaud 2008-06-15 08:33:44 +00:00
  • 4af7089ab8 * Added a generalized eigen solver for the selfadjoint case. (as new members to SelfAdjointEigenSolver) The QR module now depends on Cholesky. * Fix Transpose to correctly preserve the *TriangularBit. Gael Guennebaud 2008-06-14 19:42:12 +00:00
  • f07f907810 Add QR and Cholesky module instantiations in the lib. To try it with the unit tests set the cmake variable TEST_LIB to ON. Gael Guennebaud 2008-06-14 13:02:41 +00:00
  • 53289a8b64 * even though the _Flags default to the corrected value, still correct them in the ei_traits, so that they're guaranteed even if the user specified his own non-default flags (like before). Benoit Jacob 2008-06-13 08:09:48 +00:00
  • c90c77051f * make the _Flags template parameter of Matrix default to the corrected flags. This ensures that unless explicitly messed up otherwise, a Matrix type is equal to its own Eval type. This seriously reduces the number of types instantiated. Measured +13% compile speed, -7% binary size. Benoit Jacob 2008-06-13 07:53:45 +00:00
  • e3fac69f19 Added a Hessenberg decomposition class for both real and complex matrices. This is the first step towards a non-selfadjoint eigen solver. Notes: - We might consider merging Tridiagonalization and Hessenberg toghether ? - Or we could factorize some code into a Householder class (could also be shared with QR) Gael Guennebaud 2008-06-08 15:03:23 +00:00
  • 4dd57b585d * rewrite of the QR decomposition: - works for complex - allows direct access to the matrix R * removed the scale by the matrix dimensions in MatrixBase::isMuchSmallerThan(scalar) Gael Guennebaud 2008-06-07 22:47:11 +00:00
  • eb7b7b2cfc * remove Cross product expression: MatrixBase::cross() now returns a temporary which is even better optimized by the compiler. * Quaternion no longer inherits MatrixBase. Instead it stores the coefficients using a Matrix<> and provides only relevant methods. Gael Guennebaud 2008-06-07 13:18:29 +00:00
  • 6998037930 * move some compile time "if" to their respective unroller (assign and dot) * fix a couple of compilation issues when unrolling is disabled * reduce default unrolling limit to a more reasonable value Gael Guennebaud 2008-06-07 01:07:48 +00:00
  • a172385720 Updated fuzzy comparisons to use L2 norm as all my experiments tends to show L2 norm works very well here. (the legacy implementation is still available via a preprocessor token to allow further experiments if needed...) Gael Guennebaud 2008-06-06 18:37:53 +00:00
  • 8769bfd9aa fix a compilation issue in non debug mode Gael Guennebaud 2008-06-06 14:11:26 +00:00
  • 869394ee8b fix some compile errors with gcc 4.3, some warnings, some documentation Benoit Jacob 2008-06-06 13:10:00 +00:00
  • 2126baf9dc add an optimized path for the tridiagonalization of a 3x3 matrix. (useful for plane fitting, and covariance analysis of 3D data) Gael Guennebaud 2008-06-04 13:41:32 +00:00
  • 48262b9734 added a static assertion mechanism (see notes in Core/util/StaticAssert.h for details) Gael Guennebaud 2008-06-04 11:16:11 +00:00
  • 60726f91a9 hack to to make the nomalloc unit test compiles with -pedantic Gael Guennebaud 2008-06-04 10:15:48 +00:00
  • 42ad9c4352 update of the eigeinsolver unit test to check complex Gael Guennebaud 2008-06-03 18:04:36 +00:00
  • a0cff1a295 fix eigenvectors computations :) Gael Guennebaud 2008-06-03 18:03:55 +00:00
  • 915587d03d * add CommaInitializer::finished to allow the use of (Matrix3() << v0, v1, v2).finished() as an argument of a function. Other possibilities for the name could be "end" or "matrix" ?? * various update in Quaternion, in particular I added a lot of FIXME about the API options, these have to be discussed and fixed. Gael Guennebaud 2008-06-03 15:50:09 +00:00
  • 196f38f5db improved Quaternion class: - Euler angles and angle axis conversions, - stable spherical interpolation - documentation - update the respective unit test Gael Guennebaud 2008-06-03 13:43:29 +00:00
  • bcb32839c2 fix building of examples Gael Guennebaud 2008-06-03 09:43:59 +00:00
  • a9cf229e15 add a geometry unit test and fix a couple of typo in Quaternion.h Gael Guennebaud 2008-06-03 07:32:12 +00:00
  • 8de4d92b70 - get the doc of the enums in MatrixBase right - get the doc of the flags in Constants right - finally give up with SEPARATE_MEMBER_PAGES: it triggers too big Doxygen bugs, and produces too many small pages. So we have one huge page for MatrixBase at currently 300kb and going up, so the solution especially for users with low bandwidth will be to provide an archive of the html documentation. Benoit Jacob 2008-06-03 02:06:18 +00:00
  • 366971bea4 * start of the Geometry module with a cross product and quaternion expressions (haven't tried them yet) * applied the meta selector rule to MatrixBase::swap() Gael Guennebaud 2008-06-02 22:58:36 +00:00
  • 75de41a00b big changes in Doxygen configuration; work around bug with doxygen parsing of initialized enum values showing the last word the initializer instead of the actual enum value's name; add some more docs. Benoit Jacob 2008-06-02 20:08:37 +00:00
  • ac88feebb7 work around Doxygen bug triggered by r814874, which caused many classes to disappear from the docs. Benoit Jacob 2008-06-02 19:29:23 +00:00
  • 6209bbe286 doc improvements: fix linking in Mainpage.dox, improved Doxyfile.in Benoit Jacob 2008-06-02 18:27:37 +00:00
  • f2ebbee274 added unit tests for sizeof and dynamic memory allocation Gael Guennebaud 2008-06-02 14:54:52 +00:00
  • 54ae2ac7e8 Added the computation of eigen vectors in the symmetric eigen solver. However the eigen vectors are not correct yet, but I really cannot find the problem. Gael Guennebaud 2008-06-02 12:52:08 +00:00
  • 3b0523041a since m*m.adjoint() is positive, so are its eigenvalues, so no need for cwiseAbs() Benoit Jacob 2008-06-02 04:45:02 +00:00
  • 0444e3601a - add MatrixBase::eigenvalues() convenience method - add MatrixBase::matrixNorm(); in the non-selfadjoint case, we reduce to the selfadjoint case by using the "C*-identity" a.k.a. norm of x = sqrt(norm of x * x.adjoint()) Benoit Jacob 2008-06-02 04:42:45 +00:00
  • 92b7e2d6a1 fix a couple of issues making the eigensolver test compile and run without aborting on an assert. Had to fix a stupid bug in Block -- very strange we hadn't hit it before. Benoit Jacob 2008-06-02 02:06:33 +00:00
  • 001b01a290 Rewrite from scratch of the eigen solver for symmetric matrices which now supports selfadjoint matrix. The implementation follows Golub's famous book. Gael Guennebaud 2008-06-02 00:30:26 +00:00
  • 06752b2b77 * added a Tridiagonalization class for selfadjoint matrices * added MatrixBase::real() * added the ability to extract a selfadjoint matrix from the lower or upper part of a matrix, e.g.: m.extract<Upper|SelfAdjoint>() will ignore the strict lower part and return a selfadjoint. This is compatible with ZeroDiag and UnitDiag. Gael Guennebaud 2008-06-01 17:20:18 +00:00
  • dc5fd8dfff meagre outcome for so much time spent! * fix inverse() bug discovered by Gael's test * fix warnings introduced by the new Diagonal stuff * update Doxyfile to v1.5.6 Benoit Jacob 2008-06-01 03:36:49 +00:00
  • d5cbb1d002 added a unit test for Inverse discovered that m1 == m1.inverse().inverse() fails though m1.inverse()*m1 == I Gael Guennebaud 2008-05-31 23:44:11 +00:00
  • 64169389ed added an *optional* Eigen2 dynamic library. it allows the possiblity to save some compilation time by linking to it *and* defining the token EIGEN_EXTERN_INSTANCIATIONS Gael Guennebaud 2008-05-31 23:21:49 +00:00
  • fcf4457b78 added optimized matrix times diagonal matrix product via Diagonal flag shortcut. Gael Guennebaud 2008-05-31 21:35:11 +00:00
  • 310f7aa096 moved purely "array" related stuff to a new module Array. This include: - cwise Pow,Sin,Cos,Exp... - cwise Greater and other comparison operators - .any(), .all() and partial reduction - random Gael Guennebaud 2008-05-31 18:11:48 +00:00
  • a2f71f9d7e updated EigenSolver to use .coeff / .coeffRef Gael Guennebaud 2008-05-31 16:31:10 +00:00
  • c9fb248c36 simply a bit the basic product moving dynamic loops to the corresponding special case of the unrollers. the latter ones are therefore re-named *product_impl. Gael Guennebaud 2008-05-31 15:06:26 +00:00
  • f5e599e489 * replace compile-time-if by meta-selector in Assign.h as it speed up compilation. * fix minor typo introduced in the previous commit Gael Guennebaud 2008-05-31 14:42:07 +00:00
  • e2ac5d244e Added ArrayBit to get the ability to manipulate a Matrix like a simple scalar. In particular this flag changes the behavior of operator* to a coeff wise product. Gael Guennebaud 2008-05-29 22:33:07 +00:00
  • b501e08d81 now the unit-tests (hence all of Eigen) don't depend on Qt at all anymore. Benoit Jacob 2008-05-29 03:37:16 +00:00
  • 486fdb26a1 many small fixes and documentation improvements, this should be alpha5. Benoit Jacob 2008-05-29 03:12:30 +00:00
  • c1559d3079 * updated the assignement operator macro so that overloads in MatrixBase work * removed product_selector and cleaned Product.h a bit * cleaned Assign.h a bit Gael Guennebaud 2008-05-28 22:56:19 +00:00
  • 8711e26c8a * change Flagged to take into account NestByValue only * bugfix in Assign and cache friendly product (weird that worked before) * improved argument evaluation in Product Gael Guennebaud 2008-05-28 22:11:47 +00:00
  • 73084dc754 * added _*coeffRef members in NestedByValue * added ConjugateReturnType and AdjointReturnType that are type-defined to Derived& and Transpose<Derived> if the scalar type is not complex: this avoids abusive copies in the cache friendly Product Gael Guennebaud 2008-05-28 09:09:18 +00:00
  • f54760c889 hehe, the complicated nesting scheme in Flagged in the previous commit was a sign that we were doing something wrong. In fact, having NestByValue as a special case of Flagged was wrong, and the previous commit, while not buggy, was inefficient because then when the resulting NestByValue xpr was nested -- hence copied -- the original xpr which was already nested by value was copied again; hence instead of 1 copy we got 3 copies. The solution was to ressuscitate the old Temporary.h (renamed NestByValue.h) as it was the right approach. Benoit Jacob 2008-05-28 05:14:16 +00:00
  • aebecae510 * find the proper way of nesting the expression in Flagged: finally that's more subtle than just using ei_nested, because when flagging with NestByValueBit we want to store the expression by value already, regardless of whether it already had the NestByValueBit set. * rename temporary() ----> nestByValue() * move the old Product.h to disabled/, replace by what was ProductWIP.h * tweak -O and -g flags for tests and examples * reorder the tests -- basic things go first * simplifications, e.g. in many methoeds return derived() and count on implicit casting to the actual return type. * strip some not-really-useful stuff from the heaviest tests Benoit Jacob 2008-05-28 04:38:16 +00:00
  • 559233c73e * fix the QR module to use extract/part instead of the previous triangular stuff * added qr and eigensolver tests * fix a compilation warning in Matrix copy constructor Gael Guennebaud 2008-05-27 09:16:27 +00:00
  • 5aa00f6870 part 2 of big change: rename Triangular.h -> Extract.h (svn required to commit that separately) Benoit Jacob 2008-05-27 05:50:36 +00:00
  • 953efdbfe7 - introduce Part and Extract classes, splitting and extending the former Triangular class - full meta-unrolling in Part - move inverseProduct() to MatrixBase - compilation fix in ProductWIP: introduce a meta-selector to only do direct access on types that support it. - phase out the old Product, remove the WIP_DIRTY stuff. - misc renaming and fixes Benoit Jacob 2008-05-27 05:47:30 +00:00
  • 8f1fc80a77 some documentation fixes (Cwise* and Cholesky) Gael Guennebaud 2008-05-22 16:31:00 +00:00
  • 94e1629a1b * improved product performance: - fallback to normal product for small dynamic matrices - overloaded "c += (a * b).lazy()" to avoid the expensive and useless temporary and setZero() in such very common cases. * fix a couple of issues with the flags Gael Guennebaud 2008-05-22 14:51:25 +00:00
  • 106a0c1bef restored the product test Gael Guennebaud 2008-05-22 12:35:09 +00:00
  • 9ab6e186eb remove Like1DArrayBit in Transpose Gael Guennebaud 2008-05-22 12:25:11 +00:00
  • 522e24f2d7 update of the testing framework: replaced the QTestLib framework my custom macros and a (optional) custom script to run the tests from ctest. Gael Guennebaud 2008-05-22 12:18:55 +00:00
  • c6789a279c Fix compilation issues with MSVC and NVCC. Added a few typedef of complex return types in MatrixBase (Needed by MSVC) Gael Guennebaud 2008-05-15 09:40:11 +00:00
  • 5da60897ab Introduce generic Flagged xpr, remove already Lazy.h and Temporary.h Rename DefaultLostFlagMask --> HerediraryBits Benoit Jacob 2008-05-14 08:20:15 +00:00
  • fd2e9e5c3c * Clean a bit the eigenvalue solver: if the matrix is known to be selfadjoint at compile time, then it returns real eigenvalues. * Fix a couple of bugs with the new product. Gael Guennebaud 2008-05-13 07:40:25 +00:00
  • 3eccfd1a78 -fix certain #includes -fix CMakeLists, public headers weren't getting installed Benoit Jacob 2008-05-12 21:15:17 +00:00
  • 4317fad869 * Added several cast to int of the enums (needed for some compilers) * Fix a mistake in CwiseNullary. * Added a CoreDeclarions header that declares only the forward declarations and related basic stuffs. Gael Guennebaud 2008-05-12 18:09:30 +00:00
  • 678f18fce4 put inline keywords everywhere appropriate. So we don't need anymore to pass -finline-limit=1000 to gcc to get good performance. By the way some cleanup. Benoit Jacob 2008-05-12 17:34:46 +00:00
  • f0eb3d2d3b updated product test to carefully test all scalar types and fix an issue in the triangular test Gael Guennebaud 2008-05-12 10:26:10 +00:00
  • 45cda6704a * Draft of a eigenvalues solver (does not support complex and does not re-use the QR decomposition) Gael Guennebaud 2008-05-12 10:23:09 +00:00
  • dca416cace move arch-specific code to arch/SSE and arch/AltiVec subdirs. rename the noarch PacketMath.h to DummyPacketMath.h Benoit Jacob 2008-05-12 08:30:42 +00:00
  • 3562b01105 * Give Konstantinos a copyright line * Fix compilation of Inverse.h with vectorisation * Introduce EIGEN_GNUC_AT_LEAST(x,y) macro doing future-proof (e.g. gcc v5.0) check * Only use ProductWIP if vectorisation is enabled * rename EIGEN_ALWAYS_INLINE -> EIGEN_INLINE with fall-back to inline keyword * some cleanup/indentation Benoit Jacob 2008-05-12 08:12:40 +00:00
  • 4f6d7abc87 only include SSE3 headers if compiling with SSE3 support Benoit Jacob 2008-05-08 09:15:16 +00:00
  • 4754fa4868 removed "sort brief" in doxygen documentation Gael Guennebaud 2008-05-08 08:13:38 +00:00
  • bf5326c3ca * Added ReferencableBit flag to known if coeffRef is available. (needed by the new product implementation) * Make the packet* members template to support aligned and unaligned access. This makes Block vectorizable. Combined with ReferencableBit, we should be able to determine at runtime (in some specific cases) if an aligned vectorization is possible or not. * Improved the new product implementation to robustly handle all cases, it now passes all the tests. * Renamed the packet version ei_predux to ei_preduxp to avoid name collision. Gael Guennebaud 2008-05-08 08:12:52 +00:00
  • 64c49de7ba * split PacketMath.h to SSE and Altivec specific files * improved the flexibility of the new product implementation, now all sizes seems to be properly handled. Gael Guennebaud 2008-05-05 17:19:47 +00:00
  • 46fa4c713f * Started support for unaligned vectorization. * Introduce a new highly optimized matrix-matrix product for large matrices. The code is still highly experimental and it is activated only if you define EIGEN_WIP_PRODUCT at compile time. Currently the third dimension of the product must be a factor of the packet size (x4 for floats) and the right handed side matrix must be column major. Moreover, currently c = a*b; actually computes c += a*b !! Therefore, the code is provided for experimentation purpose only ! These limitations will be fixed soon or later to become the default product implementation. Gael Guennebaud 2008-05-05 10:23:29 +00:00
  • 8c6007f80e * Patch by Konstantinos Margaritis: AltiVec vectorization. * Fix several warnings, temporarily disable determinant test. Benoit Jacob 2008-05-03 12:21:23 +00:00
  • 0545df2149 slighly improved the cache friendly product to use mul-add only Gael Guennebaud 2008-05-03 10:01:30 +00:00
  • a6655dd91a added packet mul-add function (ei_pmad) and updated Product to use it. this change nothing for current SSE architecture but might be helpful for altivec/cell and up comming AMD processors. Gael Guennebaud 2008-05-03 00:45:08 +00:00
  • 102e029dad Removed ei_pload1, use posix_memalign to allocate aligned memory, and make Product ok when only one side is vectorizable (and the product is still vectorized) Gael Guennebaud 2008-05-02 13:30:12 +00:00
  • e19f9bc523 added a test for triangular matrices Gael Guennebaud 2008-05-02 11:35:59 +00:00
  • 890a8de962 Make products always eval into expressions. Improves performance in benchmark. Still not as fasts as explicit eval(), strangely. Benoit Jacob 2008-05-02 08:53:23 +00:00
  • ef5b20bc50 fix flag and cost computations for nested expressions Gael Guennebaud 2008-05-01 18:58:30 +00:00
  • 5588def0cf nullary xpr are now vectorized Gael Guennebaud 2008-05-01 14:28:53 +00:00
  • 02f1615d2a Enable vectorization of product with dynamic matrices, extended cache optimal product to work in any row/column major situations, and a few bugfixes (forgot to add the Cholesky header, vectorization of CwiseBinary) Gael Guennebaud 2008-05-01 13:53:05 +00:00
  • 6486991ac3 some cleaning in Cholesky and removed evil ei_sqrt of complex Gael Guennebaud 2008-04-27 18:57:28 +00:00
  • 64bacf1c3f * added ei_sqrt for complex * updated Cholesky to support complex * correct result_type for abs and abs2 functors Gael Guennebaud 2008-04-27 14:05:40 +00:00
  • 4ffffa670e added Cholesky module Gael Guennebaud 2008-04-27 10:57:32 +00:00
  • 1ec2d21ca5 Fixed a couple of issues introduced in previous commits. Added a test for Triangular. Gael Guennebaud 2008-04-26 20:28:27 +00:00
  • b4c974d059 Added triangular assignement, e.g.: m.upper() = a+b; only updates the upper triangular part of m. Note that: m = (a+b).upper(); updates all coefficients of m (but half of the additions will be skiped) Gael Guennebaud 2008-04-26 19:20:26 +00:00
  • 4c92150676 Added Triangular expression to extract upper or lower (strictly or not) part of a matrix. Triangular also provide an optimised method for forward and backward substitution. Further optimizations regarding assignments and products might come later. Gael Guennebaud 2008-04-26 18:26:05 +00:00
  • 62bf0bbd59 fix a bug in determinant of 4x4 matrices and a small type issue in Inverse Gael Guennebaud 2008-04-26 08:56:52 +00:00
  • 173e582e3c added a tough test to check the determinant that currently fails Gael Guennebaud 2008-04-25 23:13:20 +00:00
  • 6f2c72fb53 Various fixes in: - vector to vector assign - PartialRedux - Vectorization criteria of Product - returned type of normalized - SSE integer mul Gael Guennebaud 2008-04-25 23:10:37 +00:00
  • a451835bce Make the explicit vectorization much more flexible: - support dynamic sizes - support arbitrary matrix size when the matrix can be seen as a 1D array (except for fixed size matrices where the size in Bytes must be a factor of 16, this is to allow compact storage of a vector of matrices) Note that the explict vectorization is still experimental and far to be completely tested. Gael Guennebaud 2008-04-25 15:46:18 +00:00
  • 30d47b5250 forgot to add a file in the previous commit Gael Guennebaud 2008-04-24 20:25:55 +00:00
  • 9385793f71 Fix a couple of issue with the vectorization. In particular, default ei_p* functions are provided to handle not suported types seemlessly. Gael Guennebaud 2008-04-24 18:35:39 +00:00
  • 6ae037dfb5 give up on OpenMP... for now Benoit Jacob 2008-04-18 07:57:46 +00:00
  • acfd6f3bda - add _packetCoeff() to Inverse, allowing vectorization. - let Inverse take template parameter MatrixType instead of ExpressionType, in order to reduce executable code size when taking inverses of xpr's. - introduce ei_corrected_matrix_flags : the flags template parameter to the Matrix class is only a suggestion. This is also useful in ei_eval. Benoit Jacob 2008-04-16 07:18:27 +00:00
  • 43e2bc14fe +5% optimization in 4x4 inverse: -only evaluate block expressions for which that is beneficial -don't check for invertibility unless requested Benoit Jacob 2008-04-15 20:39:27 +00:00
  • 6747b45ae7 for 4x4 matrices implement the special algorithm that Markos proposed, falling back to the general algorithm in the bad case. Benoit Jacob 2008-04-15 20:15:36 +00:00
  • 2a86f052a5 - optimized determinant calculations for small matrices (size <= 4) (only 30 muls for size 4) - rework the matrix inversion: now using cofactor technique for size<=3, so the ugly unrolling is only used for size 4 anymore, and even there I'm looking to get rid of it. Benoit Jacob 2008-04-14 17:07:12 +00:00
  • 9789c04467 when evaluating an xpr, the result can now be vectorizable even if the xpr itself wasn't vectorizable. Benoit Jacob 2008-04-14 08:55:12 +00:00
  • ea3ccb1e8c * Start of the LU module, with matrix inversion already there and fully optimized. * Even if LargeBit is set, only parallelize for large enough objects (controlled by EIGEN_PARALLELIZATION_TRESHOLD). Benoit Jacob 2008-04-14 08:20:24 +00:00