GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-14 04:35:57 +08:00
Code Issues Packages Projects Releases Wiki Activity

clean a bit Hessenberg and make sure the subdiagonal coeff is real even for 2x2 matrices

Browse Source
This commit is contained in:
Gael Guennebaud 2009-12-23 12:30:05 +01:00
parent 791bac25f2
commit c90ccd089a
1 changed files with 8 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
Eigen/src/Eigenvalues/HessenbergDecomposition.h
View File

@ -1,7 +1,7 @@
// This file is part of Eigen, a lightweight C++ template library // This file is part of Eigen, a lightweight C++ template library
// for linear algebra. // for linear algebra.
// //
// Copyright (C) 2008 Gael Guennebaud <g.gael@free.fr> // Copyright (C) 2008-2009 Gael Guennebaud <g.gael@free.fr>
// //
// Eigen is free software; you can redistribute it and/or // Eigen is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public // modify it under the terms of the GNU Lesser General Public
@ -55,25 +55,21 @@ template<typename _MatrixType> class HessenbergDecomposition
}; };
typedef Matrix<Scalar, SizeMinusOne, 1> CoeffVectorType; typedef Matrix<Scalar, SizeMinusOne, 1> CoeffVectorType;
typedef Matrix<RealScalar, Size, 1> DiagonalType;
typedef Matrix<RealScalar, SizeMinusOne, 1> SubDiagonalType;
typedef typename Diagonal<MatrixType,0>::RealReturnType DiagonalReturnType;
typedef typename Diagonal<
Block<MatrixType,SizeMinusOne,SizeMinusOne>,0 >::RealReturnType SubDiagonalReturnType;
/** This constructor initializes a HessenbergDecomposition object for /** This constructor initializes a HessenbergDecomposition object for
* further use with HessenbergDecomposition::compute() * further use with HessenbergDecomposition::compute()
*/ */
HessenbergDecomposition(int size = Size==Dynamic ? 2 : Size) HessenbergDecomposition(int size = Size==Dynamic ? 2 : Size)
: m_matrix(size,size), m_hCoeffs(size-1) : m_matrix(size,size)
{} {
if(size>1)
m_hCoeffs.resize(size-1);
}
HessenbergDecomposition(const MatrixType& matrix) HessenbergDecomposition(const MatrixType& matrix)
: m_matrix(matrix) : m_matrix(matrix)
{ {
if(matrix.rows()<=2) if(matrix.rows()<2)
return; return;
m_hCoeffs.resize(matrix.rows()-1,1); m_hCoeffs.resize(matrix.rows()-1,1);
_compute(m_matrix, m_hCoeffs); _compute(m_matrix, m_hCoeffs);
@ -86,7 +82,7 @@ template<typename _MatrixType> class HessenbergDecomposition
void compute(const MatrixType& matrix) void compute(const MatrixType& matrix)
{ {
m_matrix = matrix; m_matrix = matrix;
if(matrix.rows()<=2) if(matrix.rows()<2)
return; return;
m_hCoeffs.resize(matrix.rows()-1,1); m_hCoeffs.resize(matrix.rows()-1,1);
_compute(m_matrix, m_hCoeffs); _compute(m_matrix, m_hCoeffs);