60 #include "../../../../Eigen/src/Core/util/NonMPL2.h"
62 #ifndef EIGEN_ITERATION_CONTROLLER_H
63 #define EIGEN_ITERATION_CONTROLLER_H
82 double m_resminreach, m_resadd;
92 m_resminreach = 1E50; m_resadd = 0.0;
99 void operator ++(
int) {
m_nit++; m_written =
false; m_resadd +=
m_res; }
100 void operator ++() { (*this)++; }
102 bool first() {
return m_nit == 0; }
105 int noiseLevel()
const {
return m_noise; }
106 void setNoiseLevel(
int n) {
m_noise = n; }
109 double maxResidual()
const {
return m_resmax; }
110 void setMaxResidual(
double r) {
m_resmax = r; }
112 double residual()
const {
return m_res; }
115 void setCallback(
void (*t)(
const IterationController&))
120 size_t iteration()
const {
return m_nit; }
121 void setIteration(
size_t i) {
m_nit = i; }
123 size_t maxIterarions()
const {
return m_maxiter; }
124 void setMaxIterations(
size_t i) {
m_maxiter = i; }
126 double rhsNorm()
const {
return m_rhsn; }
127 void setRhsNorm(
double r) {
m_rhsn = r; }
130 bool converged(
double nr)
132 m_res = internal::abs(nr);
133 m_resminreach = (std::min)(m_resminreach,
m_res);
136 template<
typename VectorType>
bool converged(
const VectorType &v)
137 {
return converged(v.squaredNorm()); }
139 bool finished(
double nr)
141 if (m_callback) m_callback(*
this);
149 template <
typename VectorType>
150 bool finished(
const MatrixBase<VectorType> &v)
151 {
return finished(
double(v.squaredNorm())); }
157 #endif // EIGEN_ITERATION_CONTROLLER_H