Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
model-checker : cosmetics in log message
[simgrid.git] / src / mc / mc_dpor.c
index e51d986..ac5e674 100644 (file)
-/* Copyright (c) 2008 Martin Quinson, Cristian Rosa.
  All rights reserved.                                          */
+/* Copyright (c) 2008-2014. The SimGrid Team.
* All rights reserved.                                                     */
 
 /* This program is free software; you can redistribute it and/or modify it
  * under the terms of the license (GNU LGPL) which comes with this package. */
 
-#include "private.h"
+#include "mc_private.h"
+
+#include "xbt/mmalloc/mmprivate.h"
 
 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(mc_dpor, mc,
                                 "Logging specific to MC DPOR exploration");
 
-/**
- *  \brief Initialize the DPOR exploration algorithm
- */
-void MC_dpor_init()
-{
-  mc_state_t initial_state = NULL;
-  smx_process_t process;
-  
-  /* Create the initial state and push it into the exploration stack */
-  MC_SET_RAW_MEM;
-  initial_state = MC_state_new();
-  xbt_fifo_unshift(mc_stack, initial_state);
-  MC_UNSET_RAW_MEM;
+/********** Global variables **********/
 
-  XBT_DEBUG("**************************************************");
-  XBT_DEBUG("Initial state");
+xbt_dynar_t visited_states;
+xbt_dict_t first_enabled_state;
+xbt_dynar_t initial_communications_pattern;
+xbt_dynar_t incomplete_communications_pattern;
+xbt_dynar_t communications_pattern;
+int nb_comm_pattern;
 
-  /* Wait for requests (schedules processes) */
-  MC_wait_for_requests();
+/********** Static functions ***********/
 
-  MC_SET_RAW_MEM;
-  /* Get an enabled process and insert it in the interleave set of the initial state */
-  xbt_swag_foreach(process, simix_global->process_list){
-    if(MC_process_is_enabled(process)){
-      MC_state_interleave_process(initial_state, process);
-      break;
-    }
-  }
-  MC_UNSET_RAW_MEM;
-    
-  /* FIXME: Update Statistics 
-  mc_stats->state_size +=
-      xbt_setset_set_size(initial_state->enabled_transitions); */
+static void comm_pattern_free(mc_comm_pattern_t p){
+  xbt_free(p->rdv);
+  xbt_free(p->data);
+  xbt_free(p);
+  p = NULL;
 }
 
+static void comm_pattern_free_voidp( void *p){
+  comm_pattern_free((mc_comm_pattern_t) * (void **)p);
+}
 
-/**
- *     \brief Perform the model-checking operation using a depth-first search exploration
- *         with Dynamic Partial Order Reductions
- */
-void MC_dpor(void)
-{
-  char *req_str;
-  int value;
-  smx_req_t req = NULL, prev_req = NULL;
-  mc_state_t state = NULL, prev_state = NULL, next_state = NULL;
-  smx_process_t process = NULL;
-  xbt_fifo_item_t item = NULL;
-
-  while (xbt_fifo_size(mc_stack) > 0) {
-
-    /* Get current state */
-    state = (mc_state_t) 
-      xbt_fifo_get_item_content(xbt_fifo_get_first_item(mc_stack));
-
-    XBT_DEBUG("**************************************************");
-    XBT_DEBUG("Exploration depth=%d (state=%p)(%u interleave)",
-           xbt_fifo_size(mc_stack), state,
-           MC_state_interleave_size(state));
-
-    /* Update statistics */
-    mc_stats->visited_states++;
-
-    /* If there are processes to interleave and the maximum depth has not been reached
-       then perform one step of the exploration algorithm */
-    if (xbt_fifo_size(mc_stack) < MAX_DEPTH &&
-        (req = MC_state_get_request(state, &value))) {
-
-      /* Debug information */
-      if(XBT_LOG_ISENABLED(mc_dpor, xbt_log_priority_debug)){
-        req_str = MC_request_to_string(req, value);
-        XBT_DEBUG("Execute: %s", req_str);
-        xbt_free(req_str);
-      }
-
-      MC_state_set_executed_request(state, req, value);
-      mc_stats->executed_transitions++;
-
-      /* Answer the request */
-      SIMIX_request_pre(req, value); /* After this call req is no longer usefull */
-
-      /* Wait for requests (schedules processes) */
-      MC_wait_for_requests();
-
-      /* Create the new expanded state */
-      MC_SET_RAW_MEM;
-      next_state = MC_state_new();
-      xbt_fifo_unshift(mc_stack, next_state);
-
-      /* Get an enabled process and insert it in the interleave set of the next state */
-      xbt_swag_foreach(process, simix_global->process_list){
-        if(MC_process_is_enabled(process)){
-          MC_state_interleave_process(next_state, process);
-          break;
-        }
-      }
-      MC_UNSET_RAW_MEM;
-
-      /* FIXME: Update Statistics
-      mc_stats->state_size +=
-          xbt_setset_set_size(next_state->enabled_transitions);*/
-
-      /* Let's loop again */
-
-      /* The interleave set is empty or the maximum depth is reached, let's back-track */
-    } else {
-      XBT_DEBUG("There are no more processes to interleave.");
-
-      /* Trash the current state, no longer needed */
-      MC_SET_RAW_MEM;
-      xbt_fifo_shift(mc_stack);
-      MC_state_delete(state);
-      MC_UNSET_RAW_MEM;
-
-      /* Check for deadlocks */
-      if(MC_deadlock_check()){
-        MC_show_deadlock(NULL);
-        return;
-      }
-
-      MC_SET_RAW_MEM;
-      /* Traverse the stack backwards until a state with a non empty interleave
-         set is found, deleting all the states that have it empty in the way.
-         For each deleted state, check if the request that has generated it 
-         (from it's predecesor state), depends on any other previous request 
-         executed before it. If it does then add it to the interleave set of the
-         state that executed that previous request. */
-      while ((state = xbt_fifo_shift(mc_stack)) != NULL) {
-        req = MC_state_get_internal_request(state);
-        xbt_fifo_foreach(mc_stack, item, prev_state, mc_state_t) {
-          if(MC_request_depend(req, MC_state_get_internal_request(prev_state))){
-            if(XBT_LOG_ISENABLED(mc_dpor, xbt_log_priority_debug)){
-              XBT_DEBUG("Dependent Transitions:");
-              prev_req = MC_state_get_executed_request(prev_state, &value);
-              req_str = MC_request_to_string(prev_req, value);
-              XBT_DEBUG("%s (state=%p)", req_str, prev_state);
-              xbt_free(req_str);
-              prev_req = MC_state_get_executed_request(state, &value);
-              req_str = MC_request_to_string(prev_req, value);
-              XBT_DEBUG("%s (state=%p)", req_str, state);
-              xbt_free(req_str);              
-            }
+static mc_comm_pattern_t get_comm_pattern_from_idx(xbt_dynar_t pattern, unsigned int *idx, e_smx_comm_type_t type, unsigned long proc){
+  mc_comm_pattern_t current_comm;
+  while(*idx < xbt_dynar_length(pattern)){
+    current_comm = (mc_comm_pattern_t)xbt_dynar_get_as(pattern, *idx, mc_comm_pattern_t);
+    if(current_comm->type == type && type == SIMIX_COMM_SEND){
+      if(current_comm->src_proc == proc)
+        return current_comm;
+    }else if(current_comm->type == type && type == SIMIX_COMM_RECEIVE){
+      if(current_comm->dst_proc == proc)
+        return current_comm;
+    }
+    (*idx)++;
+  }
+  return NULL;
+}
 
-            if(!MC_state_process_is_done(prev_state, req->issuer))
-              MC_state_interleave_process(prev_state, req->issuer);
-            else
-              XBT_DEBUG("Process %p is in done set", req->issuer);
+static int compare_comm_pattern(mc_comm_pattern_t comm1, mc_comm_pattern_t comm2){
+  if(strcmp(comm1->rdv, comm2->rdv) != 0)
+    return 1;
+  if(comm1->src_proc != comm2->src_proc)
+    return 1;
+  if(comm1->dst_proc != comm2->dst_proc)
+    return 1;
+  if(comm1->data_size != comm2->data_size)
+    return 1;
+  if(memcmp(comm1->data, comm2->data, comm1->data_size) != 0)
+    return 1;
+  return 0;
+}
 
-            break;
-          }
+static void deterministic_pattern(xbt_dynar_t initial_pattern, xbt_dynar_t pattern){
+
+  if(!xbt_dynar_is_empty(incomplete_communications_pattern))
+    xbt_die("Damn ! Some communications are incomplete that means one or several simcalls are not handle ... ");
+
+  unsigned int cursor = 0, send_index = 0, recv_index = 0;
+  mc_comm_pattern_t comm1, comm2;
+  int comm_comparison = 0;
+  int current_process = 0;
+  while(current_process < simix_process_maxpid){
+    while(cursor < xbt_dynar_length(initial_pattern)){
+      comm1 = (mc_comm_pattern_t)xbt_dynar_get_as(initial_pattern, cursor, mc_comm_pattern_t);
+      if(comm1->type == SIMIX_COMM_SEND && comm1->src_proc == current_process){
+        comm2 = get_comm_pattern_from_idx(pattern, &send_index, comm1->type, current_process);
+        comm_comparison = compare_comm_pattern(comm1, comm2);
+        if(comm_comparison == 1){
+          initial_state_safety->send_deterministic = 0;
+          initial_state_safety->comm_deterministic = 0;
+          return;
         }
-        if (MC_state_interleave_size(state)) {
-          /* We found a back-tracking point, let's loop */
-          xbt_fifo_unshift(mc_stack, state);
-          XBT_DEBUG("Back-tracking to depth %d", xbt_fifo_size(mc_stack));
-          MC_UNSET_RAW_MEM;
-          MC_replay(mc_stack);
-          break;
-        } else {
-          MC_state_delete(state);
+        send_index++;
+      }else if(comm1->type == SIMIX_COMM_RECEIVE && comm1->dst_proc == current_process){
+        comm2 = get_comm_pattern_from_idx(pattern, &recv_index, comm1->type, current_process);
+        comm_comparison = compare_comm_pattern(comm1, comm2);
+        if(comm_comparison == 1){
+          initial_state_safety->comm_deterministic = 0;
+          if(!_sg_mc_send_determinism)
+            return;
         }
+        recv_index++;
       }
-      MC_UNSET_RAW_MEM;
+      cursor++;
     }
+    cursor = 0;
+    send_index = 0;
+    recv_index = 0;
+    current_process++;
   }
-  MC_UNSET_RAW_MEM;
-  return;
 }
 
-
-/********************* DPOR without replay *********************/
-
-mc_state_ws_t new_state_ws(mc_snapshot_t s, mc_state_t gs){
-  mc_state_ws_t sws = NULL;
-  sws = xbt_new0(s_mc_state_ws_t, 1);
-  sws->system_state = s;
-  sws->graph_state = gs;
-  return sws;
+void complete_comm_pattern(xbt_dynar_t list, smx_action_t comm){
+  mc_comm_pattern_t current_pattern;
+  unsigned int cursor = 0;
+  int index;
+  int completed = 0;
+  void *addr_pointed;
+  xbt_dynar_foreach(incomplete_communications_pattern, cursor, index){
+    current_pattern = (mc_comm_pattern_t)xbt_dynar_get_as(list, index, mc_comm_pattern_t);
+    if(current_pattern->comm == comm){
+      current_pattern->src_proc = comm->comm.src_proc->pid;
+      current_pattern->dst_proc = comm->comm.dst_proc->pid;
+      current_pattern->src_host = simcall_host_get_name(comm->comm.src_proc->smx_host);
+      current_pattern->dst_host = simcall_host_get_name(comm->comm.dst_proc->smx_host);
+      if(current_pattern->data_size == -1){
+        current_pattern->data_size = *(comm->comm.dst_buff_size);
+        current_pattern->data = xbt_malloc0(current_pattern->data_size);
+        addr_pointed = *(void **)comm->comm.src_buff;
+        if(addr_pointed > std_heap && addr_pointed < ((xbt_mheap_t)std_heap)->breakval)
+          memcpy(current_pattern->data, addr_pointed, current_pattern->data_size);
+        else
+          memcpy(current_pattern->data, comm->comm.src_buff, current_pattern->data_size);
+      }
+      xbt_dynar_remove_at(incomplete_communications_pattern, cursor, NULL);
+      completed++;
+      if(completed == 2)
+        return;
+      cursor--;
+    }
+  }
 }
 
-void MC_dpor_stateful_init(){
-
-  XBT_DEBUG("**************************************************");
-  XBT_DEBUG("DPOR (stateful) init");
-  XBT_DEBUG("**************************************************");
-
-  mc_state_t initial_graph_state;
-  smx_process_t process; 
-  mc_snapshot_t initial_system_snapshot;
-  mc_state_ws_t initial_state ;
+void get_comm_pattern(xbt_dynar_t list, smx_simcall_t request, int call){
+  mc_comm_pattern_t pattern = NULL;
+  pattern = xbt_new0(s_mc_comm_pattern_t, 1);
+  pattern->num = ++nb_comm_pattern;
+  pattern->data_size = -1;
+  void * addr_pointed;
+  if(call == 1){ // ISEND
+    pattern->comm = simcall_comm_isend__get__result(request);
+    pattern->type = SIMIX_COMM_SEND;
+    pattern->src_proc = pattern->comm->comm.src_proc->pid;
+    pattern->src_host = simcall_host_get_name(request->issuer->smx_host);
+    pattern->data_size = pattern->comm->comm.src_buff_size;
+    pattern->data = xbt_malloc0(pattern->data_size);
+    addr_pointed = *(void **)pattern->comm->comm.src_buff;
+    if(addr_pointed > std_heap && addr_pointed < ((xbt_mheap_t)std_heap)->breakval)
+      memcpy(pattern->data, addr_pointed, pattern->data_size);
+    else
+      memcpy(pattern->data, pattern->comm->comm.src_buff, pattern->data_size);
+  }else{ // IRECV
+    pattern->comm = simcall_comm_irecv__get__result(request);
+    pattern->type = SIMIX_COMM_RECEIVE;
+    pattern->dst_proc = pattern->comm->comm.dst_proc->pid;
+    pattern->dst_host = simcall_host_get_name(pattern->comm->comm.dst_proc->smx_host);
+  }
   
-  MC_wait_for_requests();
+  if(pattern->comm->comm.rdv != NULL)
+    pattern->rdv = strdup(pattern->comm->comm.rdv->name);
+  else
+    pattern->rdv = strdup(pattern->comm->comm.rdv_cpy->name);
+  
+  xbt_dynar_push(list, &pattern);
 
-  MC_SET_RAW_MEM;
+  xbt_dynar_push_as(incomplete_communications_pattern, int, xbt_dynar_length(list) - 1);
 
-  initial_system_snapshot = xbt_new0(s_mc_snapshot_t, 1);
+}
 
-  initial_graph_state = MC_state_new();
-  xbt_swag_foreach(process, simix_global->process_list){
-    if(MC_process_is_enabled(process)){
-      MC_state_interleave_process(initial_graph_state, process);
-      break;
-    }
+static void print_communications_pattern(xbt_dynar_t comms_pattern){
+  unsigned int cursor = 0;
+  mc_comm_pattern_t current_comm;
+  xbt_dynar_foreach(comms_pattern, cursor, current_comm){
+    if(current_comm->type == SIMIX_COMM_SEND)
+      XBT_INFO("[(%lu) %s -> (%lu) %s] %s ", current_comm->src_proc, current_comm->src_host, current_comm->dst_proc, current_comm->dst_host, "iSend");
+    else
+      XBT_INFO("[(%lu) %s <- (%lu) %s] %s ", current_comm->dst_proc, current_comm->dst_host, current_comm->src_proc, current_comm->src_host, "iRecv");
   }
+}
 
-  MC_take_snapshot(initial_system_snapshot);
+static void visited_state_free(mc_visited_state_t state){
+  if(state){
+    MC_free_snapshot(state->system_state);
+    xbt_free(state);
+  }
+}
 
-  initial_state = new_state_ws(initial_system_snapshot, initial_graph_state);
-  xbt_fifo_unshift(mc_snapshot_stack, initial_state);
+static void visited_state_free_voidp(void *s){
+  visited_state_free((mc_visited_state_t) * (void **) s);
+}
+
+/** \brief Save the current state
+ *
+ *  \return Snapshot of the current state.
+ */
+static mc_visited_state_t visited_state_new(){
 
-  MC_UNSET_RAW_MEM;
+  mc_visited_state_t new_state = NULL;
+  new_state = xbt_new0(s_mc_visited_state_t, 1);
+  new_state->heap_bytes_used = mmalloc_get_bytes_used(std_heap);
+  new_state->nb_processes = xbt_swag_size(simix_global->process_list);
+  new_state->system_state = MC_take_snapshot(mc_stats->expanded_states);
+  new_state->num = mc_stats->expanded_states;
+  new_state->other_num = -1;
 
+  return new_state;
+  
 }
 
-void MC_dpor_stateful(){
+/** \brief Find a suitable subrange of candidate duplicates for a given state
+ *
+ *  \param all_ pairs dynamic array of states with candidate duplicates of the current state;
+ *  \param pair current state;
+ *  \param min (output) index of the beginning of the the subrange
+ *  \param max (output) index of the enf of the subrange
+ *
+ *  Given a suitably ordered array of state, this function extracts a subrange
+ *  (with index *min <= i <= *max) with candidate duplicates of the given state.
+ *  This function uses only fast discriminating criterions and does not use the
+ *  full state comparison algorithms.
+ *
+ *  The states in all_pairs MUST be ordered using a (given) weak order
+ *  (based on nb_processes and heap_bytes_used).
+ *  The subrange is the subrange of "equivalence" of the given state.
+ */
+static int get_search_interval(xbt_dynar_t all_states, mc_visited_state_t state, int *min, int *max){
+  XBT_VERB("Searching interval for state %i: nd_processes=%zu heap_bytes_used=%zu",
+    state->num, (size_t)state->nb_processes, (size_t)state->heap_bytes_used);
 
-  smx_process_t process = NULL;
-  
-  if(xbt_fifo_size(mc_snapshot_stack) == 0)
-    return;
+  int raw_mem_set = (mmalloc_get_current_heap() == mc_heap);
 
-  int value;
-  mc_state_t next_graph_state = NULL;
-  smx_req_t req = NULL, prev_req = NULL;
-  char *req_str;
-  xbt_fifo_item_t item = NULL;
+  MC_SET_MC_HEAP;
 
-  mc_snapshot_t next_snapshot;
-  mc_state_ws_t current_state;
-  mc_state_ws_t prev_state;
-  mc_state_ws_t next_state;
-  while(xbt_fifo_size(mc_snapshot_stack) > 0){
+  int cursor = 0, previous_cursor, next_cursor;
+  mc_visited_state_t state_test;
+  int start = 0;
+  int end = xbt_dynar_length(all_states) - 1;
+  
+  while(start <= end){
+    cursor = (start + end) / 2;
+    state_test = (mc_visited_state_t)xbt_dynar_get_as(all_states, cursor, mc_visited_state_t);
+    if(state_test->nb_processes < state->nb_processes){
+      start = cursor + 1;
+    }else if(state_test->nb_processes > state->nb_processes){
+      end = cursor - 1;
+    }else{
+      if(state_test->heap_bytes_used < state->heap_bytes_used){
+        start = cursor +1;
+      }else if(state_test->heap_bytes_used > state->heap_bytes_used){
+        end = cursor - 1;
+      }else{
+        *min = *max = cursor;
+        previous_cursor = cursor - 1;
+        while(previous_cursor >= 0){
+          state_test = (mc_visited_state_t)xbt_dynar_get_as(all_states, previous_cursor, mc_visited_state_t);
+          if(state_test->nb_processes != state->nb_processes || state_test->heap_bytes_used != state->heap_bytes_used)
+            break;
+          *min = previous_cursor;
+          previous_cursor--;
+        }
+        next_cursor = cursor + 1;
+        while(next_cursor < xbt_dynar_length(all_states)){
+          state_test = (mc_visited_state_t)xbt_dynar_get_as(all_states, next_cursor, mc_visited_state_t);
+          if(state_test->nb_processes != state->nb_processes || state_test->heap_bytes_used != state->heap_bytes_used)
+            break;
+          *max = next_cursor;
+          next_cursor++;
+        }
+        if(!raw_mem_set)
+          MC_SET_STD_HEAP;
+        return -1;
+      }
+     }
+  }
 
-    current_state = (mc_state_ws_t)xbt_fifo_get_item_content(xbt_fifo_get_first_item(mc_snapshot_stack));
+  if(!raw_mem_set)
+    MC_SET_STD_HEAP;
 
-    
-    XBT_DEBUG("**************************************************");
-    XBT_DEBUG("Depth : %d, State : %p , %u interleave", xbt_fifo_size(mc_snapshot_stack),current_state, MC_state_interleave_size(current_state->graph_state));
+  return cursor;
+}
 
-    mc_stats->visited_states++;
-    
+/** \brief Take a snapshot the current state and process it.
+ *
+ *  \return number of the duplicate state or -1 (not visited)
+ */
+static int is_visited_state(){
 
-    if((xbt_fifo_size(mc_snapshot_stack) < MAX_DEPTH) && (req = MC_state_get_request(current_state->graph_state, &value))){
+  if(_sg_mc_visited == 0)
+    return -1;
 
-      /* Debug information */
-      if(XBT_LOG_ISENABLED(mc_dpor, xbt_log_priority_debug)){
-       req_str = MC_request_to_string(req, value);
-       XBT_DEBUG("Execute: %s", req_str);
-       xbt_free(req_str);
-      }
+  int raw_mem_set = (mmalloc_get_current_heap() == mc_heap);
 
-      MC_state_set_executed_request(current_state->graph_state, req, value);
-      mc_stats->executed_transitions++;
-      
-      /* Answer the request */
-      SIMIX_request_pre(req, value);
-      
-      /* Wait for requests (schedules processes) */
-      MC_wait_for_requests();
-      
-      /* Create the new expanded graph_state */
-      MC_SET_RAW_MEM;
-      
-      next_graph_state = MC_state_new();
-      
-      /* Get an enabled process and insert it in the interleave set of the next graph_state */
-      xbt_swag_foreach(process, simix_global->process_list){
-       if(MC_process_is_enabled(process)){
-         MC_state_interleave_process(next_graph_state, process);
-         break;
-       }
-      }
+  MC_SET_MC_HEAP;
 
-      next_snapshot = xbt_new0(s_mc_snapshot_t, 1);
-      MC_take_snapshot(next_snapshot);
+  mc_visited_state_t new_state = visited_state_new();
+  
+  if(xbt_dynar_is_empty(visited_states)){
+
+    xbt_dynar_push(visited_states, &new_state); 
+
+    if(!raw_mem_set)
+      MC_SET_STD_HEAP;
+
+    return -1;
+
+  }else{
+
+    int min = -1, max = -1, index;
+    //int res;
+    mc_visited_state_t state_test;
+    int cursor;
+
+    index = get_search_interval(visited_states, new_state, &min, &max);
+
+    if(min != -1 && max != -1){
+
+      // Parallell implementation
+      /*res = xbt_parmap_mc_apply(parmap, snapshot_compare, xbt_dynar_get_ptr(visited_states, min), (max-min)+1, new_state);
+      if(res != -1){
+        state_test = (mc_visited_state_t)xbt_dynar_get_as(visited_states, (min+res)-1, mc_visited_state_t);
+        if(state_test->other_num == -1)
+          new_state->other_num = state_test->num;
+        else
+          new_state->other_num = state_test->other_num;
+        if(dot_output == NULL)
+          XBT_DEBUG("State %d already visited ! (equal to state %d)", new_state->num, state_test->num);
+        else
+          XBT_DEBUG("State %d already visited ! (equal to state %d (state %d in dot_output))", new_state->num, state_test->num, new_state->other_num);
+        xbt_dynar_remove_at(visited_states, (min + res) - 1, NULL);
+        xbt_dynar_insert_at(visited_states, (min+res) - 1, &new_state);
+        if(!raw_mem_set)
+          MC_SET_STD_HEAP;
+        return new_state->other_num;
+        }*/
+
+      cursor = min;
+      while(cursor <= max){
+        state_test = (mc_visited_state_t)xbt_dynar_get_as(visited_states, cursor, mc_visited_state_t);
+        if(snapshot_compare(state_test, new_state) == 0){
+          // The state has been visited:
+
+          if(state_test->other_num == -1)
+            new_state->other_num = state_test->num;
+          else
+            new_state->other_num = state_test->other_num;
+          if(dot_output == NULL)
+            XBT_DEBUG("State %d already visited ! (equal to state %d)", new_state->num, state_test->num);
+          else
+            XBT_DEBUG("State %d already visited ! (equal to state %d (state %d in dot_output))", new_state->num, state_test->num, new_state->other_num);
+
+          // Replace the old state with the new one (why?):
+          xbt_dynar_remove_at(visited_states, cursor, NULL);
+          xbt_dynar_insert_at(visited_states, cursor, &new_state);
+
+          if(!raw_mem_set)
+            MC_SET_STD_HEAP;
+          return new_state->other_num;
+        }
+        cursor++;
+      }
 
-      next_state = new_state_ws(next_snapshot, next_graph_state);
-      xbt_fifo_unshift(mc_snapshot_stack, next_state);
-      
-      MC_UNSET_RAW_MEM;
+      // The state has not been visited, add it to the list:
+      xbt_dynar_insert_at(visited_states, min, &new_state);
 
     }else{
-      XBT_DEBUG("There are no more processes to interleave.");
-      
-      /* Trash the current state, no longer needed */
-      MC_SET_RAW_MEM;
-      xbt_fifo_shift(mc_snapshot_stack);
-      MC_UNSET_RAW_MEM;
 
-      /* Check for deadlocks */
-      if(MC_deadlock_check()){
-        MC_show_deadlock_stateful(NULL);
-        return;
+      // The state has not been visited: insert the state in the dynamic array.
+      state_test = (mc_visited_state_t)xbt_dynar_get_as(visited_states, index, mc_visited_state_t);
+      if(state_test->nb_processes < new_state->nb_processes){
+        xbt_dynar_insert_at(visited_states, index+1, &new_state);
+      }else{
+        if(state_test->heap_bytes_used < new_state->heap_bytes_used)
+          xbt_dynar_insert_at(visited_states, index + 1, &new_state);
+        else
+          xbt_dynar_insert_at(visited_states, index, &new_state);
       }
 
-      MC_SET_RAW_MEM;
-      while((current_state = xbt_fifo_shift(mc_snapshot_stack)) != NULL){
-       req = MC_state_get_internal_request(current_state->graph_state);
-       xbt_fifo_foreach(mc_snapshot_stack, item, prev_state, mc_state_ws_t) {
-          if(MC_request_depend(req, MC_state_get_internal_request(prev_state->graph_state))){
-            if(XBT_LOG_ISENABLED(mc_dpor, xbt_log_priority_debug)){
-              XBT_DEBUG("Dependent Transitions:");
-              prev_req = MC_state_get_executed_request(prev_state->graph_state, &value);
-              req_str = MC_request_to_string(prev_req, value);
-              XBT_DEBUG("%s (state=%p)", req_str, prev_state->graph_state);
-              xbt_free(req_str);
-              prev_req = MC_state_get_executed_request(current_state->graph_state, &value);
-              req_str = MC_request_to_string(prev_req, value);
-              XBT_DEBUG("%s (state=%p)", req_str, current_state->graph_state);
-              xbt_free(req_str);              
-            }
-
-            if(!MC_state_process_is_done(prev_state->graph_state, req->issuer)){
-              MC_state_interleave_process(prev_state->graph_state, req->issuer);
-             
-           } else {
-              XBT_DEBUG("Process %p is in done set", req->issuer);
-           }
+    }
 
-            break;
-          }
+    // We have reached the maximum number of stored states;
+    if(xbt_dynar_length(visited_states) > _sg_mc_visited){
+
+      // Find the (index of the) older state:
+      int min2 = mc_stats->expanded_states;
+      unsigned int cursor2 = 0;
+      unsigned int index2 = 0;
+      xbt_dynar_foreach(visited_states, cursor2, state_test){
+        if(state_test->num < min2){
+          index2 = cursor2;
+          min2 = state_test->num;
         }
-
-       if(MC_state_interleave_size(current_state->graph_state)){
-         MC_restore_snapshot(current_state->system_state);
-         xbt_fifo_unshift(mc_snapshot_stack, current_state);
-         XBT_DEBUG("Back-tracking to depth %d", xbt_fifo_size(mc_snapshot_stack));
-         MC_UNSET_RAW_MEM;
-         break;
-       }
       }
 
-      MC_UNSET_RAW_MEM;
+      // and drop it:
+      xbt_dynar_remove_at(visited_states, index2, NULL);
+    }
 
-    } 
+    if(!raw_mem_set)
+      MC_SET_STD_HEAP;
+    
+    return -1;
+    
   }
-  MC_UNSET_RAW_MEM;
-  return;
 }
 
-
-
-/************ DPOR 2 (invisible and independant transitions) ************/
-
-xbt_dynar_t reached_pairs_prop;
-xbt_dynar_t visible_transitions;
-xbt_dynar_t enabled_processes;
-
-mc_prop_ato_t new_proposition(char* id, int value){
-  mc_prop_ato_t prop = NULL;
-  prop = xbt_new0(s_mc_prop_ato_t, 1);
-  prop->id = strdup(id);
-  prop->value = value;
-  return prop;
-}
-
-mc_pair_prop_t new_pair_prop(mc_snapshot_t sn, mc_state_t sg, xbt_state_t st){
-  mc_pair_prop_t pair = NULL;
-  pair = xbt_new0(s_mc_pair_prop_t, 1);
-  pair->system_state = sn;
-  pair->automaton_state = st;
-  pair->graph_state = sg;
-  pair->propositions = xbt_dynar_new(sizeof(mc_prop_ato_t), NULL);
-  pair->fully_expanded = 0;
-  pair->interleave = 0;
-  mc_stats_pair->expanded_pairs++;
-  return pair;
-}
-
-int reached_prop(mc_pair_prop_t pair){
-
-  char* hash_reached = malloc(sizeof(char)*160);
-  unsigned int c= 0;
-
-  MC_SET_RAW_MEM;
-  char *hash = malloc(sizeof(char)*160);
-  xbt_sha((char *)&pair, hash);
-  xbt_dynar_foreach(reached_pairs_prop, c, hash_reached){
-    if(strcmp(hash, hash_reached) == 0){
-      MC_UNSET_RAW_MEM;
-      return 1;
-    }
-  }
+/**
+ *  \brief Initialize the DPOR exploration algorithm
+ */
+void MC_dpor_init()
+{
   
-  MC_UNSET_RAW_MEM;
-  return 0;
-}
+  int raw_mem_set = (mmalloc_get_current_heap() == mc_heap);
 
-void set_pair_prop_reached(mc_pair_prop_t pair){
+  mc_state_t initial_state = NULL;
+  smx_process_t process;
+  
+  /* Create the initial state and push it into the exploration stack */
+  MC_SET_MC_HEAP;
 
-  if(reached_prop(pair) == 0){
-    MC_SET_RAW_MEM;
-    char *hash = malloc(sizeof(char)*160) ;
-    xbt_sha((char *)&pair, hash);
-    xbt_dynar_push(reached_pairs_prop, &hash); 
-    MC_UNSET_RAW_MEM;
-  }
+  if(_sg_mc_visited > 0)
+    visited_states = xbt_dynar_new(sizeof(mc_visited_state_t), visited_state_free_voidp);
 
-}
+  if(mc_reduce_kind == e_mc_reduce_dpor)
+    first_enabled_state = xbt_dict_new_homogeneous(&xbt_free_f);
 
-int invisible(mc_pair_prop_t p, mc_pair_prop_t np){
-  mc_prop_ato_t prop1;
-  mc_prop_ato_t prop2;
-  int i;
-  for(i=0;i<xbt_dynar_length(p->propositions); i++){
-    prop1 = xbt_dynar_get_as(p->propositions, i, mc_prop_ato_t);
-    prop2 = xbt_dynar_get_as(np->propositions, i, mc_prop_ato_t);
-    //XBT_DEBUG("test invisible : prop 1 = %s : %d / prop 2 = %s : %d", prop1->id, prop1->value, prop2->id, prop2->value);
-    if(prop1->value != prop2->value)
-         return 0;
+  if(_sg_mc_comms_determinism || _sg_mc_send_determinism){
+    initial_communications_pattern = xbt_dynar_new(sizeof(mc_comm_pattern_t), comm_pattern_free_voidp);
+    communications_pattern = xbt_dynar_new(sizeof(mc_comm_pattern_t), comm_pattern_free_voidp);
+    incomplete_communications_pattern = xbt_dynar_new(sizeof(int), NULL);
+    nb_comm_pattern = 0;
   }
-  return 1; 
 
-}
-
-void set_fully_expanded(mc_pair_prop_t pair){
-  pair->fully_expanded = 1;
-}
+  initial_state = MC_state_new();
 
-void MC_dpor2_init(xbt_automaton_t a)
-{
-  mc_pair_prop_t initial_pair = NULL;
-  mc_state_t initial_graph_state = NULL;
-  mc_snapshot_t initial_system_state = NULL;
-  xbt_state_t initial_automaton_state = NULL;
+  MC_SET_STD_HEAP;
 
-  MC_SET_RAW_MEM;
-  initial_graph_state = MC_state_new();
-  MC_UNSET_RAW_MEM;
+  XBT_DEBUG("**************************************************");
+  XBT_DEBUG("Initial state");
 
   /* Wait for requests (schedules processes) */
   MC_wait_for_requests();
 
-  unsigned int cursor = 0;
-  unsigned int cursor2 = 0;
-  xbt_state_t automaton_state = NULL;
-  int res;
-  xbt_transition_t transition_succ;
-
-  xbt_dynar_foreach(a->states, cursor, automaton_state){
-    if(automaton_state->type == -1){
-      xbt_dynar_foreach(automaton_state->out, cursor2, transition_succ){
-       res = MC_automaton_evaluate_label(a, transition_succ->label);
-       if((res == 1) || (res == 2)){
-         initial_automaton_state = transition_succ->dst;
-         break;
-       }
-      }
-    }
+  MC_ignore_heap(simix_global->process_to_run->data, 0);
+  MC_ignore_heap(simix_global->process_that_ran->data, 0);
 
-    if(xbt_fifo_size(mc_snapshot_stack) > 0)
-      break;
+  MC_SET_MC_HEAP;
+  /* Get an enabled process and insert it in the interleave set of the initial state */
+  xbt_swag_foreach(process, simix_global->process_list){
+    if(MC_process_is_enabled(process)){
+      MC_state_interleave_process(initial_state, process);
+      if(mc_reduce_kind != e_mc_reduce_none)
+        break;
+    }
   }
 
-  if(xbt_fifo_size(mc_snapshot_stack) == 0){
-    cursor = 0;
-    xbt_dynar_foreach(a->states, cursor, automaton_state){
-      if(automaton_state->type == -1){
-       initial_automaton_state = automaton_state;
-       break;
+  xbt_fifo_unshift(mc_stack_safety, initial_state);
+
+  if(mc_reduce_kind == e_mc_reduce_dpor){
+    /* To ensure the soundness of DPOR, we have to keep a list of 
+       processes which are still enabled at each step of the exploration. 
+       If max depth is reached, we interleave them in the state in which they have 
+       been enabled for the first time. */
+    xbt_swag_foreach(process, simix_global->process_list){
+      if(MC_process_is_enabled(process)){
+        char *key = bprintf("%lu", process->pid);
+        char *data = bprintf("%d", xbt_fifo_size(mc_stack_safety));
+        xbt_dict_set(first_enabled_state, key, data, NULL);
+        xbt_free(key);
       }
     }
   }
 
-  reached_pairs_prop = xbt_dynar_new(sizeof(char *), NULL); 
-  visible_transitions = xbt_dynar_new(sizeof(s_smx_req_t), NULL);
-  enabled_processes = xbt_dynar_new(sizeof(smx_process_t), NULL);
-
-  MC_SET_RAW_MEM;
-  initial_system_state = xbt_new0(s_mc_snapshot_t, 1);
-  MC_take_snapshot(initial_system_state);
-  initial_pair = new_pair_prop(initial_system_state, initial_graph_state, initial_automaton_state);
-  cursor = 0;
-  xbt_propositional_symbol_t ps;
-  xbt_dynar_foreach(a->propositional_symbols, cursor, ps){
-    int (*f)() = ps->function;
-    int val = (*f)();
-    mc_prop_ato_t pa = new_proposition(ps->pred, val);
-    xbt_dynar_push(initial_pair->propositions, &pa);
-  } 
-  xbt_fifo_unshift(mc_snapshot_stack, initial_pair);
-  MC_UNSET_RAW_MEM;
-
-
-  XBT_DEBUG("**************************************************");
-  XBT_DEBUG("Initial pair");
+  MC_SET_STD_HEAP;
 
-  MC_dpor2(a, 0);
-    
+  if(raw_mem_set)
+    MC_SET_MC_HEAP;
+  else
+    MC_SET_STD_HEAP;
+  
 }
 
 
-void MC_dpor2(xbt_automaton_t a, int search_cycle)
+/** \brief Model-check the application using a DFS exploration
+ *         with DPOR (Dynamic Partial Order Reductions)
+ */
+void MC_dpor(void)
 {
-  char *req_str;
+
+  char *req_str = NULL;
   int value;
-  smx_req_t req = NULL, prev_req = NULL;
-  mc_state_t next_graph_state = NULL;
-  mc_snapshot_t next_system_state = NULL;
-  xbt_state_t next_automaton_state = NULL;
-  xbt_transition_t transition_succ;
-  unsigned int cursor;
-  int res;
-  mc_pair_prop_t pair = NULL, next_pair = NULL, prev_pair = NULL;
+  smx_simcall_t req = NULL, prev_req = NULL;
+  mc_state_t state = NULL, prev_state = NULL, next_state = NULL, restored_state=NULL;
   smx_process_t process = NULL;
   xbt_fifo_item_t item = NULL;
-  int d;
+  mc_state_t state_test = NULL;
+  int pos;
+  int visited_state = -1;
+  int enabled = 0;
+  int interleave_size = 0;
+  int comm_pattern = 0;
+  smx_action_t current_comm;
 
+  while (xbt_fifo_size(mc_stack_safety) > 0) {
 
-  while (xbt_fifo_size(mc_snapshot_stack) > 0) {
-
-    /* Get current pair */
-    pair = (mc_pair_prop_t) 
-      xbt_fifo_get_item_content(xbt_fifo_get_first_item(mc_snapshot_stack));
+    /* Get current state */
+    state = (mc_state_t) 
+      xbt_fifo_get_item_content(xbt_fifo_get_first_item(mc_stack_safety));
 
     XBT_DEBUG("**************************************************");
-    XBT_DEBUG("Exploration depth=%d (pair=%p) (%d interleave)",
-             xbt_fifo_size(mc_snapshot_stack), pair, MC_state_interleave_size(pair->graph_state));
-
-    if(MC_state_interleave_size(pair->graph_state) == 0){
-    
-      xbt_dynar_reset(enabled_processes);
-
-      MC_SET_RAW_MEM;
-    
-      xbt_swag_foreach(process, simix_global->process_list){
-       if(MC_process_is_enabled(process)){ 
-         xbt_dynar_push(enabled_processes, &process);
-         //XBT_DEBUG("Process : pid=%lu",process->pid);
-       }
-      }
-
-      //XBT_DEBUG("Enabled processes before pop : %lu", xbt_dynar_length(enabled_processes));
-    
-      //XBT_DEBUG("Processes already interleaved : %d", pair->interleave);
-    
-      if(xbt_dynar_length(enabled_processes) > 0){
-       for(d=0;d<pair->interleave;d++){
-         xbt_dynar_shift(enabled_processes, NULL);
-       }
-      }
-
-      //XBT_DEBUG("Enabled processes after pop : %lu", xbt_dynar_length(enabled_processes));
-    
-      if(pair->fully_expanded == 0){
-       if(xbt_dynar_length(enabled_processes) > 0){
-         MC_state_interleave_process(pair->graph_state, xbt_dynar_get_as(enabled_processes, 0, smx_process_t));
-         //XBT_DEBUG("Interleave process");
-       }
-      }
-    
-      MC_UNSET_RAW_MEM;
-
-    }
-
+    XBT_DEBUG("Exploration depth=%d (state=%p, num %d)(%u interleave, user_max_depth %d, first_enabled_state size : %d)",
+              xbt_fifo_size(mc_stack_safety), state, state->num,
+              MC_state_interleave_size(state), user_max_depth_reached, xbt_dict_size(first_enabled_state));
+      
+    interleave_size = MC_state_interleave_size(state);
 
     /* Update statistics */
-    mc_stats_pair->visited_pairs++;
-    //sleep(1);
-
-    /*cursor = 0;
-    mc_prop_ato_t pato;
-    xbt_dynar_foreach(pair->propositions, cursor, pato){
-      XBT_DEBUG("%s : %d", pato->id, pato->value);
-      }*/
-
-    /* Test acceptance pair and acceptance cycle*/
-
-    if(pair->automaton_state->type == 1){
-      if(search_cycle == 0){
-       set_pair_prop_reached(pair);
-       XBT_DEBUG("Acceptance pair : graph=%p, automaton=%p(%s)", pair->graph_state, pair->automaton_state, pair->automaton_state->id);
-      }else{
-       if(reached_prop(pair) == 1){
-         XBT_INFO("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*");
-         XBT_INFO("|             ACCEPTANCE CYCLE            |");
-         XBT_INFO("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*");
-         XBT_INFO("Counter-example that violates formula :");
-         MC_show_snapshot_stack(mc_snapshot_stack);
-         MC_dump_snapshot_stack(mc_snapshot_stack);
-         MC_print_statistics_pairs(mc_stats_pair);
-         exit(0);
-       }
-      }
-    }
+    mc_stats->visited_states++;
 
     /* If there are processes to interleave and the maximum depth has not been reached
        then perform one step of the exploration algorithm */
-    if (xbt_fifo_size(mc_snapshot_stack) < MAX_DEPTH &&
-        (req = MC_state_get_request(pair->graph_state, &value))) {
+    if (xbt_fifo_size(mc_stack_safety) <= _sg_mc_max_depth && !user_max_depth_reached &&
+        (req = MC_state_get_request(state, &value)) && visited_state == -1) {
 
       /* Debug information */
       if(XBT_LOG_ISENABLED(mc_dpor, xbt_log_priority_debug)){
@@ -608,565 +531,305 @@ void MC_dpor2(xbt_automaton_t a, int search_cycle)
         XBT_DEBUG("Execute: %s", req_str);
         xbt_free(req_str);
       }
+      
+      MC_SET_MC_HEAP;
+      if(dot_output != NULL)
+        req_str = MC_request_get_dot_output(req, value);
+      MC_SET_STD_HEAP;
 
-      MC_state_set_executed_request(pair->graph_state, req, value);
-      //mc_stats_pairs->executed_transitions++;
-
-      /* Answer the request */
-      SIMIX_request_pre(req, value); /* After this call req is no longer usefull */
-
-      /* Wait for requests (schedules processes) */
-      MC_wait_for_requests();
-
-      /* Create the new expanded state */
-      MC_SET_RAW_MEM;
-
-      pair->interleave++;
-      //XBT_DEBUG("Process interleaved in pair %p : %u", pair, MC_state_interleave_size(pair->graph_state));
-      //xbt_dynar_pop(enabled_processes, NULL);
-
-      next_graph_state = MC_state_new();
-
-      next_system_state = xbt_new0(s_mc_snapshot_t, 1);
-      MC_take_snapshot(next_system_state);
-      MC_UNSET_RAW_MEM;
+      MC_state_set_executed_request(state, req, value);
+      mc_stats->executed_transitions++;
 
-      cursor = 0;
-      xbt_dynar_foreach(pair->automaton_state->out, cursor, transition_succ){
-       res = MC_automaton_evaluate_label(a, transition_succ->label);   
-       if((res == 1) || (res == 2)){ // enabled transition in automaton
-         next_automaton_state = transition_succ->dst;
-         break;
-       }
+      if(mc_reduce_kind ==  e_mc_reduce_dpor){
+        MC_SET_MC_HEAP;
+        char *key = bprintf("%lu", req->issuer->pid);
+        xbt_dict_remove(first_enabled_state, key); 
+        xbt_free(key);
+        MC_SET_STD_HEAP;
       }
 
-      if(next_automaton_state == NULL){
-       next_automaton_state = pair->automaton_state;
+      /* TODO : handle test and testany simcalls */
+      if(_sg_mc_comms_determinism || _sg_mc_send_determinism){
+        if(req->call == SIMCALL_COMM_ISEND)
+          comm_pattern = 1;
+        else if(req->call == SIMCALL_COMM_IRECV)
+          comm_pattern = 2;
+        else if(req->call == SIMCALL_COMM_WAIT)
+          comm_pattern = 3;
+        else if(req->call == SIMCALL_COMM_WAITANY)
+          comm_pattern = 4;
       }
 
-      MC_SET_RAW_MEM;
-      
-      next_pair = new_pair_prop(next_system_state, next_graph_state, next_automaton_state);
-      cursor = 0;
-      xbt_propositional_symbol_t ps;
-      xbt_dynar_foreach(a->propositional_symbols, cursor, ps){
-       int (*f)() = ps->function;
-       int val = (*f)();
-       mc_prop_ato_t pa = new_proposition(ps->pred, val);
-       xbt_dynar_push(next_pair->propositions, &pa);
-       //XBT_DEBUG("%s : %d", pa->id, pa->value); 
-      } 
-      xbt_fifo_unshift(mc_snapshot_stack, next_pair);
-      
-      cursor = 0;
-      if((invisible(pair, next_pair) == 0) && (pair->fully_expanded == 0)){
-       set_fully_expanded(pair);
-       if(xbt_dynar_length(enabled_processes) > 1){ // Si 1 seul process activé, déjà exécuté juste avant donc état complètement étudié
-         xbt_dynar_foreach(enabled_processes, cursor, process){
-           MC_state_interleave_process(pair->graph_state, process);
-         }
-       }
-       XBT_DEBUG("Pair %p fully expanded (%u interleave)", pair, MC_state_interleave_size(pair->graph_state));
+      /* Answer the request */
+      SIMIX_simcall_pre(req, value); /* After this call req is no longer usefull */
+
+      if(_sg_mc_comms_determinism || _sg_mc_send_determinism){
+        MC_SET_MC_HEAP;
+        if(comm_pattern == 1 || comm_pattern == 2){
+          if(!initial_state_safety->initial_communications_pattern_done)
+            get_comm_pattern(initial_communications_pattern, req, comm_pattern);
+          else
+            get_comm_pattern(communications_pattern, req, comm_pattern);
+        }else if(comm_pattern == 3){
+          current_comm = simcall_comm_wait__get__comm(req);
+          if(current_comm->comm.refcount == 1){ /* First wait only must be considered */
+            if(!initial_state_safety->initial_communications_pattern_done)
+              complete_comm_pattern(initial_communications_pattern, current_comm);
+            else
+              complete_comm_pattern(communications_pattern, current_comm);
+          }
+        }else if(comm_pattern == 4){
+          current_comm = xbt_dynar_get_as(simcall_comm_waitany__get__comms(req), value, smx_action_t);
+          if(current_comm->comm.refcount == 1){ /* First wait only must be considered */
+            if(!initial_state_safety->initial_communications_pattern_done)
+              complete_comm_pattern(initial_communications_pattern, current_comm);
+            else
+              complete_comm_pattern(communications_pattern, current_comm);
+          }
+        }
+        MC_SET_STD_HEAP;
+        comm_pattern = 0;
       }
 
-      MC_UNSET_RAW_MEM;
-
-     
-      /* Let's loop again */
-
-      /* The interleave set is empty or the maximum depth is reached, let's back-track */
-    } else {
-      XBT_DEBUG("There are no more processes to interleave.");
-
-      /* Trash the current state, no longer needed */
-      MC_SET_RAW_MEM;
-      xbt_fifo_shift(mc_snapshot_stack);
-      //MC_state_delete(state);
-      MC_UNSET_RAW_MEM;
+      /* Wait for requests (schedules processes) */
+      MC_wait_for_requests();
 
-      /* Check for deadlocks */
-      if(MC_deadlock_check()){
-        MC_show_deadlock(NULL);
-        return;
-      }
+      /* Create the new expanded state */
+      MC_SET_MC_HEAP;
 
-      MC_SET_RAW_MEM;
-      /* Traverse the stack backwards until a state with a non empty interleave
-         set is found, deleting all the states that have it empty in the way.
-         For each deleted state, check if the request that has generated it 
-         (from it's predecesor state), depends on any other previous request 
-         executed before it. If it does then add it to the interleave set of the
-         state that executed that previous request. */
-      while ((pair = xbt_fifo_shift(mc_snapshot_stack)) != NULL) {
-        req = MC_state_get_internal_request(pair->graph_state);
-        xbt_fifo_foreach(mc_snapshot_stack, item, prev_pair, mc_pair_prop_t) {
-          if(MC_request_depend(req, MC_state_get_internal_request(prev_pair->graph_state))){
-            if(XBT_LOG_ISENABLED(mc_dpor, xbt_log_priority_debug)){
-              XBT_DEBUG("Dependent Transitions:");
-              prev_req = MC_state_get_executed_request(prev_pair->graph_state, &value);
-              req_str = MC_request_to_string(prev_req, value);
-              XBT_DEBUG("%s (pair=%p)", req_str, prev_pair);
-              xbt_free(req_str);
-              prev_req = MC_state_get_executed_request(pair->graph_state, &value);
-              req_str = MC_request_to_string(prev_req, value);
-              XBT_DEBUG("%s (pair=%p)", req_str, pair);
-              xbt_free(req_str);              
-            }
+      next_state = MC_state_new();
 
-           if(prev_pair->fully_expanded == 0){
-             if(!MC_state_process_is_done(prev_pair->graph_state, req->issuer)){
-               MC_state_interleave_process(prev_pair->graph_state, req->issuer);
-               XBT_DEBUG("Process %p (%lu) interleaved in pair %p", req->issuer, req->issuer->pid, prev_pair);
-             }else{
-               XBT_DEBUG("Process %p (%lu) is in done set", req->issuer, req->issuer->pid);
-             }
-           }
-      
-            break;
+      if((visited_state = is_visited_state()) == -1){
+     
+        /* Get an enabled process and insert it in the interleave set of the next state */
+        xbt_swag_foreach(process, simix_global->process_list){
+          if(MC_process_is_enabled(process)){
+            MC_state_interleave_process(next_state, process);
+            if(mc_reduce_kind != e_mc_reduce_none)
+              break;
           }
         }
 
-       
-        if (MC_state_interleave_size(pair->graph_state) > 0) {
-          /* We found a back-tracking point, let's loop */
-         MC_restore_snapshot(pair->system_state);
-          xbt_fifo_unshift(mc_snapshot_stack, pair);
-          XBT_DEBUG("Back-tracking to depth %d", xbt_fifo_size(mc_snapshot_stack));
-          MC_UNSET_RAW_MEM;
-          break;
-        } else {
-          //MC_state_delete(state);
+        if(_sg_mc_checkpoint && ((xbt_fifo_size(mc_stack_safety) + 1) % _sg_mc_checkpoint == 0)){
+          next_state->system_state = MC_take_snapshot(next_state->num);
         }
-      }
-      MC_UNSET_RAW_MEM;
-    }
-  }
-  MC_UNSET_RAW_MEM;
-  return;
-}
-
-/************ DPOR 3 (invisible and independant transitions with coloration of pairs) ************/
-
-int expanded;
-xbt_dynar_t reached_pairs_prop_col;
-xbt_dynar_t visited_pairs_prop_col;
-
-
-
-mc_pair_prop_col_t new_pair_prop_col(mc_snapshot_t sn, mc_state_t sg, xbt_state_t st){
-  
-  mc_pair_prop_col_t pair = NULL;
-  pair = xbt_new0(s_mc_pair_prop_col_t, 1);
-  pair->system_state = sn;
-  pair->automaton_state = st;
-  pair->graph_state = sg;
-  pair->propositions = xbt_dynar_new(sizeof(mc_prop_ato_t), NULL);
-  pair->fully_expanded = 0;
-  pair->interleave = 0;
-  pair->color=ORANGE;
-  mc_stats_pair->expanded_pairs++;
-  //XBT_DEBUG("New pair %p : automaton=%p", pair, st);
-  return pair;
-}
-
-void set_expanded(mc_pair_prop_col_t pair){
-  pair->expanded = expanded;
-}
-
-int reached_prop_col(mc_pair_prop_col_t pair){
-
-  char* hash_reached = malloc(sizeof(char)*160);
-  unsigned int c= 0;
-
-  MC_SET_RAW_MEM;
-  char *hash = malloc(sizeof(char)*160);
-  xbt_sha((char *)&pair, hash);
-  xbt_dynar_foreach(reached_pairs_prop_col, c, hash_reached){
-    if(strcmp(hash, hash_reached) == 0){
-      MC_UNSET_RAW_MEM;
-      return 1;
-    }
-  }
-  
-  MC_UNSET_RAW_MEM;
-  return 0;
-}
-
-void set_pair_prop_col_reached(mc_pair_prop_col_t pair){
-
-  if(reached_prop_col(pair) == 0){
-    MC_SET_RAW_MEM;
-    char *hash = malloc(sizeof(char)*160) ;
-    xbt_sha((char *)&pair, hash);
-    xbt_dynar_push(reached_pairs_prop_col, &hash); 
-    MC_UNSET_RAW_MEM;
-  }
-
-}
-
-int invisible_col(mc_pair_prop_col_t p, mc_pair_prop_col_t np){
-  mc_prop_ato_t prop1;
-  mc_prop_ato_t prop2;
-  int i;
-  for(i=0;i<xbt_dynar_length(p->propositions); i++){
-    prop1 = xbt_dynar_get_as(p->propositions, i, mc_prop_ato_t);
-    prop2 = xbt_dynar_get_as(np->propositions, i, mc_prop_ato_t);
-    //XBT_DEBUG("test invisible : prop 1 = %s : %d / prop 2 = %s : %d", prop1->id, prop1->value, prop2->id, prop2->value);
-    if(prop1->value != prop2->value)
-         return 0;
-  }
-  return 1; 
-
-}
-
-void set_fully_expanded_col(mc_pair_prop_col_t pair){
-  pair->fully_expanded = 1;
-  pair->color = GREEN;
-}
-
-void set_pair_prop_col_visited(mc_pair_prop_col_t pair, int sc){
-
-  MC_SET_RAW_MEM;
-  mc_visited_pair_prop_col_t p = NULL;
-  p = xbt_new0(s_mc_visited_pair_prop_col_t, 1);
-  p->pair = pair;
-  p->search_cycle = sc;
-  char *hash = malloc(sizeof(char)*160);
-  xbt_sha((char *)&p, hash);
-  xbt_dynar_push(visited_pairs_prop_col, &hash); 
-  MC_UNSET_RAW_MEM;    
 
-}
-
-int visited_pair_prop_col(mc_pair_prop_col_t pair, int sc){
-
-  char* hash_visited = malloc(sizeof(char)*160);
-  unsigned int c= 0;
-
-  MC_SET_RAW_MEM;
-  mc_visited_pair_prop_col_t p = NULL;
-  p = xbt_new0(s_mc_visited_pair_prop_col_t, 1);
-  p->pair = pair;
-  p->search_cycle = sc;
-  char *hash = malloc(sizeof(char)*160);
-  xbt_sha((char *)&p, hash);
-  xbt_dynar_foreach(visited_pairs_prop_col, c, hash_visited){
-    if(strcmp(hash, hash_visited) == 0){
-      MC_UNSET_RAW_MEM;
-      return 1;
-    }
-  }
-  
-  MC_UNSET_RAW_MEM;
-  return 0;
-
-}
-
-void MC_dpor3_init(xbt_automaton_t a)
-{
-  mc_pair_prop_col_t initial_pair = NULL;
-  mc_state_t initial_graph_state = NULL;
-  mc_snapshot_t initial_system_state = NULL;
-  xbt_state_t initial_automaton_state = NULL;
-
-  MC_SET_RAW_MEM;
-  initial_graph_state = MC_state_new();
-  MC_UNSET_RAW_MEM;
-
-  /* Wait for requests (schedules processes) */
-  MC_wait_for_requests();
-
-  unsigned int cursor = 0;
-  unsigned int cursor2 = 0;
-  xbt_state_t automaton_state = NULL;
-  int res;
-  xbt_transition_t transition_succ;
-
-  xbt_dynar_foreach(a->states, cursor, automaton_state){
-    if(automaton_state->type == -1){
-      xbt_dynar_foreach(automaton_state->out, cursor2, transition_succ){
-       res = MC_automaton_evaluate_label(a, transition_succ->label);
-       if((res == 1) || (res == 2)){
-         initial_automaton_state = transition_succ->dst;
-         break;
-       }
-      }
-    }
-
-    if(initial_automaton_state != NULL)
-      break;
-  }
-
-  if(initial_automaton_state == NULL){
-    cursor = 0;
-    xbt_dynar_foreach(a->states, cursor, automaton_state){
-      if(automaton_state->type == -1){
-       initial_automaton_state = automaton_state;
-       break;
-      }
-    }
-  }
-
-  reached_pairs_prop_col = xbt_dynar_new(sizeof(char *), NULL); 
-  visited_pairs_prop_col = xbt_dynar_new(sizeof(char *), NULL); 
-  visible_transitions = xbt_dynar_new(sizeof(s_smx_req_t), NULL);
-  enabled_processes = xbt_dynar_new(sizeof(smx_process_t), NULL);
-  expanded = 0;
-
-  MC_SET_RAW_MEM;
-  initial_system_state = xbt_new0(s_mc_snapshot_t, 1);
-  MC_take_snapshot(initial_system_state);
-  initial_pair = new_pair_prop_col(initial_system_state, initial_graph_state, initial_automaton_state);
-  cursor = 0;
-  xbt_propositional_symbol_t ps;
-  xbt_dynar_foreach(a->propositional_symbols, cursor, ps){
-    int (*f)() = ps->function;
-    int val = (*f)();
-    mc_prop_ato_t pa = new_proposition(ps->pred, val);
-    xbt_dynar_push(initial_pair->propositions, &pa);
-  } 
-  xbt_fifo_unshift(mc_snapshot_stack, initial_pair);
-  MC_UNSET_RAW_MEM;
-
-
-  XBT_DEBUG("**************************************************");
-  XBT_DEBUG("Initial pair");
-
-  MC_dpor3(a, 0);
-    
-}
-
-
-void MC_dpor3(xbt_automaton_t a, int search_cycle)
-{
-  char *req_str;
-  int value;
-  smx_req_t req = NULL, prev_req = NULL;
-  mc_state_t next_graph_state = NULL;
-  mc_snapshot_t next_system_state = NULL;
-  xbt_state_t next_automaton_state = NULL;
-  xbt_transition_t transition_succ;
-  unsigned int cursor;
-  int res;
-  mc_pair_prop_col_t pair = NULL, next_pair = NULL, prev_pair = NULL;
-  smx_process_t process = NULL;
-  xbt_fifo_item_t item = NULL;
-  int d;
-
-
-  while (xbt_fifo_size(mc_snapshot_stack) > 0) {
-
-    /* Get current pair */
-    pair = (mc_pair_prop_col_t) 
-      xbt_fifo_get_item_content(xbt_fifo_get_first_item(mc_snapshot_stack));
-
-    XBT_DEBUG("**************************************************");
-    XBT_DEBUG("Exploration depth=%d (pair=%p) (%d interleave)",
-             xbt_fifo_size(mc_snapshot_stack), pair, MC_state_interleave_size(pair->graph_state));
+        if(dot_output != NULL)
+          fprintf(dot_output, "\"%d\" -> \"%d\" [%s];\n", state->num, next_state->num, req_str);
 
-    if(MC_state_interleave_size(pair->graph_state) == 0){
-    
-      xbt_dynar_reset(enabled_processes);
-
-      MC_SET_RAW_MEM;
-    
-      xbt_swag_foreach(process, simix_global->process_list){
-       if(MC_process_is_enabled(process)){ 
-         xbt_dynar_push(enabled_processes, &process);
-         //XBT_DEBUG("Process : pid=%lu",process->pid);
-       }
-      }
+      }else{
 
-      //XBT_DEBUG("Enabled processes before pop : %lu", xbt_dynar_length(enabled_processes));
-    
-      //XBT_DEBUG("Processes already interleaved : %d", pair->interleave);
-    
-      if(xbt_dynar_length(enabled_processes) > 0){
-       for(d=0;d<pair->interleave;d++){
-         xbt_dynar_shift(enabled_processes, NULL);
-       }
-      }
+        if(dot_output != NULL)
+          fprintf(dot_output, "\"%d\" -> \"%d\" [%s];\n", state->num, visited_state, req_str);
 
-      //XBT_DEBUG("Enabled processes after pop : %lu", xbt_dynar_length(enabled_processes));
-    
-      if(pair->fully_expanded == 0){
-       if(xbt_dynar_length(enabled_processes) > 0){
-         MC_state_interleave_process(pair->graph_state, xbt_dynar_get_as(enabled_processes, 0, smx_process_t));
-         //XBT_DEBUG("Interleave process");
-       }
       }
-    
-      MC_UNSET_RAW_MEM;
 
-    }
-
-
-    /* Update statistics */
-    mc_stats_pair->visited_pairs++;
-    //sleep(1);
+      xbt_fifo_unshift(mc_stack_safety, next_state);
 
-    /* Test acceptance pair and acceptance cycle*/
-
-    if(pair->automaton_state->type == 1){
-      if(search_cycle == 0){
-       set_pair_prop_col_reached(pair);
-       XBT_DEBUG("Acceptance pair : graph=%p, automaton=%p(%s)", pair->graph_state, pair->automaton_state, pair->automaton_state->id);
-      }else{
-       if(reached_prop_col(pair) == 1){
-         XBT_INFO("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*");
-         XBT_INFO("|             ACCEPTANCE CYCLE            |");
-         XBT_INFO("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*");
-         XBT_INFO("Counter-example that violates formula :");
-         MC_show_snapshot_stack(mc_snapshot_stack);
-         MC_dump_snapshot_stack(mc_snapshot_stack);
-         MC_print_statistics_pairs(mc_stats_pair);
-         exit(0);
-       }
+      if(mc_reduce_kind ==  e_mc_reduce_dpor){
+        /* Insert in dict all enabled processes, if not included yet */
+        xbt_swag_foreach(process, simix_global->process_list){
+          if(MC_process_is_enabled(process)){
+            char *key = bprintf("%lu", process->pid);
+            if(xbt_dict_get_or_null(first_enabled_state, key) == NULL){
+              char *data = bprintf("%d", xbt_fifo_size(mc_stack_safety));
+              xbt_dict_set(first_enabled_state, key, data, NULL); 
+            }
+            xbt_free(key);
+          }
+        }
       }
-    }
-
-    /* If there are processes to interleave and the maximum depth has not been reached
-       then perform one step of the exploration algorithm */
-    if (xbt_fifo_size(mc_snapshot_stack) < MAX_DEPTH &&
-        (req = MC_state_get_request(pair->graph_state, &value))) {
-
-      set_pair_prop_col_visited(pair, search_cycle);
       
-      /* Debug information */
-      if(XBT_LOG_ISENABLED(mc_dpor, xbt_log_priority_debug)){
-        req_str = MC_request_to_string(req, value);
-        XBT_DEBUG("Execute: %s", req_str);
+      if(dot_output != NULL)
         xbt_free(req_str);
-      }
 
-      MC_state_set_executed_request(pair->graph_state, req, value);
+      MC_SET_STD_HEAP;
 
-      /* Answer the request */
-      SIMIX_request_pre(req, value); /* After this call req is no longer usefull */
-
-      /* Wait for requests (schedules processes) */
-      MC_wait_for_requests();
+      /* Let's loop again */
 
-      /* Create the new expanded state */
-      MC_SET_RAW_MEM;
+      /* The interleave set is empty or the maximum depth is reached, let's back-track */
+    } else {
 
-      pair->interleave++;
+      if((xbt_fifo_size(mc_stack_safety) > _sg_mc_max_depth) || user_max_depth_reached || visited_state != -1){  
+
+        if(user_max_depth_reached && visited_state == -1)
+          XBT_DEBUG("User max depth reached !");
+        else if(visited_state == -1)
+          XBT_WARN("/!\\ Max depth reached ! /!\\ ");
+
+        visited_state = -1;
+
+        if(mc_reduce_kind ==  e_mc_reduce_dpor){
+          /* Interleave enabled processes in the state in which they have been enabled for the first time */
+          xbt_swag_foreach(process, simix_global->process_list){
+            if(MC_process_is_enabled(process)){
+              char *key = bprintf("%lu", process->pid);
+              enabled = (int)strtoul(xbt_dict_get_or_null(first_enabled_state, key), 0, 10);
+              xbt_free(key);
+              int cursor = xbt_fifo_size(mc_stack_safety);
+              xbt_fifo_foreach(mc_stack_safety, item, state_test, mc_state_t){
+                if(cursor-- == enabled){ 
+                  if(!MC_state_process_is_done(state_test, process) && state_test->num != state->num){ 
+                    XBT_DEBUG("Interleave process %lu in state %d", process->pid, state_test->num);
+                    MC_state_interleave_process(state_test, process);
+                    break;
+                  }
+                }
+              }
+            }
+          }
+        }
 
-      next_graph_state = MC_state_new();
+      }else{
 
-      next_system_state = xbt_new0(s_mc_snapshot_t, 1);
-      MC_take_snapshot(next_system_state);
-      MC_UNSET_RAW_MEM;
+        XBT_DEBUG("There are no more processes to interleave. (depth %d)", xbt_fifo_size(mc_stack_safety) + 1);
 
-      cursor = 0;
-      xbt_dynar_foreach(pair->automaton_state->out, cursor, transition_succ){
-       res = MC_automaton_evaluate_label(a, transition_succ->label);   
-       if((res == 1) || (res == 2)){ // enabled transition in automaton
-         next_automaton_state = transition_succ->dst;
-         break;
-       }
       }
 
-      if(next_automaton_state == NULL)
-       next_automaton_state = pair->automaton_state;
-
-      MC_SET_RAW_MEM;
-      
-      next_pair = new_pair_prop_col(next_system_state, next_graph_state, next_automaton_state);
-      cursor = 0;
-      xbt_propositional_symbol_t ps;
-      xbt_dynar_foreach(a->propositional_symbols, cursor, ps){
-       int (*f)() = ps->function;
-       int val = (*f)();
-       mc_prop_ato_t pa = new_proposition(ps->pred, val);
-       xbt_dynar_push(next_pair->propositions, &pa);
-      } 
-      xbt_fifo_unshift(mc_snapshot_stack, next_pair);
-      
-      cursor = 0;
-      if((invisible_col(pair, next_pair) == 0) && (pair->fully_expanded == 0)){
-       set_fully_expanded_col(pair);
-       if(xbt_dynar_length(enabled_processes) > 1){ // Si 1 seul process activé, déjà exécuté juste avant donc état complètement étudié
-         xbt_dynar_foreach(enabled_processes, cursor, process){
-           MC_state_interleave_process(pair->graph_state, process);
-         }
-       }
-       XBT_DEBUG("Pair %p fully expanded (%u interleave)", pair, MC_state_interleave_size(pair->graph_state));
+      MC_SET_MC_HEAP;
+
+      if(_sg_mc_comms_determinism || _sg_mc_send_determinism){
+        if(initial_state_safety->initial_communications_pattern_done){
+          if(interleave_size == 0){ /* if (interleave_size > 0), process interleaved but not enabled => "incorrect" path, determinism not evaluated */
+            deterministic_pattern(initial_communications_pattern, communications_pattern);
+            if(initial_state_safety->comm_deterministic == 0 && _sg_mc_comms_determinism){
+              XBT_INFO("****************************************************");
+              XBT_INFO("***** Non-deterministic communications pattern *****");
+              XBT_INFO("****************************************************");
+              XBT_INFO("Initial communications pattern:");
+              print_communications_pattern(initial_communications_pattern);
+              XBT_INFO("Communications pattern counter-example:");
+              print_communications_pattern(communications_pattern);
+              MC_print_statistics(mc_stats);
+              return;
+            }else if(initial_state_safety->send_deterministic == 0 && _sg_mc_send_determinism){
+              XBT_INFO("*********************************************************");
+              XBT_INFO("***** Non-send-deterministic communications pattern *****");
+              XBT_INFO("*********************************************************");
+              XBT_INFO("Initial communications pattern:");
+              print_communications_pattern(initial_communications_pattern);
+              XBT_INFO("Communications pattern counter-example:");
+              print_communications_pattern(communications_pattern);
+              MC_print_statistics(mc_stats);
+              return;
+            }
+          }
+        }else{
+          initial_state_safety->initial_communications_pattern_done = 1;
+        }
       }
 
-      MC_UNSET_RAW_MEM;
-
-     
-      /* Let's loop again */
-
-      /* The interleave set is empty or the maximum depth is reached, let's back-track */
-    } else {
-      XBT_DEBUG("There are no more processes to interleave.");
-
       /* Trash the current state, no longer needed */
-      MC_SET_RAW_MEM;
-      xbt_fifo_shift(mc_snapshot_stack);
-      //MC_state_delete(state);
-      MC_UNSET_RAW_MEM;
+      xbt_fifo_shift(mc_stack_safety);
+      MC_state_delete(state);
+      XBT_DEBUG("Delete state %d at depth %d", state->num, xbt_fifo_size(mc_stack_safety) + 1);
 
+      MC_SET_STD_HEAP;
+      
       /* Check for deadlocks */
       if(MC_deadlock_check()){
         MC_show_deadlock(NULL);
         return;
       }
 
-      MC_SET_RAW_MEM;
+      MC_SET_MC_HEAP;
       /* Traverse the stack backwards until a state with a non empty interleave
          set is found, deleting all the states that have it empty in the way.
          For each deleted state, check if the request that has generated it 
          (from it's predecesor state), depends on any other previous request 
          executed before it. If it does then add it to the interleave set of the
          state that executed that previous request. */
-      while ((pair = xbt_fifo_shift(mc_snapshot_stack)) != NULL) {
-        req = MC_state_get_internal_request(pair->graph_state);
-        xbt_fifo_foreach(mc_snapshot_stack, item, prev_pair, mc_pair_prop_col_t) {
-          if(MC_request_depend(req, MC_state_get_internal_request(prev_pair->graph_state))){
-            if(XBT_LOG_ISENABLED(mc_dpor, xbt_log_priority_debug)){
-              XBT_DEBUG("Dependent Transitions:");
-              prev_req = MC_state_get_executed_request(prev_pair->graph_state, &value);
-              req_str = MC_request_to_string(prev_req, value);
-              XBT_DEBUG("%s (pair=%p)", req_str, prev_pair);
-              xbt_free(req_str);
-              prev_req = MC_state_get_executed_request(pair->graph_state, &value);
-              req_str = MC_request_to_string(prev_req, value);
-              XBT_DEBUG("%s (pair=%p)", req_str, pair);
-              xbt_free(req_str);              
-            }
-
-           if(prev_pair->fully_expanded == 0){
-             if(!MC_state_process_is_done(prev_pair->graph_state, req->issuer)){
-               MC_state_interleave_process(prev_pair->graph_state, req->issuer);
-               XBT_DEBUG("Process %p (%lu) interleaved in pair %p", req->issuer, req->issuer->pid, prev_pair);
-             }else{
-               XBT_DEBUG("Process %p (%lu) is in done set", req->issuer, req->issuer->pid);
-             }
-           }
       
-            break;
+      while ((state = xbt_fifo_shift(mc_stack_safety)) != NULL) {
+        if(mc_reduce_kind == e_mc_reduce_dpor){
+          req = MC_state_get_internal_request(state);
+          xbt_fifo_foreach(mc_stack_safety, item, prev_state, mc_state_t) {
+            if(MC_request_depend(req, MC_state_get_internal_request(prev_state))){
+              if(XBT_LOG_ISENABLED(mc_dpor, xbt_log_priority_debug)){
+                XBT_DEBUG("Dependent Transitions:");
+                prev_req = MC_state_get_executed_request(prev_state, &value);
+                req_str = MC_request_to_string(prev_req, value);
+                XBT_DEBUG("%s (state=%d)", req_str, prev_state->num);
+                xbt_free(req_str);
+                prev_req = MC_state_get_executed_request(state, &value);
+                req_str = MC_request_to_string(prev_req, value);
+                XBT_DEBUG("%s (state=%d)", req_str, state->num);
+                xbt_free(req_str);              
+              }
+
+              if(!MC_state_process_is_done(prev_state, req->issuer))
+                MC_state_interleave_process(prev_state, req->issuer);
+              else
+                XBT_DEBUG("Process %p is in done set", req->issuer);
+
+              break;
+
+            }else if(req->issuer == MC_state_get_internal_request(prev_state)->issuer){
+
+              XBT_DEBUG("Simcall %d and %d with same issuer", req->call, MC_state_get_internal_request(prev_state)->call);
+              break;
+
+            }else{
+
+              XBT_DEBUG("Simcall %d, process %lu (state %d) and simcall %d, process %lu (state %d) are independant", req->call, req->issuer->pid, state->num, MC_state_get_internal_request(prev_state)->call, MC_state_get_internal_request(prev_state)->issuer->pid, prev_state->num);
+
+            }
           }
         }
-
-       
-        if (MC_state_interleave_size(pair->graph_state) > 0) {
+             
+        if (MC_state_interleave_size(state) && xbt_fifo_size(mc_stack_safety) < _sg_mc_max_depth) {
           /* We found a back-tracking point, let's loop */
-         MC_restore_snapshot(pair->system_state);
-          xbt_fifo_unshift(mc_snapshot_stack, pair);
-          XBT_DEBUG("Back-tracking to depth %d", xbt_fifo_size(mc_snapshot_stack));
-          MC_UNSET_RAW_MEM;
+          XBT_DEBUG("Back-tracking to state %d at depth %d", state->num, xbt_fifo_size(mc_stack_safety) + 1);
+          if(_sg_mc_checkpoint){
+            if(state->system_state != NULL){
+              MC_restore_snapshot(state->system_state);
+              xbt_fifo_unshift(mc_stack_safety, state);
+              MC_SET_STD_HEAP;
+            }else{
+              pos = xbt_fifo_size(mc_stack_safety);
+              item = xbt_fifo_get_first_item(mc_stack_safety);
+              while(pos>0){
+                restored_state = (mc_state_t) xbt_fifo_get_item_content(item);
+                if(restored_state->system_state != NULL){
+                  break;
+                }else{
+                  item = xbt_fifo_get_next_item(item);
+                  pos--;
+                }
+              }
+              MC_restore_snapshot(restored_state->system_state);
+              xbt_fifo_unshift(mc_stack_safety, state);
+              MC_SET_STD_HEAP;
+              MC_replay(mc_stack_safety, pos);
+            }
+          }else{
+            xbt_fifo_unshift(mc_stack_safety, state);
+            MC_SET_STD_HEAP;
+            MC_replay(mc_stack_safety, -1);
+          }
+          XBT_DEBUG("Back-tracking to state %d at depth %d done", state->num, xbt_fifo_size(mc_stack_safety));
           break;
         } else {
-          //MC_state_delete(state);
+          /*req = MC_state_get_internal_request(state);
+          if(_sg_mc_comms_determinism || _sg_mc_send_determinism){
+            if(req->call == SIMCALL_COMM_ISEND || req->call == SIMCALL_COMM_IRECV){
+              if(!xbt_dynar_is_empty(communications_pattern))
+                xbt_dynar_remove_at(communications_pattern, xbt_dynar_length(communications_pattern) - 1, NULL);
+            }
+            }*/
+          XBT_DEBUG("Delete state %d at depth %d", state->num, xbt_fifo_size(mc_stack_safety) + 1); 
+          MC_state_delete(state);
         }
       }
-      MC_UNSET_RAW_MEM;
+      MC_SET_STD_HEAP;
     }
   }
-  MC_UNSET_RAW_MEM;
+  MC_print_statistics(mc_stats);
+  MC_SET_STD_HEAP;
+
   return;
 }
+
+
+
+