Homec4science
Diffusion GooseFEM (main) Browse

rGOOSEFEM/include/xGooseFEM/Iterate.hpp
6bb8640345c8main

Last Change
Raw File

Iterate.hpp
View Options

/* =================================================================================================
(c - GPLv3) T.W.J. de Geus (Tom) | tom@geus.me | www.geus.me | github.com/tdegeus/GooseFEM
================================================================================================= */
#ifndef XGOOSEFEM_ITERATE_CPP
#define XGOOSEFEM_ITERATE_CPP
// -------------------------------------------------------------------------------------------------
#include "Iterate.h"
// ======================================= GooseFEM::Iterate =======================================
namespace xGooseFEM {
namespace Iterate {
// -------------------------------------------------------------------------------------------------
inline StopList::StopList(size_t n)
{
// resize list
m_res.resize(n);
// set all residual to infinity
reset();
}
// -------------------------------------------------------------------------------------------------
inline void StopList::reset()
{
// set all residual to infinity
for ( auto &i : m_res )
i = std::numeric_limits<double>::infinity();
}
// -------------------------------------------------------------------------------------------------
inline void StopList::reset(size_t n)
{
// resize list
m_res.resize(n);
// set all residual to infinity
reset();
}
// -------------------------------------------------------------------------------------------------
inline bool StopList::stop(double res, double tol)
{
// move residual one place back (forgetting the first)
for ( size_t i = 1 ; i < m_res.size() ; ++i )
m_res[i-1] = m_res[i];
// add new residual to the end
m_res[m_res.size()-1] = res;
// check for convergence: all residuals should be below the tolerance
for ( size_t i = 0 ; i < m_res.size() ; ++i )
if ( m_res[i] > tol )
return false;
// check for convergence: all residuals should be decreasing
for ( size_t i = 1 ; i < m_res.size() ; ++i )
if ( m_res[i] > m_res[i-1] )
return false;
// all checks passed: signal convergence
return true;
}
// -------------------------------------------------------------------------------------------------
}} // namespace ...
// =================================================================================================
#endif
c4science · Help