FixedTimeStepping.h

Go to the documentation of this file.

 
#pragma once
 
#include <memory>
#include <vector>
 
#include "TimeStepAlgorithm.h"
 
namespace NumLib
{
using RepeatDtPair = std::tuple<std::size_t, double>;
 
class FixedTimeStepping final : public TimeStepAlgorithm
{
public:
    FixedTimeStepping(double t0, double t_end, double dt);
 
    FixedTimeStepping(double t0, double tn,
                      std::vector<RepeatDtPair> 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;
 
    static bool areRepeatDtPairsValid(
        std::vector<RepeatDtPair> const& repeat_dt_pairs);
 
private:
    std::vector<double> _dt_vector;
};
 
std::size_t findDeltatInterval(Time const& t_initial,
                               std::vector<double> const& delta_ts,
                               Time const& fixed_output_time);
}  // namespace NumLib