FixedTimeStepping.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#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);
 
    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_;

    double last_prescribed_dt_ = 0.0;

    bool reuse_of_last_dt_reported_ = false;
};
 
std::size_t findDeltatInterval(Time const& t_initial,
                               std::vector<double> const& delta_ts,
                               Time const& fixed_output_time);
}  // namespace NumLib