Logo AND Algorithmique Numérique Distribuée

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