IterationNumberBasedTimeStepping.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 <vector>
 
#include "MultiplyerInterpolationType.h"
#include "TimeStepAlgorithm.h"
 
namespace NumLib
{
class IterationNumberBasedTimeStepping final : public TimeStepAlgorithm
{
public:
    IterationNumberBasedTimeStepping(
        double const t_initial,
        double const t_end,
        double const min_dt,
        double const max_dt,
        double const initial_dt,
        MultiplyerInterpolationType const multiplier_interpolation_type,
        std::vector<int>&& iter_times_vector,
        std::vector<double>&& multiplier_vector,
        std::vector<double> const& fixed_times_for_output);
 
    ~IterationNumberBasedTimeStepping() override = default;
 
    double next(double solution_error, int number_iterations,
                NumLib::TimeStep& ts_previous,
                NumLib::TimeStep& ts_current) override;
 
    bool isSolutionErrorComputationNeeded() const override { return true; }
 
    bool canReduceTimestepSize(
        NumLib::TimeStep const& timestep_previous,
        NumLib::TimeStep const& timestep_current) const override;
 
private:
    double getNextTimeStepSize(NumLib::TimeStep const& ts_previous,
                               NumLib::TimeStep const& ts_current) const;

    const std::vector<int> _iter_times_vector;
    const std::vector<double> _multiplier_vector;
    const double _min_dt;
    const double _max_dt;
    const double _initial_dt;
    const MultiplyerInterpolationType _multiplier_interpolation_type;
    const int _max_iter;
    int _iter_times = 0;
 
    std::vector<double> const _fixed_times_for_output;
};

double findMultiplier(
    int const number_iterations, bool const current_time_step_is_accepted,
    std::vector<int> const& nonlinear_iteration_numbers,
    std::vector<double> const& multipliers,
    MultiplyerInterpolationType const multiplier_interpolation_type);
 
}  // namespace NumLib