Solver.h

#ifndef HYDROBRICKS_SOLVER_H
#define HYDROBRICKS_SOLVER_H
 
#include <memory>
 
#include "Includes.h"
#include "SettingsModel.h"
 
class Processor;
 
class Solver {
  public:
    explicit Solver();
 
    virtual ~Solver() = default;
 
    virtual bool Solve(double timeStepInDays) = 0;
 
    static std::unique_ptr<Solver> Factory(const SolverSettings& solverSettings);
 
    void Connect(Processor* processor) {
        assert(processor);
        _processor = processor;
    }
 
    void InitializeContainers();
 
  protected:
    Processor* _processor;
    axxd _stateVariableChanges;
    axxd _changeRates;
    int _nIterations;
    double _timeStepInDays;
 
    void SaveStateVariables(int col);
 
    void ComputeChangeRates(int col, bool applyConstraints = true);
 
    void ApplyConstraintsFor(int col);
 
    void ResetStateVariableChanges();
 
    void SetStateVariablesToIteration(int col);
 
    void SetStateVariablesToAvgOf(int col1, int col2);
 
    void ApplyProcesses(int col) const;
 
    void ApplyProcesses(const axd& changeRates) const;
 
    void Finalize() const;
 
  private:
};
 
#endif  // HYDROBRICKS_SOLVER_H