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

reduce epsilon in QR

Browse Source
This commit is contained in:
Gael Guennebaud 2009-01-29 16:11:46 +00:00
parent cf0857c44d
commit 8e0ec3c62b
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Eigen/src/QR/QR.h
View File

@ -185,6 +185,7 @@ void QR<MatrixType>::_compute(const MatrixType& matrix)
m_qr = matrix; m_qr = matrix;
int rows = matrix.rows(); int rows = matrix.rows();
int cols = matrix.cols(); int cols = matrix.cols();
RealScalar eps2 = precision<RealScalar>()*precision<RealScalar>();
for (int k = 0; k < cols; ++k) for (int k = 0; k < cols; ++k)
{ {
@ -209,7 +210,8 @@ void QR<MatrixType>::_compute(const MatrixType& matrix)
m_hCoeffs.coeffRef(k) = 0; m_hCoeffs.coeffRef(k) = 0;
} }
} }
else if ( (!ei_isMuchSmallerThan(beta=m_qr.col(k).end(remainingSize-1).squaredNorm(),static_cast<Scalar>(1))) ) // FIXME what about ei_imag(v0) ?? else if ((beta=m_qr.col(k).end(remainingSize-1).squaredNorm())>eps2)
// FIXME what about ei_imag(v0) ??
{ {
// form k-th Householder vector // form k-th Householder vector
beta = ei_sqrt(ei_abs2(v0)+beta); beta = ei_sqrt(ei_abs2(v0)+beta);