IQNILSAcceleration.hpp

Go to the documentation of this file.

#pragma once
 
#include <Eigen/Core>
#include <map>
#include <vector>
#include "acceleration/Acceleration.hpp"
#include "acceleration/BaseQNAcceleration.hpp"
#include "acceleration/impl/SharedPointer.hpp"
 
namespace precice {
namespace acceleration {
 
class IQNILSAcceleration : public BaseQNAcceleration {
public:
  IQNILSAcceleration(
      double                  initialRelaxation,
      bool                    forceInitialRelaxation,
      int                     maxIterationsUsed,
      int                     pastTimeWindowsReused,
      int                     filter,
      double                  singularityLimit,
      std::vector<int>        dataIDs,
      impl::PtrPreconditioner preconditioner);
 
  virtual ~IQNILSAcceleration() {}
 
  virtual void initialize(const DataMap &cplData);
 
  virtual void specializedIterationsConverged(const DataMap &cplData);
 
private:
  std::map<int, Eigen::VectorXd> _secondaryOldXTildes;
 
  // @brief Secondary data x-tilde deltas.
  //
  // Stores x-tilde deltas for data not involved in least-squares computation.
  std::map<int, Eigen::MatrixXd> _secondaryMatricesW;
  std::map<int, Eigen::MatrixXd> _secondaryMatricesWBackup;
 
  virtual void updateDifferenceMatrices(const DataMap &cplData);
 
  virtual void computeQNUpdate(const DataMap &cplData, Eigen::VectorXd &xUpdate);
 
  virtual void computeUnderrelaxationSecondaryData(const DataMap &cplData);
 
  virtual void removeMatrixColumn(int columnIndex);
};
} // namespace acceleration
} // namespace precice