Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
Merge branch 'master' of https://framagit.org/simgrid/simgrid
[simgrid.git] / src / mc / explo / DFSExplorer.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 "src/mc/explo/DFSExplorer.hpp"
7 #include "src/mc/VisitedState.hpp"
8 #include "src/mc/mc_config.hpp"
9 #include "src/mc/mc_exit.hpp"
10 #include "src/mc/mc_private.hpp"
11 #include "src/mc/mc_record.hpp"
12 #include "src/mc/transition/Transition.hpp"
13
14 #include "src/xbt/mmalloc/mmprivate.h"
15 #include "xbt/log.h"
16 #include "xbt/string.hpp"
17 #include "xbt/sysdep.h"
18
19 #include <cassert>
20 #include <cstdio>
21
22 #include <memory>
23 #include <string>
24 #include <vector>
25
26 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(mc_dfs, mc, "DFS exploration algorithm of the model-checker");
27
28 namespace simgrid::mc {
29
30 xbt::signal<void(RemoteApp&)> DFSExplorer::on_exploration_start_signal;
31 xbt::signal<void(RemoteApp&)> DFSExplorer::on_backtracking_signal;
32
33 xbt::signal<void(State*, RemoteApp&)> DFSExplorer::on_state_creation_signal;
34
35 xbt::signal<void(State*, RemoteApp&)> DFSExplorer::on_restore_system_state_signal;
36 xbt::signal<void(RemoteApp&)> DFSExplorer::on_restore_initial_state_signal;
37 xbt::signal<void(Transition*, RemoteApp&)> DFSExplorer::on_transition_replay_signal;
38 xbt::signal<void(Transition*, RemoteApp&)> DFSExplorer::on_transition_execute_signal;
39
40 xbt::signal<void(RemoteApp&)> DFSExplorer::on_log_state_signal;
41
42 void DFSExplorer::check_non_termination(const State* current_state)
43 {
44   for (auto const& state : stack_) {
45     if (state->get_system_state()->equals_to(*current_state->get_system_state(),
46                                              *get_remote_app().get_remote_process_memory())) {
47       XBT_INFO("Non-progressive cycle: state %ld -> state %ld", state->get_num(), current_state->get_num());
48       XBT_INFO("******************************************");
49       XBT_INFO("*** NON-PROGRESSIVE CYCLE DETECTED ***");
50       XBT_INFO("******************************************");
51       XBT_INFO("Counter-example execution trace:");
52       for (auto const& s : get_textual_trace())
53         XBT_INFO("  %s", s.c_str());
54       XBT_INFO("You can debug the problem (and see the whole details) by rerunning out of simgrid-mc with "
55                "--cfg=model-check/replay:'%s'",
56                get_record_trace().to_string().c_str());
57       log_state();
58
59       throw TerminationError();
60     }
61   }
62 }
63
64 RecordTrace DFSExplorer::get_record_trace() // override
65 {
66   RecordTrace res;
67   for (auto const& state : stack_)
68     res.push_back(state->get_transition());
69   return res;
70 }
71
72 std::vector<std::string> DFSExplorer::get_textual_trace() // override
73 {
74   std::vector<std::string> trace;
75   for (auto const& state : stack_) {
76     const auto* t = state->get_transition();
77     trace.push_back(xbt::string_printf("%ld: %s", t->aid_, t->to_string().c_str()));
78   }
79   return trace;
80 }
81
82 void DFSExplorer::log_state() // override
83 {
84   on_log_state_signal(get_remote_app());
85   XBT_INFO("DFS exploration ended. %ld unique states visited; %lu backtracks (%lu transition replays, %lu states "
86            "visited overall)",
87            State::get_expanded_states(), backtrack_count_, visited_states_count_,
88            Transition::get_replayed_transitions());
89   Exploration::log_state();
90 }
91
92 void DFSExplorer::run()
93 {
94   on_exploration_start_signal(get_remote_app());
95   /* This function runs the DFS algorithm the state space.
96    * We do so iteratively instead of recursively, dealing with the call stack manually.
97    * This allows one to explore the call stack at will. */
98
99   while (not stack_.empty()) {
100     /* Get current state */
101     State* state = stack_.back().get();
102
103     XBT_DEBUG("**************************************************");
104     XBT_DEBUG("Exploration depth=%zu (state:#%ld; %zu interleaves todo)", stack_.size(), state->get_num(),
105               state->count_todo());
106
107     visited_states_count_++;
108
109     // Backtrack if we reached the maximum depth
110     if (stack_.size() > (std::size_t)_sg_mc_max_depth) {
111       if (reduction_mode_ == ReductionMode::dpor) {
112         XBT_ERROR("/!\\ Max depth of %d reached! THIS WILL PROBABLY BREAK the dpor reduction /!\\",
113                   _sg_mc_max_depth.get());
114         XBT_ERROR("/!\\ If bad things happen, disable dpor with --cfg=model-check/reduction:none /!\\");
115       } else
116         XBT_WARN("/!\\ Max depth reached ! /!\\ ");
117       this->backtrack();
118       continue;
119     }
120
121     // Backtrack if we are revisiting a state we saw previously while applying state-equality reduction
122     if (visited_state_ != nullptr) {
123       XBT_DEBUG("State already visited (equal to state %ld), exploration stopped on this path.",
124                 visited_state_->original_num_ == -1 ? visited_state_->num_ : visited_state_->original_num_);
125
126       visited_state_ = nullptr;
127       this->backtrack();
128       continue;
129     }
130
131     // Search for the next transition
132     // next_transition returns a pair<aid_t, double> in case we want to consider multiple state
133     auto [next, _] = state->next_transition_guided();
134
135     if (next < 0) { // If there is no more transition in the current state, backtrack.
136       XBT_VERB("There remains %lu actors, but none to interleave (depth %zu).", state->get_actor_count(),
137                stack_.size() + 1);
138
139       if (state->get_actor_count() == 0) {
140         get_remote_app().finalize_app();
141         XBT_VERB("Execution came to an end at %s (state: %ld, depth: %zu)", get_record_trace().to_string().c_str(),
142                  state->get_num(), stack_.size());
143       }
144
145       this->backtrack();
146       continue;
147     }
148
149     if (_sg_mc_sleep_set && XBT_LOG_ISENABLED(mc_dfs, xbt_log_priority_verbose)) {
150       XBT_VERB("Sleep set actually containing:");
151       for (auto& [aid, transition] : state->get_sleep_set())
152         XBT_VERB("  <%ld,%s>", aid, transition.to_string().c_str());
153     }
154
155     /* Actually answer the request: let's execute the selected request (MCed does one step) */
156     state->execute_next(next, get_remote_app());
157     on_transition_execute_signal(state->get_transition(), get_remote_app());
158
159     // If there are processes to interleave and the maximum depth has not been
160     // reached then perform one step of the exploration algorithm.
161     XBT_VERB("Execute %ld: %.60s (stack depth: %zu, state: %ld, %zu interleaves)", state->get_transition()->aid_,
162              state->get_transition()->to_string().c_str(), stack_.size(), state->get_num(), state->count_todo());
163
164     /* Create the new expanded state (copy the state of MCed into our MCer data) */
165     std::shared_ptr<State> next_state;
166     next_state = std::make_shared<State>(get_remote_app(), state);
167     on_state_creation_signal(next_state.get(), get_remote_app());
168
169     /* Sleep set procedure:
170      * adding the taken transition to the sleep set of the original state.
171      * <!> Since the parent sleep set is used to compute the child sleep set, this need to be
172      * done after next_state creation */
173     XBT_DEBUG("Marking Transition >>%s<< of process %ld done and adding it to the sleep set",
174               state->get_transition()->to_string().c_str(), state->get_transition()->aid_);
175     state->add_sleep_set(state->get_transition()); // Actors are marked done when they are considerd in ActorState
176
177     /* DPOR persistent set procedure:
178      * for each new transition considered, check if it depends on any other previous transition executed before it
179      * on another process. If there exists one, find the more recent, and add its process to the interleave set.
180      * If the process is not enabled at this  point, then add every enabled process to the interleave */
181     if (reduction_mode_ == ReductionMode::dpor) {
182       aid_t issuer_id   = state->get_transition()->aid_;
183       stack_t tmp_stack = std::list(stack_);
184       while (not tmp_stack.empty()) {
185         State* prev_state = tmp_stack.back().get();
186         if (state->get_transition()->aid_ == prev_state->get_transition()->aid_) {
187           XBT_DEBUG("Simcall >>%s<< and >>%s<< with same issuer %ld", state->get_transition()->to_string().c_str(),
188                     prev_state->get_transition()->to_string().c_str(), issuer_id);
189           tmp_stack.pop_back();
190           continue;
191         } else if (prev_state->get_transition()->depends(state->get_transition())) {
192           XBT_VERB("Dependent Transitions:");
193           XBT_VERB("  %s (state=%ld)", prev_state->get_transition()->to_string().c_str(), prev_state->get_num());
194           XBT_VERB("  %s (state=%ld)", state->get_transition()->to_string().c_str(), state->get_num());
195
196           if (prev_state->is_actor_enabled(issuer_id)) {
197             if (not prev_state->is_actor_done(issuer_id)) {
198               prev_state->consider_one(issuer_id);
199               opened_states.emplace_back(tmp_stack);
200             } else
201               XBT_DEBUG("Actor %ld is already in done set: no need to explore it again", issuer_id);
202           } else {
203             XBT_DEBUG("Actor %ld is not enabled: DPOR may be failing. To stay sound, we are marking every enabled "
204                       "transition as todo",
205                       issuer_id);
206             prev_state->consider_all();
207             opened_states.emplace_back(tmp_stack);
208           }
209           break;
210         } else {
211           XBT_VERB("INDEPENDENT Transitions:");
212           XBT_VERB("  %s (state=%ld)", prev_state->get_transition()->to_string().c_str(), prev_state->get_num());
213           XBT_VERB("  %s (state=%ld)", state->get_transition()->to_string().c_str(), state->get_num());
214         }
215         tmp_stack.pop_back();
216       }
217     }
218
219     if (_sg_mc_termination)
220       this->check_non_termination(next_state.get());
221
222     /* Check whether we already explored next_state in the past (but only if interested in state-equality reduction) */
223     if (_sg_mc_max_visited_states > 0)
224       visited_state_ = visited_states_.addVisitedState(next_state->get_num(), next_state.get(), get_remote_app());
225
226     stack_.push_back(std::move(next_state));
227
228     /* If this is a new state (or if we don't care about state-equality reduction) */
229     if (visited_state_ == nullptr) {
230       /* Get an enabled process and insert it in the interleave set of the next state */
231       if (reduction_mode_ == ReductionMode::dpor)
232         stack_.back()->consider_best(); // Take only one transition if DPOR: others may be considered later if required
233       else {
234         stack_.back()->consider_all();
235         opened_states.emplace_back(stack_);
236       }
237       dot_output("\"%ld\" -> \"%ld\" [%s];\n", state->get_num(), stack_.back()->get_num(),
238                  state->get_transition()->dot_string().c_str());
239     } else
240       dot_output("\"%ld\" -> \"%ld\" [%s];\n", state->get_num(),
241                  visited_state_->original_num_ == -1 ? visited_state_->num_ : visited_state_->original_num_,
242                  state->get_transition()->dot_string().c_str());
243   }
244
245   log_state();
246 }
247
248 void DFSExplorer::backtrack()
249 {
250   backtrack_count_++;
251   XBT_VERB("Backtracking from %s", get_record_trace().to_string().c_str());
252   on_backtracking_signal(get_remote_app());
253   get_remote_app().check_deadlock();
254
255   // if no backtracking point, then set the stack_ to empty so we can end the exploration
256   if (opened_states.size() == 0) {
257     stack_ = std::list<std::shared_ptr<State>>();
258     return;
259   }
260
261   /* We may backtrack from somewhere either because it's leaf, or because every enabled process are in done/sleep set.
262    * In the first case, we need to remove the last transition corresponding to the Finalize */
263   if (stack_.back()->get_transition()->aid_ == 0)
264     stack_.pop_back();
265
266   stack_t backtrack;
267   double min_dist = std::numeric_limits<double>::infinity();
268   aid_t min_aid   = -1;
269   for (auto& stack : opened_states) {
270     auto [aid, dist] = stack.back()->next_transition_guided();
271     if (aid == -1)
272       continue;
273     if (dist < min_dist) {
274       min_dist  = dist;
275       min_aid   = aid;
276       backtrack = stack;
277     }
278   }
279
280   if (min_aid == -1) {
281     stack_ = std::list<std::shared_ptr<State>>();
282     return;
283   }
284
285   if (backtrack.back()->count_todo() <= 1)
286     opened_states.pop_back();
287   XBT_VERB("%lu todo actors in last states of the next backtracking point %s (remaining %lu opened stacks)",
288            backtrack.back()->count_todo(), backtrack.back()->get_transition()->to_string().c_str(),
289            opened_states.size());
290
291   /* If asked to rollback on a state that has a snapshot, restore it */
292   State* last_state = backtrack.back().get();
293   if (const auto* system_state = last_state->get_system_state()) {
294     system_state->restore(*get_remote_app().get_remote_process_memory());
295     on_restore_system_state_signal(last_state, get_remote_app());
296     return;
297   }
298
299   /* if no snapshot, we need to restore the initial state and replay the transitions */
300   get_remote_app().restore_initial_state();
301   on_restore_initial_state_signal(get_remote_app());
302
303   /* Traverse the stack from the state at position start and re-execute the transitions */
304   for (std::shared_ptr<State> const& state : backtrack) {
305     if (state == backtrack.back()) /* If we are arrived on the target state, don't replay the outgoing transition */
306       break;
307     state->get_transition()->replay(get_remote_app());
308     on_transition_replay_signal(state->get_transition(), get_remote_app());
309     visited_states_count_++;
310   }
311   stack_ = backtrack;
312   XBT_VERB(">> Backtracked to %s", get_record_trace().to_string().c_str());
313 }
314
315 // DFSExplorer::DFSExplorer(const std::vector<char*>& args, bool with_dpor) : Exploration(args, _sg_mc_termination) //
316 // UNCOMMENT TO ACTIVATE REFORKS
317 DFSExplorer::DFSExplorer(const std::vector<char*>& args, bool with_dpor)
318     : Exploration(args, true) // This version does not use reforks as it breaks
319 {
320   if (with_dpor)
321     reduction_mode_ = ReductionMode::dpor;
322   else
323     reduction_mode_ = ReductionMode::none;
324
325   if (_sg_mc_termination) {
326     if (with_dpor) {
327       XBT_INFO("Check non progressive cycles (turning DPOR off)");
328       reduction_mode_ = ReductionMode::none;
329     } else {
330       XBT_INFO("Check non progressive cycles");
331     }
332   } else
333     XBT_INFO("Start a DFS exploration. Reduction is: %s.", to_c_str(reduction_mode_));
334
335   auto initial_state = std::make_unique<State>(get_remote_app());
336
337   XBT_DEBUG("**************************************************");
338
339   stack_.push_back(std::move(initial_state));
340
341   /* Get an enabled actor and insert it in the interleave set of the initial state */
342   XBT_DEBUG("Initial state. %lu actors to consider", initial_state->get_actor_count());
343   if (reduction_mode_ == ReductionMode::dpor)
344     stack_.back()->consider_best();
345   else {
346     stack_.back()->consider_all();
347     opened_states.emplace_back(stack_);
348   }
349 }
350
351 Exploration* create_dfs_exploration(const std::vector<char*>& args, bool with_dpor)
352 {
353   return new DFSExplorer(args, with_dpor);
354 }
355
356 } // namespace simgrid::mc