Revision 571ec23aca4dfcdddb19441acb74ac1fd89acf1c authored by Thomas Fischer on 25 November 2021, 14:33:49 UTC, committed by Thomas Fischer on 26 November 2021, 09:22:04 UTC
1 parent a8ffcfa
ConvergenceCriterionResidual.h
/**
* \file
* \copyright
* Copyright (c) 2012-2021, OpenGeoSys Community (http://www.opengeosys.org)
* Distributed under a Modified BSD License.
* See accompanying file LICENSE.txt or
* http://www.opengeosys.org/project/license
*
*/
#pragma once
#include <optional>
#include "ConvergenceCriterion.h"
#include "MathLib/LinAlg/LinAlgEnums.h"
namespace NumLib
{
//! Convergence criterion applying a single absolute or relative tolerance to
//! the whole residual vector.
//!
//! A check of the solution increment is not done.
//! If both an absolute and a relative tolerance are specified, at least one of
//! them has to be satisfied.
class ConvergenceCriterionResidual final : public ConvergenceCriterion
{
public:
ConvergenceCriterionResidual(std::optional<double>&& absolute_tolerance,
std::optional<double>&& relative_tolerance,
const MathLib::VecNormType norm_type);
bool hasDeltaXCheck() const override { return true; }
bool hasResidualCheck() const override { return true; }
/// The function will only do diagnostic output and no actual check of the
/// solution increment is made
void checkDeltaX(const GlobalVector& minus_delta_x,
GlobalVector const& x) override;
void checkResidual(const GlobalVector& residual) override;
private:
const std::optional<double> _abstol;
const std::optional<double> _reltol;
double _residual_norm_0;
};
std::unique_ptr<ConvergenceCriterionResidual>
createConvergenceCriterionResidual(BaseLib::ConfigTree const& config);
} // namespace NumLib
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...