models/ide_secir/simulation.h

Go to the documentation of this file.

/* 
* Copyright (C) 2020-2026 MEmilio
*
* Authors: Martin J Kuehn, Anna Wendler, Lena Ploetzke
*
* Contact: Martin J. Kuehn <Martin.Kuehn@DLR.de>
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
*     http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
#ifndef IDE_SECIR_SIMULATION_H
#define IDE_SECIR_SIMULATION_H
 
#include "ide_secir/model.h"
#include "memilio/config.h"
#include "memilio/utils/time_series.h"
#include <memory>
#include <cstdio>
 
namespace mio
{
namespace isecir
{
 
class Simulation
{
 
public:
    Simulation(Model const& model, ScalarType dt = 0.1)
        : m_model(std::make_unique<Model>(model))
        , m_dt(dt)
    {
        assert(m_dt > 0);
        m_model->set_transitiondistributions_support_max(m_dt);
        m_model->set_transitiondistributions_derivative(m_dt);
        m_model->set_transitiondistributions_in_forceofinfection(m_dt);
    }
 
    void advance(ScalarType tmax);
 
    TimeSeries<ScalarType> get_result()
    {
        return m_model->populations;
    }
 
    const TimeSeries<ScalarType>& get_result() const
    {
        return m_model->populations;
    }
 
    TimeSeries<ScalarType> const& get_transitions()
    {
        return m_model->transitions;
    }
 
    const Model& get_model() const
    {
        return *m_model;
    }
 
    Model& get_model()
    {
        return *m_model;
    }
 
    ScalarType get_dt()
    {
        return m_dt;
    }
 
private:
    std::unique_ptr<Model> m_model; 
    ScalarType m_dt; 
};
 
TimeSeries<ScalarType> simulate(ScalarType tmax, ScalarType dt, Model const& model);
 
} // namespace isecir
} // namespace mio
 
#endif //IDE_SECIR_SIMULATION_H