123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586 |
- /**
- * @file tests/pdevs/models.hpp
- * @author The PARADEVS Development Team
- * See the AUTHORS or Authors.txt file
- */
- /*
- * PARADEVS - the multimodeling and simulation environment
- * This file is a part of the PARADEVS environment
- *
- * Copyright (C) 2013-2015 ULCO http://www.univ-litoral.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/>.
- */
- #ifndef TESTS_PDEVS_MODELS_HPP
- #define TESTS_PDEVS_MODELS_HPP 1
- #include <paradevs/common/time/DoubleTime.hpp>
- #include <paradevs/common/utils/Trace.hpp>
- #include <paradevs/kernel/pdevs/Dynamics.hpp>
- #include <chrono>
- #define DELAY 100
- namespace paradevs { namespace tests { namespace pdevs {
- void delay()
- {
- for (unsigned int i = 0; i < DELAY; ++i) {
- std::vector < int > v;
- for (unsigned int j = 1000; j > 0; --j) {
- v.push_back(j);
- }
- std::sort(v.begin(), v.end());
- }
- }
- struct data
- {
- double x;
- double y;
- data() : x(0), y(0)
- { }
- data(double _x, double _y) : x(_x), y(_y)
- { }
- };
- class A :
- public paradevs::pdevs::Dynamics < common::DoubleTime >
- {
- public:
- A(const std::string& name, const common::NoParameters& parameters) :
- paradevs::pdevs::Dynamics < common::DoubleTime >(name, parameters)
- { }
- virtual ~A()
- { }
- void dint(typename common::DoubleTime::type t)
- {
- #ifndef WITH_TRACE
- (void)t;
- #endif
- #ifdef WITH_TRACE
- common::Trace < common::DoubleTime >::trace()
- << common::TraceElement < common::DoubleTime >(
- get_name(), t,
- common::DELTA_INT);
- common::Trace < common::DoubleTime >::trace().flush();
- #endif
- delay();
- ++_value.x;
- --_value.y;
- std::cout << t << ": " << get_name() << " => dint -> "
- << _value.x << " " << _value.y << std::endl;
- if (_phase == SEND) {
- _phase = WAIT;
- }
- }
- void dext(typename common::DoubleTime::type t,
- typename common::DoubleTime::type /* e */,
- const common::Bag < common::DoubleTime >& msgs)
- {
- #ifndef WITH_TRACE
- (void)t;
- (void)msgs;
- #endif
- for (common::Bag < common::DoubleTime >::const_iterator it =
- msgs.begin(); it != msgs.end(); ++it) {
- std::cout << t << ": " << get_name()
- << " => " << it->get_content().get_content < double >()
- << std::endl;
- }
- #ifdef WITH_TRACE
- common::Trace < common::DoubleTime >::trace()
- << common::TraceElement < common::DoubleTime >(
- get_name(), t,
- common::DELTA_EXT)
- << "messages = " << msgs.to_string();
- common::Trace < common::DoubleTime >::trace().flush();
- #endif
- _phase = SEND;
- }
- void dconf(typename common::DoubleTime::type t,
- typename common::DoubleTime::type /* e */,
- const common::Bag < common::DoubleTime >& msgs)
- {
- #ifndef WITH_TRACE
- (void)t;
- (void)msgs;
- #endif
- std::cout << t << ": " << get_name() << " => " << msgs.to_string()
- << std::endl;
- #ifdef WITH_TRACE
- common::Trace < common::DoubleTime >::trace()
- << common::TraceElement < common::DoubleTime >(
- get_name(), t,
- common::DELTA_CONF)
- << "messages = " << msgs.to_string();
- common::Trace < common::DoubleTime >::trace().flush();
- #endif
- }
- typename common::DoubleTime::type start(typename common::DoubleTime::type t)
- {
- #ifndef WITH_TRACE
- (void)t;
- #endif
- #ifdef WITH_TRACE
- common::Trace < common::DoubleTime >::trace()
- << common::TraceElement < common::DoubleTime >(
- get_name(), t,
- common::START);
- common::Trace < common::DoubleTime >::trace().flush();
- #endif
- _phase = WAIT;
- return 0;
- }
- typename common::DoubleTime::type ta(
- typename common::DoubleTime::type t) const
- {
- #ifndef WITH_TRACE
- (void)t;
- #endif
- #ifdef WITH_TRACE
- common::Trace < common::DoubleTime >::trace()
- << common::TraceElement < common::DoubleTime >(
- get_name(), t,
- common::TA);
- common::Trace < common::DoubleTime >::trace().flush();
- #endif
- if (_phase == WAIT) {
- return 1;
- } else {
- return 0;
- }
- }
- common::Bag < common::DoubleTime > lambda(
- typename common::DoubleTime::type t) const
- {
- #ifndef WITH_TRACE
- (void)t;
- #endif
- common::Bag < common::DoubleTime > msgs;
- msgs.push_back(
- common::ExternalEvent < common::DoubleTime >(
- "out", common::Value(_value)));
- std::cout << t << ": " << get_name() << " => lambda" << std::endl;
- #ifdef WITH_TRACE
- common::Trace < common::DoubleTime >::trace()
- << common::TraceElement < common::DoubleTime >(
- get_name(), t,
- common::LAMBDA)
- << "messages = " << msgs.to_string();
- common::Trace < common::DoubleTime >::trace().flush();
- #endif
- return msgs;
- }
- void observation(std::ostream& /* file */) const
- { }
- private:
- enum Phase { WAIT, SEND };
- Phase _phase;
- data _value;
- };
- class B :
- public paradevs::pdevs::Dynamics < common::DoubleTime >
- {
- public:
- B(const std::string& name, const common::NoParameters& parameters) :
- paradevs::pdevs::Dynamics < common::DoubleTime >(name, parameters),
- _value(0)
- { }
- virtual ~B()
- { }
- void dint(typename common::DoubleTime::type t)
- {
- #ifndef WITH_TRACE
- (void)t;
- #endif
- #ifdef WITH_TRACE
- common::Trace < common::DoubleTime >::trace()
- << common::TraceElement < common::DoubleTime >(
- get_name(), t,
- common::DELTA_INT);
- common::Trace < common::DoubleTime >::trace().flush();
- #endif
- delay();
- ++_value;
- std::cout << t << ": " << get_name() << " => dint -> "
- << _value << std::endl;
- if (_phase == SEND) {
- _phase = WAIT;
- }
- }
- void dext(typename common::DoubleTime::type t,
- typename common::DoubleTime::type /* e */,
- const common::Bag < common::DoubleTime >& msgs)
- {
- #ifndef WITH_TRACE
- (void)t;
- (void)msgs;
- #endif
- for (common::Bag < common::DoubleTime >::const_iterator it =
- msgs.begin(); it != msgs.end(); ++it) {
- std::cout << t << ": " << get_name()
- << " => " << it->get_content().get_content < data >().x
- << " " << it->get_content().get_content < data >().y
- << std::endl;
- }
- #ifdef WITH_TRACE
- common::Trace < common::DoubleTime >::trace()
- << common::TraceElement < common::DoubleTime >(
- get_name(), t,
- common::DELTA_EXT)
- << "messages = " << msgs.to_string();
- common::Trace < common::DoubleTime >::trace().flush();
- #endif
- _phase = SEND;
- }
- void dconf(typename common::DoubleTime::type t,
- typename common::DoubleTime::type /* e */,
- const common::Bag < common::DoubleTime >& msgs)
- {
- #ifndef WITH_TRACE
- (void)t;
- (void)msgs;
- #endif
- #ifdef WITH_TRACE
- common::Trace < common::DoubleTime >::trace()
- << common::TraceElement < common::DoubleTime >(
- get_name(), t,
- common::DELTA_CONF)
- << "messages = " << msgs.to_string();
- common::Trace < common::DoubleTime >::trace().flush();
- #endif
- }
- typename common::DoubleTime::type start(typename common::DoubleTime::type t)
- {
- #ifndef WITH_TRACE
- (void)t;
- #endif
- #ifdef WITH_TRACE
- common::Trace < common::DoubleTime >::trace()
- << common::TraceElement < common::DoubleTime >(
- get_name(), t,
- common::START);
- common::Trace < common::DoubleTime >::trace().flush();
- #endif
- _phase = WAIT;
- return std::numeric_limits < double >::max();
- }
- typename common::DoubleTime::type ta(
- typename common::DoubleTime::type t) const
- {
- #ifndef WITH_TRACE
- (void)t;
- #endif
- #ifdef WITH_TRACE
- common::Trace < common::DoubleTime >::trace()
- << common::TraceElement < common::DoubleTime >(
- get_name(), t,
- common::TA);
- common::Trace < common::DoubleTime >::trace().flush();
- #endif
- if (_phase == WAIT) {
- return std::numeric_limits < double >::max();
- } else {
- return 0;
- }
- }
- common::Bag < common::DoubleTime > lambda(
- typename common::DoubleTime::type t) const
- {
- #ifndef WITH_TRACE
- (void)t;
- #endif
- common::Bag < common::DoubleTime > msgs;
- std::cout << t << ": " << get_name() << " => lambda" << std::endl;
- msgs.push_back(common::ExternalEvent <
- common::DoubleTime >(
- "out", common::Value(_value)));
- #ifdef WITH_TRACE
- common::Trace < common::DoubleTime >::trace()
- << common::TraceElement < common::DoubleTime >(
- get_name(), t,
- common::LAMBDA)
- << "messages = " << msgs.to_string();
- common::Trace < common::DoubleTime >::trace().flush();
- #endif
- return msgs;
- }
- void observation(std::ostream& /* file */) const
- { }
- private:
- enum Phase { WAIT, SEND };
- Phase _phase;
- double _value;
- };
- class TwoStateModel :
- public paradevs::pdevs::Dynamics < common::DoubleTime >
- {
- public:
- TwoStateModel(const std::string& name,
- const common::NoParameters& parameters) :
- paradevs::pdevs::Dynamics < common::DoubleTime >(name, parameters)
- { }
- virtual ~TwoStateModel()
- { }
- void dint(typename common::DoubleTime::type t)
- {
- if (_phase == S1) {
- _phase = S2;
- } else if (_phase == S2) {
- _phase = S1;
- }
- _last_time = t;
- }
- typename common::DoubleTime::type start(
- typename common::DoubleTime::type t)
- {
- _phase = S1;
- _last_time = t;
- return ta(t);
- }
- typename common::DoubleTime::type ta(
- typename common::DoubleTime::type /* t */) const
- {
- if (_phase == S1) {
- return 5;
- } else {
- return 6;
- }
- }
- common::Bag < common::DoubleTime > lambda(
- typename common::DoubleTime::type t) const
- {
- std::cout << (t - _last_time) << std::endl;
- return common::Bag < common::DoubleTime >();
- }
- private:
- enum Phase { S1, S2 };
- Phase _phase;
- typename common::DoubleTime::type _last_time;
- };
- class ThreeStateModel :
- public paradevs::pdevs::Dynamics < common::DoubleTime >
- {
- public:
- ThreeStateModel(const std::string& name,
- const common::NoParameters& parameters) :
- paradevs::pdevs::Dynamics < common::DoubleTime >(name, parameters)
- { }
- virtual ~ThreeStateModel()
- { }
- void compute()
- {
- for (unsigned int i = 0; i < heights.size(); ++i) {
- if (heights[i] != -1 and heights[i] < 10) {
- heights[i] += speeds[i] * scales[i];
- }
- }
- }
- void display() const
- {
- for (std::vector < double >::const_iterator it = heights.begin();
- it != heights.end(); ++it) {
- std::cout << *it << " ";
- }
- std::cout << std::endl;
- }
- void display_full() const
- {
- unsigned int i = 1;
- for (std::vector < double >::const_iterator it = heights.begin();
- it != heights.end(); ++it, ++i) {
- if (*it > 10) {
- std::cout << "S" << i;
- }
- }
- std::cout << std::endl;
- }
- bool full() const
- {
- unsigned int n = 0;
- for (std::vector < double >::const_iterator it = heights.begin();
- it != heights.end(); ++it) {
- if (*it > 10) {
- ++n;
- }
- }
- return n > 0;
- }
- bool full_N() const
- {
- unsigned int n = 0;
- for (std::vector < double >::const_iterator it = heights.begin();
- it != heights.end(); ++it) {
- if (*it == -1) {
- ++n;
- }
- }
- return n >= 2;
- }
- void mark_full(typename common::DoubleTime::type t)
- {
- for (std::vector < double >::iterator it = heights.begin();
- it != heights.end(); ++it) {
- if (*it > 10) {
- *it = -1;
- _last_time = t;
- }
- }
- }
- void raz()
- {
- for (std::vector < double >::iterator it = heights.begin();
- it != heights.end(); ++it) {
- if (*it == -1) {
- *it = 0;
- }
- }
- }
- void dint(typename common::DoubleTime::type t)
- {
- mark_full(t);
- if (full_N()) {
- raz();
- }
- compute();
- }
- typename common::DoubleTime::type start(
- typename common::DoubleTime::type t)
- {
- heights = { 0, 0, 0, 0, 0 };
- speeds = { 0.21, 0.3, 0.7, 0.56, 0.14 };
- scales = { 1, 1, 1, 1, 1 };
- _last_time = t;
- return 0;
- }
- typename common::DoubleTime::type ta(
- typename common::DoubleTime::type /* t */) const
- { return 1; }
- common::Bag < common::DoubleTime > lambda(
- typename common::DoubleTime::type t) const
- {
- if (full()) {
- std::cout << (t - _last_time) << " ";
- // display_full();
- }
- return common::Bag < common::DoubleTime >();
- }
- private:
- std::vector < double > heights;
- std::vector < double > speeds;
- std::vector < double > scales;
- typename common::DoubleTime::type _last_time;
- };
- } } } // namespace paradevs tests pdevs
- #endif
|