src/surf/surf_c_bindings.cpp

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