diff options
| author |  Benoit Jacob <jacob.benoit.1@gmail.com> | 2009-06-19 17:01:32 +0200 |
|---|---|---|
| committer | Benoit Jacob <jacob.benoit.1@gmail.com> | 2009-06-19 17:01:32 +0200 |
| commit | a57325e971c749a42ccfc200772054a72099414c (patch) | |
| tree | b0dc2282dacac4233129a58745b88b56359cdb8d /Eigen/src/Cholesky | |
| parent | c6e81869d0a4ba18eac387199fd9a5c9ac26e050 (diff) | |
fix #14: make llt::solve() and also ldlt::solve() work with uninitialized result
Diffstat (limited to 'Eigen/src/Cholesky')
| -rw-r--r-- | Eigen/src/Cholesky/LDLT.h | 8 | ||||
| -rw-r--r-- | Eigen/src/Cholesky/LLT.h | 8 |
2 files changed, 8 insertions, 8 deletions
diff --git a/Eigen/src/Cholesky/LDLT.h b/Eigen/src/Cholesky/LDLT.h index 039e542bc..94660245e 100644 --- a/Eigen/src/Cholesky/LDLT.h +++ b/Eigen/src/Cholesky/LDLT.h @@ -117,8 +117,8 @@ template<typename MatrixType> class LDLT return m_sign == -1; } - template<typename RhsDerived, typename ResDerived> - bool solve(const MatrixBase<RhsDerived> &b, MatrixBase<ResDerived> *result) const; + template<typename RhsDerived, typename ResultType> + bool solve(const MatrixBase<RhsDerived> &b, ResultType *result) const; template<typename Derived> bool solveInPlace(MatrixBase<Derived> &bAndX) const; @@ -248,9 +248,9 @@ void LDLT<MatrixType>::compute(const MatrixType& a) * \sa LDLT::solveInPlace(), MatrixBase::ldlt() */ template<typename MatrixType> -template<typename RhsDerived, typename ResDerived> +template<typename RhsDerived, typename ResultType> bool LDLT<MatrixType> -::solve(const MatrixBase<RhsDerived> &b, MatrixBase<ResDerived> *result) const +::solve(const MatrixBase<RhsDerived> &b, ResultType *result) const { ei_assert(m_isInitialized && "LDLT is not initialized."); const int size = m_matrix.rows(); diff --git a/Eigen/src/Cholesky/LLT.h b/Eigen/src/Cholesky/LLT.h index cb2524e9d..215cda346 100644 --- a/Eigen/src/Cholesky/LLT.h +++ b/Eigen/src/Cholesky/LLT.h @@ -87,8 +87,8 @@ template<typename MatrixType> class LLT return m_matrix; } - template<typename RhsDerived, typename ResDerived> - bool solve(const MatrixBase<RhsDerived> &b, MatrixBase<ResDerived> *result) const; + template<typename RhsDerived, typename ResultType> + bool solve(const MatrixBase<RhsDerived> &b, ResultType *result) const; template<typename Derived> bool solveInPlace(MatrixBase<Derived> &bAndX) const; @@ -165,8 +165,8 @@ void LLT<MatrixType>::compute(const MatrixType& a) * \sa LLT::solveInPlace(), MatrixBase::llt() */ template<typename MatrixType> -template<typename RhsDerived, typename ResDerived> -bool LLT<MatrixType>::solve(const MatrixBase<RhsDerived> &b, MatrixBase<ResDerived> *result) const +template<typename RhsDerived, typename ResultType> +bool LLT<MatrixType>::solve(const MatrixBase<RhsDerived> &b, ResultType *result) const { ei_assert(m_isInitialized && "LLT is not initialized."); const int size = m_matrix.rows(); |