10#ifndef EIGEN_BASIC_PRECONDITIONERS_H
11#define EIGEN_BASIC_PRECONDITIONERS_H
35template <
typename _Scalar>
38 typedef _Scalar Scalar;
49 template<
typename MatType>
58 template<
typename MatType>
64 template<
typename MatType>
68 for(
int j=0; j<mat.outerSize(); ++j)
70 typename MatType::InnerIterator it(mat,j);
71 while(it && it.index()!=j) ++it;
72 if(it && it.index()==j && it.value()!=Scalar(0))
81 template<
typename MatType>
88 template<
typename Rhs,
typename Dest>
99 &&
"DiagonalPreconditioner::solve(): invalid number of rows of the right hand side matrix b");
127template <
typename _Scalar>
130 typedef _Scalar Scalar;
138 template<
typename MatType>
144 template<
typename MatType>
150 template<
typename MatType>
154 m_invdiag.
resize(mat.cols());
155 if(MatType::IsRowMajor)
158 for(
Index j=0; j<mat.outerSize(); ++j)
160 for(
typename MatType::InnerIterator it(mat,j); it; ++it)
163 for(
Index j=0; j<mat.cols(); ++j)
165 m_invdiag(j) = RealScalar(1)/
numext::real(m_invdiag(j));
169 for(
Index j=0; j<mat.outerSize(); ++j)
171 RealScalar sum = mat.col(j).squaredNorm();
172 if(sum>RealScalar(0))
173 m_invdiag(j) = RealScalar(1)/sum;
175 m_invdiag(j) = RealScalar(1);
182 template<
typename MatType>
206 template<
typename MatrixType>
209 template<
typename MatrixType>
212 template<
typename MatrixType>
215 template<
typename MatrixType>
218 template<
typename Rhs>
219 inline const Rhs&
solve(
const Rhs&
b)
const {
return b; }
EIGEN_DEVICE_FUNC RealReturnType real() const
Definition: CommonCwiseUnaryOps.h:100
#define EIGEN_NOEXCEPT
Definition: Macros.h:1428
#define EIGEN_CONSTEXPR
Definition: Macros.h:797
#define eigen_assert(x)
Definition: Macros.h:1047
A preconditioner based on the digonal entries.
Definition: BasicPreconditioners.h:37
bool m_isInitialized
Definition: BasicPreconditioners.h:107
Vector::StorageIndex StorageIndex
Definition: BasicPreconditioners.h:41
void _solve_impl(const Rhs &b, Dest &x) const
Definition: BasicPreconditioners.h:89
DiagonalPreconditioner & compute(const MatType &mat)
Definition: BasicPreconditioners.h:82
DiagonalPreconditioner & analyzePattern(const MatType &)
Definition: BasicPreconditioners.h:59
EIGEN_CONSTEXPR Index rows() const EIGEN_NOEXCEPT
Definition: BasicPreconditioners.h:55
ComputationInfo info()
Definition: BasicPreconditioners.h:103
const Solve< DiagonalPreconditioner, Rhs > solve(const MatrixBase< Rhs > &b) const
Definition: BasicPreconditioners.h:95
EIGEN_CONSTEXPR Index cols() const EIGEN_NOEXCEPT
Definition: BasicPreconditioners.h:56
Vector m_invdiag
Definition: BasicPreconditioners.h:106
@ MaxColsAtCompileTime
Definition: BasicPreconditioners.h:44
@ ColsAtCompileTime
Definition: BasicPreconditioners.h:43
DiagonalPreconditioner & factorize(const MatType &mat)
Definition: BasicPreconditioners.h:65
DiagonalPreconditioner(const MatType &mat)
Definition: BasicPreconditioners.h:50
DiagonalPreconditioner()
Definition: BasicPreconditioners.h:47
A naive preconditioner which approximates any matrix as the identity matrix.
Definition: BasicPreconditioners.h:201
IdentityPreconditioner & analyzePattern(const MatrixType &)
Definition: BasicPreconditioners.h:210
ComputationInfo info()
Definition: BasicPreconditioners.h:221
IdentityPreconditioner(const MatrixType &)
Definition: BasicPreconditioners.h:207
const Rhs & solve(const Rhs &b) const
Definition: BasicPreconditioners.h:219
IdentityPreconditioner & compute(const MatrixType &)
Definition: BasicPreconditioners.h:216
IdentityPreconditioner & factorize(const MatrixType &)
Definition: BasicPreconditioners.h:213
IdentityPreconditioner()
Definition: BasicPreconditioners.h:204
Jacobi preconditioner for LeastSquaresConjugateGradient.
Definition: BasicPreconditioners.h:129
LeastSquareDiagonalPreconditioner & factorize(const MatType &mat)
Definition: BasicPreconditioners.h:151
LeastSquareDiagonalPreconditioner()
Definition: BasicPreconditioners.h:136
ComputationInfo info()
Definition: BasicPreconditioners.h:188
LeastSquareDiagonalPreconditioner & analyzePattern(const MatType &)
Definition: BasicPreconditioners.h:145
LeastSquareDiagonalPreconditioner(const MatType &mat)
Definition: BasicPreconditioners.h:139
LeastSquareDiagonalPreconditioner & compute(const MatType &mat)
Definition: BasicPreconditioners.h:183
Base class for all dense matrices, vectors, and expressions.
Definition: MatrixBase.h:50
EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE void resize(Index rows, Index cols)
Resizes *this to a rows x cols matrix.
Definition: PlainObjectBase.h:271
EIGEN_DEVICE_FUNC Derived & setZero(Index size)
Resizes to the given size, and sets all coefficients in this expression to zero.
Definition: CwiseNullaryOp.h:562
Pseudo expression representing a solving operation.
Definition: Solve.h:63
ComputationInfo
Enum for reporting the status of a computation.
Definition: Constants.h:440
@ Success
Computation was successful.
Definition: Constants.h:442
EIGEN_DEVICE_FUNC bool abs2(bool x)
Definition: MathFunctions.h:1292
Namespace containing all symbols from the Eigen library.
Definition: MatrixExponential.h:16
EIGEN_DEFAULT_DENSE_INDEX_TYPE Index
The Index type as used for the API.
Definition: Meta.h:74
const int Dynamic
This value means that a positive quantity (e.g., a size) is not known at compile-time,...
Definition: Constants.h:22
Holds information about the various numeric (i.e.
Definition: NumTraits.h:233