1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This file is part of Eigen, a lightweight C++ template library 2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// for linear algebra. 3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Copyright (C) 2009 Gael Guennebaud <gael.guennebaud@inria.fr> 5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Copyright (C) 2010 Benoit Jacob <jacob.benoit.1@gmail.com> 6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This Source Code Form is subject to the terms of the Mozilla 8c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Public License v. 2.0. If a copy of the MPL was not distributed 9c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 10c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 11c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#ifndef EIGEN_HOUSEHOLDER_SEQUENCE_H 12c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#define EIGEN_HOUSEHOLDER_SEQUENCE_H 13c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 14c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathnamespace Eigen { 15c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 16c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath/** \ingroup Householder_Module 17c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \householder_module 18c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \class HouseholderSequence 19c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \brief Sequence of Householder reflections acting on subspaces with decreasing size 20c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \tparam VectorsType type of matrix containing the Householder vectors 21c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \tparam CoeffsType type of vector containing the Householder coefficients 22c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \tparam Side either OnTheLeft (the default) or OnTheRight 23c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 24c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * This class represents a product sequence of Householder reflections where the first Householder reflection 25c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * acts on the whole space, the second Householder reflection leaves the one-dimensional subspace spanned by 26c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * the first unit vector invariant, the third Householder reflection leaves the two-dimensional subspace 27c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * spanned by the first two unit vectors invariant, and so on up to the last reflection which leaves all but 28c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * one dimensions invariant and acts only on the last dimension. Such sequences of Householder reflections 29c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * are used in several algorithms to zero out certain parts of a matrix. Indeed, the methods 30c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * HessenbergDecomposition::matrixQ(), Tridiagonalization::matrixQ(), HouseholderQR::householderQ(), 31c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * and ColPivHouseholderQR::householderQ() all return a %HouseholderSequence. 32c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 33c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * More precisely, the class %HouseholderSequence represents an \f$ n \times n \f$ matrix \f$ H \f$ of the 34c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * form \f$ H = \prod_{i=0}^{n-1} H_i \f$ where the i-th Householder reflection is \f$ H_i = I - h_i v_i 35c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * v_i^* \f$. The i-th Householder coefficient \f$ h_i \f$ is a scalar and the i-th Householder vector \f$ 36c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * v_i \f$ is a vector of the form 37c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \f[ 38c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * v_i = [\underbrace{0, \ldots, 0}_{i-1\mbox{ zeros}}, 1, \underbrace{*, \ldots,*}_{n-i\mbox{ arbitrary entries}} ]. 39c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \f] 40c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * The last \f$ n-i \f$ entries of \f$ v_i \f$ are called the essential part of the Householder vector. 41c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 42c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * Typical usages are listed below, where H is a HouseholderSequence: 43c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \code 44c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * A.applyOnTheRight(H); // A = A * H 45c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * A.applyOnTheLeft(H); // A = H * A 46c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * A.applyOnTheRight(H.adjoint()); // A = A * H^* 47c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * A.applyOnTheLeft(H.adjoint()); // A = H^* * A 48c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * MatrixXd Q = H; // conversion to a dense matrix 49c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \endcode 50c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * In addition to the adjoint, you can also apply the inverse (=adjoint), the transpose, and the conjugate operators. 51c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 52c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * See the documentation for HouseholderSequence(const VectorsType&, const CoeffsType&) for an example. 53c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 54c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \sa MatrixBase::applyOnTheLeft(), MatrixBase::applyOnTheRight() 55c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 56c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 57c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathnamespace internal { 58c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 59c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename VectorsType, typename CoeffsType, int Side> 60c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathstruct traits<HouseholderSequence<VectorsType,CoeffsType,Side> > 61c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 62c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef typename VectorsType::Scalar Scalar; 632b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang typedef typename VectorsType::StorageIndex StorageIndex; 64c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef typename VectorsType::StorageKind StorageKind; 65c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath enum { 66c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath RowsAtCompileTime = Side==OnTheLeft ? traits<VectorsType>::RowsAtCompileTime 67c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath : traits<VectorsType>::ColsAtCompileTime, 68c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ColsAtCompileTime = RowsAtCompileTime, 69c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MaxRowsAtCompileTime = Side==OnTheLeft ? traits<VectorsType>::MaxRowsAtCompileTime 70c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath : traits<VectorsType>::MaxColsAtCompileTime, 71c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MaxColsAtCompileTime = MaxRowsAtCompileTime, 72c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Flags = 0 73c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath }; 74c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath}; 75c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 762b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wangstruct HouseholderSequenceShape {}; 772b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang 782b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wangtemplate<typename VectorsType, typename CoeffsType, int Side> 792b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wangstruct evaluator_traits<HouseholderSequence<VectorsType,CoeffsType,Side> > 802b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang : public evaluator_traits_base<HouseholderSequence<VectorsType,CoeffsType,Side> > 812b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang{ 822b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang typedef HouseholderSequenceShape Shape; 832b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang}; 842b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang 85c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename VectorsType, typename CoeffsType, int Side> 86c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathstruct hseq_side_dependent_impl 87c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 88c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef Block<const VectorsType, Dynamic, 1> EssentialVectorType; 89c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef HouseholderSequence<VectorsType, CoeffsType, OnTheLeft> HouseholderSequenceType; 90c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline const EssentialVectorType essentialVector(const HouseholderSequenceType& h, Index k) 91c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 92c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index start = k+1+h.m_shift; 93c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return Block<const VectorsType,Dynamic,1>(h.m_vectors, start, k, h.rows()-start, 1); 94c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 95c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath}; 96c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 97c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename VectorsType, typename CoeffsType> 98c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathstruct hseq_side_dependent_impl<VectorsType, CoeffsType, OnTheRight> 99c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 100c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef Transpose<Block<const VectorsType, 1, Dynamic> > EssentialVectorType; 101c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef HouseholderSequence<VectorsType, CoeffsType, OnTheRight> HouseholderSequenceType; 102c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline const EssentialVectorType essentialVector(const HouseholderSequenceType& h, Index k) 103c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 104c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index start = k+1+h.m_shift; 105c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return Block<const VectorsType,1,Dynamic>(h.m_vectors, k, start, 1, h.rows()-start).transpose(); 106c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 107c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath}; 108c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 109c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename OtherScalarType, typename MatrixType> struct matrix_type_times_scalar_type 110c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 1112b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang typedef typename ScalarBinaryOpTraits<OtherScalarType, typename MatrixType::Scalar>::ReturnType 112c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ResultScalar; 113c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef Matrix<ResultScalar, MatrixType::RowsAtCompileTime, MatrixType::ColsAtCompileTime, 114c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 0, MatrixType::MaxRowsAtCompileTime, MatrixType::MaxColsAtCompileTime> Type; 115c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath}; 116c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 117c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} // end namespace internal 118c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 119c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename VectorsType, typename CoeffsType, int Side> class HouseholderSequence 120c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath : public EigenBase<HouseholderSequence<VectorsType,CoeffsType,Side> > 121c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 1227faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez typedef typename internal::hseq_side_dependent_impl<VectorsType,CoeffsType,Side>::EssentialVectorType EssentialVectorType; 1237faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez 1247faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez public: 125c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath enum { 126c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath RowsAtCompileTime = internal::traits<HouseholderSequence>::RowsAtCompileTime, 127c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ColsAtCompileTime = internal::traits<HouseholderSequence>::ColsAtCompileTime, 128c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MaxRowsAtCompileTime = internal::traits<HouseholderSequence>::MaxRowsAtCompileTime, 129c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MaxColsAtCompileTime = internal::traits<HouseholderSequence>::MaxColsAtCompileTime 130c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath }; 131c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef typename internal::traits<HouseholderSequence>::Scalar Scalar; 132c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 133c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef HouseholderSequence< 1347faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez typename internal::conditional<NumTraits<Scalar>::IsComplex, 1357faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez typename internal::remove_all<typename VectorsType::ConjugateReturnType>::type, 1367faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez VectorsType>::type, 137c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename internal::conditional<NumTraits<Scalar>::IsComplex, 138c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename internal::remove_all<typename CoeffsType::ConjugateReturnType>::type, 139c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CoeffsType>::type, 140c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Side 141c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath > ConjugateReturnType; 142c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 143c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \brief Constructor. 144c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \param[in] v %Matrix containing the essential parts of the Householder vectors 145c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \param[in] h Vector containing the Householder coefficients 146c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 147c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * Constructs the Householder sequence with coefficients given by \p h and vectors given by \p v. The 148c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * i-th Householder coefficient \f$ h_i \f$ is given by \p h(i) and the essential part of the i-th 149c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * Householder vector \f$ v_i \f$ is given by \p v(k,i) with \p k > \p i (the subdiagonal part of the 150c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * i-th column). If \p v has fewer columns than rows, then the Householder sequence contains as many 151c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * Householder reflections as there are columns. 152c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 153c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \note The %HouseholderSequence object stores \p v and \p h by reference. 154c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 155c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * Example: \include HouseholderSequence_HouseholderSequence.cpp 156c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * Output: \verbinclude HouseholderSequence_HouseholderSequence.out 157c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 158c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \sa setLength(), setShift() 159c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 160c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath HouseholderSequence(const VectorsType& v, const CoeffsType& h) 161c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath : m_vectors(v), m_coeffs(h), m_trans(false), m_length(v.diagonalSize()), 162c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_shift(0) 163c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 164c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 165c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 166c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \brief Copy constructor. */ 167c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath HouseholderSequence(const HouseholderSequence& other) 168c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath : m_vectors(other.m_vectors), 169c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_coeffs(other.m_coeffs), 170c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_trans(other.m_trans), 171c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_length(other.m_length), 172c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_shift(other.m_shift) 173c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 174c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 175c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 176c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \brief Number of rows of transformation viewed as a matrix. 177c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \returns Number of rows 178c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \details This equals the dimension of the space that the transformation acts on. 179c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 180c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index rows() const { return Side==OnTheLeft ? m_vectors.rows() : m_vectors.cols(); } 181c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 182c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \brief Number of columns of transformation viewed as a matrix. 183c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \returns Number of columns 184c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \details This equals the dimension of the space that the transformation acts on. 185c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 186c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index cols() const { return rows(); } 187c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 188c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \brief Essential part of a Householder vector. 189c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \param[in] k Index of Householder reflection 190c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \returns Vector containing non-trivial entries of k-th Householder vector 191c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 192c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * This function returns the essential part of the Householder vector \f$ v_i \f$. This is a vector of 193c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * length \f$ n-i \f$ containing the last \f$ n-i \f$ entries of the vector 194c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \f[ 195c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * v_i = [\underbrace{0, \ldots, 0}_{i-1\mbox{ zeros}}, 1, \underbrace{*, \ldots,*}_{n-i\mbox{ arbitrary entries}} ]. 196c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \f] 197c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * The index \f$ i \f$ equals \p k + shift(), corresponding to the k-th column of the matrix \p v 198c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * passed to the constructor. 199c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 200c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \sa setShift(), shift() 201c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 202c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath const EssentialVectorType essentialVector(Index k) const 203c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 204c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath eigen_assert(k >= 0 && k < m_length); 205c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return internal::hseq_side_dependent_impl<VectorsType,CoeffsType,Side>::essentialVector(*this, k); 206c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 207c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 208c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \brief %Transpose of the Householder sequence. */ 209c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath HouseholderSequence transpose() const 210c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 211c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return HouseholderSequence(*this).setTrans(!m_trans); 212c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 213c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 214c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \brief Complex conjugate of the Householder sequence. */ 215c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ConjugateReturnType conjugate() const 216c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 2177faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez return ConjugateReturnType(m_vectors.conjugate(), m_coeffs.conjugate()) 218c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath .setTrans(m_trans) 219c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath .setLength(m_length) 220c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath .setShift(m_shift); 221c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 222c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 223c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \brief Adjoint (conjugate transpose) of the Householder sequence. */ 224c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ConjugateReturnType adjoint() const 225c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 226c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return conjugate().setTrans(!m_trans); 227c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 228c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 229c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \brief Inverse of the Householder sequence (equals the adjoint). */ 230c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ConjugateReturnType inverse() const { return adjoint(); } 231c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 232c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \internal */ 233c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath template<typename DestType> inline void evalTo(DestType& dst) const 234c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 235c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Matrix<Scalar, DestType::RowsAtCompileTime, 1, 236c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath AutoAlign|ColMajor, DestType::MaxRowsAtCompileTime, 1> workspace(rows()); 237c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath evalTo(dst, workspace); 238c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 239c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 240c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \internal */ 241c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath template<typename Dest, typename Workspace> 242c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath void evalTo(Dest& dst, Workspace& workspace) const 243c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 244c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath workspace.resize(rows()); 245c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index vecs = m_length; 2462b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang if(internal::is_same_dense(dst,m_vectors)) 247c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 248c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // in-place 249c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath dst.diagonal().setOnes(); 250c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath dst.template triangularView<StrictlyUpper>().setZero(); 251c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(Index k = vecs-1; k >= 0; --k) 252c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 253c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index cornerSize = rows() - k - m_shift; 254c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if(m_trans) 255c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath dst.bottomRightCorner(cornerSize, cornerSize) 256c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath .applyHouseholderOnTheRight(essentialVector(k), m_coeffs.coeff(k), workspace.data()); 257c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath else 258c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath dst.bottomRightCorner(cornerSize, cornerSize) 259c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath .applyHouseholderOnTheLeft(essentialVector(k), m_coeffs.coeff(k), workspace.data()); 260c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 261c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // clear the off diagonal vector 262c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath dst.col(k).tail(rows()-k-1).setZero(); 263c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 264c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // clear the remaining columns if needed 265c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(Index k = 0; k<cols()-vecs ; ++k) 266c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath dst.col(k).tail(rows()-k-1).setZero(); 267c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 268c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath else 269c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 270c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath dst.setIdentity(rows(), rows()); 271c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(Index k = vecs-1; k >= 0; --k) 272c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 273c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index cornerSize = rows() - k - m_shift; 274c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if(m_trans) 275c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath dst.bottomRightCorner(cornerSize, cornerSize) 276c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath .applyHouseholderOnTheRight(essentialVector(k), m_coeffs.coeff(k), &workspace.coeffRef(0)); 277c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath else 278c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath dst.bottomRightCorner(cornerSize, cornerSize) 279c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath .applyHouseholderOnTheLeft(essentialVector(k), m_coeffs.coeff(k), &workspace.coeffRef(0)); 280c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 281c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 282c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 283c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 284c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \internal */ 285c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath template<typename Dest> inline void applyThisOnTheRight(Dest& dst) const 286c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 287c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Matrix<Scalar,1,Dest::RowsAtCompileTime,RowMajor,1,Dest::MaxRowsAtCompileTime> workspace(dst.rows()); 288c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath applyThisOnTheRight(dst, workspace); 289c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 290c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 291c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \internal */ 292c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath template<typename Dest, typename Workspace> 293c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline void applyThisOnTheRight(Dest& dst, Workspace& workspace) const 294c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 295c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath workspace.resize(dst.rows()); 296c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(Index k = 0; k < m_length; ++k) 297c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 298c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index actual_k = m_trans ? m_length-k-1 : k; 299c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath dst.rightCols(rows()-m_shift-actual_k) 300c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath .applyHouseholderOnTheRight(essentialVector(actual_k), m_coeffs.coeff(actual_k), workspace.data()); 301c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 302c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 303c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 304c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \internal */ 305c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath template<typename Dest> inline void applyThisOnTheLeft(Dest& dst) const 306c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 3072b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang Matrix<Scalar,1,Dest::ColsAtCompileTime,RowMajor,1,Dest::MaxColsAtCompileTime> workspace; 308c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath applyThisOnTheLeft(dst, workspace); 309c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 310c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 311c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \internal */ 312c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath template<typename Dest, typename Workspace> 313c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline void applyThisOnTheLeft(Dest& dst, Workspace& workspace) const 314c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 3152b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang const Index BlockSize = 48; 3162b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang // if the entries are large enough, then apply the reflectors by block 3172b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang if(m_length>=BlockSize && dst.cols()>1) 318c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 3192b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang for(Index i = 0; i < m_length; i+=BlockSize) 3202b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang { 3212b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang Index end = m_trans ? (std::min)(m_length,i+BlockSize) : m_length-i; 3222b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang Index k = m_trans ? i : (std::max)(Index(0),end-BlockSize); 3232b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang Index bs = end-k; 3242b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang Index start = k + m_shift; 3252b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang 3262b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang typedef Block<typename internal::remove_all<VectorsType>::type,Dynamic,Dynamic> SubVectorsType; 3272b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang SubVectorsType sub_vecs1(m_vectors.const_cast_derived(), Side==OnTheRight ? k : start, 3282b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang Side==OnTheRight ? start : k, 3292b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang Side==OnTheRight ? bs : m_vectors.rows()-start, 3302b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang Side==OnTheRight ? m_vectors.cols()-start : bs); 3312b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang typename internal::conditional<Side==OnTheRight, Transpose<SubVectorsType>, SubVectorsType&>::type sub_vecs(sub_vecs1); 3322b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang Block<Dest,Dynamic,Dynamic> sub_dst(dst,dst.rows()-rows()+m_shift+k,0, rows()-m_shift-k,dst.cols()); 3332b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang apply_block_householder_on_the_left(sub_dst, sub_vecs, m_coeffs.segment(k, bs), !m_trans); 3342b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang } 3352b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang } 3362b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang else 3372b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang { 3382b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang workspace.resize(dst.cols()); 3392b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang for(Index k = 0; k < m_length; ++k) 3402b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang { 3412b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang Index actual_k = m_trans ? k : m_length-k-1; 3422b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang dst.bottomRows(rows()-m_shift-actual_k) 3432b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang .applyHouseholderOnTheLeft(essentialVector(actual_k), m_coeffs.coeff(actual_k), workspace.data()); 3442b8756b6f1de65d3f8bffab45be6c44ceb7411fcMiao Wang } 345c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 346c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 347c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 348c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \brief Computes the product of a Householder sequence with a matrix. 349c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \param[in] other %Matrix being multiplied. 350c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \returns Expression object representing the product. 351c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 352c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * This function computes \f$ HM \f$ where \f$ H \f$ is the Householder sequence represented by \p *this 353c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * and \f$ M \f$ is the matrix \p other. 354c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 355c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath template<typename OtherDerived> 356c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename internal::matrix_type_times_scalar_type<Scalar, OtherDerived>::Type operator*(const MatrixBase<OtherDerived>& other) const 357c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 358c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename internal::matrix_type_times_scalar_type<Scalar, OtherDerived>::Type 359c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath res(other.template cast<typename internal::matrix_type_times_scalar_type<Scalar,OtherDerived>::ResultScalar>()); 360c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath applyThisOnTheLeft(res); 361c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return res; 362c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 363c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 364c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath template<typename _VectorsType, typename _CoeffsType, int _Side> friend struct internal::hseq_side_dependent_impl; 365c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 366c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \brief Sets the length of the Householder sequence. 367c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \param [in] length New value for the length. 368c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 369c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * By default, the length \f$ n \f$ of the Householder sequence \f$ H = H_0 H_1 \ldots H_{n-1} \f$ is set 370c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * to the number of columns of the matrix \p v passed to the constructor, or the number of rows if that 371c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * is smaller. After this function is called, the length equals \p length. 372c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 373c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \sa length() 374c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 375c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath HouseholderSequence& setLength(Index length) 376c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 377c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_length = length; 378c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return *this; 379c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 380c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 381c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \brief Sets the shift of the Householder sequence. 382c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \param [in] shift New value for the shift. 383c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 384c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * By default, a %HouseholderSequence object represents \f$ H = H_0 H_1 \ldots H_{n-1} \f$ and the i-th 385c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * column of the matrix \p v passed to the constructor corresponds to the i-th Householder 386c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * reflection. After this function is called, the object represents \f$ H = H_{\mathrm{shift}} 387c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * H_{\mathrm{shift}+1} \ldots H_{n-1} \f$ and the i-th column of \p v corresponds to the (shift+i)-th 388c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * Householder reflection. 389c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 390c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \sa shift() 391c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 392c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath HouseholderSequence& setShift(Index shift) 393c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 394c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_shift = shift; 395c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return *this; 396c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 397c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 398c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index length() const { return m_length; } /**< \brief Returns the length of the Householder sequence. */ 399c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index shift() const { return m_shift; } /**< \brief Returns the shift of the Householder sequence. */ 400c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 401c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /* Necessary for .adjoint() and .conjugate() */ 402c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath template <typename VectorsType2, typename CoeffsType2, int Side2> friend class HouseholderSequence; 403c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 404c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath protected: 405c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 406c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** \brief Sets the transpose flag. 407c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \param [in] trans New value of the transpose flag. 408c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 409c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * By default, the transpose flag is not set. If the transpose flag is set, then this object represents 410c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \f$ H^T = H_{n-1}^T \ldots H_1^T H_0^T \f$ instead of \f$ H = H_0 H_1 \ldots H_{n-1} \f$. 411c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 412c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \sa trans() 413c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 414c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath HouseholderSequence& setTrans(bool trans) 415c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 416c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_trans = trans; 417c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return *this; 418c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 419c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 420c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath bool trans() const { return m_trans; } /**< \brief Returns the transpose flag. */ 421c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 422c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename VectorsType::Nested m_vectors; 423c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename CoeffsType::Nested m_coeffs; 424c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath bool m_trans; 425c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index m_length; 426c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Index m_shift; 427c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath}; 428c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 429c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath/** \brief Computes the product of a matrix with a Householder sequence. 430c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \param[in] other %Matrix being multiplied. 431c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \param[in] h %HouseholderSequence being multiplied. 432c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \returns Expression object representing the product. 433c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 434c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * This function computes \f$ MH \f$ where \f$ M \f$ is the matrix \p other and \f$ H \f$ is the 435c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * Householder sequence represented by \p h. 436c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 437c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename OtherDerived, typename VectorsType, typename CoeffsType, int Side> 438c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtypename internal::matrix_type_times_scalar_type<typename VectorsType::Scalar,OtherDerived>::Type operator*(const MatrixBase<OtherDerived>& other, const HouseholderSequence<VectorsType,CoeffsType,Side>& h) 439c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 440c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename internal::matrix_type_times_scalar_type<typename VectorsType::Scalar,OtherDerived>::Type 441c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath res(other.template cast<typename internal::matrix_type_times_scalar_type<typename VectorsType::Scalar,OtherDerived>::ResultScalar>()); 442c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath h.applyThisOnTheRight(res); 443c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return res; 444c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 445c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 446c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath/** \ingroup Householder_Module \householder_module 447c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \brief Convenience function for constructing a Householder sequence. 448c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \returns A HouseholderSequence constructed from the specified arguments. 449c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 450c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename VectorsType, typename CoeffsType> 451c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan KamathHouseholderSequence<VectorsType,CoeffsType> householderSequence(const VectorsType& v, const CoeffsType& h) 452c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 453c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return HouseholderSequence<VectorsType,CoeffsType,OnTheLeft>(v, h); 454c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 455c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 456c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath/** \ingroup Householder_Module \householder_module 457c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \brief Convenience function for constructing a Householder sequence. 458c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \returns A HouseholderSequence constructed from the specified arguments. 459c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \details This function differs from householderSequence() in that the template argument \p OnTheSide of 460c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * the constructed HouseholderSequence is set to OnTheRight, instead of the default OnTheLeft. 461c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 462c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename VectorsType, typename CoeffsType> 463c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan KamathHouseholderSequence<VectorsType,CoeffsType,OnTheRight> rightHouseholderSequence(const VectorsType& v, const CoeffsType& h) 464c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 465c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return HouseholderSequence<VectorsType,CoeffsType,OnTheRight>(v, h); 466c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 467c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 468c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} // end namespace Eigen 469c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 470c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif // EIGEN_HOUSEHOLDER_SEQUENCE_H 471