AbsoluteConvergenceMeasure.hpp

Go to the documentation of this file.

#pragma once
 
#include <Eigen/Core>
#include <iomanip>
#include <ostream>
#include <string>
#include "ConvergenceMeasure.hpp"
#include "logging/Logger.hpp"
#include "utils/IntraComm.hpp"
 
namespace precice {
namespace cplscheme {
namespace tests {
class AbsoluteConvergenceMeasureTest;
}
} // namespace cplscheme
} // namespace precice
 
namespace precice {
namespace cplscheme {
namespace impl {
 
class AbsoluteConvergenceMeasure : public ConvergenceMeasure {
public:
  explicit AbsoluteConvergenceMeasure(double convergenceLimit);
 
  virtual ~AbsoluteConvergenceMeasure(){};
 
  virtual void newMeasurementSeries()
  {
    _isConvergence = false;
  }
 
  virtual void measure(
      const Eigen::VectorXd &oldValues,
      const Eigen::VectorXd &newValues)
  {
    _normDiff      = utils::IntraComm::l2norm(newValues - oldValues);
    _isConvergence = _normDiff <= _convergenceLimit;
  }
 
  virtual bool isConvergence() const
  {
    return _isConvergence;
  }
 
  virtual std::string printState(const std::string &dataName)
  {
    std::ostringstream os;
    os << "absolute convergence measure: ";
    os << "two-norm diff of data \"" << dataName << "\" = ";
    os << std::scientific << std::setprecision(2) << _normDiff;
    os << ", limit = " << _convergenceLimit;
    os << ", conv = ";
    if (_isConvergence)
      os << "true";
    else
      os << "false";
    return os.str();
  }
 
  virtual double getNormResidual()
  {
    return _normDiff;
  }
 
  virtual std::string getAbbreviation() const
  {
    return "Abs";
  }
 
private:
  logging::Logger _log{"cplscheme::AbsoluteConvergenceMeasure"};
 
  double _convergenceLimit;
 
  double _normDiff = 0;
 
  bool _isConvergence = false;
};
} // namespace impl
} // namespace cplscheme
} // namespace precice