src/mc/mc_diff.cpp

   1 /* mc_diff - Memory snapshooting and comparison                             */
   2
   3 /* Copyright (c) 2008-2014. The SimGrid Team.
   4  * All rights reserved.                                                     */
   5
   6 /* This program is free software; you can redistribute it and/or modify it
   7  * under the terms of the license (GNU LGPL) which comes with this package. */
   8
   9 #include "xbt/ex_interface.h"   /* internals of backtrace setup */
  10 #include "xbt/str.h"
  11 #include "mc/mc.h"
  12 #include "xbt/mmalloc.h"
  13 #include "mc_object_info.h"
  14 #include "mc/datatypes.h"
  15 #include "mc/mc_private.h"
  16 #include "mc/mc_snapshot.h"
  17
  18 using simgrid::mc::remote;
  19
  20 extern "C" {
  21
  22 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(mc_diff, xbt,
  23                                 "Logging specific to mc_diff in mc");
  24
  25 xbt_dynar_t mc_heap_comparison_ignore;
  26 xbt_dynar_t stacks_areas;
  27
  28
  29
  30 /********************************* Backtrace ***********************************/
  31 /******************************************************************************/
  32
  33 static void mmalloc_backtrace_block_display(void *heapinfo, int block)
  34 {
  35
  36   /* xbt_ex_t e; */
  37
  38   /* if (((malloc_info *)heapinfo)[block].busy_block.bt_size == 0) { */
  39   /*   fprintf(stderr, "No backtrace available for that block, sorry.\n"); */
  40   /*   return; */
  41   /* } */
  42
  43   /* memcpy(&e.bt,&(((malloc_info *)heapinfo)[block].busy_block.bt),sizeof(void*)*XBT_BACKTRACE_SIZE); */
  44   /* e.used = ((malloc_info *)heapinfo)[block].busy_block.bt_size; */
  45
  46   /* xbt_ex_setup_backtrace(&e); */
  47   /* if (e.used == 0) { */
  48   /*   fprintf(stderr, "(backtrace not set)\n"); */
  49   /* } else if (e.bt_strings == NULL) { */
  50   /*   fprintf(stderr, "(backtrace not ready to be computed. %s)\n",xbt_binary_name?"Dunno why":"xbt_binary_name not setup yet"); */
  51   /* } else { */
  52   /*   int i; */
  53
  54   /*   fprintf(stderr, "Backtrace of where the block %d was malloced (%d frames):\n", block ,e.used); */
  55   /*   for (i = 0; i < e.used; i++)       /\* no need to display "xbt_backtrace_display" *\/{ */
  56   /*     fprintf(stderr, "%d ---> %s\n",i, e.bt_strings[i] + 4); */
  57   /*   } */
  58   /* } */
  59 }
  60
  61 static void mmalloc_backtrace_fragment_display(void *heapinfo, int block,
  62                                                int frag)
  63 {
  64
  65   /* xbt_ex_t e; */
  66
  67   /* memcpy(&e.bt,&(((malloc_info *)heapinfo)[block].busy_frag.bt[frag]),sizeof(void*)*XBT_BACKTRACE_SIZE); */
  68   /* e.used = XBT_BACKTRACE_SIZE; */
  69
  70   /* xbt_ex_setup_backtrace(&e); */
  71   /* if (e.used == 0) { */
  72   /*   fprintf(stderr, "(backtrace not set)\n"); */
  73   /* } else if (e.bt_strings == NULL) { */
  74   /*   fprintf(stderr, "(backtrace not ready to be computed. %s)\n",xbt_binary_name?"Dunno why":"xbt_binary_name not setup yet"); */
  75   /* } else { */
  76   /*   int i; */
  77
  78   /*   fprintf(stderr, "Backtrace of where the fragment %d in block %d was malloced (%d frames):\n", frag, block ,e.used); */
  79   /*   for (i = 0; i < e.used; i++)       /\* no need to display "xbt_backtrace_display" *\/{ */
  80   /*     fprintf(stderr, "%d ---> %s\n",i, e.bt_strings[i] + 4); */
  81   /*   } */
  82   /* } */
  83
  84 }
  85
  86 static void mmalloc_backtrace_display(void *addr)
  87 {
  88
  89   /* size_t block, frag_nb; */
  90   /* int type; */
  91
  92   /* block = (((char*) (addr) - (char*) heap -> heapbase) / BLOCKSIZE + 1); */
  93
  94   /* type = heap->heapinfo[block].type; */
  95
  96   /* switch(type){ */
  97   /* case MMALLOC_TYPE_HEAPINFO :  */
  98   /* case MMALLOC_TYPE_FREE : /\* Free block *\/ */
  99   /*   fprintf(stderr, "Asked to display the backtrace of a block that is free. I'm puzzled\n"); */
 100   /*   xbt_abort(); */
 101   /*   break;  */
 102   /* case 0: /\* Large block *\/ */
 103   /*   mmalloc_backtrace_block_display(heap->heapinfo, block); */
 104   /*   break; */
 105   /* default: /\* Fragmented block *\/ */
 106   /*   frag_nb = RESIDUAL(addr, BLOCKSIZE) >> type; */
 107   /*   if(heap->heapinfo[block].busy_frag.frag_size[frag_nb] == -1){ */
 108   /*     fprintf(stderr , "Asked to display the backtrace of a fragment that is free. I'm puzzled\n"); */
 109   /*     xbt_abort(); */
 110   /*   } */
 111   /*   mmalloc_backtrace_fragment_display(heap->heapinfo, block, frag_nb); */
 112   /*   break; */
 113   /* } */
 114 }
 115
 116
 117 static int compare_backtrace(int b1, int f1, int b2, int f2)
 118 {
 119   /*int i = 0;
 120      if(f1 != -1){
 121      for(i=0; i< XBT_BACKTRACE_SIZE; i++){
 122      if(heapinfo1[b1].busy_frag.bt[f1][i] != heapinfo2[b2].busy_frag.bt[f2][i]){
 123      //mmalloc_backtrace_fragment_display((void*)heapinfo1, b1, f1);
 124      //mmalloc_backtrace_fragment_display((void*)heapinfo2, b2, f2);
 125      return 1;
 126      }
 127      }
 128      }else{
 129      for(i=0; i< heapinfo1[b1].busy_block.bt_size; i++){
 130      if(heapinfo1[b1].busy_block.bt[i] != heapinfo2[b2].busy_block.bt[i]){
 131      //mmalloc_backtrace_block_display((void*)heapinfo1, b1);
 132      //mmalloc_backtrace_block_display((void*)heapinfo2, b2);
 133      return 1;
 134      }
 135      }
 136      } */
 137   return 0;
 138 }
 139
 140
 141 /*********************************** Heap comparison ***********************************/
 142 /***************************************************************************************/
 143
 144 typedef char *type_name;
 145
 146 struct s_mc_diff {
 147   s_xbt_mheap_t std_heap_copy;
 148   size_t heaplimit;
 149   // Number of blocks in the heaps:
 150   size_t heapsize1, heapsize2;
 151   std::vector<s_mc_heap_ignore_region_t>* to_ignore1;
 152   std::vector<s_mc_heap_ignore_region_t>* to_ignore2;
 153   s_heap_area_t *equals_to1, *equals_to2;
 154   mc_type_t *types1, *types2;
 155   size_t available;
 156 };
 157
 158 #define equals_to1_(i,j) equals_to1[ MAX_FRAGMENT_PER_BLOCK*(i) + (j)]
 159 #define equals_to2_(i,j) equals_to2[ MAX_FRAGMENT_PER_BLOCK*(i) + (j)]
 160 #define types1_(i,j) types1[ MAX_FRAGMENT_PER_BLOCK*(i) + (j)]
 161 #define types2_(i,j) types2[ MAX_FRAGMENT_PER_BLOCK*(i) + (j)]
 162
 163 __thread struct s_mc_diff *mc_diff_info = NULL;
 164
 165 /*********************************** Free functions ************************************/
 166
 167 static void heap_area_pair_free(heap_area_pair_t pair)
 168 {
 169   xbt_free(pair);
 170   pair = NULL;
 171 }
 172
 173 static void heap_area_pair_free_voidp(void *d)
 174 {
 175   heap_area_pair_free((heap_area_pair_t) * (void **) d);
 176 }
 177
 178 static void heap_area_free(heap_area_t area)
 179 {
 180   xbt_free(area);
 181   area = NULL;
 182 }
 183
 184 /************************************************************************************/
 185
 186 static s_heap_area_t make_heap_area(int block, int fragment)
 187 {
 188   s_heap_area_t area;
 189   area.valid = 1;
 190   area.block = block;
 191   area.fragment = fragment;
 192   return area;
 193 }
 194
 195
 196 static int is_new_heap_area_pair(xbt_dynar_t list, int block1, int fragment1,
 197                                  int block2, int fragment2)
 198 {
 199
 200   unsigned int cursor = 0;
 201   heap_area_pair_t current_pair;
 202
 203   xbt_dynar_foreach(list, cursor, current_pair) {
 204     if (current_pair->block1 == block1 && current_pair->block2 == block2
 205         && current_pair->fragment1 == fragment1
 206         && current_pair->fragment2 == fragment2)
 207       return 0;
 208   }
 209
 210   return 1;
 211 }
 212
 213 static int add_heap_area_pair(xbt_dynar_t list, int block1, int fragment1,
 214                               int block2, int fragment2)
 215 {
 216
 217   if (is_new_heap_area_pair(list, block1, fragment1, block2, fragment2)) {
 218     heap_area_pair_t pair = NULL;
 219     pair = xbt_new0(s_heap_area_pair_t, 1);
 220     pair->block1 = block1;
 221     pair->fragment1 = fragment1;
 222     pair->block2 = block2;
 223     pair->fragment2 = fragment2;
 224
 225     xbt_dynar_push(list, &pair);
 226
 227     return 1;
 228   }
 229
 230   return 0;
 231 }
 232
 233 static ssize_t heap_comparison_ignore_size(std::vector<s_mc_heap_ignore_region_t>* ignore_list,
 234                                            const void *address)
 235 {
 236   int start = 0;
 237   int end = ignore_list->size() - 1;
 238
 239   while (start <= end) {
 240     unsigned int cursor = (start + end) / 2;
 241     s_mc_heap_ignore_region_t region = (*ignore_list)[cursor];
 242     if (region.address == address)
 243       return region.size;
 244     if (region.address < address)
 245       start = cursor + 1;
 246     if (region.address > address)
 247       end = cursor - 1;
 248   }
 249
 250   return -1;
 251 }
 252
 253 static int is_stack(const void *address)
 254 {
 255   unsigned int cursor = 0;
 256   stack_region_t stack;
 257
 258   xbt_dynar_foreach(stacks_areas, cursor, stack) {
 259     if (address == stack->address)
 260       return 1;
 261   }
 262
 263   return 0;
 264 }
 265
 266 // TODO, this should depend on the snapshot?
 267 static int is_block_stack(int block)
 268 {
 269   unsigned int cursor = 0;
 270   stack_region_t stack;
 271
 272   xbt_dynar_foreach(stacks_areas, cursor, stack) {
 273     if (block == stack->block)
 274       return 1;
 275   }
 276
 277   return 0;
 278 }
 279
 280 static void match_equals(struct s_mc_diff *state, xbt_dynar_t list)
 281 {
 282
 283   unsigned int cursor = 0;
 284   heap_area_pair_t current_pair;
 285
 286   xbt_dynar_foreach(list, cursor, current_pair) {
 287
 288     if (current_pair->fragment1 != -1) {
 289
 290       state->equals_to1_(current_pair->block1, current_pair->fragment1) =
 291           make_heap_area(current_pair->block2, current_pair->fragment2);
 292       state->equals_to2_(current_pair->block2, current_pair->fragment2) =
 293           make_heap_area(current_pair->block1, current_pair->fragment1);
 294
 295     } else {
 296
 297       state->equals_to1_(current_pair->block1, 0) =
 298           make_heap_area(current_pair->block2, current_pair->fragment2);
 299       state->equals_to2_(current_pair->block2, 0) =
 300           make_heap_area(current_pair->block1, current_pair->fragment1);
 301
 302     }
 303
 304   }
 305 }
 306
 307 /** Check whether two blocks are known to be matching
 308  *
 309  *  @param state  State used
 310  *  @param b1     Block of state 1
 311  *  @param b2     Block of state 2
 312  *  @return       if the blocks are known to be matching
 313  */
 314 static int equal_blocks(struct s_mc_diff *state, int b1, int b2)
 315 {
 316
 317   if (state->equals_to1_(b1, 0).block == b2
 318       && state->equals_to2_(b2, 0).block == b1)
 319     return 1;
 320
 321   return 0;
 322 }
 323
 324 /** Check whether two fragments are known to be matching
 325  *
 326  *  @param state  State used
 327  *  @param b1     Block of state 1
 328  *  @param f1     Fragment of state 1
 329  *  @param b2     Block of state 2
 330  *  @param f2     Fragment of state 2
 331  *  @return       if the fragments are known to be matching
 332  */
 333 static int equal_fragments(struct s_mc_diff *state, int b1, int f1, int b2,
 334                            int f2)
 335 {
 336
 337   if (state->equals_to1_(b1, f1).block == b2
 338       && state->equals_to1_(b1, f1).fragment == f2
 339       && state->equals_to2_(b2, f2).block == b1
 340       && state->equals_to2_(b2, f2).fragment == f1)
 341     return 1;
 342
 343   return 0;
 344 }
 345
 346 }
 347
 348 int init_heap_information(xbt_mheap_t heap1, xbt_mheap_t heap2,
 349                           std::vector<s_mc_heap_ignore_region_t>* i1,
 350                           std::vector<s_mc_heap_ignore_region_t>* i2)
 351 {
 352   if (mc_diff_info == NULL) {
 353     mc_diff_info = xbt_new0(struct s_mc_diff, 1);
 354     mc_diff_info->equals_to1 = NULL;
 355     mc_diff_info->equals_to2 = NULL;
 356     mc_diff_info->types1 = NULL;
 357     mc_diff_info->types2 = NULL;
 358   }
 359   struct s_mc_diff *state = mc_diff_info;
 360
 361   if ((((struct mdesc *) heap1)->heaplimit !=
 362        ((struct mdesc *) heap2)->heaplimit)
 363       ||
 364       ((((struct mdesc *) heap1)->heapsize !=
 365         ((struct mdesc *) heap2)->heapsize)))
 366     return -1;
 367
 368   state->heaplimit = ((struct mdesc *) heap1)->heaplimit;
 369
 370   state->std_heap_copy = *mc_model_checker->process().get_heap();
 371
 372   state->heapsize1 = heap1->heapsize;
 373   state->heapsize2 = heap2->heapsize;
 374
 375   state->to_ignore1 = i1;
 376   state->to_ignore2 = i2;
 377
 378   if (state->heaplimit > state->available) {
 379     state->equals_to1 = (s_heap_area_t*)
 380         realloc(state->equals_to1,
 381                 state->heaplimit * MAX_FRAGMENT_PER_BLOCK *
 382                 sizeof(s_heap_area_t));
 383     state->types1 = (simgrid::mc::Type**)
 384         realloc(state->types1,
 385                 state->heaplimit * MAX_FRAGMENT_PER_BLOCK *
 386                 sizeof(type_name *));
 387     state->equals_to2 = (s_heap_area_t*)
 388         realloc(state->equals_to2,
 389                 state->heaplimit * MAX_FRAGMENT_PER_BLOCK *
 390                 sizeof(s_heap_area_t));
 391     state->types2 = (simgrid::mc::Type**)
 392         realloc(state->types2,
 393                 state->heaplimit * MAX_FRAGMENT_PER_BLOCK *
 394                 sizeof(type_name *));
 395     state->available = state->heaplimit;
 396   }
 397
 398   memset(state->equals_to1, 0,
 399          state->heaplimit * MAX_FRAGMENT_PER_BLOCK * sizeof(s_heap_area_t));
 400   memset(state->equals_to2, 0,
 401          state->heaplimit * MAX_FRAGMENT_PER_BLOCK * sizeof(s_heap_area_t));
 402   memset(state->types1, 0,
 403          state->heaplimit * MAX_FRAGMENT_PER_BLOCK * sizeof(type_name *));
 404   memset(state->types2, 0,
 405          state->heaplimit * MAX_FRAGMENT_PER_BLOCK * sizeof(type_name *));
 406
 407   return 0;
 408
 409 }
 410
 411 extern "C" {
 412
 413 void reset_heap_information()
 414 {
 415
 416 }
 417
 418 // TODO, have a robust way to find it in O(1)
 419 static inline
 420 mc_mem_region_t MC_get_heap_region(mc_snapshot_t snapshot)
 421 {
 422   size_t n = snapshot->snapshot_regions.size();
 423   for (size_t i=0; i!=n; ++i) {
 424     mc_mem_region_t region = snapshot->snapshot_regions[i].get();
 425     if (region->region_type() == simgrid::mc::RegionType::Heap)
 426       return region;
 427   }
 428   xbt_die("No heap region");
 429 }
 430
 431 int mmalloc_compare_heap(mc_snapshot_t snapshot1, mc_snapshot_t snapshot2)
 432 {
 433   mc_process_t process = &mc_model_checker->process();
 434   struct s_mc_diff *state = mc_diff_info;
 435
 436   /* Start comparison */
 437   size_t i1, i2, j1, j2, k;
 438   void *addr_block1, *addr_block2, *addr_frag1, *addr_frag2;
 439   int nb_diff1 = 0, nb_diff2 = 0;
 440
 441   int equal, res_compare = 0;
 442
 443   /* Check busy blocks */
 444
 445   i1 = 1;
 446
 447   malloc_info heapinfo_temp1, heapinfo_temp2;
 448   malloc_info heapinfo_temp2b;
 449
 450   mc_mem_region_t heap_region1 = MC_get_heap_region(snapshot1);
 451   mc_mem_region_t heap_region2 = MC_get_heap_region(snapshot2);
 452
 453   // This is the address of std_heap->heapinfo in the application process:
 454   void* heapinfo_address = &((xbt_mheap_t) process->heap_address)->heapinfo;
 455
 456   // This is in snapshot do not use them directly:
 457   const malloc_info* heapinfos1 = snapshot1->read<malloc_info*>(
 458     (std::uint64_t)heapinfo_address, simgrid::mc::ProcessIndexMissing);
 459   const malloc_info* heapinfos2 = snapshot2->read<malloc_info*>(
 460     (std::uint64_t)heapinfo_address, simgrid::mc::ProcessIndexMissing);
 461
 462   while (i1 <= state->heaplimit) {
 463
 464     const malloc_info* heapinfo1 = (const malloc_info*) MC_region_read(heap_region1, &heapinfo_temp1, &heapinfos1[i1], sizeof(malloc_info));
 465     const malloc_info* heapinfo2 = (const malloc_info*) MC_region_read(heap_region2, &heapinfo_temp2, &heapinfos2[i1], sizeof(malloc_info));
 466
 467     if (heapinfo1->type == MMALLOC_TYPE_FREE || heapinfo1->type == MMALLOC_TYPE_HEAPINFO) {      /* Free block */
 468       i1 ++;
 469       continue;
 470     }
 471
 472     if (heapinfo1->type < 0) {
 473       fprintf(stderr, "Unkown mmalloc block type.\n");
 474       abort();
 475     }
 476
 477     addr_block1 =
 478         ((void *) (((ADDR2UINT(i1)) - 1) * BLOCKSIZE +
 479                    (char *) state->std_heap_copy.heapbase));
 480
 481     if (heapinfo1->type == MMALLOC_TYPE_UNFRAGMENTED) {       /* Large block */
 482
 483       if (is_stack(addr_block1)) {
 484         for (k = 0; k < heapinfo1->busy_block.size; k++)
 485           state->equals_to1_(i1 + k, 0) = make_heap_area(i1, -1);
 486         for (k = 0; k < heapinfo2->busy_block.size; k++)
 487           state->equals_to2_(i1 + k, 0) = make_heap_area(i1, -1);
 488         i1 += heapinfo1->busy_block.size;
 489         continue;
 490       }
 491
 492       if (state->equals_to1_(i1, 0).valid) {
 493         i1++;
 494         continue;
 495       }
 496
 497       i2 = 1;
 498       equal = 0;
 499       res_compare = 0;
 500
 501       /* Try first to associate to same block in the other heap */
 502       if (heapinfo2->type == heapinfo1->type) {
 503
 504         if (state->equals_to2_(i1, 0).valid == 0) {
 505
 506           addr_block2 = (ADDR2UINT(i1) - 1) * BLOCKSIZE +
 507                          (char *) state->std_heap_copy.heapbase;
 508
 509           res_compare =
 510               compare_heap_area(simgrid::mc::ProcessIndexMissing, addr_block1, addr_block2, snapshot1, snapshot2,
 511                                 NULL, NULL, 0);
 512
 513           if (res_compare != 1) {
 514             for (k = 1; k < heapinfo2->busy_block.size; k++)
 515               state->equals_to2_(i1 + k, 0) = make_heap_area(i1, -1);
 516             for (k = 1; k < heapinfo1->busy_block.size; k++)
 517               state->equals_to1_(i1 + k, 0) = make_heap_area(i1, -1);
 518             equal = 1;
 519             i1 += heapinfo1->busy_block.size;
 520           }
 521
 522         }
 523
 524       }
 525
 526       while (i2 <= state->heaplimit && !equal) {
 527
 528         addr_block2 = (ADDR2UINT(i2) - 1) * BLOCKSIZE +
 529                        (char *) state->std_heap_copy.heapbase;
 530
 531         if (i2 == i1) {
 532           i2++;
 533           continue;
 534         }
 535
 536         const malloc_info* heapinfo2b = (const malloc_info*) MC_region_read(heap_region2, &heapinfo_temp2b, &heapinfos2[i2], sizeof(malloc_info));
 537
 538         if (heapinfo2b->type != MMALLOC_TYPE_UNFRAGMENTED) {
 539           i2++;
 540           continue;
 541         }
 542
 543         if (state->equals_to2_(i2, 0).valid) {
 544           i2++;
 545           continue;
 546         }
 547
 548         res_compare =
 549             compare_heap_area(simgrid::mc::ProcessIndexMissing, addr_block1, addr_block2, snapshot1, snapshot2,
 550                               NULL, NULL, 0);
 551
 552         if (res_compare != 1) {
 553           for (k = 1; k < heapinfo2b->busy_block.size; k++)
 554             state->equals_to2_(i2 + k, 0) = make_heap_area(i1, -1);
 555           for (k = 1; k < heapinfo1->busy_block.size; k++)
 556             state->equals_to1_(i1 + k, 0) = make_heap_area(i2, -1);
 557           equal = 1;
 558           i1 += heapinfo1->busy_block.size;
 559         }
 560
 561         i2++;
 562
 563       }
 564
 565       if (!equal) {
 566         XBT_DEBUG("Block %zu not found (size_used = %zu, addr = %p)", i1,
 567                   heapinfo1->busy_block.busy_size, addr_block1);
 568         i1 = state->heaplimit + 1;
 569         nb_diff1++;
 570         //i1++;
 571       }
 572
 573     } else {                    /* Fragmented block */
 574
 575       for (j1 = 0; j1 < (size_t) (BLOCKSIZE >> heapinfo1->type); j1++) {
 576
 577         if (heapinfo1->busy_frag.frag_size[j1] == -1) /* Free fragment */
 578           continue;
 579
 580         if (state->equals_to1_(i1, j1).valid)
 581           continue;
 582
 583         addr_frag1 =
 584             (void *) ((char *) addr_block1 + (j1 << heapinfo1->type));
 585
 586         i2 = 1;
 587         equal = 0;
 588
 589         /* Try first to associate to same fragment in the other heap */
 590         if (heapinfo2->type == heapinfo1->type) {
 591
 592           if (state->equals_to2_(i1, j1).valid == 0) {
 593
 594             addr_block2 = (ADDR2UINT(i1) - 1) * BLOCKSIZE +
 595                            (char *) state->std_heap_copy.heapbase;
 596             addr_frag2 =
 597                 (void *) ((char *) addr_block2 +
 598                           (j1 << heapinfo2->type));
 599
 600             res_compare =
 601                 compare_heap_area(simgrid::mc::ProcessIndexMissing, addr_frag1, addr_frag2, snapshot1, snapshot2,
 602                                   NULL, NULL, 0);
 603
 604             if (res_compare != 1)
 605               equal = 1;
 606
 607           }
 608
 609         }
 610
 611         while (i2 <= state->heaplimit && !equal) {
 612
 613           const malloc_info* heapinfo2b = (const malloc_info*) MC_region_read(
 614             heap_region2, &heapinfo_temp2b, &heapinfos2[i2],
 615             sizeof(malloc_info));
 616
 617           if (heapinfo2b->type == MMALLOC_TYPE_FREE || heapinfo2b->type == MMALLOC_TYPE_HEAPINFO) {
 618             i2 ++;
 619             continue;
 620           }
 621
 622           // We currently do not match fragments with unfragmented blocks (maybe we should).
 623           if (heapinfo2b->type == MMALLOC_TYPE_UNFRAGMENTED) {
 624             i2++;
 625             continue;
 626           }
 627
 628           if (heapinfo2b->type < 0) {
 629             fprintf(stderr, "Unkown mmalloc block type.\n");
 630             abort();
 631           }
 632
 633           for (j2 = 0; j2 < (size_t) (BLOCKSIZE >> heapinfo2b->type);
 634                j2++) {
 635
 636             if (i2 == i1 && j2 == j1)
 637               continue;
 638
 639             if (state->equals_to2_(i2, j2).valid)
 640               continue;
 641
 642             addr_block2 = (ADDR2UINT(i2) - 1) * BLOCKSIZE +
 643                            (char *) state->std_heap_copy.heapbase;
 644             addr_frag2 =
 645                 (void *) ((char *) addr_block2 +
 646                           (j2 << heapinfo2b->type));
 647
 648             res_compare =
 649                 compare_heap_area(simgrid::mc::ProcessIndexMissing, addr_frag1, addr_frag2, snapshot2, snapshot2,
 650                                   NULL, NULL, 0);
 651
 652             if (res_compare != 1) {
 653               equal = 1;
 654               break;
 655             }
 656
 657           }
 658
 659           i2++;
 660
 661         }
 662
 663         if (!equal) {
 664           XBT_DEBUG
 665               ("Block %zu, fragment %zu not found (size_used = %zd, address = %p)\n",
 666                i1, j1, heapinfo1->busy_frag.frag_size[j1],
 667                addr_frag1);
 668           i2 = state->heaplimit + 1;
 669           i1 = state->heaplimit + 1;
 670           nb_diff1++;
 671           break;
 672         }
 673
 674       }
 675
 676       i1++;
 677
 678     }
 679
 680   }
 681
 682   /* All blocks/fragments are equal to another block/fragment ? */
 683   size_t i = 1, j = 0;
 684
 685   for(i = 1; i <= state->heaplimit; i++) {
 686     const malloc_info* heapinfo1 = (const malloc_info*) MC_region_read(
 687       heap_region1, &heapinfo_temp1, &heapinfos1[i], sizeof(malloc_info));
 688     if (heapinfo1->type == MMALLOC_TYPE_UNFRAGMENTED) {
 689       if (i1 == state->heaplimit) {
 690         if (heapinfo1->busy_block.busy_size > 0) {
 691           if (state->equals_to1_(i, 0).valid == 0) {
 692             if (XBT_LOG_ISENABLED(mc_diff, xbt_log_priority_debug)) {
 693               // TODO, add address
 694               XBT_DEBUG("Block %zu not found (size used = %zu)", i,
 695                         heapinfo1->busy_block.busy_size);
 696               //mmalloc_backtrace_block_display((void*)heapinfo1, i);
 697             }
 698             nb_diff1++;
 699           }
 700         }
 701       }
 702     }
 703     if (heapinfo1->type > 0) {
 704       for (j = 0; j < (size_t) (BLOCKSIZE >> heapinfo1->type); j++) {
 705         if (i1 == state->heaplimit) {
 706           if (heapinfo1->busy_frag.frag_size[j] > 0) {
 707             if (state->equals_to1_(i, j).valid == 0) {
 708               if (XBT_LOG_ISENABLED(mc_diff, xbt_log_priority_debug)) {
 709                 // TODO, print fragment address
 710                 XBT_DEBUG
 711                     ("Block %zu, Fragment %zu not found (size used = %zd)",
 712                      i, j,
 713                      heapinfo1->busy_frag.frag_size[j]);
 714                 //mmalloc_backtrace_fragment_display((void*)heapinfo1, i, j);
 715               }
 716               nb_diff1++;
 717             }
 718           }
 719         }
 720       }
 721     }
 722   }
 723
 724   if (i1 == state->heaplimit)
 725     XBT_DEBUG("Number of blocks/fragments not found in heap1 : %d", nb_diff1);
 726
 727   for (i=1; i <= state->heaplimit; i++) {
 728     const malloc_info* heapinfo2 = (const malloc_info*) MC_region_read(
 729       heap_region2, &heapinfo_temp2, &heapinfos2[i], sizeof(malloc_info));
 730     if (heapinfo2->type == MMALLOC_TYPE_UNFRAGMENTED) {
 731       if (i1 == state->heaplimit) {
 732         if (heapinfo2->busy_block.busy_size > 0) {
 733           if (state->equals_to2_(i, 0).valid == 0) {
 734             if (XBT_LOG_ISENABLED(mc_diff, xbt_log_priority_debug)) {
 735               // TODO, print address of the block
 736               XBT_DEBUG("Block %zu not found (size used = %zu)", i,
 737                         heapinfo2->busy_block.busy_size);
 738               //mmalloc_backtrace_block_display((void*)heapinfo2, i);
 739             }
 740             nb_diff2++;
 741           }
 742         }
 743       }
 744     }
 745     if (heapinfo2->type > 0) {
 746       for (j = 0; j < (size_t) (BLOCKSIZE >> heapinfo2->type); j++) {
 747         if (i1 == state->heaplimit) {
 748           if (heapinfo2->busy_frag.frag_size[j] > 0) {
 749             if (state->equals_to2_(i, j).valid == 0) {
 750               if (XBT_LOG_ISENABLED(mc_diff, xbt_log_priority_debug)) {
 751                 // TODO, print address of the block
 752                 XBT_DEBUG
 753                     ("Block %zu, Fragment %zu not found (size used = %zd)",
 754                      i, j,
 755                      heapinfo2->busy_frag.frag_size[j]);
 756                 //mmalloc_backtrace_fragment_display((void*)heapinfo2, i, j);
 757               }
 758               nb_diff2++;
 759             }
 760           }
 761         }
 762       }
 763     }
 764   }
 765
 766   if (i1 == state->heaplimit)
 767     XBT_DEBUG("Number of blocks/fragments not found in heap2 : %d", nb_diff2);
 768
 769   return ((nb_diff1 > 0) || (nb_diff2 > 0));
 770 }
 771
 772 /**
 773  *
 774  * @param state
 775  * @param real_area1     Process address for state 1
 776  * @param real_area2     Process address for state 2
 777  * @param snapshot1      Snapshot of state 1
 778  * @param snapshot2      Snapshot of state 2
 779  * @param previous
 780  * @param size
 781  * @param check_ignore
 782  */
 783 static int compare_heap_area_without_type(struct s_mc_diff *state, int process_index,
 784                                           const void *real_area1, const void *real_area2,
 785                                           mc_snapshot_t snapshot1,
 786                                           mc_snapshot_t snapshot2,
 787                                           xbt_dynar_t previous, int size,
 788                                           int check_ignore)
 789 {
 790   mc_process_t process = &mc_model_checker->process();
 791
 792   int i = 0;
 793   const void *addr_pointed1, *addr_pointed2;
 794   int pointer_align, res_compare;
 795   ssize_t ignore1, ignore2;
 796
 797   mc_mem_region_t heap_region1 = MC_get_heap_region(snapshot1);
 798   mc_mem_region_t heap_region2 = MC_get_heap_region(snapshot2);
 799
 800   while (i < size) {
 801
 802     if (check_ignore > 0) {
 803       if ((ignore1 =
 804            heap_comparison_ignore_size(state->to_ignore1,
 805                                        (char *) real_area1 + i)) != -1) {
 806         if ((ignore2 =
 807              heap_comparison_ignore_size(state->to_ignore2,
 808                                          (char *) real_area2 + i)) == ignore1) {
 809           if (ignore1 == 0) {
 810             check_ignore--;
 811             return 0;
 812           } else {
 813             i = i + ignore2;
 814             check_ignore--;
 815             continue;
 816           }
 817         }
 818       }
 819     }
 820
 821     if (MC_snapshot_region_memcmp(((char *) real_area1) + i, heap_region1, ((char *) real_area2) + i, heap_region2, 1) != 0) {
 822
 823       pointer_align = (i / sizeof(void *)) * sizeof(void *);
 824       addr_pointed1 = snapshot1->read(
 825         remote((void**)((char *) real_area1 + pointer_align)), process_index);
 826       addr_pointed2 = snapshot2->read(
 827         remote((void**)((char *) real_area2 + pointer_align)), process_index);
 828
 829       if (process->in_maestro_stack(remote(addr_pointed1))
 830         && process->in_maestro_stack(remote(addr_pointed2))) {
 831         i = pointer_align + sizeof(void *);
 832         continue;
 833       } else if (addr_pointed1 > state->std_heap_copy.heapbase
 834                  && addr_pointed1 < mc_snapshot_get_heap_end(snapshot1)
 835                  && addr_pointed2 > state->std_heap_copy.heapbase
 836                  && addr_pointed2 < mc_snapshot_get_heap_end(snapshot2)) {
 837         // Both addreses are in the heap:
 838         res_compare =
 839             compare_heap_area(process_index, addr_pointed1, addr_pointed2, snapshot1,
 840                               snapshot2, previous, NULL, 0);
 841         if (res_compare == 1) {
 842           return res_compare;
 843         }
 844         i = pointer_align + sizeof(void *);
 845         continue;
 846       } else {
 847         return 1;
 848       }
 849
 850     }
 851
 852     i++;
 853
 854   }
 855
 856   return 0;
 857
 858 }
 859
 860 /**
 861  *
 862  * @param state
 863  * @param real_area1     Process address for state 1
 864  * @param real_area2     Process address for state 2
 865  * @param snapshot1      Snapshot of state 1
 866  * @param snapshot2      Snapshot of state 2
 867  * @param previous
 868  * @param type_id
 869  * @param area_size      either a byte_size or an elements_count (?)
 870  * @param check_ignore
 871  * @param pointer_level
 872  * @return               0 (same), 1 (different), -1 (unknown)
 873  */
 874 static int compare_heap_area_with_type(struct s_mc_diff *state, int process_index,
 875                                        const void *real_area1, const void *real_area2,
 876                                        mc_snapshot_t snapshot1,
 877                                        mc_snapshot_t snapshot2,
 878                                        xbt_dynar_t previous, mc_type_t type,
 879                                        int area_size, int check_ignore,
 880                                        int pointer_level)
 881 {
 882 top:
 883   if (is_stack(real_area1) && is_stack(real_area2))
 884     return 0;
 885
 886   ssize_t ignore1, ignore2;
 887
 888   if ((check_ignore > 0)
 889       && ((ignore1 = heap_comparison_ignore_size(state->to_ignore1, real_area1))
 890           > 0)
 891       && ((ignore2 = heap_comparison_ignore_size(state->to_ignore2, real_area2))
 892           == ignore1)) {
 893     return 0;
 894   }
 895
 896   mc_type_t subtype, subsubtype;
 897   int res, elm_size;
 898   const void *addr_pointed1, *addr_pointed2;
 899
 900   mc_mem_region_t heap_region1 = MC_get_heap_region(snapshot1);
 901   mc_mem_region_t heap_region2 = MC_get_heap_region(snapshot2);
 902
 903   switch (type->type) {
 904   case DW_TAG_unspecified_type:
 905     return 1;
 906
 907   case DW_TAG_base_type:
 908     if (!type->name.empty() && type->name == "char") {        /* String, hence random (arbitrary ?) size */
 909       if (real_area1 == real_area2)
 910         return -1;
 911       else
 912         return (MC_snapshot_region_memcmp(real_area1, heap_region1, real_area2, heap_region2, area_size) != 0);
 913     } else {
 914       if (area_size != -1 && type->byte_size != area_size)
 915         return -1;
 916       else {
 917         return (MC_snapshot_region_memcmp(real_area1, heap_region1, real_area2, heap_region2, type->byte_size) != 0);
 918       }
 919     }
 920     break;
 921   case DW_TAG_enumeration_type:
 922     if (area_size != -1 && type->byte_size != area_size)
 923       return -1;
 924     else
 925       return (MC_snapshot_region_memcmp(real_area1, heap_region1, real_area2, heap_region2, type->byte_size) != 0);
 926     break;
 927   case DW_TAG_typedef:
 928   case DW_TAG_const_type:
 929   case DW_TAG_volatile_type:
 930     // Poor man's TCO:
 931     type = type->subtype;
 932     goto top;
 933     break;
 934   case DW_TAG_array_type:
 935     subtype = type->subtype;
 936     switch (subtype->type) {
 937     case DW_TAG_unspecified_type:
 938       return 1;
 939
 940     case DW_TAG_base_type:
 941     case DW_TAG_enumeration_type:
 942     case DW_TAG_pointer_type:
 943     case DW_TAG_reference_type:
 944     case DW_TAG_rvalue_reference_type:
 945     case DW_TAG_structure_type:
 946     case DW_TAG_class_type:
 947     case DW_TAG_union_type:
 948       if (subtype->full_type)
 949         subtype = subtype->full_type;
 950       elm_size = subtype->byte_size;
 951       break;
 952       // TODO, just remove the type indirection?
 953     case DW_TAG_const_type:
 954     case DW_TAG_typedef:
 955     case DW_TAG_volatile_type:
 956       subsubtype = subtype->subtype;
 957       if (subsubtype->full_type)
 958         subsubtype = subsubtype->full_type;
 959       elm_size = subsubtype->byte_size;
 960       break;
 961     default:
 962       return 0;
 963       break;
 964     }
 965     for (int i = 0; i < type->element_count; i++) {
 966       // TODO, add support for variable stride (DW_AT_byte_stride)
 967       res =
 968           compare_heap_area_with_type(state, process_index,
 969                                       (char *) real_area1 + (i * elm_size),
 970                                       (char *) real_area2 + (i * elm_size),
 971                                       snapshot1, snapshot2, previous,
 972                                       type->subtype, subtype->byte_size,
 973                                       check_ignore, pointer_level);
 974       if (res == 1)
 975         return res;
 976     }
 977     break;
 978   case DW_TAG_reference_type:
 979   case DW_TAG_rvalue_reference_type:
 980   case DW_TAG_pointer_type:
 981     if (type->subtype && type->subtype->type == DW_TAG_subroutine_type) {
 982       addr_pointed1 = snapshot1->read(remote((void**)real_area1), process_index);
 983       addr_pointed2 = snapshot2->read(remote((void**)real_area2), process_index);
 984       return (addr_pointed1 != addr_pointed2);;
 985     } else {
 986       pointer_level++;
 987       if (pointer_level > 1) {  /* Array of pointers */
 988         for (size_t i = 0; i < (area_size / sizeof(void *)); i++) {
 989           addr_pointed1 = snapshot1->read(
 990             remote((void**)((char*) real_area1 + i * sizeof(void *))),
 991             process_index);
 992           addr_pointed2 = snapshot2->read(
 993             remote((void**)((char*) real_area2 + i * sizeof(void *))),
 994             process_index);
 995           if (addr_pointed1 > state->std_heap_copy.heapbase
 996               && addr_pointed1 < mc_snapshot_get_heap_end(snapshot1)
 997               && addr_pointed2 > state->std_heap_copy.heapbase
 998               && addr_pointed2 < mc_snapshot_get_heap_end(snapshot2))
 999             res =
1000                 compare_heap_area(process_index, addr_pointed1, addr_pointed2, snapshot1,
1001                                   snapshot2, previous, type->subtype,
1002                                   pointer_level);
1003           else
1004             res = (addr_pointed1 != addr_pointed2);
1005           if (res == 1)
1006             return res;
1007         }
1008       } else {
1009         addr_pointed1 = snapshot1->read(remote((void**)real_area1), process_index);
1010         addr_pointed2 = snapshot2->read(remote((void**)real_area2), process_index);
1011         if (addr_pointed1 > state->std_heap_copy.heapbase
1012             && addr_pointed1 < mc_snapshot_get_heap_end(snapshot1)
1013             && addr_pointed2 > state->std_heap_copy.heapbase
1014             && addr_pointed2 < mc_snapshot_get_heap_end(snapshot2))
1015           return compare_heap_area(process_index, addr_pointed1, addr_pointed2, snapshot1,
1016                                    snapshot2, previous, type->subtype,
1017                                    pointer_level);
1018         else
1019           return (addr_pointed1 != addr_pointed2);
1020       }
1021     }
1022     break;
1023   case DW_TAG_structure_type:
1024   case DW_TAG_class_type:
1025     if (type->full_type)
1026       type = type->full_type;
1027     if (area_size != -1 && type->byte_size != area_size) {
1028       if (area_size > type->byte_size && area_size % type->byte_size == 0) {
1029         for (size_t i = 0; i < (size_t)(area_size / type->byte_size); i++) {
1030           res =
1031               compare_heap_area_with_type(state, process_index,
1032                                           (char *) real_area1 + i * type->byte_size,
1033                                           (char *) real_area2 + i * type->byte_size,
1034                                           snapshot1, snapshot2, previous, type, -1,
1035                                           check_ignore, 0);
1036           if (res == 1)
1037             return res;
1038         }
1039       } else {
1040         return -1;
1041       }
1042     } else {
1043       for(simgrid::mc::Type& member : type->members) {
1044         // TODO, optimize this? (for the offset case)
1045         void *real_member1 =
1046             mc_member_resolve(real_area1, type, &member, (mc_address_space_t) snapshot1, process_index);
1047         void *real_member2 =
1048             mc_member_resolve(real_area2, type, &member, (mc_address_space_t) snapshot2, process_index);
1049         res =
1050             compare_heap_area_with_type(state, process_index, real_member1, real_member2,
1051                                         snapshot1, snapshot2,
1052                                         previous, member.subtype, -1,
1053                                         check_ignore, 0);
1054         if (res == 1) {
1055           return res;
1056         }
1057       }
1058     }
1059     break;
1060   case DW_TAG_union_type:
1061     return compare_heap_area_without_type(state, process_index, real_area1, real_area2,
1062                                           snapshot1, snapshot2, previous,
1063                                           type->byte_size, check_ignore);
1064     break;
1065   default:
1066     break;
1067   }
1068
1069   return 0;
1070
1071 }
1072
1073 /** Infer the type of a part of the block from the type of the block
1074  *
1075  * TODO, handle DW_TAG_array_type as well as arrays of the object ((*p)[5], p[5])
1076  *
1077  * TODO, handle subfields ((*p).bar.foo, (*p)[5].bar…)
1078  *
1079  * @param  type_id            DWARF type ID of the root address
1080  * @param  area_size
1081  * @return                    DWARF type ID for given offset
1082  */
1083 static mc_type_t get_offset_type(void *real_base_address, mc_type_t type,
1084                                  int offset, int area_size,
1085                                  mc_snapshot_t snapshot, int process_index)
1086 {
1087
1088   // Beginning of the block, the infered variable type if the type of the block:
1089   if (offset == 0)
1090     return type;
1091
1092   switch (type->type) {
1093   case DW_TAG_structure_type:
1094   case DW_TAG_class_type:
1095     if (type->full_type)
1096       type = type->full_type;
1097
1098     if (area_size != -1 && type->byte_size != area_size) {
1099       if (area_size > type->byte_size && area_size % type->byte_size == 0)
1100         return type;
1101       else
1102         return NULL;
1103     } else {
1104       for(simgrid::mc::Type& member : type->members) {
1105
1106         if (member.has_offset_location()) {
1107           // We have the offset, use it directly (shortcut):
1108           if (member.offset() == offset)
1109             return member.subtype;
1110         } else {
1111           void *real_member =
1112             mc_member_resolve(real_base_address, type, &member,
1113               snapshot, process_index);
1114           if ((char*) real_member - (char *) real_base_address == offset)
1115             return member.subtype;
1116         }
1117
1118       }
1119       return NULL;
1120     }
1121     break;
1122   default:
1123     /* FIXME : other cases ? */
1124     return NULL;
1125     break;
1126   }
1127 }
1128
1129 /**
1130  *
1131  * @param area1          Process address for state 1
1132  * @param area2          Process address for state 2
1133  * @param snapshot1      Snapshot of state 1
1134  * @param snapshot2      Snapshot of state 2
1135  * @param previous       Pairs of blocks already compared on the current path (or NULL)
1136  * @param type_id        Type of variable
1137  * @param pointer_level
1138  * @return 0 (same), 1 (different), -1
1139  */
1140 int compare_heap_area(int process_index, const void *area1, const void *area2, mc_snapshot_t snapshot1,
1141                       mc_snapshot_t snapshot2, xbt_dynar_t previous,
1142                       mc_type_t type, int pointer_level)
1143 {
1144   mc_process_t process = &mc_model_checker->process();
1145
1146   struct s_mc_diff *state = mc_diff_info;
1147
1148   int res_compare;
1149   ssize_t block1, frag1, block2, frag2;
1150   ssize_t size;
1151   int check_ignore = 0;
1152
1153   void *real_addr_block1, *real_addr_block2, *real_addr_frag1, *real_addr_frag2;
1154   int type_size = -1;
1155   int offset1 = 0, offset2 = 0;
1156   int new_size1 = -1, new_size2 = -1;
1157   mc_type_t new_type1 = NULL, new_type2 = NULL;
1158
1159   int match_pairs = 0;
1160
1161   // This is the address of std_heap->heapinfo in the application process:
1162   void* heapinfo_address = &((xbt_mheap_t) process->heap_address)->heapinfo;
1163
1164   const malloc_info* heapinfos1 = snapshot1->read(
1165     remote((const malloc_info**)heapinfo_address), process_index);
1166   const malloc_info* heapinfos2 = snapshot2->read(
1167     remote((const malloc_info**)heapinfo_address), process_index);
1168
1169   malloc_info heapinfo_temp1, heapinfo_temp2;
1170
1171   if (previous == NULL) {
1172     previous =
1173         xbt_dynar_new(sizeof(heap_area_pair_t), heap_area_pair_free_voidp);
1174     match_pairs = 1;
1175   }
1176   // Get block number:
1177   block1 =
1178       ((char *) area1 -
1179        (char *) state->std_heap_copy.heapbase) / BLOCKSIZE + 1;
1180   block2 =
1181       ((char *) area2 -
1182        (char *) state->std_heap_copy.heapbase) / BLOCKSIZE + 1;
1183
1184   // If either block is a stack block:
1185   if (is_block_stack((int) block1) && is_block_stack((int) block2)) {
1186     add_heap_area_pair(previous, block1, -1, block2, -1);
1187     if (match_pairs) {
1188       match_equals(state, previous);
1189       xbt_dynar_free(&previous);
1190     }
1191     return 0;
1192   }
1193   // If either block is not in the expected area of memory:
1194   if (((char *) area1 < (char *) state->std_heap_copy.heapbase)
1195       || (block1 > (ssize_t) state->heapsize1) || (block1 < 1)
1196       || ((char *) area2 < (char *) state->std_heap_copy.heapbase)
1197       || (block2 > (ssize_t) state->heapsize2) || (block2 < 1)) {
1198     if (match_pairs) {
1199       xbt_dynar_free(&previous);
1200     }
1201     return 1;
1202   }
1203
1204   // Process address of the block:
1205   real_addr_block1 = (ADDR2UINT(block1) - 1) * BLOCKSIZE +
1206                  (char *) state->std_heap_copy.heapbase;
1207   real_addr_block2 = (ADDR2UINT(block2) - 1) * BLOCKSIZE +
1208                  (char *) state->std_heap_copy.heapbase;
1209
1210   if (type) {
1211
1212     if (type->full_type)
1213       type = type->full_type;
1214
1215     // This assume that for "boring" types (volatile ...) byte_size is absent:
1216     while (type->byte_size == 0 && type->subtype != NULL)
1217       type = type->subtype;
1218
1219     // Find type_size:
1220     if ((type->type == DW_TAG_pointer_type)
1221         || ((type->type == DW_TAG_base_type) && !type->name.empty()
1222             && type->name == "char"))
1223       type_size = -1;
1224     else
1225       type_size = type->byte_size;
1226
1227   }
1228
1229   mc_mem_region_t heap_region1 = MC_get_heap_region(snapshot1);
1230   mc_mem_region_t heap_region2 = MC_get_heap_region(snapshot2);
1231
1232   const malloc_info* heapinfo1 = (const malloc_info*) MC_region_read(
1233     heap_region1, &heapinfo_temp1, &heapinfos1[block1], sizeof(malloc_info));
1234   const malloc_info* heapinfo2 = (const malloc_info*) MC_region_read(
1235     heap_region2, &heapinfo_temp2, &heapinfos2[block2], sizeof(malloc_info));
1236
1237   if ((heapinfo1->type == MMALLOC_TYPE_FREE || heapinfo1->type==MMALLOC_TYPE_HEAPINFO)
1238     && (heapinfo2->type == MMALLOC_TYPE_FREE || heapinfo2->type ==MMALLOC_TYPE_HEAPINFO)) {
1239
1240     /* Free block */
1241     if (match_pairs) {
1242       match_equals(state, previous);
1243       xbt_dynar_free(&previous);
1244     }
1245     return 0;
1246
1247   } else if (heapinfo1->type == MMALLOC_TYPE_UNFRAGMENTED
1248     && heapinfo2->type == MMALLOC_TYPE_UNFRAGMENTED) {
1249     /* Complete block */
1250
1251     // TODO, lookup variable type from block type as done for fragmented blocks
1252
1253     offset1 = (char *) area1 - (char *) real_addr_block1;
1254     offset2 = (char *) area2 - (char *) real_addr_block2;
1255
1256     if (state->equals_to1_(block1, 0).valid
1257         && state->equals_to2_(block2, 0).valid) {
1258       if (equal_blocks(state, block1, block2)) {
1259         if (match_pairs) {
1260           match_equals(state, previous);
1261           xbt_dynar_free(&previous);
1262         }
1263         return 0;
1264       }
1265     }
1266
1267     if (type_size != -1) {
1268       if (type_size != (ssize_t) heapinfo1->busy_block.busy_size
1269           && type_size != (ssize_t)   heapinfo2->busy_block.busy_size
1270           && (type->name.empty() || type->name == "struct s_smx_context")) {
1271         if (match_pairs) {
1272           match_equals(state, previous);
1273           xbt_dynar_free(&previous);
1274         }
1275         return -1;
1276       }
1277     }
1278
1279     if (heapinfo1->busy_block.size !=
1280         heapinfo2->busy_block.size) {
1281       if (match_pairs) {
1282         xbt_dynar_free(&previous);
1283       }
1284       return 1;
1285     }
1286
1287     if (heapinfo1->busy_block.busy_size !=
1288         heapinfo2->busy_block.busy_size) {
1289       if (match_pairs) {
1290         xbt_dynar_free(&previous);
1291       }
1292       return 1;
1293     }
1294
1295     if (!add_heap_area_pair(previous, block1, -1, block2, -1)) {
1296       if (match_pairs) {
1297         match_equals(state, previous);
1298         xbt_dynar_free(&previous);
1299       }
1300       return 0;
1301     }
1302
1303     size = heapinfo1->busy_block.busy_size;
1304
1305     // Remember (basic) type inference.
1306     // The current data structure only allows us to do this for the whole block.
1307     if (type != NULL && area1 == real_addr_block1) {
1308       state->types1_(block1, 0) = type;
1309     }
1310     if (type != NULL && area2 == real_addr_block2) {
1311       state->types2_(block2, 0) = type;
1312     }
1313
1314     if (size <= 0) {
1315       if (match_pairs) {
1316         match_equals(state, previous);
1317         xbt_dynar_free(&previous);
1318       }
1319       return 0;
1320     }
1321
1322     frag1 = -1;
1323     frag2 = -1;
1324
1325     if ((heapinfo1->busy_block.ignore > 0)
1326         && (heapinfo2->busy_block.ignore ==
1327             heapinfo1->busy_block.ignore))
1328       check_ignore = heapinfo1->busy_block.ignore;
1329
1330   } else if ((heapinfo1->type > 0) && (heapinfo2->type > 0)) {      /* Fragmented block */
1331
1332     // Fragment number:
1333     frag1 =
1334         ((uintptr_t) (ADDR2UINT(area1) % (BLOCKSIZE))) >> heapinfo1->type;
1335     frag2 =
1336         ((uintptr_t) (ADDR2UINT(area2) % (BLOCKSIZE))) >> heapinfo2->type;
1337
1338     // Process address of the fragment:
1339     real_addr_frag1 =
1340         (void *) ((char *) real_addr_block1 +
1341                   (frag1 << heapinfo1->type));
1342     real_addr_frag2 =
1343         (void *) ((char *) real_addr_block2 +
1344                   (frag2 << heapinfo2->type));
1345
1346     // Check the size of the fragments against the size of the type:
1347     if (type_size != -1) {
1348       if (heapinfo1->busy_frag.frag_size[frag1] == -1
1349           || heapinfo2->busy_frag.frag_size[frag2] == -1) {
1350         if (match_pairs) {
1351           match_equals(state, previous);
1352           xbt_dynar_free(&previous);
1353         }
1354         return -1;
1355       }
1356       // ?
1357       if (type_size != heapinfo1->busy_frag.frag_size[frag1]
1358           || type_size != heapinfo2->busy_frag.frag_size[frag2]) {
1359         if (match_pairs) {
1360           match_equals(state, previous);
1361           xbt_dynar_free(&previous);
1362         }
1363         return -1;
1364       }
1365     }
1366
1367     // Check if the blocks are already matched together:
1368     if (state->equals_to1_(block1, frag1).valid
1369         && state->equals_to2_(block2, frag2).valid) {
1370       if (offset1==offset2 && equal_fragments(state, block1, frag1, block2, frag2)) {
1371         if (match_pairs) {
1372           match_equals(state, previous);
1373           xbt_dynar_free(&previous);
1374         }
1375         return 0;
1376       }
1377     }
1378     // Compare the size of both fragments:
1379     if (heapinfo1->busy_frag.frag_size[frag1] !=
1380         heapinfo2->busy_frag.frag_size[frag2]) {
1381       if (type_size == -1) {
1382         if (match_pairs) {
1383           match_equals(state, previous);
1384           xbt_dynar_free(&previous);
1385         }
1386         return -1;
1387       } else {
1388         if (match_pairs) {
1389           xbt_dynar_free(&previous);
1390         }
1391         return 1;
1392       }
1393     }
1394
1395     // Size of the fragment:
1396     size = heapinfo1->busy_frag.frag_size[frag1];
1397
1398     // Remember (basic) type inference.
1399     // The current data structure only allows us to do this for the whole fragment.
1400     if (type != NULL && area1 == real_addr_frag1) {
1401       state->types1_(block1, frag1) = type;
1402     }
1403     if (type != NULL && area2 == real_addr_frag2) {
1404       state->types2_(block2, frag2) = type;
1405     }
1406     // The type of the variable is already known:
1407     if (type) {
1408       new_type1 = type;
1409       new_type2 = type;
1410     }
1411     // Type inference from the block type.
1412     else if (state->types1_(block1, frag1) != NULL
1413              || state->types2_(block2, frag2) != NULL) {
1414
1415       offset1 = (char *) area1 - (char *) real_addr_frag1;
1416       offset2 = (char *) area2 - (char *) real_addr_frag2;
1417
1418       if (state->types1_(block1, frag1) != NULL
1419           && state->types2_(block2, frag2) != NULL) {
1420         new_type1 =
1421             get_offset_type(real_addr_frag1, state->types1_(block1, frag1),
1422                             offset1, size, snapshot1, process_index);
1423         new_type2 =
1424             get_offset_type(real_addr_frag2, state->types2_(block2, frag2),
1425                             offset1, size, snapshot2, process_index);
1426       } else if (state->types1_(block1, frag1) != NULL) {
1427         new_type1 =
1428             get_offset_type(real_addr_frag1, state->types1_(block1, frag1),
1429                             offset1, size, snapshot1, process_index);
1430         new_type2 =
1431             get_offset_type(real_addr_frag2, state->types1_(block1, frag1),
1432                             offset2, size, snapshot2, process_index);
1433       } else if (state->types2_(block2, frag2) != NULL) {
1434         new_type1 =
1435             get_offset_type(real_addr_frag1, state->types2_(block2, frag2),
1436                             offset1, size, snapshot1, process_index);
1437         new_type2 =
1438             get_offset_type(real_addr_frag2, state->types2_(block2, frag2),
1439                             offset2, size, snapshot2, process_index);
1440       } else {
1441         if (match_pairs) {
1442           match_equals(state, previous);
1443           xbt_dynar_free(&previous);
1444         }
1445         return -1;
1446       }
1447
1448       if (new_type1 != NULL && new_type2 != NULL && new_type1 != new_type2) {
1449
1450         type = new_type1;
1451         while (type->byte_size == 0 && type->subtype != NULL)
1452           type = type->subtype;
1453         new_size1 = type->byte_size;
1454
1455         type = new_type2;
1456         while (type->byte_size == 0 && type->subtype != NULL)
1457           type = type->subtype;
1458         new_size2 = type->byte_size;
1459
1460       } else {
1461         if (match_pairs) {
1462           match_equals(state, previous);
1463           xbt_dynar_free(&previous);
1464         }
1465         return -1;
1466       }
1467     }
1468
1469     if (new_size1 > 0 && new_size1 == new_size2) {
1470       type = new_type1;
1471       size = new_size1;
1472     }
1473
1474     if (offset1 == 0 && offset2 == 0) {
1475       if (!add_heap_area_pair(previous, block1, frag1, block2, frag2)) {
1476         if (match_pairs) {
1477           match_equals(state, previous);
1478           xbt_dynar_free(&previous);
1479         }
1480         return 0;
1481       }
1482     }
1483
1484     if (size <= 0) {
1485       if (match_pairs) {
1486         match_equals(state, previous);
1487         xbt_dynar_free(&previous);
1488       }
1489       return 0;
1490     }
1491
1492     if ((heapinfo1->busy_frag.ignore[frag1] > 0)
1493         && (heapinfo2->busy_frag.ignore[frag2] ==
1494             heapinfo1->busy_frag.ignore[frag1]))
1495       check_ignore = heapinfo1->busy_frag.ignore[frag1];
1496
1497   } else {
1498
1499     if (match_pairs) {
1500       xbt_dynar_free(&previous);
1501     }
1502     return 1;
1503
1504   }
1505
1506
1507   /* Start comparison */
1508   if (type) {
1509     res_compare =
1510         compare_heap_area_with_type(state, process_index, area1, area2, snapshot1, snapshot2,
1511                                     previous, type, size, check_ignore,
1512                                     pointer_level);
1513   } else {
1514     res_compare =
1515         compare_heap_area_without_type(state, process_index, area1, area2, snapshot1, snapshot2,
1516                                        previous, size, check_ignore);
1517   }
1518   if (res_compare == 1) {
1519     if (match_pairs)
1520       xbt_dynar_free(&previous);
1521     return res_compare;
1522   }
1523
1524   if (match_pairs) {
1525     match_equals(state, previous);
1526     xbt_dynar_free(&previous);
1527   }
1528
1529   return 0;
1530 }
1531
1532 /*********************************************** Miscellaneous ***************************************************/
1533 /****************************************************************************************************************/
1534
1535 // Not used and broken code:
1536 # if 0
1537
1538 // Not used:
1539 static int get_pointed_area_size(void *area, int heap)
1540 {
1541
1542   struct s_mc_diff *state = mc_diff_info;
1543
1544   int block, frag;
1545   malloc_info *heapinfo;
1546
1547   if (heap == 1)
1548     heapinfo = state->heapinfo1;
1549   else
1550     heapinfo = state->heapinfo2;
1551
1552   block =
1553       ((char *) area -
1554        (char *) state->std_heap_copy.heapbase) / BLOCKSIZE + 1;
1555
1556   if (((char *) area < (char *) state->std_heap_copy.heapbase)
1557       || (block > state->heapsize1) || (block < 1))
1558     return -1;
1559
1560   if (heapinfo[block].type == MMALLOC_TYPE_FREE || heapinfo[block].type == MMALLOC_TYPE_HEAPINFO) {     /* Free block */
1561     return -1;
1562   } else if (heapinfo[block].type == MMALLOC_TYPE_UNFRAGMENTED) {       /* Complete block */
1563     return (int) heapinfo[block].busy_block.busy_size;
1564   } else {
1565     frag =
1566         ((uintptr_t) (ADDR2UINT(area) % (BLOCKSIZE))) >> heapinfo[block].type;
1567     return (int) heapinfo[block].busy_frag.frag_size[frag];
1568   }
1569 }
1570
1571 #ifndef max
1572 #define max( a, b ) ( ((a) > (b)) ? (a) : (b) )
1573 #endif
1574
1575 // Not used:
1576 int mmalloc_linear_compare_heap(xbt_mheap_t heap1, xbt_mheap_t heap2)
1577 {
1578
1579   struct s_mc_diff *state = mc_diff_info;
1580
1581   if (heap1 == NULL && heap1 == NULL) {
1582     XBT_DEBUG("Malloc descriptors null");
1583     return 0;
1584   }
1585
1586   if (heap1->heaplimit != heap2->heaplimit) {
1587     XBT_DEBUG("Different limit of valid info table indices");
1588     return 1;
1589   }
1590
1591   /* Heap information */
1592   state->heaplimit = ((struct mdesc *) heap1)->heaplimit;
1593
1594   state->std_heap_copy = *mc_model_checker->process().get_heap();
1595
1596   state->heapbase1 = (char *) heap1 + BLOCKSIZE;
1597   state->heapbase2 = (char *) heap2 + BLOCKSIZE;
1598
1599   state->heapinfo1 =
1600       (malloc_info *) ((char *) heap1 +
1601                        ((uintptr_t)
1602                         ((char *) heap1->heapinfo - (char *) state->s_heap)));
1603   state->heapinfo2 =
1604       (malloc_info *) ((char *) heap2 +
1605                        ((uintptr_t)
1606                         ((char *) heap2->heapinfo - (char *) state->s_heap)));
1607
1608   state->heapsize1 = heap1->heapsize;
1609   state->heapsize2 = heap2->heapsize;
1610
1611   /* Start comparison */
1612   size_t i, j, k;
1613   void *addr_block1, *addr_block2, *addr_frag1, *addr_frag2;
1614
1615   int distance = 0;
1616
1617   /* Check busy blocks */
1618
1619   i = 1;
1620
1621   while (i <= state->heaplimit) {
1622
1623     addr_block1 =
1624         ((void *) (((ADDR2UINT(i)) - 1) * BLOCKSIZE +
1625                    (char *) state->heapbase1));
1626     addr_block2 =
1627         ((void *) (((ADDR2UINT(i)) - 1) * BLOCKSIZE +
1628                    (char *) state->heapbase2));
1629
1630     if (state->heapinfo1[i].type != state->heapinfo2[i].type) {
1631
1632       distance += BLOCKSIZE;
1633       XBT_DEBUG("Different type of blocks (%zu) : %d - %d -> distance = %d", i,
1634                 state->heapinfo1[i].type, state->heapinfo2[i].type, distance);
1635       i++;
1636
1637     } else {
1638
1639       if (state->heapinfo1[i].type == MMALLOC_TYPE_FREE
1640         || state->heapinfo1[i].type == MMALLOC_TYPE_HAPINFO) {     /* Free block */
1641         i++;
1642         continue;
1643       }
1644
1645       if (state->heapinfo1[i].type == MMALLOC_TYPE_UNFRAGMENTED) {      /* Large block */
1646
1647         if (state->heapinfo1[i].busy_block.size !=
1648             state->heapinfo2[i].busy_block.size) {
1649           distance +=
1650               BLOCKSIZE * max(state->heapinfo1[i].busy_block.size,
1651                               state->heapinfo2[i].busy_block.size);
1652           i += max(state->heapinfo1[i].busy_block.size,
1653                    state->heapinfo2[i].busy_block.size);
1654           XBT_DEBUG
1655               ("Different larger of cluster at block %zu : %zu - %zu -> distance = %d",
1656                i, state->heapinfo1[i].busy_block.size,
1657                state->heapinfo2[i].busy_block.size, distance);
1658           continue;
1659         }
1660
1661         /*if(heapinfo1[i].busy_block.busy_size != heapinfo2[i].busy_block.busy_size){
1662            distance += max(heapinfo1[i].busy_block.busy_size, heapinfo2[i].busy_block.busy_size);
1663            i += max(heapinfo1[i].busy_block.size, heapinfo2[i].busy_block.size);
1664            XBT_DEBUG("Different size used oin large cluster at block %zu : %zu - %zu -> distance = %d", i, heapinfo1[i].busy_block.busy_size, heapinfo2[i].busy_block.busy_size, distance);
1665            continue;
1666            } */
1667
1668         k = 0;
1669
1670         //while(k < (heapinfo1[i].busy_block.busy_size)){
1671         while (k < state->heapinfo1[i].busy_block.size * BLOCKSIZE) {
1672           if (memcmp((char *) addr_block1 + k, (char *) addr_block2 + k, 1) !=
1673               0) {
1674             distance++;
1675           }
1676           k++;
1677         }
1678
1679         i++;
1680
1681       } else {                  /* Fragmented block */
1682
1683         for (j = 0; j < (size_t) (BLOCKSIZE >> state->heapinfo1[i].type); j++) {
1684
1685           addr_frag1 =
1686               (void *) ((char *) addr_block1 + (j << state->heapinfo1[i].type));
1687           addr_frag2 =
1688               (void *) ((char *) addr_block2 + (j << state->heapinfo2[i].type));
1689
1690           if (state->heapinfo1[i].busy_frag.frag_size[j] == 0
1691               && state->heapinfo2[i].busy_frag.frag_size[j] == 0) {
1692             continue;
1693           }
1694
1695
1696           /*if(heapinfo1[i].busy_frag.frag_size[j] != heapinfo2[i].busy_frag.frag_size[j]){
1697              distance += max(heapinfo1[i].busy_frag.frag_size[j], heapinfo2[i].busy_frag.frag_size[j]);
1698              XBT_DEBUG("Different size used in fragment %zu in block %zu : %d - %d -> distance = %d", j, i, heapinfo1[i].busy_frag.frag_size[j], heapinfo2[i].busy_frag.frag_size[j], distance);
1699              continue;
1700              } */
1701
1702           k = 0;
1703
1704           //while(k < max(heapinfo1[i].busy_frag.frag_size[j], heapinfo2[i].busy_frag.frag_size[j])){
1705           while (k < (BLOCKSIZE / (BLOCKSIZE >> state->heapinfo1[i].type))) {
1706             if (memcmp((char *) addr_frag1 + k, (char *) addr_frag2 + k, 1) !=
1707                 0) {
1708               distance++;
1709             }
1710             k++;
1711           }
1712
1713         }
1714
1715         i++;
1716
1717       }
1718
1719     }
1720
1721   }
1722
1723   return distance;
1724
1725 }
1726 #endif
1727
1728 }