56 Commits

Author	SHA1	Message	Date
Steve Bronder	e7b8643d70	Revert "Revert "Adds EIGEN_CONSTEXPR and EIGEN_NOEXCEPT to rows(), cols(), innerStride(), outerStride(), and size()"" ... This reverts commit 5f0b4a4010.	2021-03-24 18:14:56 +00:00
David Tellenbach	5f0b4a4010	Revert "Adds EIGEN_CONSTEXPR and EIGEN_NOEXCEPT to rows(), cols(), innerStride(), outerStride(), and size()" ... This reverts commit 6cbb3038ac because it breaks clang-10 builds on x86 and aarch64 when C++11 is enabled.	2021-03-05 13:16:43 +01:00
Steve Bronder	6cbb3038ac	Adds EIGEN_CONSTEXPR and EIGEN_NOEXCEPT to rows(), cols(), innerStride(), outerStride(), and size()	2021-03-04 18:58:08 +00:00
Patrick Peltzer	15e53d5d93	PR 567: makes all dense solvers inherit SoverBase (LU,Cholesky,QR,SVD). ... This changeset also includes: * add HouseholderSequence::conjugateIf * define int as the StorageIndex type for all dense solvers * dedicated unit tests, including assertion checking * _check_solve_assertion(): this method can be implemented in derived solver classes to implement custom checks * CompleteOrthogonalDecompositions: add applyZOnTheLeftInPlace, fix scalar type in applyZAdjointOnTheLeftInPlace(), add missing assertions * Cholesky: add missing assertions * FullPivHouseholderQR: Corrected Scalar type in _solve_impl() * BDCSVD: Unambiguous return type for ternary operator * SVDBase: Corrected Scalar type in _solve_impl()	2019-01-17 01:17:39 +01:00
Gael Guennebaud	5539587b1f	Some warning fixes	2018-07-17 10:29:12 +02:00
Gael Guennebaud	8a5955a052	Optimize the product of a householder-sequence with the identity, and optimize the evaluation of a HouseholderSequence to a dense matrix using faster blocked product.	2018-07-11 17:16:50 +02:00
Andrea Bocci	f7124b3e46	Extend CUDA support to matrix inversion and selfadjointeigensolver	2018-06-11 18:33:24 +02:00
Christoph Hertzberg	50633d1a83	Renamed .trans() et al. to .reverseFlag() et at. Adapted documentation of .setReverseFlag()	2018-04-17 11:30:27 +02:00
Christoph Hertzberg	42715533f1	bug #1493: Make representation of HouseholderSequence consistent and working for complex numbers. Made corresponding unit test actually test that. Also simplify implementation of QR decompositions	2018-04-15 10:15:28 +02:00
Gael Guennebaud	66e99ab6a1	Relax mixing-type constraints for binary coefficient-wise operators: ... - Replace internal::scalar_product_traits<A,B> by Eigen::ScalarBinaryOpTraits<A,B,OP> - Remove the "functor_is_product_like" helper (was pretty ugly) - Currently, OP is not used, but it is available to the user for fine grained tuning - Currently, only the following operators have been generalized: *,/,+,-,=,*=,/=,+=,-= - TODO: generalize all other binray operators (comparisons,pow,etc.) - TODO: handle "scalar op array" operators (currently only * is handled) - TODO: move the handling of the "void" scalar type to ScalarBinaryOpTraits	2016-06-06 15:11:41 +02:00
Gael Guennebaud	969b8959a0	Fix compilation: Matrix does not indirectly live in the internal namespace anymore!	2016-06-03 07:44:58 +02:00
Gael Guennebaud	30d97c03ce	Defer the allocation of the working space: ... - it is not always needed, - and this fixes a long-to-float conversion warning	2016-05-26 17:39:42 +02:00
Gael Guennebaud	6f960b83ff	Make use of is_same_dense helper instead of extract_data to detect input/outputs are the same.	2016-04-13 18:47:12 +02:00
Gael Guennebaud	b7716c0328	Fix incomplete previous patch on matrix comparision.	2016-04-13 18:32:56 +02:00
Gael Guennebaud	fc202bab39	Index refactoring: StorageIndex must be used for storage only (and locally when it make sense). In all other cases use the global Index type.	2015-02-13 18:57:41 +01:00
Christoph Hertzberg	e8cdbedefb	bug #877, bug #572: Introduce a global Index typedef. Rename Sparse*::Index to StorageIndex, make Dense*::StorageIndex an alias to DenseIndex. Overall this commit gets rid of all Index conversion warnings.	2014-12-04 22:48:53 +01:00
Gael Guennebaud	0ca43f7e9a	Remove deprecated code not used by evaluators	2014-09-18 15:15:27 +02:00
Gael Guennebaud	af9c9f7706	Fix comparison to block size	2014-09-14 17:33:39 +02:00
Gael Guennebaud	b23556bbbd	Oops, a block size of 1 is not very useful, set it to 48 as in HouseholderQR	2014-09-05 08:50:50 +02:00
Gael Guennebaud	8846aa6d1b	Optimization: enable cache-efficient application of HouseholderSequence.	2014-09-04 09:15:59 +02:00
Gael Guennebaud	b868bfb84a	Make operator=(EigenBase<>) uses the new assignment mechanism and introduce a generic EigenBase to EigenBase assignment kind based on the previous evalTo mechanism.	2014-06-25 17:23:52 +02:00
Gael Guennebaud	55365566b2	Fix HouseholderSequence::conjugate() and ::adjoint() and add respective unit tests.	2013-06-17 00:14:42 +02:00
Gael Guennebaud	2b6528effc	HouseholderSequence should expose standard enums (Rows/Cols, etc.))	2013-06-10 11:42:14 +02:00
Benoit Jacob	69124cfca2	Automatic relicensing to MPL2 using Keirs script. Manual fixup follows.	2012-07-13 14:42:47 -04:00
Jitse Niesen	3c412183b2	Get rid of include directives inside namespace blocks (bug #339).	2012-04-15 11:06:28 +01:00
Adolfo Rodriguez Tsourouksdissian	7bf0e8cd82	bug #206 - part 2: For HouseholderSequence objects, added non-allocating versions of evalTo() and applyThisOnTheRight/Left that take additional working vector parameters.	2011-10-30 23:55:16 -04:00
Gael Guennebaud	af712e80e6	fix bug #73: weird compilation error in HouseholderSequence where double and float were mixed. Hopefuly this also solve bug #91...	2011-01-28 12:35:32 +01:00
Jitse Niesen	9111d73017	Fix compilation error in HouseholderSequence introduced in my previous commit.	2011-01-07 13:46:23 +00:00
Jitse Niesen	2cc75f4922	Make HouseholderSequence::setTrans() protected (cf. bug #50). ... Users can call .transpose() instead.	2011-01-06 11:30:19 +00:00
Jitse Niesen	47a9d2ed54	Document HouseholderSequence. ... Incomplete: I did not explain the difference between OnTheLeft and OnTheRight, and there is only one example.	2011-01-02 16:59:44 +00:00
Benoit Jacob	dbd9c5fd50	fix HouseholderSequence API, bug #50: ... * remove ctors taking more than 2 ints * rename actualVectors to length * add length/shift/trans accessors/mutators	2010-12-30 04:18:40 -05:00
Benoit Jacob	75b7d98665	bug #54 - really fix const correctness except in Sparse	2010-12-22 17:45:37 -05:00
Gael Guennebaud	5a3a229550	fix return type of rightHouseholderSequence()	2010-11-15 11:11:22 +01:00
Hauke Heibel	c738cd56eb	Renamed cleantype to remove_all since it is close to remove_{const|pointer|reference}.	2010-10-26 16:47:01 +02:00
Hauke Heibel	7bc8e3ac09	Initial fixes for bug #85. ... Renamed meta_{true|false} to {true|false}_type, meta_if to conditional, is_same_type to is_same, un{ref|pointer|const} to remove_{reference|pointer|const} and makeconst to add_const. Changed boolean type 'ret' member to 'value'. Changed 'ret' members refering to types to 'type'. Adapted all code occurences.	2010-10-25 22:13:49 +02:00
Benoit Jacob	4716040703	bug #86 : use internal:: namespace instead of ei_ prefix	2010-10-25 10:15:22 -04:00
Gael Guennebaud	28e64b0da3	email change	2010-06-24 23:21:58 +02:00
Gael Guennebaud	469382407c	* Make HouseholderSequence::evalTo works in place ... * Clean a bit the Triadiagonalization making sure it the inplace function really works inplace ;), and that only the lower triangular part of the matrix is referenced. * Remove the Tridiagonalization member object of SelfAdjointEigenSolver exploiting the in place capability of HouseholdeSequence. * Update unit test to check SelfAdjointEigenSolver only consider the lower triangular part.	2010-06-10 16:39:46 +02:00
Gael Guennebaud	8350ab9fb8	* remove ei_index, and let ei_traits propagate the index types ... * add an Index type template parapeter to sparse objects	2010-06-03 08:41:11 +02:00
Benoit Jacob	aaaade4b3d	the Index types change. ... As discussed on the list (too long to explain here).	2010-05-30 16:00:58 -04:00
Jitse Niesen	e7dc772554	Respect MaxRowsAtCompileTime in HouseholderSequence::evalTo(). ... This fixes the failing test nomalloc_4. Also remove a print inserted for debugging in schur_real test.	2010-05-25 12:40:42 +01:00
Benoit Jacob	9962c59b56	* implement the corner() API change: new methods topLeftCorner() etc ... * get rid of BlockReturnType: it was not needed, and code was not always using it consistently anyway * add topRows(), leftCols(), bottomRows(), rightCols() * add corners unit-test covering all of that * adapt docs, expand "porting from eigen 2 to 3" * adapt Eigen2Support	2010-04-22 14:11:18 -04:00
Hauke Heibel	f0c8dcf1e2	Renamed AnyMatrixBase to EigenBase.	2010-02-20 15:26:02 +01:00
Benoit Jacob	bfe6fdde24	allow to multiply a householder sequence and a matrix when one is real and one is complex. ... This is especially important as in bidiagonalization, the band matrix is real.	2010-01-15 00:35:26 -05:00
Benoit Jacob	f1d1756cdd	Introduce third template parameter to HouseholderSequence: int Side. ... When it's OnTheRight, we read householder vectors as rows above the diagonal. With unit test. The use case will be bidiagonalization.	2010-01-14 19:16:49 -05:00
Benoit Jacob	24a09ceae8	* Fix a bug in HouseholderQR with mixed fixed/dynamic size: must use EIGEN_SIZE_MIN instead of EIGEN_ENUM_MIN, and there are many other occurences throughout Eigen! ... * HouseholderSequence: - add shift parameter - add essentialVector() method to start abstracting the direction - add unit test in householder.cpp	2010-01-11 08:48:39 -05:00
Benoit Jacob	b05f59ee07	Backed out changeset 58fb27cd56	2010-01-07 22:00:45 -05:00
Benoit Jacob	58fb27cd56	undo	2010-01-07 21:53:52 -05:00
Benoit Jacob	39ac57fa6d	Big renaming: ... start ---> head end ---> tail Much frustration with sed syntax. Need to learn perl some day.	2010-01-04 21:24:43 -05:00
Hauke Heibel	832045d363	Warning 4512 (assignment operators could not be generated) is now simply disabled. ... All unimplemented assignment operators have been removed.	2009-12-14 10:32:43 +01:00