Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
fix dot output with file descriptor checkpoint/restore
[simgrid.git] / src / mc / mc_liveness.c
index 980d7a4..1b8c3dd 100644 (file)
-#include "private.h"
-#include "unistd.h"
+/* Copyright (c) 2011-2014. The SimGrid Team.
+ * All rights reserved.                                                     */
 
-XBT_LOG_NEW_DEFAULT_SUBCATEGORY(mc_liveness, mc,
-                                "Logging specific to algorithms for liveness properties verification");
+/* 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. */
 
-xbt_dynar_t initial_pairs = NULL;
-int max_pair = 0;
-xbt_dynar_t visited_pairs;
-xbt_dynar_t reached_pairs;
-int expanded = 0;
-
-mc_pairs_t new_pair(mc_snapshot_t sn, mc_state_t sg, xbt_state_t st){
-  mc_pairs_t p = NULL;
-  p = xbt_new0(s_mc_pairs_t, 1);
-  p->system_state = sn;
-  p->automaton_state = st;
-  p->graph_state = sg;
-  p->num = max_pair;
-  max_pair++;
-  mc_stats_pair->expanded_pairs++;
-  return p;
-    
-}
+#include <unistd.h>
+#include <sys/wait.h>
 
-void set_pair_visited(mc_state_t gs, xbt_state_t as, int sc){
+#include <xbt/dynar.h>
+#include <xbt/automaton.h>
 
-  MC_SET_RAW_MEM;
-  mc_visited_pairs_t p = NULL;
-  p = xbt_new0(s_mc_visited_pairs_t, 1);
-  p->graph_state = gs;
-  p->automaton_state = as;
-  p->search_cycle = sc;
-  char *hash = malloc(sizeof(char)*128);
-  xbt_sha((char *)&p, hash);
-  xbt_dynar_push(visited_pairs, &hash); 
-  MC_UNSET_RAW_MEM;    
+#include "mc_request.h"
+#include "mc_liveness.h"
+#include "mc_private.h"
+#include "mc_record.h"
 
-}
+XBT_LOG_NEW_DEFAULT_SUBCATEGORY(mc_liveness, mc,
+                                "Logging specific to algorithms for liveness properties verification");
 
-int reached(mc_state_t gs, xbt_state_t as ){
-
-  /*mc_reached_pairs_t rp = NULL;
-  unsigned int c= 0;
-  unsigned int i;
-  int different_process;
-
-  MC_SET_RAW_MEM;
-  xbt_dynar_foreach(reached_pairs, c, rp){
-    if(rp->automaton_state == as){
-      different_process=0;
-      for(i=0; i < gs->max_pid; i++){
-       if(gs->proc_status[i].state != rp->graph_state->proc_status[i].state){
-         different_process++;
-       }
-      }
-      if(different_process==0){
-       MC_UNSET_RAW_MEM;
-       XBT_DEBUG("Pair (graph=%p, automaton=%p(%s)) already reached (graph=%p)!", gs, as, as->id, rp->graph_state);
-       return 1;
-      }
-    }
-  }
-  MC_UNSET_RAW_MEM;
-  return 0;*/
-  //XBT_DEBUG("Search acceptance pair already reach !"); 
-
-  char* hash_reached = malloc(sizeof(char)*128);
-  unsigned int c= 0;
-
-  MC_SET_RAW_MEM;
-  mc_reached_pairs_t p = NULL;
-  p = xbt_new0(s_mc_reached_pairs_t, 1);
-  p->graph_state = gs;
-  p->automaton_state = as;
-  char *hash = malloc(sizeof(char)*128);
-  xbt_sha((char *)&p, hash);
-  //XBT_DEBUG("Hash : %s", hash);
-  xbt_dynar_foreach(reached_pairs, c, hash_reached){
-    if(strcmp(hash, hash_reached) == 0){
-      MC_UNSET_RAW_MEM;
-      //XBT_DEBUG("Pair (graph=%p, automaton=%p(%s)) already reached !", gs, as, as->id);
-      return 1;
-    }
-  }
-  
-  MC_UNSET_RAW_MEM;
-  return 0;
-}
+/********* Global variables *********/
 
-void set_pair_reached(mc_state_t gs, xbt_state_t as){
-
-  //XBT_DEBUG("Set pair (graph=%p, automaton=%p) reached ", gs, as);
-  if(reached(gs, as) == 0){
-    MC_SET_RAW_MEM;
-    mc_reached_pairs_t p = NULL;
-    p = xbt_new0(s_mc_reached_pairs_t, 1);
-    p->graph_state = gs;
-    p->automaton_state = as;
-    char *hash = malloc(sizeof(char)*128) ;
-    xbt_sha((char *)&p, hash);
-    xbt_dynar_push(reached_pairs, &hash); 
-    //XBT_DEBUG("New reached pair : graph=%p, automaton=%p(%s)", gs, as, as->id);
-    MC_UNSET_RAW_MEM;
-  }
+xbt_dynar_t acceptance_pairs;
+xbt_dynar_t successors;
+xbt_parmap_t parmap;
 
-}
+/********* Static functions *********/
 
