src/surf/surf_c_bindings.cpp

   1 /* Copyright (c) 2013-2021. 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/s4u/Engine.hpp"
   7 #include "src/include/surf/surf.hpp"
   8 #include "src/instr/instr_private.hpp"
   9 #include "src/kernel/EngineImpl.hpp"
  10 #include "src/kernel/resource/DiskImpl.hpp"
  11 #include "src/kernel/resource/profile/FutureEvtSet.hpp"
  12 #include "src/plugins/vm/VirtualMachineImpl.hpp"
  13
  14 #include <algorithm>
  15
  16 XBT_LOG_EXTERNAL_DEFAULT_CATEGORY(surf_kernel);
  17
  18 /*********
  19  * TOOLS *
  20  *********/
  21
  22 extern double NOW;
  23
  24 void surf_presolve()
  25 {
  26   XBT_DEBUG("Consume all trace events occurring before the starting time.");
  27   double next_event_date;
  28   while ((next_event_date = simgrid::kernel::profile::future_evt_set.next_date()) != -1.0) {
  29     if (next_event_date > NOW)
  30       break;
  31
  32     double value                                  = -1.0;
  33     simgrid::kernel::resource::Resource* resource = nullptr;
  34     while (auto* event = simgrid::kernel::profile::future_evt_set.pop_leq(next_event_date, &value, &resource)) {
  35       if (value >= 0)
  36         resource->apply_event(event, value);
  37     }
  38   }
  39
  40   XBT_DEBUG("Set every models in the right state by updating them to 0.");
  41   for (auto const& model : simgrid::kernel::EngineImpl::get_instance()->get_all_models())
  42     model->update_actions_state(NOW, 0.0);
  43 }
  44
  45 double surf_solve(double max_date)
  46 {
  47   double time_delta                             = -1.0; /* duration */
  48   double value                                  = -1.0;
  49   simgrid::kernel::resource::Resource* resource = nullptr;
  50
  51   if (max_date != -1.0) {
  52     xbt_assert(max_date >= NOW, "You asked to simulate up to %f, but that's in the past already", max_date);
  53
  54     time_delta = max_date - NOW;
  55   }
  56
  57   XBT_DEBUG("Looking for next event in all models");
  58   auto engine = simgrid::kernel::EngineImpl::get_instance();
  59   for (auto model : engine->get_all_models()) {
  60     if (not model->next_occurring_event_is_idempotent()) {
  61       continue;
  62     }
  63     double next_event = model->next_occurring_event(NOW);
  64     if ((time_delta < 0.0 || next_event < time_delta) && next_event >= 0.0) {
  65       time_delta = next_event;
  66     }
  67   }
  68
  69   XBT_DEBUG("Min for resources (remember that NS3 don't update that value): %f", time_delta);
  70
  71   XBT_DEBUG("Looking for next trace event");
  72
  73   while (true) { // Handle next occurring events until none remains
  74     double next_event_date = simgrid::kernel::profile::future_evt_set.next_date();
  75     XBT_DEBUG("Next TRACE event: %f", next_event_date);
  76
  77     for (auto model : engine->get_all_models()) {
  78       /* Skip all idempotent models, they were already treated above
  79        * NS3 is the one to handled here */
  80       if (model->next_occurring_event_is_idempotent())
  81         continue;
  82
  83       if (next_event_date != -1.0) {
  84         time_delta = std::min(next_event_date - NOW, time_delta);
  85       } else {
  86         time_delta = std::max(next_event_date - NOW, time_delta); // Get the positive component
  87       }
  88
  89       XBT_DEBUG("Run the NS3 network at most %fs", time_delta);
  90       // run until min or next flow
  91       double model_next_action_end = model->next_occurring_event(time_delta);
  92
  93       XBT_DEBUG("Min for network : %f", model_next_action_end);
  94       if (model_next_action_end >= 0.0)
  95         time_delta = model_next_action_end;
  96     }
  97
  98     if (next_event_date < 0.0 || (next_event_date > NOW + time_delta)) {
  99       // next event may have already occurred or will after the next resource change, then bail out
 100       XBT_DEBUG("no next usable TRACE event. Stop searching for it");
 101       break;
 102     }
 103
 104     XBT_DEBUG("Updating models (min = %g, NOW = %g, next_event_date = %g)", time_delta, NOW, next_event_date);
 105
 106     while (auto* event = simgrid::kernel::profile::future_evt_set.pop_leq(next_event_date, &value, &resource)) {
 107       if (resource->is_used() || (watched_hosts().find(resource->get_cname()) != watched_hosts().end())) {
 108         time_delta = next_event_date - NOW;
 109         XBT_DEBUG("This event invalidates the next_occurring_event() computation of models. Next event set to %f",
 110                   time_delta);
 111       }
 112       // FIXME: I'm too lame to update NOW live, so I change it and restore it so that the real update with surf_min
 113       // will work
 114       double round_start = NOW;
 115       NOW                = next_event_date;
 116       /* update state of the corresponding resource to the new value. Does not touch lmm.
 117          It will be modified if needed when updating actions */
 118       XBT_DEBUG("Calling update_resource_state for resource %s", resource->get_cname());
 119       resource->apply_event(event, value);
 120       NOW = round_start;
 121     }
 122   }
 123
 124   /* FIXME: Moved this test to here to avoid stopping simulation if there are actions running on cpus and all cpus are
 125    * with availability = 0. This may cause an infinite loop if one cpu has a trace with periodicity = 0 and the other a
 126    * trace with periodicity > 0.
 127    * The options are: all traces with same periodicity(0 or >0) or we need to change the way how the events are managed
 128    */
 129   if (time_delta < 0) {
 130     XBT_DEBUG("No next event at all. Bail out now.");
 131     return -1.0;
 132   }
 133
 134   XBT_DEBUG("Duration set to %f", time_delta);
 135
 136   // Bump the time: jump into the future
 137   NOW = NOW + time_delta;
 138
 139   // Inform the models of the date change
 140   for (auto const& model : simgrid::kernel::EngineImpl::get_instance()->get_all_models())
 141     model->update_actions_state(NOW, time_delta);
 142
 143   simgrid::s4u::Engine::on_time_advance(time_delta);
 144
 145   return time_delta;
 146 }