SolverSequential.h

#ifndef HYDROBRICKS_SOLVER_SEQUENTIAL_H
#define HYDROBRICKS_SOLVER_SEQUENTIAL_H
 
#include "Includes.h"
#include "Solver.h"
 
class Brick;
 
class SolverSequential : public Solver {
  public:
    void InitializeContainers() override;
 
    bool Solve(double timeStepInDays) final;
 
  protected:
    axd _rates;
 
    virtual void ComputeBrickRates(Brick* brick, double content, double inflow, double timeStepInDays,
                                   int iRateStart) = 0;
 
    static double TotalRateAt(Brick* brick, double* contentDelta, double offset);
 
    static double StoreRatesAndTotalAt(Brick* brick, double* contentDelta, double offset, vecDouble& rates);
 
    static bool SumAffineResponse(Brick* brick, double& rate, double& offset);
};
 
#endif  // HYDROBRICKS_SOLVER_SEQUENTIAL_H