Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
[mc] mc_comm_pattern.c containing all communication pattern functions
[simgrid.git] / src / mc / mc_liveness.c
index 9ea615f..4d7c048 100644 (file)
-/* Copyright (c) 2008-2012 Da SimGrid Team. All rights reserved.            */
+/* Copyright (c) 2011-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 "mc_private.h"
 #include <unistd.h>
+#include <sys/wait.h>
 
-XBT_LOG_NEW_DEFAULT_SUBCATEGORY(mc_liveness, mc,
-                                "Logging specific to algorithms for liveness properties verification");
-
-xbt_dynar_t reached_pairs;
-xbt_dynar_t reached_pairs_hash;
-xbt_dynar_t visited_pairs;
-xbt_dynar_t visited_pairs_hash;
-xbt_dynar_t successors;
-
-xbt_dynar_t hosts_table;
-
-
-/* fast implementation of djb2 algorithm */
-unsigned int hash_region(char *str, int str_len){
-
-  int c;
-  register unsigned int hash = 5381;
-
-  while (str_len--) {
-    c = *str++;
-    hash = ((hash << 5) + hash) + c;    /* hash * 33 + c */
-  }
+#include <xbt/dynar.h>
+#include <xbt/automaton.h>
 
-  return hash;
-
-}
-
-const char* get_memory_map_addr(void *addr){
+#include "mc_request.h"
+#include "mc_liveness.h"
+#include "mc_private.h"
+#include "mc_record.h"
+#include "mc_smx.h"
 
-  FILE *fp;                     /* File pointer to process's proc maps file */
-  char *line = NULL;            /* Temporal storage for each line that is readed */
-  ssize_t read;                 /* Number of bytes readed */
-  size_t n = 0;                 /* Amount of bytes to read by getline */
+XBT_LOG_NEW_DEFAULT_SUBCATEGORY(mc_liveness, mc,
+                                "Logging specific to algorithms for liveness properties verification");
 
-  fp = fopen("/proc/self/maps", "r");
-  
-  if(fp == NULL)
-    perror("fopen failed");
+/********* Global variables *********/
 
-  if(addr == NULL){
-    fclose(fp);
-    return "nil";
-  }
+xbt_dynar_t acceptance_pairs;
+xbt_parmap_t parmap;
 
-  xbt_dynar_t lfields = NULL;
-  xbt_dynar_t start_end  = NULL;
-  void *start_addr;
-  void *end_addr;
-
-  while ((read = getline(&line, &n, fp)) != -1) {
-
-    xbt_str_trim(line, NULL);
-    xbt_str_strip_spaces(line);
-    lfields = xbt_str_split(line,NULL);
-
-    start_end = xbt_str_split(xbt_dynar_get_as(lfields, 0, char*), "-");
-    start_addr = (void *) strtoul(xbt_dynar_get_as(start_end, 0, char*), NULL, 16);
-    end_addr = (void *) strtoul(xbt_dynar_get_as(start_end, 1, char*), NULL, 16);
-
-    if((addr > start_addr) && ( addr < end_addr)){
-      free(line);
-      fclose(fp);
-      if(start_addr == std_heap)
-       return "std_heap";
-      if(start_addr == raw_heap)
-       return "raw_heap";
-      if(xbt_dynar_length(lfields) == 6)
-       return xbt_dynar_get_as(lfields, xbt_dynar_length(lfields) - 1, char*);
-      else
-       return "Anonymous";
-    }
+/********* Static functions *********/
 
+static xbt_dynar_t get_atomic_propositions_values()
+{
+  unsigned int cursor = 0;
+  xbt_automaton_propositional_symbol_t ps = NULL;
+  xbt_dynar_t values = xbt_dynar_new(sizeof(int), NULL);
+  xbt_dynar_foreach(_mc_property_automaton->propositional_symbols, cursor, ps) {
+    int res = xbt_automaton_propositional_symbol_evaluate(ps);
+    xbt_dynar_push_as(values, int, res);
   }
 
-  free(line);
-  fclose(fp);
-  return "Unknown area";
-
+  return values;
 }
 
-int data_program_region_compare(void *d1, void *d2, size_t size){
-  int distance = 0;
-  int pointer_align;
-  int i;
-  
-  for(i=0; i<size; i++){
-    if(memcmp(((char *)d1) + i, ((char *)d2) + i, 1) != 0){
-      fprintf(stderr,"Different byte (offset=%d) (%p - %p) in data program region\n", i, (char *)d1 + i, (char *)d2 + i);
-      distance++;
-      pointer_align = (i /sizeof(void *)) * sizeof(void *);
-      fprintf(stderr, "Pointed address : %p (in %s) - %p (in %s)\n", *((void **)((char *)d1 + pointer_align)),  get_memory_map_addr(*((void **)((char *)d1 + pointer_align))), *((void **)((char *)d2 + pointer_align)), get_memory_map_addr(*((void **)((char *)d2 + pointer_align))));
+static mc_visited_pair_t is_reached_acceptance_pair(mc_pair_t pair) {
+  xbt_mheap_t heap = mmalloc_set_current_heap(mc_heap);
+
+  mc_visited_pair_t new_pair = NULL;
+  new_pair = MC_visited_pair_new(pair->num, pair->automaton_state, pair->atomic_propositions, pair->graph_state);
+  new_pair->acceptance_pair = 1;
+
+  if (xbt_dynar_is_empty(acceptance_pairs)) {
+
+    xbt_dynar_push(acceptance_pairs, &new_pair);
+
+  } else {
+
+    int min = -1, max = -1, index;
+    //int res;
+    mc_visited_pair_t pair_test;
+    int cursor;
+
+    index = get_search_interval(acceptance_pairs, new_pair, &min, &max);
+
+    if (min != -1 && max != -1) {       // Acceptance pair with same number of processes and same heap bytes used exists
+
+      // Parallell implementation
+      /*res = xbt_parmap_mc_apply(parmap, snapshot_compare, xbt_dynar_get_ptr(acceptance_pairs, min), (max-min)+1, pair);
+         if(res != -1){
+         if(!raw_mem_set)
+         MC_SET_STD_HEAP;
+         return ((mc_pair_t)xbt_dynar_get_as(acceptance_pairs, (min+res)-1, mc_pair_t))->num;
+         } */
+
+      cursor = min;
+      if(pair->search_cycle == 1){
+        while (cursor <= max) {
+          pair_test = (mc_visited_pair_t) xbt_dynar_get_as(acceptance_pairs, cursor, mc_visited_pair_t);
+          if (xbt_automaton_state_compare(pair_test->automaton_state, new_pair->automaton_state) == 0) {
+            if (xbt_automaton_propositional_symbols_compare_value(pair_test->atomic_propositions, new_pair->atomic_propositions) == 0) {
+              if (snapshot_compare(pair_test, new_pair) == 0) {
+                XBT_INFO("Pair %d already reached (equal to pair %d) !", new_pair->num, pair_test->num);
+                xbt_fifo_shift(mc_stack);
+                if (dot_output != NULL)
+                  fprintf(dot_output, "\"%d\" -> \"%d\" [%s];\n", initial_global_state->prev_pair, pair_test->num, initial_global_state->prev_req);
+                mmalloc_set_current_heap(heap);
+                return NULL;
+              }
+            }
+          }
+          cursor++;
+        }
+      }
+      xbt_dynar_insert_at(acceptance_pairs, min, &new_pair);
+    } else {
+      pair_test = (mc_visited_pair_t) xbt_dynar_get_as(acceptance_pairs, index, mc_visited_pair_t);
+      if (pair_test->nb_processes < new_pair->nb_processes) {
+        xbt_dynar_insert_at(acceptance_pairs, index + 1, &new_pair);
+      } else {
+        if (pair_test->heap_bytes_used < new_pair->heap_bytes_used)
+          xbt_dynar_insert_at(acceptance_pairs, index + 1, &new_pair);
+        else
+          xbt_dynar_insert_at(acceptance_pairs, index, &new_pair);
+      }
     }
-  }
-  
-  fprintf(stderr, "Hamming distance between data program regions : %d\n", distance);
-
-  return distance;
-}
 
-int data_libsimgrid_region_compare(void *d1, void *d2, size_t size){
-  int distance = 0;
-  int pointer_align;
-  int i;
-  
-  for(i=0; i<size; i++){
-    if(memcmp(((char *)d1) + i, ((char *)d2) + i, 1) != 0){
-      fprintf(stderr, "Different byte (offset=%d) (%p - %p) in data libsimgrid region\n", i, (char *)d1 + i, (char *)d2 + i);
-      distance++;
-      pointer_align = (i /sizeof(void *)) * sizeof(void *);
-      fprintf(stderr, "Pointed address : %p (in %s) - %p (in %s)\n", *((void **)((char *)d1 + pointer_align)), get_memory_map_addr(*((void **)((char *)d1 + pointer_align))), *((void **)((char *)d2 + pointer_align)), get_memory_map_addr(*((void **)((char *)d2 + pointer_align))));
-    }
   }
-  
-  fprintf(stderr, "Hamming distance between data libsimgrid regions : %d\n", distance);
-  
-  return distance;
+  mmalloc_set_current_heap(heap);
+  return new_pair;
 }
 
-int snapshot_compare(mc_snapshot_t s1, mc_snapshot_t s2, void* s_heap, void* r_heap){
-
-  
-  if(s1->num_reg != s2->num_reg){
-    XBT_INFO("Different num_reg (s1 = %u, s2 = %u)", s1->num_reg, s2->num_reg);
-    return 1;
-  }
-
-  int i;
-  int errors = 0;
-
-  for(i=0 ; i< s1->num_reg ; i++){
+static void remove_acceptance_pair(int pair_num)
+{
+  xbt_mheap_t heap = mmalloc_set_current_heap(mc_heap);
 
-    if(s1->regions[i]->type != s2->regions[i]->type){
-       XBT_INFO("Different type of region");
-       errors++;
-    }
+  unsigned int cursor = 0;
+  mc_visited_pair_t pair_test = NULL;
+  int pair_found = 0;
 
-    switch(s1->regions[i]->type){
-      case 0:
-      if(s1->regions[i]->size != s2->regions[i]->size){
-       XBT_INFO("Different size of heap (s1 = %zu, s2 = %zu)", s1->regions[i]->size, s2->regions[i]->size);
-       errors++;
-      }
-      if(s1->regions[i]->start_addr != s2->regions[i]->start_addr){
-       XBT_INFO("Different start addr of heap (s1 = %p, s2 = %p)", s1->regions[i]->start_addr, s2->regions[i]->start_addr);
-       errors++;
-      }
-      if(mmalloc_compare_heap(s1->regions[i]->data, s2->regions[i]->data, s_heap, r_heap)){
-       XBT_INFO("Different heap (mmalloc_compare)");
-       errors++; 
-      }
-      break;
-    case 1 :
-      if(s1->regions[i]->size != s2->regions[i]->size){
-       XBT_INFO("Different size of libsimgrid (s1 = %zu, s2 = %zu)", s1->regions[i]->size, s2->regions[i]->size);
-       errors++;
-      }
-      if(s1->regions[i]->start_addr != s2->regions[i]->start_addr){
-       XBT_INFO("Different start addr of libsimgrid (s1 = %p, s2 = %p)", s1->regions[i]->start_addr, s2->regions[i]->start_addr);
-       errors++;
-      }
-      if(data_libsimgrid_region_compare(s1->regions[i]->data, s2->regions[i]->data, s1->regions[i]->size) != 0){
-       XBT_INFO("Different memcmp for data in libsimgrid");
-       errors++;
-      }
-      break;
-    case 2 :
-      if(s1->regions[i]->size != s2->regions[i]->size){
-       XBT_INFO("Different size of data program (s1 = %zu, s2 = %zu)", s1->regions[i]->size, s2->regions[i]->size);
-       errors++;
-      }
-      if(s1->regions[i]->start_addr != s2->regions[i]->start_addr){
-       XBT_INFO("Different start addr of data program (s1 = %p, s2 = %p)", s1->regions[i]->start_addr, s2->regions[i]->start_addr);
-       errors++;
-      }
-      if(data_program_region_compare(s1->regions[i]->data, s2->regions[i]->data, s1->regions[i]->size) != 0){
-       XBT_INFO("Different memcmp for data in program");
-       errors++;
-      }
-      break;
-    default:
+  xbt_dynar_foreach(acceptance_pairs, cursor, pair_test) {
+    if (pair_test->num == pair_num) {
+       pair_found = 1;
       break;
     }
   }
 
-  return (errors > 0);
-  
-}
-
-int reached(xbt_state_t st){
-
+  if(pair_found == 1) {
+    xbt_dynar_remove_at(acceptance_pairs, cursor, &pair_test);
 
-  if(xbt_dynar_is_empty(reached_pairs)){
+    pair_test->acceptance_removed = 1;
 
-    return 0;
+    if (_sg_mc_visited == 0 || pair_test->visited_removed == 1) 
+      MC_visited_pair_delete(pair_test);
 
-  }else{
-
-    MC_SET_RAW_MEM;
-
-    mc_snapshot_t sn = xbt_new0(s_mc_snapshot_t, 1);
-    MC_take_snapshot_liveness(sn);    
-    
-    xbt_dynar_t prop_ato = xbt_dynar_new(sizeof(int), NULL);
-    int res;
-    int_f_void_t f;
-
-    /* Get values of propositional symbols */
-    unsigned int cursor = 0;
-    xbt_propositional_symbol_t ps = NULL;
-    xbt_dynar_foreach(automaton->propositional_symbols, cursor, ps){
-      f = (int_f_void_t)ps->function;
-      res = (*f)();
-      xbt_dynar_push_as(prop_ato, int, res);
-    }
-
-    cursor = 0;
-    mc_pair_reached_t pair_test = NULL;
-
-    //xbt_dict_t current_rdv_points = SIMIX_get_rdv_points();
-     
-    xbt_dynar_foreach(reached_pairs, cursor, pair_test){
-      XBT_INFO("Pair reached #%u", cursor+1);
-      if(automaton_state_compare(pair_test->automaton_state, st) == 0){
-       if(propositional_symbols_compare_value(pair_test->prop_ato, prop_ato) == 0){
-         //XBT_INFO("Rdv points size %d - %d", xbt_dict_length(pair_test->rdv_points), xbt_dict_length(current_rdv_points));
-         //if(xbt_dict_length(pair_test->rdv_points) == xbt_dict_length(current_rdv_points)){
-         //if(rdv_points_compare(pair_test->rdv_points, current_rdv_points) == 0){
-         if(snapshot_compare(pair_test->system_state, sn, std_heap, raw_heap) == 0){
-           MC_free_snapshot(sn);
-           xbt_dynar_reset(prop_ato);
-           xbt_free(prop_ato);
-           MC_UNSET_RAW_MEM;
-           return 1;
-         }
-         /* }
-         }else{
-           XBT_INFO("Different size of rdv points (%d - %d)",xbt_dict_length(pair_test->rdv_points), xbt_dict_length(current_rdv_points) );
-           }*/
-       }else{
-         XBT_INFO("Different values of propositional symbols");
-       }
-      }else{
-       XBT_INFO("Different automaton state");
-      }
-    }
-
-    MC_free_snapshot(sn);
-    xbt_dynar_reset(prop_ato);
-    xbt_free(prop_ato);
-    MC_UNSET_RAW_MEM;
-    return 0;
-    
   }
-}
 
-int rdv_points_compare(xbt_dict_t d1, xbt_dict_t d2){ /* d1 = pair_test, d2 = current_pair */ 
-  
-   xbt_dict_cursor_t cursor_dict = NULL;
-   char *key;
-   char *data;
-   smx_rdv_t rdv1, rdv2;
-   xbt_fifo_item_t item1, item2;
-   smx_action_t action1, action2;
-   xbt_fifo_item_t item_req1, item_req2;
-   smx_simcall_t req1, req2;
-   int i=0;
-   int j=0;
-
-   xbt_dict_foreach(d1, cursor_dict, key, data){
-     rdv1 = (smx_rdv_t)data;
-     rdv2 = xbt_dict_get_or_null(d2, rdv1->name);
-     if(rdv2 == NULL){
-       XBT_INFO("Rdv point unknown");
-       return 1;
-     }else{
-       if(xbt_fifo_size(rdv1->comm_fifo) != xbt_fifo_size(rdv2->comm_fifo)){
-        XBT_INFO("Different total of actions in mailbox \"%s\" (%d - %d)", rdv1->name, xbt_fifo_size(rdv1->comm_fifo),xbt_fifo_size(rdv2->comm_fifo) );
-        return 1;
-       }else{
-        
-        XBT_INFO("Total of actions in mailbox \"%s\" : %d", rdv1->name, xbt_fifo_size(rdv1->comm_fifo)); 
-        
-        item1 = xbt_fifo_get_first_item(rdv1->comm_fifo);      
-        item2 = xbt_fifo_get_first_item(rdv2->comm_fifo);
-
-        while(i<xbt_fifo_size(rdv1->comm_fifo)){
-          action1 = (smx_action_t) xbt_fifo_get_item_content(item1);
-          action2 = (smx_action_t) xbt_fifo_get_item_content(item2);
-
-          if(action1->type != action2->type){
-            XBT_INFO("Different type of action");
-            return 1;
-          }
-
-          if(action1->state != action2->state){
-            XBT_INFO("Different state of action");
-            return 1;
-          }
-
-          if(xbt_fifo_size(action1->simcalls) != xbt_fifo_size(action2->simcalls)){
-            XBT_INFO("Different size of simcall list (%d - %d", xbt_fifo_size(action1->simcalls), xbt_fifo_size(action2->simcalls));
-            return 1;
-          }else{
-
-            item_req1 = xbt_fifo_get_first_item(action1->simcalls);    
-            item_req2 = xbt_fifo_get_first_item(action2->simcalls);
-
-            while(j<xbt_fifo_size(action1->simcalls)){
-
-              req1 = (smx_simcall_t) xbt_fifo_get_item_content(item_req1);
-              req2 = (smx_simcall_t) xbt_fifo_get_item_content(item_req2);
-              
-              if(req1->call != req2->call){
-                XBT_INFO("Different simcall call in simcalls of action (%d - %d)", (int)req1->call, (int)req2->call);
-                return 1;
-              }
-              if(req1->issuer->pid != req2->issuer->pid){
-                XBT_INFO("Different simcall issuer in simcalls of action (%lu- %lu)", req1->issuer->pid, req2->issuer->pid);
-                return 1;
-              }
-
-              item_req1 = xbt_fifo_get_next_item(item_req1);   
-              item_req2 = xbt_fifo_get_next_item(item_req2);
-              j++;
-              
-            }
-          }
-
-          switch(action1->type){
-          case 0: /* execution */
-          case 1: /* parallel execution */
-            if(strcmp(action1->execution.host->name, action2->execution.host->name) != 0)
-              return 1;
-            break;
-          case 2: /* comm */
-            if(action1->comm.type != action2->comm.type)
-              return 1;
-            //XBT_INFO("Type of comm : %d", action1->comm.type);
-            
-            switch(action1->comm.type){
-            case 0: /* SEND */
-              if(action1->comm.src_proc->pid != action2->comm.src_proc->pid)
-                return 1;
-              if(strcmp(action1->comm.src_proc->smx_host->name, action2->comm.src_proc->smx_host->name) != 0)
-                return 1;
-              break;
-            case 1: /* RECEIVE */
-              if(action1->comm.dst_proc->pid != action2->comm.dst_proc->pid)
-                return 1;
-              if(strcmp(action1->comm.dst_proc->smx_host->name, action2->comm.dst_proc->smx_host->name) != 0)
-                return 1;
-              break;
-            case 2: /* READY */
-              if(action1->comm.src_proc->pid != action2->comm.src_proc->pid)
-                return 1;
-              if(strcmp(action1->comm.src_proc->smx_host->name, action2->comm.src_proc->smx_host->name) != 0)
-                return 1;
-              if(action1->comm.dst_proc->pid != action2->comm.dst_proc->pid)
-                return 1;
-              if(strcmp(action1->comm.dst_proc->smx_host->name, action2->comm.dst_proc->smx_host->name) != 0)
-                return 1;
-              break;
-            case 3: /* DONE */
-              if(action1->comm.src_proc->pid != action2->comm.src_proc->pid)
-                return 1;
-              if(strcmp(action1->comm.src_proc->smx_host->name, action2->comm.src_proc->smx_host->name) != 0)
-                return 1;
-              if(action1->comm.dst_proc->pid != action2->comm.dst_proc->pid)
-                return 1;
-              if(strcmp(action1->comm.dst_proc->smx_host->name, action2->comm.dst_proc->smx_host->name) != 0)
-                return 1;
-              break;
-              
-            } /* end of switch on type of comm */
-            
-            if(action1->comm.refcount != action2->comm.refcount)
-              return 1;
-            if(action1->comm.detached != action2->comm.detached)
-              return 1;
-            if(action1->comm.rate != action2->comm.rate)
-              return 1;
-            if(action1->comm.task_size != action2->comm.task_size)
-              return 1;
-            if(action1->comm.src_buff != action2->comm.src_buff)
-              return 1;
-            if(action1->comm.dst_buff != action2->comm.dst_buff)
-              return 1;
-            if(action1->comm.src_data != action2->comm.src_data)
-              return 1;
-            if(action1->comm.dst_data != action2->comm.dst_data)
-              return 1;
-            
-            break;
-          case 3: /* sleep */
-            if(strcmp(action1->sleep.host->name, action2->sleep.host->name) != 0)
-              return 1;
-            break;
-          case 4: /* synchro */
-            
-            break;
-          default:
-            break;
-          }
-
-          item1 = xbt_fifo_get_next_item(item1);       
-          item2 = xbt_fifo_get_next_item(item2);
-          i++;
-        }
-       }
-     }
-   }
-
-   return 0;
-   
+  mmalloc_set_current_heap(heap);
 }
 
-void set_pair_reached(xbt_state_t st){
 
-  MC_SET_RAW_MEM;
-  
-  mc_pair_reached_t pair = NULL;
-  pair = xbt_new0(s_mc_pair_reached_t, 1);
-  pair->automaton_state = st;
-  pair->prop_ato = xbt_dynar_new(sizeof(int), NULL);
-  pair->system_state = xbt_new0(s_mc_snapshot_t, 1);
-  //pair->rdv_points = xbt_dict_new();  
-  MC_take_snapshot_liveness(pair->system_state);
-
-  /* Get values of propositional symbols */
-  unsigned int cursor = 0;
-  xbt_propositional_symbol_t ps = NULL;
-  int res;
-  int_f_void_t f;
-
-  xbt_dynar_foreach(automaton->propositional_symbols, cursor, ps){
-    f = (int_f_void_t)ps->function;
-    res = (*f)();
-    xbt_dynar_push_as(pair->prop_ato, int, res);
-  }
+static int MC_automaton_evaluate_label(xbt_automaton_exp_label_t l,
+                                       xbt_dynar_t atomic_propositions_values)
+{
 
-  /*xbt_dict_t rdv_points = SIMIX_get_rdv_points();
-
-  xbt_dict_cursor_t cursor_dict = NULL;
-  char *key;
-  char *data;
-  xbt_fifo_item_t item;
-  smx_action_t action;
-
-  xbt_dict_foreach(rdv_points, cursor_dict, key, data){
-    smx_rdv_t new_rdv = xbt_new0(s_smx_rvpoint_t, 1);
-    new_rdv->name = strdup(((smx_rdv_t)data)->name);
-    new_rdv->comm_fifo = xbt_fifo_new();
-    xbt_fifo_foreach(((smx_rdv_t)data)->comm_fifo, item, action, smx_action_t) {
-      smx_action_t a = xbt_new0(s_smx_action_t, 1);
-      memcpy(a, action, sizeof(s_smx_action_t));
-      xbt_fifo_push(new_rdv->comm_fifo, a);
-      XBT_INFO("New action (type = %d, state = %d) in mailbox \"%s\"", action->type, action->state, key);
-      if(action->type==2)
-       XBT_INFO("Type of communication : %d, Ref count = %d", action->comm.type, action->comm.refcount);
+  switch (l->type) {
+  case 0:{
+      int left_res = MC_automaton_evaluate_label(l->u.or_and.left_exp, atomic_propositions_values);
+      int right_res = MC_automaton_evaluate_label(l->u.or_and.right_exp, atomic_propositions_values);
+      return (left_res || right_res);
     }
-    //new_rdv->comm_fifo = xbt_fifo_copy(((smx_rdv_t)data)->comm_fifo);
-    xbt_dict_set(pair->rdv_points, new_rdv->name, new_rdv, NULL);
-    }*/
-  xbt_dynar_push(reached_pairs, &pair); 
-
-  MC_UNSET_RAW_MEM;
-  
-}
-
-
-int reached_hash(xbt_state_t st){
-
-
-  if(xbt_dynar_is_empty(reached_pairs_hash)){
-
-    return 0;
-
-  }else{
-
-    MC_SET_RAW_MEM;
-
-    mc_snapshot_t sn = xbt_new0(s_mc_snapshot_t, 1);
-    MC_take_snapshot_liveness(sn);
-
-    int j;
-    unsigned int hash_regions[sn->num_reg];
-    for(j=0; j<sn->num_reg; j++){
-      hash_regions[j] = hash_region(sn->regions[j]->data, sn->regions[j]->size);
+  case 1:{
+      int left_res = MC_automaton_evaluate_label(l->u.or_and.left_exp, atomic_propositions_values);
+      int right_res = MC_automaton_evaluate_label(l->u.or_and.right_exp, atomic_propositions_values);
+      return (left_res && right_res);
     }
-
-
-    /* Get values of propositional symbols */
-    xbt_dynar_t prop_ato = xbt_dynar_new(sizeof(int), NULL);
-    unsigned int cursor = 0;
-    xbt_propositional_symbol_t ps = NULL;
-    int res;
-    int_f_void_t f;
-
-    xbt_dynar_foreach(automaton->propositional_symbols, cursor, ps){
-      f = (int_f_void_t)ps->function;
-      res = (*f)();
-      xbt_dynar_push_as(prop_ato, int, res);
+  case 2:{
+      int res = MC_automaton_evaluate_label(l->u.exp_not, atomic_propositions_values);
+      return (!res);
     }
-
-    mc_pair_reached_hash_t pair_test = NULL;
-
-    int region_diff = 0;
-
-    cursor = 0;
-
-    xbt_dynar_foreach(reached_pairs_hash, cursor, pair_test){
-
-      if(automaton_state_compare(pair_test->automaton_state, st) == 0){
-       if(propositional_symbols_compare_value(pair_test->prop_ato, prop_ato) == 0){
-         for(j=0 ; j< sn->num_reg ; j++){
-           if(hash_regions[j] != pair_test->hash_regions[j]){
-             region_diff++;
-           }
-         }
-         if(region_diff == 0){
-           MC_free_snapshot(sn);
-           xbt_dynar_reset(prop_ato);
-           xbt_free(prop_ato);
-           MC_UNSET_RAW_MEM;
-           return 1;
-         }else{
-           XBT_INFO("Different snapshot");
-         }
-       }else{
-         XBT_INFO("Different values of propositional symbols");
-       }
-      }else{
-       XBT_INFO("Different automaton state");
+  case 3:{
+      unsigned int cursor = 0;
+      xbt_automaton_propositional_symbol_t p = NULL;
+      xbt_dynar_foreach(_mc_property_automaton->propositional_symbols, cursor, p) {
+        if (strcmp(xbt_automaton_propositional_symbol_get_name(p), l->u.predicat) == 0)
+          return (int) xbt_dynar_get_as(atomic_propositions_values, cursor, int);
       }
-
-      region_diff = 0;
+      return -1;
     }
-    
-    MC_free_snapshot(sn);
-    xbt_dynar_reset(prop_ato);
-    xbt_free(prop_ato);
-    MC_UNSET_RAW_MEM;
-    return 0;
-    
-  }
-}
-
-void set_pair_reached_hash(xbt_state_t st){
-  MC_SET_RAW_MEM;
-
-  mc_snapshot_t sn = xbt_new0(s_mc_snapshot_t, 1);
-  MC_take_snapshot_liveness(sn);
-  mc_pair_reached_hash_t pair = NULL;
-  pair = xbt_new0(s_mc_pair_reached_hash_t, 1);
-  pair->automaton_state = st;
-  pair->prop_ato = xbt_dynar_new(sizeof(int), NULL);
-  pair->hash_regions = malloc(sizeof(unsigned int) * sn->num_reg);
-  
-  int i;
-
-  for(i=0 ; i< sn->num_reg ; i++){
-    pair->hash_regions[i] = hash_region(sn->regions[i]->data, sn->regions[i]->size);
-  }
-  
-  /* Get values of propositional symbols */
-  unsigned int cursor = 0;
-  xbt_propositional_symbol_t ps = NULL;
-  int res;
-  int_f_void_t f;
-
-  xbt_dynar_foreach(automaton->propositional_symbols, cursor, ps){
-    f = (int_f_void_t)ps->function;
-    res = (*f)();
-    xbt_dynar_push_as(pair->prop_ato, int, res);
-  }
-  
-  xbt_dynar_push(reached_pairs_hash, &pair);
-
-  MC_free_snapshot(sn);
-  
-  MC_UNSET_RAW_MEM;
-    
-}
-
-
-int visited(xbt_state_t st, int sc){
-
-
-  if(xbt_dynar_is_empty(visited_pairs)){
-
-    return 0;
-
-  }else{
-
-    MC_SET_RAW_MEM;
-
-    mc_snapshot_t sn = xbt_new0(s_mc_snapshot_t, 1);
-    MC_take_snapshot_liveness(sn);
-
-    xbt_dynar_t prop_ato = xbt_dynar_new(sizeof(int), NULL);
-
-    /* Get values of propositional symbols */
-    unsigned int cursor = 0;
-    xbt_propositional_symbol_t ps = NULL;
-    int res;
-    int_f_void_t f;
-
-    xbt_dynar_foreach(automaton->propositional_symbols, cursor, ps){
-      f = (int_f_void_t)ps->function;
-      res = (*f)();
-      xbt_dynar_push_as(prop_ato, int, res);
-    }
-
-    cursor = 0;
-    mc_pair_visited_t pair_test = NULL;
-
-    xbt_dynar_foreach(visited_pairs, cursor, pair_test){
-      if(pair_test->search_cycle == sc) {
-       if(automaton_state_compare(pair_test->automaton_state, st) == 0){
-         if(propositional_symbols_compare_value(pair_test->prop_ato, prop_ato) == 0){
-           if(snapshot_compare(pair_test->system_state, sn, std_heap, raw_heap) == 0){
-           
-             MC_free_snapshot(sn);
-             xbt_dynar_reset(prop_ato);
-             xbt_free(prop_ato);
-             MC_UNSET_RAW_MEM;
-               
-             return 1;
-       
-           }else{
-             XBT_INFO("Different snapshot");
-           }
-         }else{
-           XBT_INFO("Different values of propositional symbols"); 
-         }
-       }else{
-         XBT_INFO("Different automaton state");
-       }
-      }else{
-       XBT_INFO("Different value of search_cycle");
-      }
+  case 4:{
+      return 2;
     }
-
-
-    MC_free_snapshot(sn);
-    xbt_dynar_reset(prop_ato);
-    xbt_free(prop_ato);
-    MC_UNSET_RAW_MEM;
-    return 0;
-    
+  default:
+    return -1;
   }
 }
 
+void MC_pre_modelcheck_liveness(void) {
 
-int visited_hash(xbt_state_t st, int sc){
-
-
-  if(xbt_dynar_is_empty(visited_pairs_hash)){
-
-    return 0;
-
-  }else{
-
-    MC_SET_RAW_MEM;
+  initial_global_state->raw_mem_set = (mmalloc_get_current_heap() == mc_heap);
 
-    mc_snapshot_t sn = xbt_new0(s_mc_snapshot_t, 1);
-    MC_take_snapshot_liveness(sn);
+  mc_pair_t initial_pair = NULL;
+  smx_process_t process;
 
-    int j;
-    unsigned int hash_regions[sn->num_reg];
-    for(j=0; j<sn->num_reg; j++){
-      hash_regions[j] = hash_region(sn->regions[j]->data, sn->regions[j]->size);
-    }
+  MC_wait_for_requests();
 
-    
-    /* Get values of propositional symbols */
-    xbt_dynar_t prop_ato = xbt_dynar_new(sizeof(int), NULL);
-    unsigned int cursor = 0;
-    xbt_propositional_symbol_t ps = NULL;
-    int res;
-    int_f_void_t f;
-
-    xbt_dynar_foreach(automaton->propositional_symbols, cursor, ps){
-      f = (int_f_void_t)ps->function;
-      res = (*f)();
-      xbt_dynar_push_as(prop_ato, int, res);
-    }
+  MC_SET_MC_HEAP;
 
-    mc_pair_visited_hash_t pair_test = NULL;
+  acceptance_pairs = xbt_dynar_new(sizeof(mc_visited_pair_t), NULL);
+  if(_sg_mc_visited > 0)
+    visited_pairs = xbt_dynar_new(sizeof(mc_visited_pair_t), NULL);
 
-    int region_diff = 0;
-    cursor = 0;
+  initial_global_state->snapshot = MC_take_snapshot(0);
+  initial_global_state->prev_pair = 0;
 
-    xbt_dynar_foreach(visited_pairs_hash, cursor, pair_test){
+  unsigned int cursor = 0;
+  xbt_automaton_state_t automaton_state;
   
-      if(pair_test->search_cycle == sc) {
-       if(automaton_state_compare(pair_test->automaton_state, st) == 0){
-         if(propositional_symbols_compare_value(pair_test->prop_ato, prop_ato) == 0){
-           for(j=0 ; j< sn->num_reg ; j++){
-             if(hash_regions[j] != pair_test->hash_regions[j]){
-               region_diff++;
-             }
-           }
-           if(region_diff == 0){
-             MC_free_snapshot(sn);
-             xbt_dynar_reset(prop_ato);
-             xbt_free(prop_ato);
-             MC_UNSET_RAW_MEM;
-             return 1;
-           }else{
-             //XBT_INFO("Different snapshot");
-           }
-         }else{
-           //XBT_INFO("Different values of propositional symbols"); 
-         }
-       }else{
-         //XBT_INFO("Different automaton state");
-       }
-      }else{
-       //XBT_INFO("Different value of search_cycle");
-      }
-      
-      region_diff = 0;
+  xbt_dynar_foreach(_mc_property_automaton->states, cursor, automaton_state) {
+    if (automaton_state->type == -1) {  /* Initial automaton state */
+
+      initial_pair = MC_pair_new();
+      initial_pair->automaton_state = automaton_state;
+      initial_pair->graph_state = MC_state_new();
+      initial_pair->atomic_propositions = get_atomic_propositions_values();
+      initial_pair->depth = 1;
+
+      /* Get enabled processes and insert them in the interleave set of the graph_state */
+      MC_EACH_SIMIX_PROCESS(process,
+        if (MC_process_is_enabled(process)) {
+          MC_state_interleave_process(initial_pair->graph_state, process);
+        }
+      );
+
+      initial_pair->requests = MC_state_interleave_size(initial_pair->graph_state);
+      initial_pair->search_cycle = 0;
+
+      xbt_fifo_unshift(mc_stack, initial_pair);
     }
-    
-    MC_free_snapshot(sn);
-    xbt_dynar_reset(prop_ato);
-    xbt_free(prop_ato);
-    MC_UNSET_RAW_MEM;
-    return 0;
-    
-  }
-}
-
-void set_pair_visited_hash(xbt_state_t st, int sc){
-  MC_SET_RAW_MEM;
-
-  mc_snapshot_t sn = xbt_new0(s_mc_snapshot_t, 1);
-  MC_take_snapshot_liveness(sn);
-  mc_pair_visited_hash_t pair = NULL;
-  pair = xbt_new0(s_mc_pair_visited_hash_t, 1);
-  pair->automaton_state = st;
-  pair->prop_ato = xbt_dynar_new(sizeof(int), NULL);
-  pair->search_cycle = sc;
-  pair->hash_regions = malloc(sizeof(unsigned int) * sn->num_reg);
-  
-  int i;
-
-  for(i=0 ; i< sn->num_reg ; i++){
-    pair->hash_regions[i] = hash_region(sn->regions[i]->data, sn->regions[i]->size);
-  }
-  
-  /* Get values of propositional symbols */
-  unsigned int cursor = 0;
-  xbt_propositional_symbol_t ps = NULL;
-  int res;
-  int_f_void_t f;
-
-  xbt_dynar_foreach(automaton->propositional_symbols, cursor, ps){
-    f = (int_f_void_t)ps->function;
-    res = (*f)();
-    xbt_dynar_push_as(pair->prop_ato, int, res);
   }
-  
-  xbt_dynar_push(visited_pairs_hash, &pair);
 
-  MC_free_snapshot(sn);
+  MC_SET_STD_HEAP;
   
-  MC_UNSET_RAW_MEM;
-    
-}
-
-void set_pair_visited(xbt_state_t st, int sc){
+  MC_modelcheck_liveness();
 
-  MC_SET_RAW_MEM;
-  mc_pair_visited_t pair = NULL;
-  pair = xbt_new0(s_mc_pair_visited_t, 1);
-  pair->automaton_state = st;
-  pair->prop_ato = xbt_dynar_new(sizeof(int), NULL);
-  pair->search_cycle = sc;
-  pair->system_state = xbt_new0(s_mc_snapshot_t, 1);
-  MC_take_snapshot_liveness(pair->system_state);
-
-
-  /* Get values of propositional symbols */
-  unsigned int cursor = 0;
-  xbt_propositional_symbol_t ps = NULL;
-  int res;
-  int_f_void_t f;
-
-  xbt_dynar_foreach(automaton->propositional_symbols, cursor, ps){
-    f = (int_f_void_t)ps->function;
-    res = (*f)();
-    xbt_dynar_push_as(pair->prop_ato, int, res);
-  }
-  
-  xbt_dynar_push(visited_pairs, &pair);
-  
-  MC_UNSET_RAW_MEM;
-  
-  
-}
-
-void MC_pair_delete(mc_pair_t pair){
-  xbt_free(pair->graph_state->proc_status);
-  xbt_free(pair->graph_state);
-  xbt_free(pair);
+  if (initial_global_state->raw_mem_set)
+    MC_SET_MC_HEAP;
 }
 
-
-
-int MC_automaton_evaluate_label(xbt_exp_label_t l){
+void MC_modelcheck_liveness()
+{
+  smx_process_t process = NULL;
+  mc_pair_t current_pair = NULL;
+  int value, res, visited_num = -1;
+  smx_simcall_t req = NULL;
+  xbt_automaton_transition_t transition_succ = NULL;
+  int cursor = 0;
+  mc_pair_t next_pair = NULL;
+  xbt_dynar_t prop_values = NULL;
+  mc_visited_pair_t reached_pair = NULL;
   
-  switch(l->type){
-  case 0 : {
-    int left_res = MC_automaton_evaluate_label(l->u.or_and.left_exp);
-    int right_res = MC_automaton_evaluate_label(l->u.or_and.right_exp);
-    return (left_res || right_res);
-  }
-  case 1 : {
-    int left_res = MC_automaton_evaluate_label(l->u.or_and.left_exp);
-    int right_res = MC_automaton_evaluate_label(l->u.or_and.right_exp);
-    return (left_res && right_res);
-  }
-  case 2 : {
-    int res = MC_automaton_evaluate_label(l->u.exp_not);
-    return (!res);
-  }
-  case 3 : { 
-    unsigned int cursor = 0;
-    xbt_propositional_symbol_t p = NULL;
-    int_f_void_t f;
-    xbt_dynar_foreach(automaton->propositional_symbols, cursor, p){
-      if(strcmp(p->pred, l->u.predicat) == 0){
-       f = (int_f_void_t)p->function;
-       return (*f)();
+  while(xbt_fifo_size(mc_stack) > 0){
+
+    /* Get current pair */
+    current_pair = (mc_pair_t) xbt_fifo_get_item_content(xbt_fifo_get_first_item(mc_stack));
+
+    /* Update current state in buchi automaton */
+    _mc_property_automaton->current_state = current_pair->automaton_state;
+
+    XBT_DEBUG("********************* ( Depth = %d, search_cycle = %d, interleave size %d, pair_num %d)",
+       current_pair->depth, current_pair->search_cycle,
+       MC_state_interleave_size(current_pair->graph_state), current_pair->num);
+
+    if (current_pair->requests > 0) {
+
+      if (current_pair->automaton_state->type == 1 && current_pair->exploration_started == 0) {
+        /* If new acceptance pair, return new pair */
+        if ((reached_pair = is_reached_acceptance_pair(current_pair)) == NULL) {
+          int counter_example_depth = current_pair->depth;
+          XBT_INFO("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*");
+          XBT_INFO("|             ACCEPTANCE CYCLE            |");
+          XBT_INFO("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*");
+          XBT_INFO("Counter-example that violates formula :");
+          MC_show_stack_liveness(mc_stack);
+          MC_dump_stack_liveness(mc_stack);
+          MC_print_statistics(mc_stats);
+          XBT_INFO("Counter-example depth : %d", counter_example_depth);
+          xbt_abort();
+        }
       }
-    }
-    return -1;
-  }
-  case 4 : {
-    return 2;
-  }
-  default : 
-    return -1;
-  }
-}
-
-
-/********************* Double-DFS stateless *******************/
-
-void MC_pair_stateless_delete(mc_pair_stateless_t pair){
-  xbt_free(pair->graph_state->proc_status);
-  xbt_free(pair->graph_state);
-  xbt_free(pair);
-}
-
-mc_pair_stateless_t new_pair_stateless(mc_state_t sg, xbt_state_t st, int r){
-  mc_pair_stateless_t p = NULL;
-  p = xbt_new0(s_mc_pair_stateless_t, 1);
-  p->automaton_state = st;
-  p->graph_state = sg;
-  p->requests = r;
-  mc_stats_pair->expanded_pairs++;
-  return p;
-}
-
-void MC_ddfs_init(void){
-
-  XBT_INFO("**************************************************");
-  XBT_INFO("Double-DFS init");
-  XBT_INFO("**************************************************");
-
-  mc_pair_stateless_t mc_initial_pair = NULL;
-  mc_state_t initial_graph_state = NULL;
-  smx_process_t process; 
 
-  MC_wait_for_requests();
-
-  MC_SET_RAW_MEM;
-
-  initial_graph_state = MC_state_pair_new();
-  xbt_swag_foreach(process, simix_global->process_list){
-    if(MC_process_is_enabled(process)){
-      MC_state_interleave_process(initial_graph_state, process);
-    }
-  }
-
-  reached_pairs = xbt_dynar_new(sizeof(mc_pair_reached_t), NULL);
-  //reached_pairs_hash = xbt_dynar_new(sizeof(mc_pair_reached_hash_t), NULL);
-  //visited_pairs = xbt_dynar_new(sizeof(mc_pair_visited_t), NULL);
-  visited_pairs_hash = xbt_dynar_new(sizeof(mc_pair_visited_hash_t), NULL);
-  successors = xbt_dynar_new(sizeof(mc_pair_stateless_t), NULL);
-
-  /* Save the initial state */
-  initial_snapshot_liveness = xbt_new0(s_mc_snapshot_t, 1);
-  MC_take_snapshot_liveness(initial_snapshot_liveness);
-
-  MC_UNSET_RAW_MEM; 
-
-  unsigned int cursor = 0;
-  xbt_state_t state;
-
-  xbt_dynar_foreach(automaton->states, cursor, state){
-    if(state->type == -1){
-      
-      MC_SET_RAW_MEM;
-      mc_initial_pair = new_pair_stateless(initial_graph_state, state, MC_state_interleave_size(initial_graph_state));
-      xbt_fifo_unshift(mc_stack_liveness, mc_initial_pair);
-      MC_UNSET_RAW_MEM;
-      
-      if(cursor != 0){
-       MC_restore_snapshot(initial_snapshot_liveness);
-       MC_UNSET_RAW_MEM;
-      }
-
-      MC_ddfs(0);
+      /* Pair already visited ? stop the exploration on the current path */
+      if ((current_pair->exploration_started == 0) && (visited_num = is_visited_pair(reached_pair, current_pair)) != -1) {
+
+        if (dot_output != NULL){
+          MC_SET_MC_HEAP;
+          fprintf(dot_output, "\"%d\" -> \"%d\" [%s];\n", initial_global_state->prev_pair, visited_num, initial_global_state->prev_req);
+          fflush(dot_output);
+          MC_SET_STD_HEAP;
+        }
+
+        XBT_DEBUG("Pair already visited (equal to pair %d), exploration on the current path stopped.", visited_num);
+        MC_SET_MC_HEAP;
+        current_pair->requests = 0;
+        MC_SET_STD_HEAP;
+        goto backtracking;
+        
+      }else{
 
-    }else{
-      if(state->type == 2){
+        req = MC_state_get_request(current_pair->graph_state, &value);
+
+         if (dot_output != NULL) {
+           MC_SET_MC_HEAP;
+           if (initial_global_state->prev_pair != 0 && initial_global_state->prev_pair != current_pair->num) {
+             fprintf(dot_output, "\"%d\" -> \"%d\" [%s];\n", initial_global_state->prev_pair, current_pair->num, initial_global_state->prev_req);
+             xbt_free(initial_global_state->prev_req);
+           }
+           initial_global_state->prev_pair = current_pair->num;
+           initial_global_state->prev_req = MC_request_get_dot_output(req, value);
+           if (current_pair->search_cycle)
+             fprintf(dot_output, "%d [shape=doublecircle];\n", current_pair->num);
+           fflush(dot_output);
+           MC_SET_STD_HEAP;
+         }
+
+         char* req_str = MC_request_to_string(req, value);  
+         XBT_DEBUG("Execute: %s", req_str);                 
+         xbt_free(req_str);
+
+         /* Set request as executed */
+         MC_state_set_executed_request(current_pair->graph_state, req, value);
+
+         /* Update mc_stats */
+         mc_stats->executed_transitions++;
+         if(current_pair->exploration_started == 0)
+           mc_stats->visited_pairs++;
+
+         /* Answer the request */
+         MC_simcall_handle(req, value);
+         
+         /* Wait for requests (schedules processes) */
+         MC_wait_for_requests();
+
+         MC_SET_MC_HEAP;
+
+         current_pair->requests--;
+         current_pair->exploration_started = 1;
+
+         /* Get values of atomic propositions (variables used in the property formula) */ 
+         prop_values = get_atomic_propositions_values();
+
+         /* Evaluate enabled/true transitions in automaton according to atomic propositions values and create new pairs */
+         cursor = xbt_dynar_length(current_pair->automaton_state->out) - 1;
+         while (cursor >= 0) {
+           transition_succ = (xbt_automaton_transition_t)xbt_dynar_get_as(current_pair->automaton_state->out, cursor, xbt_automaton_transition_t);
+           res = MC_automaton_evaluate_label(transition_succ->label, prop_values);
+           if (res == 1 || res == 2) { /* 1 = True transition (always enabled), 2 = enabled transition according to atomic prop values */
+              next_pair = MC_pair_new();
+              next_pair->graph_state = MC_state_new();
+              next_pair->automaton_state = transition_succ->dst;
+              next_pair->atomic_propositions = get_atomic_propositions_values();
+              next_pair->depth = current_pair->depth + 1;
+              /* Get enabled processes and insert them in the interleave set of the next graph_state */
+              MC_EACH_SIMIX_PROCESS(process,
+                if (MC_process_is_enabled(process)) {
+                  MC_state_interleave_process(next_pair->graph_state, process);
+                }
+              );
+
+              next_pair->requests = MC_state_interleave_size(next_pair->graph_state);
+
+              /* FIXME : get search_cycle value for each acceptant state */
+              if (next_pair->automaton_state->type == 1 || current_pair->search_cycle)
+                next_pair->search_cycle = 1;
+
+              /* Add new pair to the exploration stack */
+              xbt_fifo_unshift(mc_stack, next_pair);
+
+           }
+           cursor--;
+         }
+         
+         MC_SET_STD_HEAP;
+        
+      } /* End of visited_pair test */
       
-       MC_SET_RAW_MEM;
-       mc_initial_pair = new_pair_stateless(initial_graph_state, state, MC_state_interleave_size(initial_graph_state));
-       xbt_fifo_unshift(mc_stack_liveness, mc_initial_pair);
-       MC_UNSET_RAW_MEM;
-
-       set_pair_reached(state);
-       //set_pair_reached_hash(state);
-
-       if(cursor != 0){
-         MC_restore_snapshot(initial_snapshot_liveness);
-         MC_UNSET_RAW_MEM;
-       }
-       
-       MC_ddfs(1);
-       
-      }
-    }
-  } 
-
-}
-
-
-void MC_ddfs(int search_cycle){
-
-  smx_process_t process;
-  mc_pair_stateless_t current_pair = NULL;
-
-  if(xbt_fifo_size(mc_stack_liveness) == 0)
-    return;
-
-
-  /* Get current pair */
-  current_pair = (mc_pair_stateless_t)xbt_fifo_get_item_content(xbt_fifo_get_first_item(mc_stack_liveness));
-
-  /* Update current state in buchi automaton */
-  automaton->current_state = current_pair->automaton_state;
-
-  XBT_INFO("********************* ( Depth = %d, search_cycle = %d )", xbt_fifo_size(mc_stack_liveness), search_cycle);
-  XBT_INFO("Pair : graph=%p, automaton=%p(%s), %u interleave", current_pair->graph_state, current_pair->automaton_state, current_pair->automaton_state->id, MC_state_interleave_size(current_pair->graph_state));
-  mc_stats_pair->visited_pairs++;
-
-  //sleep(1);
-
-  int value;
-  mc_state_t next_graph_state = NULL;
-  smx_simcall_t req = NULL;
-  char *req_str;
-
-  xbt_transition_t transition_succ;
-  unsigned int cursor = 0;
-  int res;
-
-  mc_pair_stateless_t next_pair = NULL;
-  mc_pair_stateless_t pair_succ;
-
-  if(xbt_fifo_size(mc_stack_liveness) < MAX_DEPTH_LIVENESS){
+    } else {
 
-    //set_pair_visited(current_pair->automaton_state, search_cycle);
-    set_pair_visited_hash(current_pair->automaton_state, search_cycle);
-
-    //XBT_INFO("Visited pairs : %lu", xbt_dynar_length(visited_pairs));
-    XBT_INFO("Visited pairs : %lu", xbt_dynar_length(visited_pairs_hash));
-
-    if(current_pair->requests > 0){
-
-      while((req = MC_state_get_request(current_pair->graph_state, &value)) != NULL){
-   
-       /* Debug information */
-
-       req_str = MC_request_to_string(req, value);
-       XBT_INFO("Execute: %s", req_str);
-       xbt_free(req_str);
-
-       MC_state_set_executed_request(current_pair->graph_state, req, value);   
+    backtracking:
+      if(visited_num == -1)
+        XBT_DEBUG("No more request to execute. Looking for backtracking point.");
     
-       /* Answer the request */
-       SIMIX_simcall_pre(req, value);
-
-       /* Wait for requests (schedules processes) */
-       MC_wait_for_requests();
-
-
-       MC_SET_RAW_MEM;
-
-       /* Create the new expanded graph_state */
-       next_graph_state = MC_state_pair_new();
-
-       /* Get 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);
-         }
-       }
-
-       xbt_dynar_reset(successors);
-
-       MC_UNSET_RAW_MEM;
-
-
-       cursor= 0;
-       xbt_dynar_foreach(current_pair->automaton_state->out, cursor, transition_succ){
-
-         res = MC_automaton_evaluate_label(transition_succ->label);
-
-         if(res == 1){ // enabled transition in automaton
-           MC_SET_RAW_MEM;
-           next_pair = new_pair_stateless(next_graph_state, transition_succ->dst, MC_state_interleave_size(next_graph_state));
-           xbt_dynar_push(successors, &next_pair);
-           MC_UNSET_RAW_MEM;
-         }
-
-       }
-
-       cursor = 0;
-   
-       xbt_dynar_foreach(current_pair->automaton_state->out, cursor, transition_succ){
-      
-         res = MC_automaton_evaluate_label(transition_succ->label);
-       
-         if(res == 2){ // true transition in automaton
-           MC_SET_RAW_MEM;
-           next_pair = new_pair_stateless(next_graph_state, transition_succ->dst, MC_state_interleave_size(next_graph_state));
-           xbt_dynar_push(successors, &next_pair);
-           MC_UNSET_RAW_MEM;
-         }
-
-       }
-
-       cursor = 0; 
-       
-       xbt_dynar_foreach(successors, cursor, pair_succ){
-
-         if(search_cycle == 1){
-
-           if((pair_succ->automaton_state->type == 1) || (pair_succ->automaton_state->type == 2)){ 
-                     
-             if(reached(pair_succ->automaton_state)){
-               //if(reached_hash(pair_succ->automaton_state)){
-             
-               XBT_INFO("Next pair (depth = %d, %u interleave) already reached !", xbt_fifo_size(mc_stack_liveness) + 1, MC_state_interleave_size(pair_succ->graph_state));
-
-               XBT_INFO("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*");
-               XBT_INFO("|             ACCEPTANCE CYCLE            |");
-               XBT_INFO("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*");
-               XBT_INFO("Counter-example that violates formula :");
-               MC_show_stack_liveness(mc_stack_liveness);
-               MC_dump_stack_liveness(mc_stack_liveness);
-               MC_print_statistics_pairs(mc_stats_pair);
-               exit(0);
-
-             }else{
-
-               XBT_INFO("Next pair (depth =%d) -> Acceptance pair : graph=%p, automaton=%p(%s)", xbt_fifo_size(mc_stack_liveness) + 1, pair_succ->graph_state, pair_succ->automaton_state, pair_succ->automaton_state->id);
-             
-               set_pair_reached(pair_succ->automaton_state);
-               //set_pair_reached_hash(pair_succ->automaton_state);
-
-               XBT_INFO("Reached pairs : %lu", xbt_dynar_length(reached_pairs));
-               //XBT_INFO("Reached pairs : %lu", xbt_dynar_length(reached_pairs_hash));
-
-               MC_SET_RAW_MEM;
-               xbt_fifo_unshift(mc_stack_liveness, pair_succ);
-               MC_UNSET_RAW_MEM;
-               
-               MC_ddfs(search_cycle);
-
-             }
-
-           }else{
-
-             if(!visited_hash(pair_succ->automaton_state, search_cycle)){
-               //if(!visited(pair_succ->automaton_state, search_cycle)){
-               
-               MC_SET_RAW_MEM;
-               xbt_fifo_unshift(mc_stack_liveness, pair_succ);
-               MC_UNSET_RAW_MEM;
-               
-               MC_ddfs(search_cycle);
-               
-             }else{
-
-               XBT_INFO("Next pair already visited ! ");
-
-             }
-             
-           }
-
-         }else{
-         
-           if(((pair_succ->automaton_state->type == 1) || (pair_succ->automaton_state->type == 2))){
-
-             XBT_INFO("Next pair (depth =%d) -> Acceptance pair : graph=%p, automaton=%p(%s)", xbt_fifo_size(mc_stack_liveness) + 1, pair_succ->graph_state, pair_succ->automaton_state, pair_succ->automaton_state->id);
-           
-             set_pair_reached(pair_succ->automaton_state); 
-             //set_pair_reached_hash(pair_succ->automaton_state);
-
-             search_cycle = 1;
-
-             XBT_INFO("Reached pairs : %lu", xbt_dynar_length(reached_pairs));
-             //XBT_INFO("Reached pairs : %lu", xbt_dynar_length(reached_pairs_hash));
-
-           }
-
-           if(!visited_hash(pair_succ->automaton_state, search_cycle)){
-             //if(!visited(pair_succ->automaton_state, search_cycle)){
-             
-             MC_SET_RAW_MEM;
-             xbt_fifo_unshift(mc_stack_liveness, pair_succ);
-             MC_UNSET_RAW_MEM;
-             
-             MC_ddfs(search_cycle);
-             
-           }else{
-
-             XBT_INFO("Next pair already visited ! ");
-
-           }
-
-         }
-
-        
-         /* Restore system before checking others successors */
-         if(cursor != (xbt_dynar_length(successors) - 1))
-           MC_replay_liveness(mc_stack_liveness, 1);
-       
-         
-       }
-
-       if(MC_state_interleave_size(current_pair->graph_state) > 0){
-         XBT_INFO("Backtracking to depth %d", xbt_fifo_size(mc_stack_liveness));
-         MC_replay_liveness(mc_stack_liveness, 0);
-       }
-      }
-
-    }else{  /*No request to execute, search evolution in Büchi automaton */
-
-      MC_SET_RAW_MEM;
-
-      /* Create the new expanded graph_state */
-      next_graph_state = MC_state_pair_new();
-
-      xbt_dynar_reset(successors);
-
-      MC_UNSET_RAW_MEM;
-
-
-      cursor= 0;
-      xbt_dynar_foreach(current_pair->automaton_state->out, cursor, transition_succ){
-
-       res = MC_automaton_evaluate_label(transition_succ->label);
-
-       if(res == 1){ // enabled transition in automaton
-         MC_SET_RAW_MEM;
-         next_pair = new_pair_stateless(next_graph_state, transition_succ->dst, MC_state_interleave_size(next_graph_state));
-         xbt_dynar_push(successors, &next_pair);
-         MC_UNSET_RAW_MEM;
-       }
-
-      }
-
-      cursor = 0;
-   
-      xbt_dynar_foreach(current_pair->automaton_state->out, cursor, transition_succ){
-      
-       res = MC_automaton_evaluate_label(transition_succ->label);
-       
-       if(res == 2){ // true transition in automaton
-         MC_SET_RAW_MEM;
-         next_pair = new_pair_stateless(next_graph_state, transition_succ->dst, MC_state_interleave_size(next_graph_state));
-         xbt_dynar_push(successors, &next_pair);
-         MC_UNSET_RAW_MEM;
-       }
-
-      }
-
-      cursor = 0; 
-     
-      xbt_dynar_foreach(successors, cursor, pair_succ){
-
-       if(search_cycle == 1){
-
-         if((pair_succ->automaton_state->type == 1) || (pair_succ->automaton_state->type == 2)){ 
-
-           if(reached(pair_succ->automaton_state)){
-             //if(reached_hash(pair_succ->automaton_state)){
-
-             XBT_INFO("Next pair (depth = %d) already reached !", xbt_fifo_size(mc_stack_liveness) + 1);
-
-             XBT_INFO("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*");
-             XBT_INFO("|             ACCEPTANCE CYCLE            |");
-             XBT_INFO("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*");
-             XBT_INFO("Counter-example that violates formula :");
-             MC_show_stack_liveness(mc_stack_liveness);
-             MC_dump_stack_liveness(mc_stack_liveness);
-             MC_print_statistics_pairs(mc_stats_pair);
-             exit(0);
-
-           }else{
-
-             XBT_INFO("Next pair (depth = %d) -> Acceptance pair : graph=%p, automaton=%p(%s)", xbt_fifo_size(mc_stack_liveness) + 1, pair_succ->graph_state, pair_succ->automaton_state, pair_succ->automaton_state->id);
-             
-             set_pair_reached(pair_succ->automaton_state);
-             //set_pair_reached_hash(pair_succ->automaton_state);
-               
-             XBT_INFO("Reached pairs : %lu", xbt_dynar_length(reached_pairs));
-             //XBT_INFO("Reached pairs : %lu", xbt_dynar_length(reached_pairs_hash));
-
-             MC_SET_RAW_MEM;
-             xbt_fifo_unshift(mc_stack_liveness, pair_succ);
-             MC_UNSET_RAW_MEM;
-             
-             MC_ddfs(search_cycle);
-
-           }
-
-         }else{
-
-           if(!visited_hash(pair_succ->automaton_state, search_cycle)){
-             //if(!visited(pair_succ->automaton_state, search_cycle)){
-
-             MC_SET_RAW_MEM;
-             xbt_fifo_unshift(mc_stack_liveness, pair_succ);
-             MC_UNSET_RAW_MEM;
-             
-             MC_ddfs(search_cycle);
-             
-           }else{
-
-             XBT_INFO("Next pair already visited ! ");
-
-           }
-         }
-           
-
-       }else{
-           
-         if(((pair_succ->automaton_state->type == 1) || (pair_succ->automaton_state->type == 2))){
-
-           set_pair_reached(pair_succ->automaton_state);
-           //set_pair_reached_hash(pair_succ->automaton_state);
-                   
-           search_cycle = 1;
-
-           XBT_INFO("Reached pairs : %lu", xbt_dynar_length(reached_pairs));
-           //XBT_INFO("Reached pairs : %lu", xbt_dynar_length(reached_pairs_hash));
-
-         }
-
-         if(!visited_hash(pair_succ->automaton_state, search_cycle)){
-           //if(!visited(pair_succ->automaton_state, search_cycle)){
-
-           MC_SET_RAW_MEM;
-           xbt_fifo_unshift(mc_stack_liveness, pair_succ);
-           MC_UNSET_RAW_MEM;
-           
-           MC_ddfs(search_cycle);
-           
-         }else{
-
-           XBT_INFO("Next pair already visited ! ");
-
-         }
-
-       }
-
-       /* Restore system before checking others successors */
-       if(cursor != xbt_dynar_length(successors) - 1)
-         MC_replay_liveness(mc_stack_liveness, 1);
-
-        
+      MC_SET_MC_HEAP;
+    
+      xbt_dynar_free(&prop_values);
+    
+      /* Traverse the stack backwards until a pair with a non empty interleave
+         set is found, deleting all the pairs that have it empty in the way. */
+      while ((current_pair = xbt_fifo_shift(mc_stack)) != NULL) {
+        if (current_pair->requests > 0) {
+          /* We found a backtracking point */
+          XBT_DEBUG("Backtracking to depth %d", current_pair->depth);
+          xbt_fifo_unshift(mc_stack, current_pair);
+          MC_SET_STD_HEAP;
+          MC_replay_liveness(mc_stack);
+          XBT_DEBUG("Backtracking done");
+          break;
+        }else{
+          /* Delete pair */
+          XBT_DEBUG("Delete pair %d at depth %d", current_pair->num, current_pair->depth);
+          if (current_pair->automaton_state->type == 1) 
+            remove_acceptance_pair(current_pair->num);
+          MC_pair_delete(current_pair);
+        }
       }
-           
-    }
     
-  }else{
+      MC_SET_STD_HEAP;
+    } /* End of if (current_pair->requests > 0) else ... */
     
-    XBT_INFO("Max depth reached");
-
-  }
-
-  if(xbt_fifo_size(mc_stack_liveness) == MAX_DEPTH_LIVENESS ){
-    XBT_INFO("Pair (graph=%p, automaton =%p, search_cycle = %d, depth = %d) shifted in stack, maximum depth reached", current_pair->graph_state, current_pair->automaton_state, search_cycle, xbt_fifo_size(mc_stack_liveness) );
-  }else{
-    XBT_INFO("Pair (graph=%p, automaton =%p, search_cycle = %d, depth = %d) shifted in stack", current_pair->graph_state, current_pair->automaton_state, search_cycle, xbt_fifo_size(mc_stack_liveness) );
-  }
-
+  } /* End of while(xbt_fifo_size(mc_stack) > 0) */
   
-  MC_SET_RAW_MEM;
-  xbt_fifo_shift(mc_stack_liveness);
-  if((current_pair->automaton_state->type == 1) || (current_pair->automaton_state->type == 2)){
-    xbt_dynar_pop(reached_pairs, NULL);
-    //xbt_dynar_pop(reached_pairs_hash, NULL);
-  }
-  MC_UNSET_RAW_MEM;
-  
-  
-
 }