integrating_controller.h

Go to the documentation of this file.

#ifndef WXM_INTEGRATOR_CONTROLLER_H
#define WXM_INTEGRATOR_CONTROLLER_H
 
#include "hostactions/time_step/timestep_decision.h"
#include "time_stepper_base.h"
 
#include "hostactions/wmhostaction.h"
 
#include "dfem/tools/integrator.h"
 
#include <memory>
#include <vector>
#include <array>
#include <numeric>
 
namespace wxm
{
namespace dt_calc
{
class integrating_controller : public time_stepper_base
{
    // holder for child subdomain integrators
    std::vector<std::unique_ptr<WmHostAction>> _processors;
 
    // integrators to run to compute the error estimate
    std::vector<wxm::dfem::dg_integrator*> _integrators;
 
    // pair of variable/component
    typedef std::vector<std::pair<wxm::variable::distributed_variable_t*, size_t>>
        variables_type;
 
    // pair of distributed variable and component index
    std::vector<variables_type> _variables;
 
    bool add_variables(const WxCryptSet& wxc, const std::string& name);
 
public:
    double max_dt;
    double max_growth;
 
    double kappa;
 
    double error_tol;
 
    double max_error;
 
    double alpha;
 
    bool stability_limit;
 
    integrating_controller(const WmConstOpaquePtr& parent_solver)
        : time_stepper_base(parent_solver)
    {
    }
 
    void setup(const WxCryptSet& wxc) override;
 
    bool allow_retry() const override
    {
        return true;
    }
 
    TimestepDecision test_dt(const WxStepperStatus& stab) override;
};
} // namespace dt_calc
} // namespace wxm
 
#endif