EvolutionaryPIDcontroller.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
{
class EvolutionaryPIDcontroller final : public TimeStepAlgorithm
{
public:
    EvolutionaryPIDcontroller(const double t0, const double t_end,
                              const double h0, const double h_min,
                              const double h_max, const double rel_h_min,
                              const double rel_h_max, const double tol,
                              std::vector<double> const& fixed_times_for_output)
        : TimeStepAlgorithm(t0, t_end),
          _h0(h0),
          _h_min(h_min),
          _h_max(h_max),
          _rel_h_min(rel_h_min),
          _rel_h_max(rel_h_max),
          _tol(tol),
          _e_n_minus1(0.),
          _e_n_minus2(0.),
          _fixed_times_for_output(fixed_times_for_output)
    {
    }
 
    double next(double solution_error,
                int number_iterations,
                NumLib::TimeStep& timestep_previous,
                NumLib::TimeStep& timestep_current) override;
 
    bool isSolutionErrorComputationNeeded() const override { return true; }
 
    virtual bool canReduceTimestepSize(
        NumLib::TimeStep const& timestep_previous,
        NumLib::TimeStep const& timestep_current) const override;
 
private:
    const double _kP = 0.075;  
    const double _kI = 0.175;  
    const double _kD = 0.01;   
 
    const double _h0;     
    const double _h_min;  
    const double _h_max;  

    const double _rel_h_min;
    const double _rel_h_max;
 
    const double _tol;
 
    double _e_n_minus1;  
    double _e_n_minus2;  
 
    std::vector<double> const _fixed_times_for_output;

    double limitStepSize(const double h_new,
                         const bool previous_step_accepted,
                         NumLib::TimeStep const& timestep_current) const;
};
 
}  // end of namespace NumLib