simu.cpp

Go to the documentation of this file.

#include <mecacell/mecacell.h>
#include "scenarios/Scenario.hpp"
#include <fstream>
#include <numeric>
 
using namespace std;
using scenario_t = Scenario;
 
 
class Simu {
 
    scenario_t scenario;
    MecaCell::Config::random_engine_t local_random_engine;
    int seed = -1;
    void setStaticVariables() {
        MecaCell::Config::setCurrentEngine(&local_random_engine);
    }
public:
 
    CellPlugin<PrimoCell<CellBody>>* plugins;
 
    Simu(string params) : scenario() {
        setStaticVariables();
        seed = MecaCell::Config::globalRand().default_seed;
        nlohmann::json configuration;
        configuration = nlohmann::json::parse(params);
        scenario.init(configuration);
        plugins = &(scenario.getWorld().cellPlugin);
    }
    Simu(string params, int _seed, bool benchmark) : scenario() {
        setStaticVariables();
        seed = _seed;
        Rand::initRandomSeed(seed);
        nlohmann::json configuration;
        configuration = nlohmann::json::parse(params);
        scenario.getWorld().benchmark = benchmark;
        scenario.init(configuration);
        plugins = &(scenario.getWorld().cellPlugin);
    }
 
    void update() {
        setStaticVariables();
        scenario.loop();
    }
 
    void update(int nbSteps) {
        setStaticVariables();
        for (int i = 0; i < nbSteps; i++)
            scenario.loop();
    }
 
    int getCurrentStep() {
        return scenario.getWorld().getNbUpdates();
    }
 
    double getTime(char format) {
        if (format == 'd') {
            return (double)scenario.getWorld().getNbUpdates() * scenario.getWorld().getDtInH() / 24.0;
        } else if (format == 'h') {
            return (double)scenario.getWorld().getNbUpdates() * scenario.getWorld().getDtInH();
        } else if (format == 'm') {
            return (double)scenario.getWorld().getNbUpdates() * scenario.getWorld().getDtInMn();
        } else if (format == 's') {
            return (double)scenario.getWorld().getNbUpdates() * scenario.getWorld().getDtInS();
        } else {
            throw std::invalid_argument("received invalid format ['d': day, 'h': hours, 'm': minutes, 's': seconds]");
        }
    }
 
 
    int getSeed() const {
        return seed;
    }
 
    vector<PrimoCell<CellBody>*>* _getCells() {
        setStaticVariables();
        return &(scenario.getWorld().cells);
    }
    void showBenchmarksSummaries() const {
        scenario.showBenchmarksSummaries();
    }
 
    string _getBenchmarksJSON() const {
        return scenario.getBenchmarksJSON();
    }
 
    string __repr__() {
        std::ostringstream os;
        os << "<Simu step=" << getCurrentStep() << " nbCells=" << _getCells()->size() << " seed=" << getSeed() << " benchmark=" << scenario.getWorld().benchmark << ">";
        return os.str();
    }
 
    /*isigen:eval@SimuUtilityDeclarations'''
        from string import Template
        cellAttributeGettersTemplate = Template("""
        vector<$type> get$nameUpper(){
        \tvector<$type> ret;
        \tret.resize(scenario.getWorld().cells.size());
        \ttransform(scenario.getWorld().cells.begin(), scenario.getWorld().cells.end(), ret.begin(), [](auto* c){return $name;});
        \treturn ret;
        }
        """)
        attributesToIgnore = ['type', 'state']
        validType = ['bool','double','int','void']
        validType = [*validType,*[ 'vector<'+vt+'>' for vt in validType]]
        validType = [*validType,*[ 'vector<'+vt+'>' for vt in validType]]
        def checkValidType(type):
            return type.replace(' ','') in validType
        cellAttributeGetters = [a for a in builderData["cell"]["attributes"] if checkValidType(a['type']) and a['value'] not in attributesToIgnore]
        if "enums" in builderData:
            cellAttributeGettersEnum = [a for a in builderData["cell"]["attributes"] if a['type'] in builderData["enums"] and a['value'] not in attributesToIgnore]
        '''*/
 
    /*isigen:insert@functionsScenarioWrapper'''
        import re
        scenarioFunWrapper = ''
        for fun in builderData["scenario"]["functions"].values():
            res = re.search(r"(?P<ret>.+) +(?P<fun>.+) *\‍((?P<params>.*)\‍)",fun)
            if checkValidType(res.group('ret')):
                scenarioFunWrapper+="{ret} {funName}({params}){{\n\t{content}\n}}\n".format(ret=res.group('ret'),funName=res.group('fun'),params=res.group('params'),content=('return 'if res.group('ret')!='void' else '')+'scenario.'+res.group('fun')+'('+ ','.join([ s.split(' ')[-1] for s in res.group('params').split(',')])+');')
        print(scenarioFunWrapper)
        '''*/
};