-int visited(mc_state_t gs, xbt_state_t as, int sc){
-  /*unsigned int c = 0;
-  mc_visited_pairs_t p = NULL;
-  unsigned int i;
-  int different_process;
-
-  //XBT_DEBUG("Visited pair length : %lu", xbt_dynar_length(visited_pairs));
-
-  MC_SET_RAW_MEM;
-  xbt_dynar_foreach(visited_pairs, c, p){
-    //XBT_DEBUG("Test pair visited");
-    //sleep(1);
-    if((p->automaton_state == as) && (p->search_cycle == sc)){
-      different_process=0;
-      for(i=0; i < gs->max_pid; i++){
-       if(gs->proc_status[i].state != p->graph_state->proc_status[i].state){
-         different_process++;
-       }
-      }
-      if(different_process==0){
-       MC_UNSET_RAW_MEM;
-       return 1;
-      }
-    }
-  }
-  MC_UNSET_RAW_MEM;
-  return 0;*/
-
-  //XBT_DEBUG("Test visited pair");
-
-  char* hash_visited = malloc(sizeof(char)*128);
-  unsigned int c= 0;
-
-  MC_SET_RAW_MEM;
-  mc_visited_pairs_t p = NULL;
-  p = xbt_new0(s_mc_visited_pairs_t, 1);
-  p->graph_state = gs;
-  p->automaton_state = as;
-  char *hash = malloc(sizeof(char)*128);
-  xbt_sha((char *)&p, hash);
-  //XBT_DEBUG("Hash : %s", hash);
-  xbt_dynar_foreach(visited_pairs, c, hash_visited){
-    if(strcmp(hash, hash_visited) == 0){
-      MC_UNSET_RAW_MEM;
-      return 1;
-    }
+static xbt_dynar_t get_atomic_propositions_values()
+{
+  int res;
+  int_f_void_t f;
+  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) {
+    f = (int_f_void_t) ps->function;
+    res = f();
+    xbt_dynar_push_as(values, int, res);
   }
-  
-  MC_UNSET_RAW_MEM;
-  //XBT_DEBUG("End test visited pair");
-  return 0;
 
+  return values;
 }
 
-void MC_dfs_init(xbt_automaton_t a){
-
-  XBT_DEBUG("**************************************************");
-  XBT_DEBUG("DFS init");
-  XBT_DEBUG("**************************************************");
-  mc_pairs_t mc_initial_pair;
-  mc_state_t initial_graph_state;
-  smx_process_t process; 
-  mc_snapshot_t init_snapshot;
-  MC_wait_for_requests();
-
-  MC_SET_RAW_MEM;
-
-  init_snapshot = xbt_new0(s_mc_snapshot_t, 1);
-
-  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);
-    }
-  }
-
-  visited_pairs = xbt_dynar_new(sizeof(char *), NULL); 
-  reached_pairs = xbt_dynar_new(sizeof(char *), NULL); 
 
