FixedTimeStepping.h

Go to the documentation of this file.

 
#pragma once
 
#include <memory>
#include <vector>
 
#include "TimeStepAlgorithm.h"
 
namespace NumLib
{
class FixedTimeStepping final : public TimeStepAlgorithm
{
public:
    FixedTimeStepping(double t0, double t_end, double dt);
 
    FixedTimeStepping(
        double t0, double tn,
        std::vector<std::pair<std::size_t, double>> const& repeat_dt_pairs,
        std::vector<double> const& fixed_times_for_output);
 
    std::tuple<bool, double> next(double solution_error, int number_iterations,
                                  NumLib::TimeStep& ts_previous,
                                  NumLib::TimeStep& ts_current) override;
 
    void resetCurrentTimeStep(const double dt, TimeStep& /*ts_previous*/,
                              TimeStep& /*ts_current*/) override
    {
        _dt_vector.push_back(dt);
    }
 
private:
    std::vector<double> _dt_vector;
};
 
std::size_t findDeltatInterval(double const t_initial,
                               std::vector<double> const& delta_ts,
                               double const fixed_output_time);
}  // namespace NumLib