CreateEvolutionaryPIDcontroller.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "CreateEvolutionaryPIDcontroller.h"
 
#include "BaseLib/Algorithm.h"
#include "BaseLib/ConfigTree.h"
#include "EvolutionaryPIDcontroller.h"
#include "NumLib/TimeStepping/Time.h"
#include "TimeStepAlgorithm.h"
 
namespace NumLib
{
class TimeStepAlgorithm;
 
EvolutionaryPIDcontrollerParameters parseEvolutionaryPIDcontroller(
    BaseLib::ConfigTree const& config)
{
    config.checkConfigParameter("type", "EvolutionaryPIDcontroller");

    auto const t0 = config.getConfigParameter<double>("t_initial");
    auto const t_end = config.getConfigParameter<double>("t_end");

    auto const h0 = config.getConfigParameter<double>("dt_guess");

    auto const h_min = config.getConfigParameter<double>("dt_min");
    auto const h_max = config.getConfigParameter<double>("dt_max");
    auto const rel_h_min = config.getConfigParameter<double>("rel_dt_min");
    auto const rel_h_max = config.getConfigParameter<double>("rel_dt_max");

    auto const tol = config.getConfigParameter<double>("tol");
 
    return {t0, t_end, h0, h_min, h_max, rel_h_min, rel_h_max, tol};
}
 
std::unique_ptr<TimeStepAlgorithm> createEvolutionaryPIDcontroller(
    EvolutionaryPIDcontrollerParameters const& config,
    std::vector<double> const& fixed_times_for_output)
{
    if (Time(config.t_end) <= Time(config.t0))
    {
        OGS_FATAL(
            "Evolutionary PID controller timestepping: end time ({}) must be "
            "larger than initial time ({})",
            config.t_end,
            config.t0);
    }
 
    return std::make_unique<EvolutionaryPIDcontroller>(config.t0,
                                                       config.t_end,
                                                       config.h0,
                                                       config.h_min,
                                                       config.h_max,
                                                       config.rel_h_min,
                                                       config.rel_h_max,
                                                       config.tol,
                                                       fixed_times_for_output);
}
}  // end of namespace NumLib