-  MC_take_snapshot(init_snapshot);
+static mc_visited_pair_t is_reached_acceptance_pair(int pair_num,
+                                                    xbt_automaton_state_t
+                                                    automaton_state,
+                                                    xbt_dynar_t
+                                                    atomic_propositions)
+{
 
-  MC_UNSET_RAW_MEM; 
+  int raw_mem_set = (mmalloc_get_current_heap() == mc_heap);
 
-  /* unsigned int cursor = 0; */
-  /* xbt_state_t state = NULL; */
-  /* int nb_init_state = 0; */
+  MC_SET_MC_HEAP;
 
-  /* xbt_dynar_foreach(a->states, cursor, state){ */
-  /*   if(state->type == -1){ */
+  mc_visited_pair_t pair = NULL;
+  pair = MC_visited_pair_new(pair_num, automaton_state, atomic_propositions);
+  pair->acceptance_pair = 1;
 
-  /*     MC_SET_RAW_MEM; */
-  /*     mc_initial_pair = new_pair(init_snapshot, initial_graph_state, state); */
+  if (xbt_dynar_is_empty(acceptance_pairs)) {
 
-  /*     xbt_fifo_unshift(mc_snapshot_stack, mc_initial_pair); */
+    xbt_dynar_push(acceptance_pairs, &pair);
 
-  /*     XBT_DEBUG("**************************************************"); */
-  /*     XBT_DEBUG("Initial state=%p ", mc_initial_pair); */
+  } else {
 
-  /*     MC_UNSET_RAW_MEM; */
+    int min = -1, max = -1, index;
+    //int res;
+    mc_visited_pair_t pair_test;
+    int cursor;
 
-  /*     set_pair_visited(mc_initial_pair->graph_state, mc_initial_pair->automaton_state, 0); */
-      
-  /*     if(nb_init_state == 0) */
-  /*   MC_dfs(a, 0, 0); */
-  /*     else */
-  /*   MC_dfs(a, 0, 1); */
-      
-  /*     nb_init_state++; */
-  /*   } */
-  /* } */
+    index = get_search_interval(acceptance_pairs, pair, &min, &max);
 
-  /* regarder dans l'automate toutes les transitions activables grâce à l'état initial du système 
-    -> donnera les états initiaux de la propriété consistants avec l'état initial du système */
+    if (min != -1 && max != -1) {       // Acceptance pair with same number of processes and same heap bytes used exists
 
-  unsigned int cursor = 0;
-  unsigned int cursor2 = 0;
-  xbt_state_t state = NULL;
-  int res;
-  xbt_transition_t transition_succ;
-  xbt_dynar_t successors = xbt_dynar_new(sizeof(mc_pairs_t), NULL);
-  mc_pairs_t pair_succ;
+      // 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;
+         } */
 
-  xbt_dynar_foreach(a->states, cursor, state){
-    if(state->type == -1){
-      xbt_dynar_foreach(state->out, cursor2, transition_succ){
-       res = MC_automaton_evaluate_label(a, transition_succ->label);
-       
-       if((res == 1) || (res == 2)){
-        
-         MC_SET_RAW_MEM;
+      cursor = min;
+      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, pair->automaton_state) == 0) {
+          if (xbt_automaton_propositional_symbols_compare_value
+              (pair_test->atomic_propositions,
+               pair->atomic_propositions) == 0) {
+            if (snapshot_compare(pair_test, pair) == 0) {
+              XBT_INFO("Pair %d already reached (equal to pair %d) !",
+                       pair->num, pair_test->num);
 
-         mc_initial_pair = new_pair(init_snapshot, initial_graph_state, transition_succ->dst);
-         xbt_dynar_push(successors, &mc_initial_pair);
+              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);
 
-         MC_UNSET_RAW_MEM;
+              if (!raw_mem_set)
+                MC_SET_STD_HEAP;
 
-       }
+              return NULL;
+            }
+          }
+        }
+        cursor++;
       }
-    }
-  }
-
-  cursor = 0;
-
-  if(xbt_dynar_length(successors) == 0){
-    xbt_dynar_foreach(a->states, cursor, state){
-      if(state->type == -1){
-       MC_SET_RAW_MEM;
-
-       mc_initial_pair = new_pair(init_snapshot, initial_graph_state, transition_succ->dst);
-       xbt_dynar_push(successors, &mc_initial_pair);
-       
-       MC_UNSET_RAW_MEM;
+      xbt_dynar_insert_at(acceptance_pairs, min, &pair);
+    } else {
+      pair_test =
+          (mc_visited_pair_t) xbt_dynar_get_as(acceptance_pairs, index,
+                                               mc_visited_pair_t);
+      if (pair_test->nb_processes < pair->nb_processes) {
+        xbt_dynar_insert_at(acceptance_pairs, index + 1, &pair);
+      } else {
+        if (pair_test->heap_bytes_used < pair->heap_bytes_used)
+          xbt_dynar_insert_at(acceptance_pairs, index + 1, &pair);
+        else
+          xbt_dynar_insert_at(acceptance_pairs, index, &pair);
       }
     }
-  }
-
-  cursor = 0;
-  xbt_dynar_foreach(successors, cursor, pair_succ){
-    MC_SET_RAW_MEM;
-
-    xbt_fifo_unshift(mc_snapshot_stack, pair_succ);
-    set_pair_visited(pair_succ->graph_state, pair_succ->automaton_state, 0);
 
-    MC_UNSET_RAW_MEM;
-
-    if(cursor == 0){
-      MC_dfs(a, 0, 0);
-    }else{
-      MC_dfs(a, 0, 1);
-    }
   }
 
+  if (!raw_mem_set)
+    MC_SET_STD_HEAP;
 
-  
-}
+  return pair;
 
+}
 
