123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301 |
- /**
- * @file pdevs/tests.cpp
- * @author The ARTIS Development Team
- * See the AUTHORS or Authors.txt file
- */
- /*
- * ARTIS - the multimodeling and simulation environment
- * This file is a part of the ARTIS environment
- *
- * Copyright (C) 2013-2019 ULCO http://www.univ-littoral.fr
- *
- * This program is free software: you can redistribute it and/or modify
- * it under the terms of the GNU General Public License as published by
- * the Free Software Foundation, either version 3 of the License, or
- * (at your option) any later version.
- *
- * This program is distributed in the hope that it will be useful,
- * but WITHOUT ANY WARRANTY; without even the implied warranty of
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
- * GNU General Public License for more details.
- *
- * You should have received a copy of the GNU General Public License
- * along with this program. If not, see <http://www.gnu.org/licenses/>.
- */
- #include <tests/qss/graph_manager.hpp>
- #include <artis-star/common/RootCoordinator.hpp>
- #include <artis-star/common/observer/Iterator.hpp>
- #include <artis-star/common/observer/Output.hpp>
- #include <fstream>
- #include <iostream>
- #include <boost/archive/binary_oarchive.hpp>
- using namespace artis::tests::qss;
- void test_parabola()
- {
- artis::qss::QSSParameters<ParabolaParameters> parameters = {{0.5},
- {true, true, 0.001, 3},
- {0.2}};
- class View : public artis::observer::View<artis::common::DoubleTime>
- {
- public:
- View()
- {
- selector("Value",
- {OnlyOneParabolaGraphManager::A, ParabolaGraphManager::submodel::S_Integrator,
- artis::qss::Integrator<artis::common::DoubleTime>::var::VALUE});
- }
- };
- artis::common::context::Context<artis::common::DoubleTime> context(0, 5);
- artis::common::RootCoordinator<
- artis::common::DoubleTime, artis::pdevs::Coordinator<
- artis::common::DoubleTime,
- OnlyOneParabolaGraphManager,
- artis::qss::QSSParameters<ParabolaParameters>>
- > rc(context, "root", parameters, artis::common::NoParameters());
- rc.attachView("Value", new View());
- rc.run(context);
- const View::Values &values = rc.observer().view("Value").get("Value");
- for (const auto &value: values) {
- double v;
- value.second(v);
- std::cout << value.first << ": " << v << std::endl;
- }
- }
- class PredatorPreyView : public artis::observer::View<artis::common::DoubleTime>
- {
- public:
- PredatorPreyView()
- {
- selector("PredatorView",
- {PreyPredatorGraphManager::PREDATOR, PredatorGraphManager::submodel::S_Integrator,
- artis::qss::Integrator<artis::common::DoubleTime>::var::VALUE});
- selector("PreyView",
- {PreyPredatorGraphManager::PREY, PreyGraphManager::submodel::S_Integrator,
- artis::qss::Integrator<artis::common::DoubleTime>::var::VALUE});
- }
- };
- void run_predator_prey(artis::common::context::Context<artis::common::DoubleTime> &context,
- const PreyPredatorGraphManagerParameters ¶meters)
- {
- artis::common::RootCoordinator<
- artis::common::DoubleTime, artis::pdevs::Coordinator<
- artis::common::DoubleTime,
- PreyPredatorGraphManager,
- PreyPredatorGraphManagerParameters>
- > rc(context, "root", parameters, artis::common::NoParameters());
- rc.attachView("Value1", new PredatorPreyView());
- rc.run(context);
- artis::observer::Output<artis::common::DoubleTime,
- artis::observer::TimedIterator<artis::common::DoubleTime>>
- output(rc.observer());
- output(context.begin(), context.end(), {context.begin(), 0.1});
- // rc.save(context);
- //
- // std::ofstream os("state");
- // boost::archive::binary_oarchive oa(os);
- //
- // oa << context;
- //
- // {
- // artis::observer::DiscreteTimeIterator<artis::common::DoubleTime> it(
- // rc.observer().view("Value").get("PredatorView"), context.begin(), 0.1);
- //
- // while (it.has_next()) {
- // double v;
- //
- // (*it).second(v);
- // std::cout << (*it).first << ";" << v << std::endl;
- // ++it;
- // }
- // }
- // {
- // artis::observer::DiscreteTimeIterator<artis::common::DoubleTime> it(
- // rc.observer().view("Value").get("PreyView"), context.begin(), 0.1);
- //
- // while (it.has_next()) {
- // double v;
- //
- // (*it).second(v);
- // std::cout << (*it).first << ";" << v << std::endl;
- // ++it;
- // }
- // }
- }
- void test_predator_prey()
- {
- PreyPredatorGraphManagerParameters parameters = {{{45.},
- {true, true, 0.1, 3},
- {0.5, 0.01, 0.01, 0.2}},
- {{5000.},
- {true, true, 1, 3},
- {0.5, 0.01, 0.01, 0.2}}};
- artis::common::context::Context<artis::common::DoubleTime> context(0, 100);
- run_predator_prey(context, parameters);
- // artis::common::context::Context<artis::common::DoubleTime> new_context(context);
- //
- // new_context.end(200);
- // run_predator_prey(new_context, parameters);
- }
- class PredatorPreySmartGardenerView : public artis::observer::View<artis::common::DoubleTime>
- {
- public:
- PredatorPreySmartGardenerView()
- {
- selector("PredatorView",
- {PreyPredatorSmartGardenerGraphManager::PREY_PREDATOR,
- PreyPredatorGraphManager::PREDATOR, PredatorGraphManager::submodel::S_Integrator,
- artis::qss::Integrator<artis::common::DoubleTime>::var::VALUE});
- selector("PreyView",
- {PreyPredatorSmartGardenerGraphManager::PREY_PREDATOR,
- PreyPredatorGraphManager::PREY, PreyGraphManager::submodel::S_Integrator,
- artis::qss::Integrator<artis::common::DoubleTime>::var::VALUE});
- }
- };
- void test_predator_prey_smart_gardener()
- {
- PreyPredatorSmartGardenerGraphManagerParameters parameters = {{{{45.},
- {true, true, 0.1, 3},
- {0.5, 0.01, 0.01, 0.2}},
- {{5000.},
- {true, true, 1, 3},
- {0.5, 0.01, 0.01, 0.2}}},
- {2000, 0.25, 0.75, 5}};
- artis::common::context::Context<artis::common::DoubleTime> context(0, 100);
- artis::common::RootCoordinator<
- artis::common::DoubleTime, artis::pdevs::Coordinator<
- artis::common::DoubleTime,
- PreyPredatorSmartGardenerGraphManager,
- PreyPredatorSmartGardenerGraphManagerParameters>
- > rc(context, "root", parameters, artis::common::NoParameters());
- rc.attachView("Value2", new PredatorPreySmartGardenerView());
- rc.run(context);
- artis::observer::Output<artis::common::DoubleTime,
- artis::observer::TimedIterator<artis::common::DoubleTime>>
- output(rc.observer());
- output(context.begin(), context.end(), {context.begin(), 0.1});
- }
- class MixedPredatorPreyView : public artis::observer::View<artis::common::DoubleTime>
- {
- public:
- MixedPredatorPreyView()
- {
- selector("PredatorView",
- {MixedPreyPredatorGraphManager::PREDATOR, DiscretePredatorGraphManager::PREDATOR,
- DiscretePredator::VALUE});
- selector("PreyView",
- {MixedPreyPredatorGraphManager::PREY, PreyGraphManager::submodel::S_Integrator,
- artis::qss::Integrator<artis::common::DoubleTime>::var::VALUE});
- }
- };
- void test_mixed_predator_prey()
- {
- MixedPreyPredatorGraphManagerParameters parameters = {{0.0001, 45., 0.5, 0.01, 0.01, 0.2},
- {{5000.},
- {true, true, 1, 3},
- {0.5, 0.01, 0.01, 0.2}}};
- artis::common::context::Context<artis::common::DoubleTime> context(0, 100);
- artis::common::RootCoordinator<
- artis::common::DoubleTime, artis::pdevs::Coordinator<
- artis::common::DoubleTime,
- MixedPreyPredatorGraphManager,
- MixedPreyPredatorGraphManagerParameters>
- > rc(context, "root", parameters, artis::common::NoParameters());
- rc.attachView("Value3", new MixedPredatorPreyView());
- rc.run(context);
- artis::observer::Output<artis::common::DoubleTime,
- artis::observer::TimedIterator<artis::common::DoubleTime>>
- output(rc.observer());
- output(context.begin(), context.end(), {context.begin(), 0.1});
- }
- class MultiPredatorPreyView : public artis::observer::View<artis::common::DoubleTime>
- {
- public:
- MultiPredatorPreyView()
- {
- selector("PredatorView",
- {MultiPreyPredatorGraphManager::submodel::S_Integrator,
- artis::qss::Integrator<artis::common::DoubleTime>::var::VALUE});
- selector("PreyView",
- {MultiPreyPredatorGraphManager::submodel::S_Integrator + 1,
- artis::qss::Integrator<artis::common::DoubleTime>::var::VALUE});
- selector("PredatorDerivativeView",
- {MultiPreyPredatorGraphManager::submodel::S_Derivative,
- PreyPredator::var::VALUE});
- selector("PreyDerivativeView",
- {MultiPreyPredatorGraphManager::submodel::S_Derivative,
- PreyPredator::var::VALUE + 1});
- }
- };
- void test_multi_predator_prey()
- {
- artis::qss::MultiQSSParameters<PreyPredatorParameters> parameters = {
- {{5000.}, {45.}},
- {{true, true, 1, 3}, {true, true, 0.1, 3}},
- {0.5, 0.01, 0.01, 0.2}
- };
- artis::common::context::Context<artis::common::DoubleTime> context(0, 100);
- artis::common::RootCoordinator<
- artis::common::DoubleTime, artis::pdevs::Coordinator<
- artis::common::DoubleTime,
- MultiPreyPredatorGraphManager,
- artis::qss::MultiQSSParameters<PreyPredatorParameters>>
- > rc(context, "root", parameters, artis::common::NoParameters());
- rc.attachView("Value4", new MultiPredatorPreyView());
- rc.run(context);
- artis::observer::Output<artis::common::DoubleTime,
- artis::observer::TimedIterator<artis::common::DoubleTime>>
- output(rc.observer());
- output(context.begin(), context.end(), {context.begin(), 0.1});
- }
- int main()
- {
- test_parabola();
- test_predator_prey();
- test_predator_prey_smart_gardener();
- test_mixed_predator_prey();
- test_multi_predator_prey();
- }
|