src/kernel/EngineImpl.cpp

   1 /* Copyright (c) 2016-2023. The SimGrid Team. All rights reserved.          */
   2
   3 /* This program is free software; you can redistribute it and/or modify it
   4  * under the terms of the license (GNU LGPL) which comes with this package. */
   5
   6 #include <simgrid/Exception.hpp>
   7 #include <simgrid/kernel/Timer.hpp>
   8 #include <simgrid/kernel/routing/NetPoint.hpp>
   9 #include <simgrid/kernel/routing/NetZoneImpl.hpp>
  10 #include <simgrid/s4u/Host.hpp>
  11
  12 #include "src/kernel/EngineImpl.hpp"
  13 #include "src/kernel/resource/StandardLinkImpl.hpp"
  14 #include "src/kernel/resource/profile/Profile.hpp"
  15 #include "src/kernel/xml/platf.hpp"
  16 #include "src/mc/mc.h"
  17 #include "src/mc/mc_record.hpp"
  18 #include "src/mc/mc_replay.hpp"
  19 #include "src/simgrid/math_utils.h"
  20 #include "src/simgrid/sg_config.hpp"
  21 #include "src/smpi/include/smpi_actor.hpp"
  22 #include "src/xbt/xbt_modinter.h" /* whether initialization was already done */
  23
  24 #include "xbt/log.hpp"
  25
  26 #include <boost/algorithm/string/predicate.hpp>
  27 #include <dlfcn.h>
  28
  29 #if SIMGRID_HAVE_MC
  30 #include "src/mc/remote/AppSide.hpp"
  31 #endif
  32
  33 XBT_LOG_NEW_DEFAULT_CATEGORY(ker_engine, "Logging specific to Engine (kernel)");
  34
  35 namespace simgrid::kernel {
  36 double EngineImpl::now_           = 0.0;
  37 EngineImpl* EngineImpl::instance_ = nullptr; /* That singleton is awful too. */
  38
  39 config::Flag<double> cfg_breakpoint{"debug/breakpoint",
  40                                     "When non-negative, raise a SIGTRAP after given (simulated) time", -1.0};
  41 config::Flag<bool> cfg_verbose_exit{"debug/verbose-exit", "Display the actor status at exit", true};
  42
  43 constexpr std::initializer_list<std::pair<const char*, context::ContextFactory* (*)()>> context_factories = {
  44 #if HAVE_RAW_CONTEXTS
  45     {"raw", &context::raw_factory},
  46 #endif
  47 #if HAVE_UCONTEXT_CONTEXTS
  48     {"ucontext", &context::sysv_factory},
  49 #endif
  50 #if HAVE_BOOST_CONTEXTS
  51     {"boost", &context::boost_factory},
  52 #endif
  53     {"thread", &context::thread_factory},
  54 };
  55
  56 static_assert(context_factories.size() > 0, "No context factories are enabled for this build");
  57
  58 // Create the list of possible contexts:
  59 static inline std::string contexts_list()
  60 {
  61   std::string res;
  62   std::string sep = "";
  63   for (auto const& [factory_name, _] : context_factories) {
  64     res += sep + factory_name;
  65     sep = ", ";
  66   }
  67   return res;
  68 }
  69
  70 static config::Flag<std::string> context_factory_name("contexts/factory",
  71                                                       ("Possible values: " + contexts_list()).c_str(),
  72                                                       context_factories.begin()->first);
  73
  74 } // namespace simgrid::kernel
  75
  76 XBT_ATTRIB_NORETURN static void inthandler(int)
  77 {
  78   if (simgrid::kernel::cfg_verbose_exit) {
  79     XBT_INFO("CTRL-C pressed. The current status will be displayed before exit (disable that behavior with option "
  80              "'debug/verbose-exit').");
  81     simgrid::kernel::EngineImpl::get_instance()->display_all_actor_status();
  82   } else {
  83     XBT_INFO("CTRL-C pressed, exiting. Hiding the current process status since 'debug/verbose-exit' is set to false.");
  84   }
  85   exit(1);
  86 }
  87
  88 static void segvhandler(int signum, siginfo_t* siginfo, void* /*context*/)
  89 {
  90   if ((siginfo->si_signo == SIGSEGV && siginfo->si_code == SEGV_ACCERR) || siginfo->si_signo == SIGBUS) {
  91     fprintf(stderr,
  92             "Access violation or Bus error detected.\n"
  93             "This probably comes from a programming error in your code, or from a stack\n"
  94             "overflow. If you are certain of your code, try increasing the stack size\n"
  95             "   --cfg=contexts/stack-size:XXX (current size is %u KiB).\n"
  96             "\n"
  97             "If it does not help, this may have one of the following causes:\n"
  98             "a bug in SimGrid, a bug in the OS or a bug in a third-party libraries.\n"
  99             "Failing hardware can sometimes generate such errors too.\n"
 100             "\n"
 101             "If you think you've found a bug in SimGrid, please report it along with a\n"
 102             "Minimal Working Example (MWE) reproducing your problem and a full backtrace\n"
 103             "of the fault captured with gdb or valgrind.\n",
 104             simgrid::kernel::context::Context::stack_size / 1024);
 105   } else if (siginfo->si_signo == SIGSEGV) {
 106     fprintf(stderr, "Segmentation fault.\n");
 107 #if HAVE_SMPI
 108     if (SMPI_is_inited() && smpi_cfg_privatization() == SmpiPrivStrategies::NONE) {
 109 #if HAVE_PRIVATIZATION
 110       fprintf(stderr, "Try to enable SMPI variable privatization with --cfg=smpi/privatization:yes.\n");
 111 #else
 112       fprintf(stderr, "Sadly, your system does not support --cfg=smpi/privatization:yes (yet).\n");
 113 #endif /* HAVE_PRIVATIZATION */
 114     }
 115 #endif /* HAVE_SMPI */
 116   }
 117   std::raise(signum);
 118 }
 119
 120 static void install_signal_handlers()
 121 {
 122   /* Install signal handler for SIGINT */
 123   std::signal(SIGINT, inthandler);
 124
 125   /* Install signal handler for SIGSEGV */
 126   if (simgrid::kernel::context::Context::install_sigsegv_stack(true) == -1) {
 127     XBT_WARN("Failed to register alternate signal stack: %s", strerror(errno));
 128     return;
 129   }
 130
 131   struct sigaction action;
 132   action.sa_sigaction = &segvhandler;
 133   action.sa_flags     = SA_ONSTACK | SA_RESETHAND | SA_SIGINFO;
 134   sigemptyset(&action.sa_mask);
 135
 136   /* Linux tend to raise only SIGSEGV where other systems also raise SIGBUS on severe error */
 137   for (int sig : {SIGSEGV, SIGBUS}) {
 138     if (sigaction(sig, &action, nullptr) == -1)
 139       XBT_WARN("Failed to register signal handler for signal %d: %s", sig, strerror(errno));
 140   }
 141 }
 142
 143 static simgrid::config::Flag<bool> cfg_dbg_clean_atexit{
 144     "debug/clean-atexit", "Whether to cleanup SimGrid at exit. Disable it if your code segfaults after its end.", true};
 145 static void xbt_postexit()
 146 {
 147   if (not cfg_dbg_clean_atexit)
 148     return;
 149   xbt_log_postexit();
 150 }
 151
 152 namespace simgrid::kernel {
 153
 154 EngineImpl::~EngineImpl()
 155 {
 156   /* Also delete the other data */
 157   delete netzone_root_;
 158   for (auto const& [_, netpoint] : netpoints_)
 159     delete netpoint;
 160
 161   for (auto const& [_, mailbox] : mailboxes_)
 162     delete mailbox;
 163
 164   /* Kill all actors (but maestro) */
 165   maestro_->kill_all();
 166   run_all_actors();
 167   empty_trash();
 168
 169   delete maestro_;
 170   delete context_factory_;
 171
 172   /* clear models before freeing handle, network models can use external callback defined in the handle */
 173   models_prio_.clear();
 174 }
 175
 176 void EngineImpl::initialize(int* argc, char** argv)
 177 {
 178   xbt_assert(EngineImpl::instance_ == nullptr,
 179              "It is currently forbidden to create more than one instance of kernel::EngineImpl");
 180   EngineImpl::instance_ = this;
 181 #if SIMGRID_HAVE_MC
 182   // The communication initialization is done ASAP, as we need to get some init parameters from the MC for different
 183   // layers. But instance_ needs to be created, as we send the address of some of its fields to the MC that wants to
 184   // read them directly.
 185   simgrid::mc::AppSide::initialize();
 186 #endif
 187
 188   static bool inited = false;
 189   if (not inited) {
 190     inited = true;
 191     atexit(xbt_postexit);
 192     xbt_log_init(argc, argv);
 193
 194     simgrid::xbt::install_exception_handler();
 195
 196     for (int i = 0; i < *argc; i++)
 197       cmdline_.emplace_back(argv[i]);
 198
 199     sg_config_init(argc, argv);
 200   }
 201
 202   instance_->context_mod_init();
 203
 204   install_signal_handlers();
 205
 206   /* register a function to be called after the environment creation */
 207   s4u::Engine::on_platform_created_cb([this]() { this->presolve(); });
 208
 209   if (cfg_dbg_clean_atexit)
 210     atexit(shutdown);
 211 }
 212
 213 void EngineImpl::context_mod_init() const
 214 {
 215   xbt_assert(not instance_->has_context_factory());
 216
 217 #if HAVE_SMPI && defined(__NetBSD__)
 218   smpi_init_options_internal(false);
 219   std::string priv = config::get_value<std::string>("smpi/privatization");
 220   if (context_factory_name == "thread" && (priv == "dlopen" || priv == "yes" || priv == "default" || priv == "1")) {
 221     XBT_WARN("dlopen+thread broken on Apple and BSD. Switching to raw contexts.");
 222     context_factory_name = "raw";
 223   }
 224 #endif
 225
 226 #if HAVE_SMPI && defined(__FreeBSD__)
 227   smpi_init_options_internal(false);
 228   if (context_factory_name == "thread" && config::get_value<std::string>("smpi/privatization") != "no") {
 229     XBT_WARN("mmap broken on FreeBSD, but dlopen+thread broken too. Switching to dlopen+raw contexts.");
 230     context_factory_name = "raw";
 231   }
 232 #endif
 233
 234   /* use the factory specified by --cfg=contexts/factory:value */
 235   for (auto const& [factory_name, factory] : context_factories)
 236     if (context_factory_name == factory_name) {
 237       instance_->set_context_factory(factory());
 238       break;
 239     }
 240
 241   if (not instance_->has_context_factory()) {
 242     XBT_ERROR("Invalid context factory specified. Valid factories on this machine:");
 243 #if HAVE_RAW_CONTEXTS
 244     XBT_ERROR("  raw: high performance context factory implemented specifically for SimGrid");
 245 #else
 246     XBT_ERROR("  (raw contexts were disabled at compilation time on this machine -- check configure logs for details)");
 247 #endif
 248 #if HAVE_UCONTEXT_CONTEXTS
 249     XBT_ERROR("  ucontext: classical system V contexts (implemented with makecontext, swapcontext and friends)");
 250 #else
 251     XBT_ERROR("  (ucontext was disabled at compilation time on this machine -- check configure logs for details)");
 252 #endif
 253 #if HAVE_BOOST_CONTEXTS
 254     XBT_ERROR("  boost: this uses the boost libraries context implementation");
 255 #else
 256     XBT_ERROR("  (boost was disabled at compilation time on this machine -- check configure logs for details. Did you "
 257               "install the libboost-context-dev package?)");
 258 #endif
 259     XBT_ERROR("  thread: slow portability layer using standard threads as provided by libstdc");
 260     xbt_die("Please use a valid factory.");
 261   }
 262 }
 263
 264 void EngineImpl::shutdown()
 265 {
 266   if (EngineImpl::instance_ == nullptr)
 267     return;
 268   XBT_DEBUG("EngineImpl::shutdown() called. Simulation's over.");
 269 #if HAVE_SMPI
 270   if (not instance_->actor_list_.empty()) {
 271     if (smpi_process() && smpi_process()->initialized()) {
 272       xbt_die("Process exited without calling MPI_Finalize - Killing simulation");
 273     } else {
 274       XBT_WARN("Process called exit when leaving - Skipping cleanups");
 275       return;
 276     }
 277   }
 278 #endif
 279
 280   if (instance_->has_actors_to_run() && simgrid_get_clock() <= 0.0) {
 281     XBT_CRITICAL("   ");
 282     XBT_CRITICAL("The time is still 0, and you still have processes ready to run.");
 283     XBT_CRITICAL("It seems that you forgot to run the simulation that you setup.");
 284     xbt_die("Bailing out to avoid that stop-before-start madness. Please fix your code.");
 285   }
 286
 287   while (not timer::kernel_timers().empty()) {
 288     delete timer::kernel_timers().top().second;
 289     timer::kernel_timers().pop();
 290   }
 291
 292   tmgr_finalize();
 293   sg_platf_parser_finalize();
 294
 295   delete instance_;
 296   instance_ = nullptr;
 297 }
 298
 299 void EngineImpl::seal_platform() const
 300 {
 301   /* Seal only once */
 302   static bool sealed = false;
 303   if (sealed)
 304     return;
 305   sealed = true;
 306
 307   /* seal netzone root, recursively seal children netzones, hosts and disks */
 308   netzone_root_->seal();
 309 }
 310
 311 void EngineImpl::load_platform(const std::string& platf)
 312 {
 313   double start = xbt_os_time();
 314   if (boost::algorithm::ends_with(platf, ".so") || boost::algorithm::ends_with(platf, ".dylib")) {
 315     void* handle = dlopen(platf.c_str(), RTLD_LAZY);
 316     xbt_assert(handle, "Impossible to open platform file: %s", platf.c_str());
 317     platf_handle_           = std::unique_ptr<void, std::function<int(void*)>>(handle, dlclose);
 318     using load_fct_t = void (*)(const simgrid::s4u::Engine&);
 319     auto callable           = (load_fct_t)dlsym(platf_handle_.get(), "load_platform");
 320     const char* dlsym_error = dlerror();
 321     xbt_assert(not dlsym_error, "Error: %s", dlsym_error);
 322     callable(*simgrid::s4u::Engine::get_instance());
 323   } else {
 324     parse_platform_file(platf);
 325   }
 326
 327   double end = xbt_os_time();
 328   XBT_DEBUG("PARSE TIME: %g", (end - start));
 329 }
 330
 331 void EngineImpl::load_deployment(const std::string& file) const
 332 {
 333   sg_platf_parser_finalize();
 334
 335   simgrid_parse_open(file);
 336   simgrid_parse();
 337   simgrid_parse_close();
 338 }
 339
 340 void EngineImpl::register_function(const std::string& name, const actor::ActorCodeFactory& code)
 341 {
 342   registered_functions[name] = code;
 343 }
 344 void EngineImpl::register_default(const actor::ActorCodeFactory& code)
 345 {
 346   default_function = code;
 347 }
 348
 349 void EngineImpl::add_model(std::shared_ptr<resource::Model> model, const std::vector<resource::Model*>& dependencies)
 350 {
 351   auto model_name = model->get_name();
 352   xbt_assert(models_prio_.find(model_name) == models_prio_.end(),
 353              "Model %s already exists, use model.set_name() to change its name", model_name.c_str());
 354
 355   for (const auto* dep : dependencies) {
 356     xbt_assert(models_prio_.find(dep->get_name()) != models_prio_.end(),
 357                "Model %s doesn't exists. Impossible to use it as dependency.", dep->get_name().c_str());
 358   }
 359   models_.push_back(model.get());
 360   models_prio_[model_name] = std::move(model);
 361 }
 362
 363 /** Wake up all actors waiting for an action to finish */
 364 void EngineImpl::handle_ended_actions() const
 365 {
 366   for (auto const& model : models_) {
 367     XBT_DEBUG("Handling the failed actions (if any)");
 368     while (auto* action = model->extract_failed_action()) {
 369       XBT_DEBUG("   Handling Action %p", action);
 370       if (action->get_activity() != nullptr) { // Skip vcpu actions
 371         // Action failures are not automatically reported when the action is started by maestro (as in SimDAG)
 372         if (action->get_activity()->get_actor() == maestro_)
 373           action->get_activity()->get_iface()->complete(s4u::Activity::State::FAILED);
 374
 375         activity::ActivityImplPtr(action->get_activity())->post();
 376       }
 377     }
 378     XBT_DEBUG("Handling the terminated actions (if any)");
 379     while (auto* action = model->extract_done_action()) {
 380       XBT_DEBUG("   Handling Action %p", action);
 381       if (action->get_activity() != nullptr) {
 382         // Action termination are not automatically reported when the action is started by maestro (as in SimDAG)
 383         action->get_activity()->set_finish_time(action->get_finish_time());
 384
 385         if (action->get_activity()->get_actor() == maestro_)
 386           action->get_activity()->get_iface()->complete(s4u::Activity::State::FINISHED);
 387
 388         activity::ActivityImplPtr(action->get_activity())->post();
 389       }
 390     }
 391   }
 392 }
 393 /**
 394  * @brief Executes the actors in actors_to_run.
 395  *
 396  * The actors in actors_to_run are run (in parallel if possible). On exit, actors_to_run is empty, and actors_that_ran
 397  * contains the list of actors that just ran.  The two lists are swapped so, be careful when using them before and after
 398  * a call to this function.
 399  */
 400 void EngineImpl::run_all_actors()
 401 {
 402   instance_->get_context_factory()->run_all(actors_to_run_);
 403
 404   for (auto const& actor : actors_to_run_)
 405     if (actor->to_be_freed())
 406       actor->cleanup_from_kernel();
 407
 408   actors_to_run_.swap(actors_that_ran_);
 409   actors_to_run_.clear();
 410 }
 411
 412 actor::ActorImpl* EngineImpl::get_actor_by_pid(aid_t pid)
 413 {
 414   auto item = actor_list_.find(pid);
 415   return item == actor_list_.end() ? nullptr : item->second;
 416 }
 417
 418 void EngineImpl::remove_daemon(actor::ActorImpl* actor)
 419 {
 420   auto it = daemons_.find(actor);
 421   xbt_assert(it != daemons_.end(), "The dying daemon is not a daemon after all. Please report that bug.");
 422   daemons_.erase(it);
 423 }
 424
 425 void EngineImpl::add_actor_to_run_list_no_check(actor::ActorImpl* actor)
 426 {
 427   XBT_DEBUG("Inserting [%p] %s(%s) in the to_run list", actor, actor->get_cname(), actor->get_host()->get_cname());
 428   actors_to_run_.push_back(actor);
 429 }
 430
 431 void EngineImpl::add_actor_to_run_list(actor::ActorImpl* actor)
 432 {
 433   if (std::find(begin(actors_to_run_), end(actors_to_run_), actor) != end(actors_to_run_)) {
 434     XBT_DEBUG("Actor %s is already in the to_run list", actor->get_cname());
 435   } else {
 436     XBT_DEBUG("Inserting [%p] %s(%s) in the to_run list", actor, actor->get_cname(), actor->get_host()->get_cname());
 437     actors_to_run_.push_back(actor);
 438   }
 439 }
 440 void EngineImpl::empty_trash()
 441 {
 442   while (not actors_to_destroy_.empty()) {
 443     actor::ActorImpl* actor = &actors_to_destroy_.front();
 444     actors_to_destroy_.pop_front();
 445     XBT_DEBUG("Getting rid of %s (refcount: %d)", actor->get_cname(), actor->get_refcount());
 446     intrusive_ptr_release(actor);
 447   }
 448 }
 449
 450 void EngineImpl::display_all_actor_status() const
 451 {
 452   XBT_INFO("%zu actors are still running, waiting for something.", actor_list_.size());
 453   /*  List the actors and their state */
 454   XBT_INFO("Legend of the following listing: \"Actor <pid> (<name>@<host>): <status>\"");
 455   for (auto const& [_, actor] : actor_list_) {
 456     if (actor->waiting_synchro_) {
 457       const char* synchro_description = "unknown";
 458
 459       if (boost::dynamic_pointer_cast<kernel::activity::ExecImpl>(actor->waiting_synchro_) != nullptr)
 460         synchro_description = "execution";
 461
 462       if (boost::dynamic_pointer_cast<kernel::activity::CommImpl>(actor->waiting_synchro_) != nullptr)
 463         synchro_description = "communication";
 464
 465       if (boost::dynamic_pointer_cast<kernel::activity::SleepImpl>(actor->waiting_synchro_) != nullptr)
 466         synchro_description = "sleeping";
 467
 468       if (boost::dynamic_pointer_cast<kernel::activity::SynchroImpl>(actor->waiting_synchro_) != nullptr)
 469         synchro_description = "synchronization";
 470
 471       if (boost::dynamic_pointer_cast<kernel::activity::IoImpl>(actor->waiting_synchro_) != nullptr)
 472         synchro_description = "I/O";
 473
 474       XBT_INFO("Actor %ld (%s@%s): waiting for %s activity %#zx (%s) in state %s to finish %s", actor->get_pid(),
 475                actor->get_cname(), actor->get_host()->get_cname(), synchro_description,
 476                (xbt_log_no_loc ? (size_t)0xDEADBEEF : (size_t)actor->waiting_synchro_.get()),
 477                actor->waiting_synchro_->get_cname(), actor->waiting_synchro_->get_state_str(),
 478                (actor->simcall_.observer_ != nullptr && not xbt_log_no_loc
 479                     ? actor->simcall_.observer_->to_string().c_str()
 480                     : ""));
 481     } else {
 482       XBT_INFO("Actor %ld (%s@%s) simcall %s", actor->get_pid(), actor->get_cname(), actor->get_host()->get_cname(),
 483                (actor->simcall_.observer_ != nullptr ? actor->simcall_.observer_->to_string().c_str()
 484                                                      : actor->simcall_.get_cname()));
 485     }
 486   }
 487 }
 488
 489 void EngineImpl::presolve() const
 490 {
 491   XBT_DEBUG("Consume all trace events occurring before the starting time.");
 492   double next_event_date;
 493   while ((next_event_date = profile::future_evt_set.next_date()) != -1.0) {
 494     if (next_event_date > now_)
 495       break;
 496
 497     double value                 = -1.0;
 498     resource::Resource* resource = nullptr;
 499     while (auto* event = profile::future_evt_set.pop_leq(next_event_date, &value, &resource)) {
 500       if (value >= 0)
 501         resource->apply_event(event, value);
 502     }
 503   }
 504
 505   XBT_DEBUG("Set every models in the right state by updating them to 0.");
 506   for (auto const& model : models_)
 507     model->update_actions_state(now_, 0.0);
 508 }
 509
 510 double EngineImpl::solve(double max_date) const
 511 {
 512   double time_delta            = -1.0; /* duration */
 513   double value                 = -1.0;
 514   resource::Resource* resource = nullptr;
 515
 516   if (max_date != -1.0) {
 517     xbt_assert(max_date >= now_, "You asked to simulate up to %f, but that's in the past already", max_date);
 518
 519     time_delta = max_date - now_;
 520   }
 521
 522   XBT_DEBUG("Looking for next event in all models");
 523   for (auto model : models_) {
 524     if (not model->next_occurring_event_is_idempotent()) {
 525       continue;
 526     }
 527     double next_event = model->next_occurring_event(now_);
 528     if ((time_delta < 0.0 || next_event < time_delta) && next_event >= 0.0) {
 529       time_delta = next_event;
 530     }
 531   }
 532
 533   XBT_DEBUG("Min for resources (remember that NS3 don't update that value): %f", time_delta);
 534
 535   XBT_DEBUG("Looking for next trace event");
 536
 537   while (true) { // Handle next occurring events until none remains
 538     double next_event_date = profile::future_evt_set.next_date();
 539     XBT_DEBUG("Next TRACE event: %f", next_event_date);
 540
 541     for (auto model : models_) {
 542       /* Skip all idempotent models, they were already treated above
 543        * NS3 is the one to handled here */
 544       if (model->next_occurring_event_is_idempotent())
 545         continue;
 546
 547       if (next_event_date != -1.0) {
 548         time_delta = std::min(next_event_date - now_, time_delta);
 549       } else {
 550         time_delta = std::max(next_event_date - now_, time_delta); // Get the positive component
 551       }
 552
 553       XBT_DEBUG("Run the NS3 network at most %fs", time_delta);
 554       // run until min or next flow
 555       double model_next_action_end = model->next_occurring_event(time_delta);
 556
 557       XBT_DEBUG("Min for network : %f", model_next_action_end);
 558       if (model_next_action_end >= 0.0)
 559         time_delta = model_next_action_end;
 560     }
 561
 562     if (next_event_date < 0.0 || (next_event_date > now_ + time_delta)) {
 563       // next event may have already occurred or will after the next resource change, then bail out
 564       XBT_DEBUG("no next usable TRACE event. Stop searching for it");
 565       break;
 566     }
 567
 568     XBT_DEBUG("Updating models (min = %g, NOW = %g, next_event_date = %g)", time_delta, now_, next_event_date);
 569
 570     while (auto* event = profile::future_evt_set.pop_leq(next_event_date, &value, &resource)) {
 571       if(value<0)
 572               continue;
 573       if (resource->is_used()) {
 574         time_delta = next_event_date - now_;
 575         XBT_DEBUG("This event invalidates the next_occurring_event() computation of models. Next event set to %f",
 576                   time_delta);
 577       }
 578       // FIXME: I'm too lame to update now_ live, so I change it and restore it so that the real update works
 579       double round_start = now_;
 580       now_               = next_event_date;
 581       /* update state of the corresponding resource to the new value. Does not touch lmm.
 582          It will be modified if needed when updating actions */
 583       XBT_DEBUG("Calling update_resource_state for resource %s", resource->get_cname());
 584       resource->apply_event(event, value);
 585       now_ = round_start;
 586     }
 587   }
 588
 589   /* FIXME: Moved this test to here to avoid stopping simulation if there are actions running on cpus and all cpus are
 590    * with availability = 0. This may cause an infinite loop if one cpu has a trace with periodicity = 0 and the other a
 591    * trace with periodicity > 0.
 592    * The options are: all traces with same periodicity(0 or >0) or we need to change the way how the events are managed
 593    */
 594   if (time_delta < 0) {
 595     XBT_DEBUG("No next event at all. Bail out now.");
 596     return -1.0;
 597   }
 598
 599   XBT_DEBUG("Duration set to %f", time_delta);
 600
 601   // Bump the time: jump into the future
 602   now_ += time_delta;
 603
 604   // Inform the models of the date change
 605   for (auto const& model : models_)
 606     model->update_actions_state(now_, time_delta);
 607
 608   s4u::Engine::on_time_advance(time_delta);
 609
 610   return time_delta;
 611 }
 612
 613 void EngineImpl::run(double max_date)
 614 {
 615   seal_platform();
 616
 617   if (MC_is_active()) {
 618 #if SIMGRID_HAVE_MC
 619     mc::AppSide::get()->main_loop();
 620 #else
 621     xbt_die("MC_is_active() is not supposed to return true in non-MC settings");
 622 #endif
 623     THROW_IMPOSSIBLE; // main_loop never returns
 624   }
 625
 626   if (MC_record_replay_is_active()) {
 627     mc::RecordTrace::replay(MC_record_path());
 628     empty_trash();
 629     return;
 630   }
 631
 632   double elapsed_time = -1;
 633   const std::set<s4u::Activity*>* vetoed_activities = s4u::Activity::get_vetoed_activities();
 634
 635   do {
 636     XBT_DEBUG("New Schedule Round; size(queue)=%zu", actors_to_run_.size());
 637
 638     if (cfg_breakpoint >= 0.0 && simgrid_get_clock() >= cfg_breakpoint) {
 639       XBT_DEBUG("Breakpoint reached (%g)", cfg_breakpoint.get());
 640       cfg_breakpoint = -1.0; // Let the simulation continue without hiting the breakpoint again and again
 641 #ifdef SIGTRAP
 642       std::raise(SIGTRAP);
 643 #else
 644       std::raise(SIGABRT);
 645 #endif
 646     }
 647
 648     while (not actors_to_run_.empty()) {
 649       XBT_DEBUG("New Sub-Schedule Round; size(queue)=%zu", actors_to_run_.size());
 650
 651       /* Run all actors that are ready to run, possibly in parallel */
 652       run_all_actors();
 653
 654       /* answer sequentially and in a fixed arbitrary order all the simcalls that were issued during that sub-round.
 655        * The order must be fixed for the simulation to be reproducible (see RR-7653). It's OK here because only maestro
 656        * changes the list. Killer actors are moved to the end to let victims finish their simcall before dying, but
 657        * the order remains reproducible (even if arbitrarily). No need to sort the vector for sake of reproducibility.
 658        */
 659       for (auto const& actor : actors_that_ran_)
 660         if (actor->simcall_.call_ != actor::Simcall::Type::NONE)
 661           actor->simcall_handle(0);
 662
 663       handle_ended_actions();
 664
 665       /* If only daemon actors remain, cancel their actions, mark them to die and reschedule them */
 666       if (actor_list_.size() == daemons_.size())
 667         for (auto const& dmon : daemons_) {
 668           XBT_DEBUG("Kill %s", dmon->get_cname());
 669           maestro_->kill(dmon);
 670         }
 671     }
 672
 673     // Compute the max_date of the next solve.
 674     // It's either when a timer occurs, or when user-specified deadline is reached, or -1 if none is given
 675     double next_time = timer::Timer::next();
 676     if (next_time < 0 && max_date > -1) {
 677       next_time = max_date;
 678     } else if (next_time > -1 && max_date > -1) { // either both <0, or both >0
 679       next_time = std::min(next_time, max_date);
 680     }
 681
 682     XBT_DEBUG("Calling solve(%g) %g", next_time, now_);
 683     elapsed_time = solve(next_time);
 684     XBT_DEBUG("Moving time ahead. NOW=%g; elapsed: %g", now_, elapsed_time);
 685
 686     // Execute timers until there isn't anything to be done:
 687     bool again = false;
 688     do {
 689       again = timer::Timer::execute_all();
 690       handle_ended_actions();
 691     } while (again);
 692
 693     /* Clean actors to destroy */
 694     empty_trash();
 695
 696     XBT_DEBUG("### elapsed time %f, #actors %zu, #to_run %zu, #vetoed %d", elapsed_time, actor_list_.size(),
 697               actors_to_run_.size(), (vetoed_activities == nullptr ? -1 : static_cast<int>(vetoed_activities->size())));
 698
 699     if (elapsed_time < 0. && actors_to_run_.empty() && not actor_list_.empty()) {
 700       if (actor_list_.size() <= daemons_.size()) {
 701         XBT_CRITICAL("Oops! Daemon actors cannot do any blocking activity (communications, synchronization, etc) "
 702                      "once the simulation is over. Please fix your on_exit() functions.");
 703       } else {
 704         XBT_CRITICAL("Oops! Deadlock detected, some activities are still around but will never complete. "
 705                      "This usually happens when the user code is not perfectly clean.");
 706       }
 707       display_all_actor_status();
 708       simgrid::s4u::Engine::on_deadlock();
 709       for (auto const& [_, actor] : actor_list_) {
 710         XBT_DEBUG("Kill %s", actor->get_cname());
 711         maestro_->kill(actor);
 712       }
 713     }
 714   } while ((vetoed_activities == nullptr || vetoed_activities->empty()) &&
 715            ((elapsed_time > -1.0 && not double_equals(max_date, now_, 0.00001)) || has_actors_to_run()));
 716
 717   if (not actor_list_.empty() && max_date < 0 && not(vetoed_activities == nullptr || vetoed_activities->empty()))
 718     THROW_IMPOSSIBLE;
 719
 720   simgrid::s4u::Engine::on_simulation_end();
 721 }
 722
 723 double EngineImpl::get_clock()
 724 {
 725   return now_;
 726 }
 727 } // namespace simgrid::kernel