-void MC_dfs(xbt_automaton_t a, int search_cycle, int restore){
+static void remove_acceptance_pair(int pair_num)
+{
 
-  smx_process_t process = NULL;
+  int raw_mem_set = (mmalloc_get_current_heap() == mc_heap);
 
+  MC_SET_MC_HEAP;
 
-  if(xbt_fifo_size(mc_snapshot_stack) == 0)
-    return;
+  unsigned int cursor = 0;
+  mc_visited_pair_t pair_test = NULL;
 
-  if(restore == 1){
-    MC_restore_snapshot(((mc_pairs_t)xbt_fifo_get_item_content(xbt_fifo_get_first_item(mc_snapshot_stack)))->system_state);
-    MC_UNSET_RAW_MEM;
+  xbt_dynar_foreach(acceptance_pairs, cursor, pair_test) {
+    if (pair_test->num == pair_num) {
+      break;
+    }
   }
 
+  xbt_dynar_remove_at(acceptance_pairs, cursor, &pair_test);
 
-  /* Get current state */
-  mc_pairs_t current_pair = (mc_pairs_t)xbt_fifo_get_item_content(xbt_fifo_get_first_item(mc_snapshot_stack));
+  pair_test->acceptance_removed = 1;
 
-  /*if(restore==1){
-    xbt_swag_foreach(process, simix_global->process_list){
-      if(MC_process_is_enabled(process)){
-       //XBT_DEBUG("Pid : %lu", process->pid);
-       MC_state_interleave_process(current_pair->graph_state, process);
-      }
-    }
-    }*/
+  if (_sg_mc_visited == 0) {
+    MC_visited_pair_delete(pair_test);
+  } else if (pair_test->visited_removed == 1) {
+    MC_visited_pair_delete(pair_test);
+  }
 
-  XBT_DEBUG("************************************************** ( search_cycle = %d )", search_cycle);
-  XBT_DEBUG("State : 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));
-  
+  if (!raw_mem_set)
+    MC_SET_STD_HEAP;
+}
 
-  mc_stats_pair->visited_pairs++;
 
-  int value;
-  mc_state_t next_graph_state = NULL;
-  smx_req_t req = NULL;
-  char *req_str;
+static int MC_automaton_evaluate_label(xbt_automaton_exp_label_t l,
+                                       xbt_dynar_t atomic_propositions_values)
+{
 
-  mc_pairs_t pair_succ;
-  xbt_transition_t transition_succ;
-  unsigned int cursor;
-  int res;
-  //int enabled_transition = 0;
-
-  xbt_dynar_t successors = xbt_dynar_new(sizeof(mc_pairs_t), NULL);
-
-  mc_pairs_t next_pair;
-  mc_snapshot_t next_snapshot;
-  mc_snapshot_t current_snapshot;
-  
-  
-  while((req = MC_state_get_request(current_pair->graph_state, &value)) != NULL){
-
-    MC_SET_RAW_MEM;
-    current_snapshot = xbt_new0(s_mc_snapshot_t, 1);
-    MC_take_snapshot(current_snapshot);
-    MC_UNSET_RAW_MEM;
-   
-    
-    /* Debug information */
-    if(XBT_LOG_ISENABLED(mc_liveness, xbt_log_priority_debug)){
-      req_str = MC_request_to_string(req, value);
-      XBT_DEBUG("Execute: %s", req_str);
-      xbt_free(req_str);
+  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);
     }
-
-    MC_state_set_executed_request(current_pair->graph_state, req, value);
-
-    /* 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_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);
+  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);
+    }
+  case 2:{
+      int res =
+          MC_automaton_evaluate_label(l->u.exp_not, atomic_propositions_values);
+      return (!res);
+    }
+  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(p->pred, l->u.predicat) == 0)
+          return (int) xbt_dynar_get_as(atomic_propositions_values, cursor,
+                                        int);
       }
+      return -1;
     }
+  case 4:{
+      return 2;
+    }
+  default:
+    return -1;
+  }
+}
 
-    next_snapshot = xbt_new0(s_mc_snapshot_t, 1);
-    MC_take_snapshot(next_snapshot);
-      
-    MC_UNSET_RAW_MEM;
-
-    xbt_dynar_reset(successors);
-    
-    cursor = 0;
-    xbt_dynar_foreach(current_pair->automaton_state->out, cursor, transition_succ){
-
-      res = MC_automaton_evaluate_label(a, transition_succ->label);
+void MC_pre_modelcheck_liveness(void)
+{
 
-      MC_SET_RAW_MEM;
-       
-      if((res == 1) || (res == 2)){ // enabled transition in automaton
-       next_pair = new_pair(next_snapshot,next_graph_state, transition_succ->dst);
-       xbt_dynar_push(successors, &next_pair);
-      }
+  initial_global_state->raw_mem_set = (mmalloc_get_current_heap() == mc_heap);
 
-      MC_UNSET_RAW_MEM;
-    }
+  mc_pair_t initial_pair = NULL;
+  smx_process_t process;
 
-   
-    if(xbt_dynar_length(successors) == 0){
+  MC_wait_for_requests();
 
-      MC_SET_RAW_MEM;
-      next_pair = new_pair(next_snapshot, next_graph_state, current_pair->automaton_state);
-      xbt_dynar_push(successors, &next_pair);
-      MC_UNSET_RAW_MEM;
-       
-    }
+  MC_SET_MC_HEAP;
 
+  acceptance_pairs = xbt_dynar_new(sizeof(mc_visited_pair_t), NULL);
+  visited_pairs = xbt_dynar_new(sizeof(mc_visited_pair_t), NULL);
+  successors = xbt_dynar_new(sizeof(mc_pair_t), NULL);
 
-    cursor = 0;
-    xbt_dynar_foreach(successors, cursor, pair_succ){
+  initial_global_state->snapshot = MC_take_snapshot(0);
+  initial_global_state->prev_pair = 0;
 
-      //XBT_DEBUG("Search visited pair : graph=%p, automaton=%p", pair_succ->graph_state, pair_succ->automaton_state);
+  MC_SET_STD_HEAP;
 
-      if((search_cycle == 1) && (reached(pair_succ->graph_state, pair_succ->automaton_state) == 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(visited(pair_succ->graph_state, pair_succ->automaton_state, search_cycle) == 0){
+  unsigned int cursor = 0;
+  xbt_automaton_state_t automaton_state;
 
-       //XBT_DEBUG("Unvisited pair !");
+  xbt_dynar_foreach(_mc_property_automaton->states, cursor, automaton_state) {
+    if (automaton_state->type == -1) {  /* Initial automaton state */
 
-       mc_stats_pair->executed_transitions++;
-       MC_SET_RAW_MEM;
-       xbt_fifo_unshift(mc_snapshot_stack, pair_succ);
-       set_pair_visited(pair_succ->graph_state, pair_succ->automaton_state, search_cycle);
-       MC_UNSET_RAW_MEM;
+      MC_SET_MC_HEAP;
 
-       MC_dfs(a, search_cycle, 0);
+      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();
 
-       if((search_cycle == 0) && (current_pair->automaton_state->type == 1)){
+      /* Get enabled processes and insert them in the interleave set of the graph_state */
+      xbt_swag_foreach(process, simix_global->process_list) {
+        if (MC_process_is_enabled(process)) {
+          MC_state_interleave_process(initial_pair->graph_state, process);
+        }
+      }
 
-         MC_restore_snapshot(current_pair->system_state);
-         MC_UNSET_RAW_MEM;
-         /*xbt_swag_foreach(process, simix_global->process_list){
-           if(MC_process_is_enabled(process)){
-             //XBT_DEBUG("Pid : %lu", process->pid);
-             MC_state_interleave_process(current_pair->graph_state, process);
-           }
-           }*/
-           
-         set_pair_reached(current_pair->graph_state, current_pair->automaton_state);
-         XBT_DEBUG("Acceptance pair : graph=%p, automaton=%p(%s)", current_pair->graph_state, current_pair->automaton_state, current_pair->automaton_state->id);
-         MC_dfs(a, 1, 1);
+      initial_pair->requests =
+          MC_state_interleave_size(initial_pair->graph_state);
+      initial_pair->search_cycle = 0;
 
-       }
-      }else{
+      xbt_fifo_unshift(mc_stack, initial_pair);
+      
+      MC_SET_STD_HEAP;
 
-       XBT_DEBUG("Pair already visited !");
+      MC_modelcheck_liveness();
 
+      if (cursor != 0) {
+        MC_restore_snapshot(initial_global_state->snapshot);
+        MC_SET_STD_HEAP;
       }
     }
-    
-    
-    if(MC_state_interleave_size(current_pair->graph_state) > 0){
-      MC_restore_snapshot(current_snapshot);
-      MC_UNSET_RAW_MEM;
-      //XBT_DEBUG("Snapshot restored");
-    }
   }
 
-  
-  MC_SET_RAW_MEM;
-  //remove_pair_reached(current_pair->graph_state);
-  xbt_fifo_shift(mc_snapshot_stack);
-  XBT_DEBUG("State (graph=%p, automaton =%p) shifted in snapshot_stack", current_pair->graph_state, current_pair->automaton_state);
-  MC_UNSET_RAW_MEM;
+  if (initial_global_state->raw_mem_set)
+    MC_SET_MC_HEAP;
+  else
+    MC_SET_STD_HEAP;
 
-}
 
-void MC_show_snapshot_stack(xbt_fifo_t stack){
-  int value;
-  mc_pairs_t pair;
-  xbt_fifo_item_t item;
-  smx_req_t req;
-  char *req_str = NULL;
-  
-  for (item = xbt_fifo_get_last_item(stack);
-       (item ? (pair = (mc_pairs_t) (xbt_fifo_get_item_content(item)))
-        : (NULL)); item = xbt_fifo_get_prev_item(item)) {
-    req = MC_state_get_executed_request(pair->graph_state, &value);
-    if(req){
-      req_str = MC_request_to_string(req, value);
-      XBT_INFO("%s", req_str);
-      xbt_free(req_str);
-    }
-  }
 }
 
-void MC_dump_snapshot_stack(xbt_fifo_t stack){
-  mc_pairs_t pair;
 
-  MC_SET_RAW_MEM;
-  while ((pair = (mc_pairs_t) xbt_fifo_pop(stack)) != NULL)
-    MC_pair_delete(pair);
-  MC_UNSET_RAW_MEM;
-}
+void MC_modelcheck_liveness()
+{
 
-void MC_pair_delete(mc_pairs_t pair){
-  xbt_free(pair->graph_state->proc_status);
-  xbt_free(pair->graph_state);
-  //xbt_free(pair->automaton_state); -> FIXME : à implémenter
-  xbt_free(pair);
-}
+  smx_process_t process;
+  mc_pair_t current_pair = NULL;
 
+  if (xbt_fifo_size(mc_stack) == 0)
+    return;
 
-int MC_automaton_evaluate_label(xbt_automaton_t a, xbt_exp_label_t l){
-  
-  switch(l->type){
-  case 0 : {
-    int left_res = MC_automaton_evaluate_label(a, l->u.or_and.left_exp);
-    int right_res = MC_automaton_evaluate_label(a, l->u.or_and.right_exp);
-    return (left_res || right_res);
-    break;
-  }
-  case 1 : {
-    int left_res = MC_automaton_evaluate_label(a, l->u.or_and.left_exp);
-    int right_res = MC_automaton_evaluate_label(a, l->u.or_and.right_exp);
-    return (left_res && right_res);
-    break;
-  }
-  case 2 : {
-    int res = MC_automaton_evaluate_label(a, l->u.exp_not);
-    return (!res);
-    break;
-  }
-  case 3 : { 
-    unsigned int cursor = 0;
-    xbt_propositional_symbol_t p = NULL;
-    xbt_dynar_foreach(a->propositional_symbols, cursor, p){
-      if(strcmp(p->pred, l->u.predicat) == 0){
-       int (*f)() = p->function;
-       return (*f)();
-      }
-    }
-    return -1;
-    break;
-  }
-  case 4 : {
-    return 2;
-    break;
-  }
-  default : 
-    return -1;
-  }
-}
+  /* 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;
 
-/********************* Double DFS using color pair *********************/
+  XBT_DEBUG
+      ("********************* ( Depth = %d, search_cycle = %d, interleave size %d, pair_num %d)",
+       xbt_fifo_size(mc_stack), current_pair->search_cycle,
+       MC_state_interleave_size(current_pair->graph_state), current_pair->num);
 
+  mc_stats->visited_pairs++;
 
+  int value;
+  smx_simcall_t req = NULL;
 
+  xbt_automaton_transition_t transition_succ;
+  unsigned int cursor = 0;
+  int res;
+  int visited_num;
+
+  mc_pair_t next_pair = NULL;
+  xbt_dynar_t prop_values = NULL;
+  mc_visited_pair_t reached_pair = NULL;
+  int counter_example_depth = 0;
+
+  if (xbt_fifo_size(mc_stack) < _sg_mc_max_depth) {
+
+    if (current_pair->requests > 0) {
+
+      if (current_pair->search_cycle) {
+
+        if ((current_pair->automaton_state->type == 1)
+            || (current_pair->automaton_state->type == 2)) {
+          if ((reached_pair =
+               is_reached_acceptance_pair(current_pair->num,
+                                          current_pair->automaton_state,
+                                          current_pair->atomic_propositions)) ==
+              NULL) {
+
+            counter_example_depth = xbt_fifo_size(mc_stack);
+            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();
 
+          }
+        }
+      }
 
+      if ((visited_num =
+           is_visited_pair(reached_pair, current_pair->num,
+                           current_pair->automaton_state,
+                           current_pair->atomic_propositions)) != -1) {
+
+        MC_SET_MC_HEAP;
+        if (dot_output != NULL)
+          fprintf(dot_output, "\"%d\" -> \"%d\" [%s];\n",
+                  initial_global_state->prev_pair, visited_num,
+                  initial_global_state->prev_req);
+        MC_SET_STD_HEAP;
+
+      } else {
+
+        while ((req =
+                MC_state_get_request(current_pair->graph_state,
+                                     &value)) != NULL) {
+
+          MC_SET_MC_HEAP;
+          if (dot_output != NULL) {
+            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;
+          }
+          MC_SET_STD_HEAP;
 
-/********************* DPOR without replay *********************/
+          MC_LOG_REQUEST(mc_liveness, req, value);
 
+          MC_SET_MC_HEAP;
+          if (dot_output != NULL) {
+            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);
+          }
+          MC_SET_STD_HEAP;
 
+          MC_state_set_executed_request(current_pair->graph_state, req, value);
+          mc_stats->executed_transitions++;
 
-void MC_dpor_with_restore_snapshot_init(xbt_automaton_t a){
+          /* Answer the request */
+          SIMIX_simcall_handle(req, value);
 
-  XBT_DEBUG("**************************************************");
-  XBT_DEBUG("DPOR (with restore snapshot) init");
-  XBT_DEBUG("**************************************************");
+          /* Wait for requests (schedules processes) */
+          MC_wait_for_requests();
 
-  mc_pairs_t initial_pair;
-  mc_state_t initial_graph_state;
-  smx_process_t process; 
-  mc_snapshot_t initial_snapshot;
-  
-  MC_wait_for_requests();
+          MC_SET_MC_HEAP;
+          prop_values = get_atomic_propositions_values();
+          MC_SET_STD_HEAP;
 
-  MC_SET_RAW_MEM;
+          int new_pair = 0;
 
-  initial_snapshot = xbt_new0(s_mc_snapshot_t, 1);
+          /* Evaluate enabled transition according to atomic propositions values */
+          cursor = 0;
+          xbt_dynar_foreach(current_pair->automaton_state->out, cursor,
+                            transition_succ) {
 
-  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);
-      break;
-    }
-  }
+            res =
+                MC_automaton_evaluate_label(transition_succ->label,
+                                            prop_values);
 
-  visited_pairs = xbt_dynar_new(sizeof(mc_visited_pairs_t), NULL); 
-  reached_pairs = xbt_dynar_new(sizeof(mc_reached_pairs_t), NULL); 
+            if (res == 1) {     // enabled transition in automaton
 
-  MC_take_snapshot(initial_snapshot);
+              if (new_pair)
+                MC_replay_liveness(mc_stack, 1);
 
-  MC_UNSET_RAW_MEM; 
+              MC_SET_MC_HEAP;
 
-  unsigned int cursor = 0;
-  unsigned int cursor2 = 0;
-  xbt_state_t state = NULL;
-  int res;
-  xbt_transition_t transition_succ;
+              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();
 
-  xbt_dynar_foreach(a->states, cursor, state){
-    if(state->type == -1){
-      xbt_dynar_foreach(state->out, cursor2, transition_succ){
-       res = MC_automaton_evaluate_label(a, transition_succ->label);
-       
-       if((res == 1) || (res == 2)){
-        
-         MC_SET_RAW_MEM;
+              /* Get enabled processes and insert them 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_pair->graph_state, process);
+                }
+              }
 
-         initial_pair = new_pair(initial_snapshot, initial_graph_state, transition_succ->dst);
-         xbt_fifo_unshift(mc_snapshot_stack, initial_pair);
+              next_pair->requests =
+                  MC_state_interleave_size(next_pair->graph_state);
 
-         MC_UNSET_RAW_MEM;
-         
-         /* Add pair to list of visited pairs*/
-         set_pair_visited(initial_pair->graph_state, initial_pair->automaton_state, 0);
+              if (next_pair->automaton_state->type == 1
+                  || next_pair->automaton_state->type == 2
+                  || current_pair->search_cycle)
+                next_pair->search_cycle = 1;
 
-         MC_dpor_with_restore_snapshot(a, 0, 0);
+              xbt_fifo_unshift(mc_stack, next_pair);
 
-         break;
+              if (mc_stats->expanded_pairs % 1000000 == 0)
+                XBT_INFO("Expanded pairs : %lu", mc_stats->expanded_pairs);
 
-       }
-      }
-    }
+              MC_SET_STD_HEAP;
 
-    if(xbt_fifo_size(mc_snapshot_stack)>0)
-      break;
-  }
+              new_pair = 1;
 
-  /* Initial state of automaton is initial state for initial pair*/
-  if(xbt_fifo_size(mc_snapshot_stack) == 0){
-
-    cursor = 0;
-    
-    xbt_dynar_foreach(a->states, cursor, state){
-      if(state->type == -1){
-    
-       MC_SET_RAW_MEM;
-
-       initial_pair = new_pair(initial_snapshot, initial_graph_state, state);
-       xbt_fifo_unshift(mc_snapshot_stack, initial_pair);
-    
-       MC_UNSET_RAW_MEM;
-    
-       /* Add pair to list of visited pairs*/
-       set_pair_visited(initial_pair->graph_state, initial_pair->automaton_state, 0);
-    
-       MC_dpor_with_restore_snapshot(a, 0, 0);
-       
-       break;
-      }
-    }
+              MC_modelcheck_liveness();
 
-  }
+            }
 
-}
+          }
 
+          /* Then, evaluate true transitions (always true, whatever atomic propositions values) */
+          cursor = 0;
+          xbt_dynar_foreach(current_pair->automaton_state->out, cursor,
+                            transition_succ) {
 
-/* Test du dpor sans rejeu des états de la pile à partir de l'état initial du système mais avec restauration au milieu de l'arbre -> s'applique aux safety properties mais pas d'intégration pour les liveness */
-void MC_dpor_with_restore_snapshot(xbt_automaton_t a, int search_cycle, int restore){
+            res =
+                MC_automaton_evaluate_label(transition_succ->label,
+                                            prop_values);
 
-  smx_process_t process = NULL;
-  
-  if(xbt_fifo_size(mc_snapshot_stack) == 0)
-    return;
+            if (res == 2) {     // true transition in automaton
 
-  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;
+              if (new_pair)
+                MC_replay_liveness(mc_stack, 1);
 
-  xbt_transition_t transition_succ;
-  unsigned int cursor;
-  int res;
+              MC_SET_MC_HEAP;
 
-  mc_pairs_t next_pair;
-  mc_snapshot_t next_snapshot;
-  mc_snapshot_t current_snapshot;
-  mc_pairs_t current_pair;
-  mc_pairs_t prev_pair;
-  int new_pair_created;
+              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();
 
-  while(xbt_fifo_size(mc_snapshot_stack) > 0){
+              /* 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_pair->graph_state, process);
+                }
+              }
 
-    current_pair = (mc_pairs_t)xbt_fifo_get_item_content(xbt_fifo_get_first_item(mc_snapshot_stack));
+              next_pair->requests =
+                  MC_state_interleave_size(next_pair->graph_state);
 
-    
-    XBT_DEBUG("**************************************************");
-    XBT_DEBUG("Depth : %d, State : g=%p, a=%p(%s), %u interleave", xbt_fifo_size(mc_snapshot_stack),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++;
+              if (next_pair->automaton_state->type == 1
+                  || next_pair->automaton_state->type == 2
+                  || current_pair->search_cycle)
+                next_pair->search_cycle = 1;
 
-    /* Add pair to list of visited pairs*/
-    set_pair_visited(current_pair->graph_state, current_pair->automaton_state, search_cycle);
+              xbt_fifo_unshift(mc_stack, next_pair);
 
-    if((xbt_fifo_size(mc_snapshot_stack) < MAX_DEPTH) && (req = MC_state_get_request(current_pair->graph_state, &value))){
+              if (mc_stats->expanded_pairs % 1000000 == 0)
+                XBT_INFO("Expanded pairs : %lu", mc_stats->expanded_pairs);
 
-      MC_SET_RAW_MEM;
-      current_snapshot = xbt_new0(s_mc_snapshot_t, 1);
-      MC_take_snapshot(current_snapshot);
-      MC_UNSET_RAW_MEM;
+              MC_SET_STD_HEAP;
 
-      /* Debug information */
-      if(XBT_LOG_ISENABLED(mc_liveness, xbt_log_priority_debug)){
-       req_str = MC_request_to_string(req, value);
-       XBT_DEBUG("Execute: %s", req_str);
-       xbt_free(req_str);
-      }
+              new_pair = 1;
 
-      MC_state_set_executed_request(current_pair->graph_state, req, value);
+              MC_modelcheck_liveness();
 
-      /* 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_pair_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;
-       }
-      }
+            }
 
-      next_snapshot = xbt_new0(s_mc_snapshot_t, 1);
-      MC_take_snapshot(next_snapshot);
-      
-      MC_UNSET_RAW_MEM;
+          }
 
-      new_pair_created = 0;
-      
-      cursor = 0;
-      xbt_dynar_foreach(current_pair->automaton_state->out, cursor, transition_succ){
-       res = MC_automaton_evaluate_label(a, transition_succ->label);
-
-       if((res == 1) || (res == 2)){
-        
-         MC_SET_RAW_MEM;
-         next_pair = new_pair(next_snapshot,next_graph_state, transition_succ->dst);
-         xbt_fifo_unshift(mc_snapshot_stack, next_pair);
-         MC_UNSET_RAW_MEM;
-         
-         new_pair_created = 1;
-
-         break;
-        
-       }
-      }
+          if (MC_state_interleave_size(current_pair->graph_state) > 0) {
+            XBT_DEBUG("Backtracking to depth %d", xbt_fifo_size(mc_stack));
+            MC_replay_liveness(mc_stack, 0);
+          }
 
-      /* Pas d'avancement possible dans l'automate */
-      if(new_pair_created == 0){ 
-         
-       MC_SET_RAW_MEM;
-       next_pair = new_pair(next_snapshot,next_graph_state,current_pair->automaton_state);
-       xbt_fifo_unshift(mc_snapshot_stack, next_pair);
-       MC_UNSET_RAW_MEM;
+        }
 
       }
-      
-    }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;
-
-      while((current_pair = xbt_fifo_shift(mc_snapshot_stack)) != NULL){
-       req = MC_state_get_internal_request(current_pair->graph_state);
-       xbt_fifo_foreach(mc_snapshot_stack, item, prev_pair, mc_pairs_t) {
-          if(MC_request_depend(req, MC_state_get_internal_request(prev_pair->graph_state))){
-            if(XBT_LOG_ISENABLED(mc_liveness, 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 (state=%p)", req_str, prev_pair->graph_state);
-              xbt_free(req_str);
-              prev_req = MC_state_get_executed_request(current_pair->graph_state, &value);
-              req_str = MC_request_to_string(prev_req, value);
-              XBT_DEBUG("%s (state=%p)", req_str, current_pair->graph_state);
-              xbt_free(req_str);              
-            }
 
-            if(!MC_state_process_is_done(prev_pair->graph_state, req->issuer)){
-              MC_state_interleave_process(prev_pair->graph_state, req->issuer);
-             
-           } else {
-              XBT_DEBUG("Process %p is in done set", req->issuer);
-           }
+    }
 
-            break;
-          }
-        }
+  } else {
 
-       if(MC_state_interleave_size(current_pair->graph_state)){
-         MC_restore_snapshot(current_pair->system_state);
-         xbt_fifo_unshift(mc_snapshot_stack, current_pair);
-         XBT_DEBUG("Back-tracking to depth %d", xbt_fifo_size(mc_snapshot_stack));
-         MC_UNSET_RAW_MEM;
-         break;
-       }
-      }
+    XBT_WARN("/!\\ Max depth reached ! /!\\ ");
+    if (MC_state_interleave_size(current_pair->graph_state) > 0) {
+      XBT_WARN
+          ("/!\\ But, there are still processes to interleave. Model-checker will not be able to ensure the soundness of the verification from now. /!\\ ");
+      if (_sg_mc_max_depth == 1000)
+        XBT_WARN
+            ("Notice : the default value of max depth is 1000 but you can change it with cfg=model-check/max_depth:value.");
+    }
 
-      MC_UNSET_RAW_MEM;
+  }
 
-    } 
+  if (xbt_fifo_size(mc_stack) == _sg_mc_max_depth) {
+    XBT_DEBUG("Pair %d (depth = %d) shifted in stack, maximum depth reached",
+              current_pair->num, xbt_fifo_size(mc_stack));
+  } else {
+    XBT_DEBUG("Pair %d (depth = %d) shifted in stack", current_pair->num,
+              xbt_fifo_size(mc_stack));
   }
-  MC_UNSET_RAW_MEM;
-  return;
+
+
+  MC_SET_MC_HEAP;
+  xbt_dynar_free(&prop_values);
+  current_pair = xbt_fifo_shift(mc_stack);
+  if (xbt_fifo_size(mc_stack) != _sg_mc_max_depth - 1
+      && current_pair->requests > 0 && current_pair->search_cycle) {
+    remove_acceptance_pair(current_pair->num);
+  }
+  MC_pair_delete(current_pair);
+
+  MC_SET_STD_HEAP;
+
 }