ResidualSumPreconditioner.hpp

Go to the documentation of this file.

#pragma once
 
#include <Eigen/Core>
#include <stddef.h>
#include <string>
#include <vector>
#include "acceleration/impl/Preconditioner.hpp"
#include "logging/Logger.hpp"
 
namespace precice::acceleration::impl {

class ResidualSumPreconditioner : public Preconditioner {
public:
  ResidualSumPreconditioner(int maxNonConstTimeWindows, bool preconditionerUpdateOnThreshold);
  ~ResidualSumPreconditioner() override = default;
 
  void initialize(std::vector<size_t> &svs) override;
 
private:
  void _update_(bool timeWindowComplete, const Eigen::VectorXd &oldValues, const Eigen::VectorXd &res) override;
 
  logging::Logger _log{"acceleration::ResidualSumPreconditioner"};
 
  std::vector<double> _residualSum;
  std::vector<double> _previousResidualSum;
  bool                _firstTimeWindow = true;
  bool                _preconditionerUpdateOnThreshold;
};
 
} // namespace precice::acceleration::impl