GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-07-07 21:51:47 +08:00
Code Issues Packages Projects Releases Wiki Activity

Added support for custom-scalars

Browse Source
This commit is contained in:
Pavel Holoborodko 2013-09-02 19:09:39 +09:00
parent ace2ed7b87
commit ebd6a7a46c
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
Eigen/src/SparseQR/SparseQR.h
View File

@ -205,7 +205,7 @@ class SparseQR
/** \brief Reports whether previous computation was successful. /** \brief Reports whether previous computation was successful.
* *
* \returns \c Success if computation was succesful, * \returns \c Success if computation was successful,
* \c NumericalIssue if the QR factorization reports a numerical problem * \c NumericalIssue if the QR factorization reports a numerical problem
* \c InvalidInput if the input matrix is invalid * \c InvalidInput if the input matrix is invalid
* *
@ -256,7 +256,7 @@ class SparseQR
/** \brief Preprocessing step of a QR factorization /** \brief Preprocessing step of a QR factorization
* *
* In this step, the fill-reducing permutation is computed and applied to the columns of A * In this step, the fill-reducing permutation is computed and applied to the columns of A
* and the column elimination tree is computed as well. Only the sparcity pattern of \a mat is exploited. * and the column elimination tree is computed as well. Only the sparsity pattern of \a mat is exploited.
* *
* \note In this step it is assumed that there is no empty row in the matrix \a mat. * \note In this step it is assumed that there is no empty row in the matrix \a mat.
*/ */
@ -292,7 +292,7 @@ void SparseQR<MatrixType,OrderingType>::analyzePattern(const MatrixType& mat)
/** \brief Performs the numerical QR factorization of the input matrix /** \brief Performs the numerical QR factorization of the input matrix
* *
* The function SparseQR::analyzePattern(const MatrixType&) must have been called beforehand with * The function SparseQR::analyzePattern(const MatrixType&) must have been called beforehand with
* a matrix having the same sparcity pattern than \a mat. * a matrix having the same sparsity pattern than \a mat.
* *
* \param mat The sparse column-major matrix * \param mat The sparse column-major matrix
*/ */
@ -445,7 +445,8 @@ void SparseQR<MatrixType,OrderingType>::factorize(const MatrixType& mat)
} }
else else
{ {
beta = std::sqrt(numext::abs2(c0) + sqrNorm); using std::sqrt;
beta = sqrt(numext::abs2(c0) + sqrNorm);
if(numext::real(c0) >= RealScalar(0)) if(numext::real(c0) >= RealScalar(0))
beta = -beta; beta = -beta;
tval(Qidx(0)) = 1; tval(Qidx(0)) = 1;