ConstantPreconditioner.cpp

Go to the documentation of this file.

#include "acceleration/impl/ConstantPreconditioner.hpp"
#include <algorithm>
#include <utility>
 
#include "logging/LogMacros.hpp"
#include "utils/assertion.hpp"
 
namespace precice::acceleration::impl {
 
ConstantPreconditioner::ConstantPreconditioner(std::vector<double> factors)
    : Preconditioner(-1),
      _factors(std::move(factors))
{
}
 
void ConstantPreconditioner::initialize(std::vector<size_t> &svs)
{
  PRECICE_TRACE();
  Preconditioner::initialize(svs);
 
  // is always constant by definition
  _frozen = true;
  PRECICE_ASSERT(_maxNonConstTimeWindows == -1, _maxNonConstTimeWindows);
 
  PRECICE_ASSERT(_factors.size() == _subVectorSizes.size());
 
  int offset = 0;
  for (size_t k = 0; k < _subVectorSizes.size(); k++) {
    for (size_t i = 0; i < _subVectorSizes[k]; i++) {
      _weights[i + offset]    = 1.0 / _factors[k];
      _invWeights[i + offset] = _factors[k];
    }
    offset += _subVectorSizes[k];
  }
}
 
void ConstantPreconditioner::_update_(bool                   timeWindowComplete,
                                      const Eigen::VectorXd &oldValues,
                                      const Eigen::VectorXd &res)
{
 
  //nothing to do here
}
 
} // namespace precice::acceleration::impl