Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
c681d8b72890e060030bd75681b80197f7420d35
[simgrid.git] / src / xbt / mmalloc / mm_diff.c
1 /* mm_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/datatypes.h"
14 #include "mc/mc_private.h"
15
16 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(mm_diff, xbt,
17                                 "Logging specific to mm_diff in mmalloc");
18
19 xbt_dynar_t mc_heap_comparison_ignore;
20 xbt_dynar_t stacks_areas;
21 void *maestro_stack_start, *maestro_stack_end;
22
23
24 /********************************* Backtrace ***********************************/
25 /******************************************************************************/
26
27 static void mmalloc_backtrace_block_display(void* heapinfo, int block){
28
29   /* xbt_ex_t e; */
30
31   /* if (((malloc_info *)heapinfo)[block].busy_block.bt_size == 0) { */
32   /*   fprintf(stderr, "No backtrace available for that block, sorry.\n"); */
33   /*   return; */
34   /* } */
35
36   /* memcpy(&e.bt,&(((malloc_info *)heapinfo)[block].busy_block.bt),sizeof(void*)*XBT_BACKTRACE_SIZE); */
37   /* e.used = ((malloc_info *)heapinfo)[block].busy_block.bt_size; */
38
39   /* xbt_ex_setup_backtrace(&e); */
40   /* if (e.used == 0) { */
41   /*   fprintf(stderr, "(backtrace not set)\n"); */
42   /* } else if (e.bt_strings == NULL) { */
43   /*   fprintf(stderr, "(backtrace not ready to be computed. %s)\n",xbt_binary_name?"Dunno why":"xbt_binary_name not setup yet"); */
44   /* } else { */
45   /*   int i; */
46
47   /*   fprintf(stderr, "Backtrace of where the block %d was malloced (%d frames):\n", block ,e.used); */
48   /*   for (i = 0; i < e.used; i++)       /\* no need to display "xbt_backtrace_display" *\/{ */
49   /*     fprintf(stderr, "%d ---> %s\n",i, e.bt_strings[i] + 4); */
50   /*   } */
51   /* } */
52 }
53
54 static void mmalloc_backtrace_fragment_display(void* heapinfo, int block, int frag){
55
56   /* xbt_ex_t e; */
57
58   /* memcpy(&e.bt,&(((malloc_info *)heapinfo)[block].busy_frag.bt[frag]),sizeof(void*)*XBT_BACKTRACE_SIZE); */
59   /* e.used = XBT_BACKTRACE_SIZE; */
60
61   /* xbt_ex_setup_backtrace(&e); */
62   /* if (e.used == 0) { */
63   /*   fprintf(stderr, "(backtrace not set)\n"); */
64   /* } else if (e.bt_strings == NULL) { */
65   /*   fprintf(stderr, "(backtrace not ready to be computed. %s)\n",xbt_binary_name?"Dunno why":"xbt_binary_name not setup yet"); */
66   /* } else { */
67   /*   int i; */
68
69   /*   fprintf(stderr, "Backtrace of where the fragment %d in block %d was malloced (%d frames):\n", frag, block ,e.used); */
70   /*   for (i = 0; i < e.used; i++)       /\* no need to display "xbt_backtrace_display" *\/{ */
71   /*     fprintf(stderr, "%d ---> %s\n",i, e.bt_strings[i] + 4); */
72   /*   } */
73   /* } */
74
75 }
76
77 static void mmalloc_backtrace_display(void *addr){
78
79   /* size_t block, frag_nb; */
80   /* int type; */
81   
82   /* xbt_mheap_t heap = __mmalloc_current_heap ?: (xbt_mheap_t) mmalloc_preinit(); */
83
84   /* block = (((char*) (addr) - (char*) heap -> heapbase) / BLOCKSIZE + 1); */
85
86   /* type = heap->heapinfo[block].type; */
87
88   /* switch(type){ */
89   /* case -1 : /\* Free block *\/ */
90   /*   fprintf(stderr, "Asked to display the backtrace of a block that is free. I'm puzzled\n"); */
91   /*   xbt_abort(); */
92   /*   break;  */
93   /* case 0: /\* Large block *\/ */
94   /*   mmalloc_backtrace_block_display(heap->heapinfo, block); */
95   /*   break; */
96   /* default: /\* Fragmented block *\/ */
97   /*   frag_nb = RESIDUAL(addr, BLOCKSIZE) >> type; */
98   /*   if(heap->heapinfo[block].busy_frag.frag_size[frag_nb] == -1){ */
99   /*     fprintf(stderr , "Asked to display the backtrace of a fragment that is free. I'm puzzled\n"); */
100   /*     xbt_abort(); */
101   /*   } */
102   /*   mmalloc_backtrace_fragment_display(heap->heapinfo, block, frag_nb); */
103   /*   break; */
104   /* } */
105 }
106
107
108 static int compare_backtrace(int b1, int f1, int b2, int f2){
109   /*int i = 0;
110   if(f1 != -1){
111     for(i=0; i< XBT_BACKTRACE_SIZE; i++){
112       if(heapinfo1[b1].busy_frag.bt[f1][i] != heapinfo2[b2].busy_frag.bt[f2][i]){
113         //mmalloc_backtrace_fragment_display((void*)heapinfo1, b1, f1);
114         //mmalloc_backtrace_fragment_display((void*)heapinfo2, b2, f2);
115         return 1;
116       }
117     }
118   }else{
119     for(i=0; i< heapinfo1[b1].busy_block.bt_size; i++){
120       if(heapinfo1[b1].busy_block.bt[i] != heapinfo2[b2].busy_block.bt[i]){
121         //mmalloc_backtrace_block_display((void*)heapinfo1, b1);
122         //mmalloc_backtrace_block_display((void*)heapinfo2, b2);
123         return 1;
124       }
125     }
126     }*/
127   return 0;
128 }
129
130
131 /*********************************** Heap comparison ***********************************/
132 /***************************************************************************************/
133
134 typedef char* type_name;
135
136 struct s_mm_diff {
137   void *s_heap, *heapbase1, *heapbase2;
138   malloc_info *heapinfo1, *heapinfo2;
139   size_t heaplimit, heapsize1, heapsize2;
140   xbt_dynar_t to_ignore1, to_ignore2;
141   heap_area_t **equals_to1, **equals_to2;
142   type_name **types1, **types2;
143 };
144
145 __thread struct s_mm_diff* mm_diff_info = NULL;
146
147 /*********************************** Free functions ************************************/
148
149 static void heap_area_pair_free(heap_area_pair_t pair){
150   xbt_free(pair);
151   pair = NULL;
152 }
153
154 static void heap_area_pair_free_voidp(void *d){
155   heap_area_pair_free((heap_area_pair_t) * (void **) d);
156 }
157
158 static void heap_area_free(heap_area_t area){
159   xbt_free(area);
160   area = NULL;
161 }
162
163 /************************************************************************************/
164
165 static heap_area_t new_heap_area(int block, int fragment){
166   heap_area_t area = NULL;
167   area = xbt_new0(s_heap_area_t, 1);
168   area->block = block;
169   area->fragment = fragment;
170   return area;
171 }
172
173  
174 static int is_new_heap_area_pair(xbt_dynar_t list, int block1, int fragment1, int block2, int fragment2){
175   
176   unsigned int cursor = 0;
177   heap_area_pair_t current_pair;
178
179   xbt_dynar_foreach(list, cursor, current_pair){
180     if(current_pair->block1 == block1 && current_pair->block2 == block2 && current_pair->fragment1 == fragment1 && current_pair->fragment2 == fragment2)
181       return 0; 
182   }
183   
184   return 1;
185 }
186
187 static int add_heap_area_pair(xbt_dynar_t list, int block1, int fragment1, int block2, int fragment2){
188
189   if(is_new_heap_area_pair(list, block1, fragment1, block2, fragment2)){
190     heap_area_pair_t pair = NULL;
191     pair = xbt_new0(s_heap_area_pair_t, 1);
192     pair->block1 = block1;
193     pair->fragment1 = fragment1;
194     pair->block2 = block2;
195     pair->fragment2 = fragment2;
196     
197     xbt_dynar_push(list, &pair); 
198
199     return 1;
200   }
201
202   return 0;
203 }
204
205 static ssize_t heap_comparison_ignore_size(xbt_dynar_t ignore_list, void *address){
206
207   unsigned int cursor = 0;
208   int start = 0;
209   int end = xbt_dynar_length(ignore_list) - 1;
210   mc_heap_ignore_region_t region;
211
212   while(start <= end){
213     cursor = (start + end) / 2;
214     region = (mc_heap_ignore_region_t)xbt_dynar_get_as(ignore_list, cursor, mc_heap_ignore_region_t);
215     if(region->address == address)
216       return region->size;
217     if(region->address < address)
218       start = cursor + 1;
219     if(region->address > address)
220       end = cursor - 1;   
221   }
222
223   return -1;
224 }
225
226 static int is_stack(void *address){
227   unsigned int cursor = 0;
228   stack_region_t stack;
229
230   xbt_dynar_foreach(stacks_areas, cursor, stack){
231     if(address == stack->address)
232       return 1;
233   }
234
235   return 0;
236 }
237
238 static int is_block_stack(int block){
239   unsigned int cursor = 0;
240   stack_region_t stack;
241
242   xbt_dynar_foreach(stacks_areas, cursor, stack){
243     if(block == stack->block)
244       return 1;
245   }
246
247   return 0;
248 }
249
250 static void match_equals(struct s_mm_diff *state, xbt_dynar_t list){
251
252   unsigned int cursor = 0;
253   heap_area_pair_t current_pair;
254   heap_area_t previous_area;
255
256   xbt_dynar_foreach(list, cursor, current_pair){
257
258     if(current_pair->fragment1 != -1){
259
260       if(state->equals_to1[current_pair->block1][current_pair->fragment1] != NULL){
261         previous_area = state->equals_to1[current_pair->block1][current_pair->fragment1];
262         heap_area_free(state->equals_to2[previous_area->block][previous_area->fragment]);
263         state->equals_to2[previous_area->block][previous_area->fragment] = NULL;
264         heap_area_free(previous_area);
265       }
266       if(state->equals_to2[current_pair->block2][current_pair->fragment2] != NULL){
267         previous_area = state->equals_to2[current_pair->block2][current_pair->fragment2];
268         heap_area_free(state->equals_to1[previous_area->block][previous_area->fragment]);
269         state->equals_to1[previous_area->block][previous_area->fragment] = NULL;
270         heap_area_free(previous_area);
271       }
272
273       state->equals_to1[current_pair->block1][current_pair->fragment1] = new_heap_area(current_pair->block2, current_pair->fragment2);
274       state->equals_to2[current_pair->block2][current_pair->fragment2] = new_heap_area(current_pair->block1, current_pair->fragment1);
275       
276     }else{
277
278       if(state->equals_to1[current_pair->block1][0] != NULL){
279         previous_area = state->equals_to1[current_pair->block1][0];
280         heap_area_free(state->equals_to2[previous_area->block][0]);
281         state->equals_to2[previous_area->block][0] = NULL;
282         heap_area_free(previous_area);
283       }
284       if(state->equals_to2[current_pair->block2][0] != NULL){
285         previous_area = state->equals_to2[current_pair->block2][0];
286         heap_area_free(state->equals_to1[previous_area->block][0]);
287         state->equals_to1[previous_area->block][0] = NULL;
288         heap_area_free(previous_area);
289       }
290
291       state->equals_to1[current_pair->block1][0] = new_heap_area(current_pair->block2, current_pair->fragment2);
292       state->equals_to2[current_pair->block2][0] = new_heap_area(current_pair->block1, current_pair->fragment1);
293
294     }
295
296   }
297 }
298
299 static int equal_blocks(struct s_mm_diff *state, int b1, int b2){
300   
301   if(state->equals_to1[b1][0]->block == b2 && state->equals_to2[b2][0]->block == b1)
302     return 1;
303
304   return 0;
305 }
306
307 static int equal_fragments(struct s_mm_diff *state, int b1, int f1, int b2, int f2){
308   
309   if(state->equals_to1[b1][f1]->block == b2
310     && state->equals_to1[b1][f1]->fragment == f2
311     && state->equals_to2[b2][f2]->block == b1
312     && state->equals_to2[b2][f2]->fragment == f1)
313     return 1;
314
315   return 0;
316 }
317
318 int init_heap_information(xbt_mheap_t heap1, xbt_mheap_t heap2, xbt_dynar_t i1, xbt_dynar_t i2){
319   if(mm_diff_info==NULL) {
320     mm_diff_info = xbt_new0(struct s_mm_diff, 1);
321   }
322   struct s_mm_diff *state = mm_diff_info;
323
324   if((((struct mdesc *)heap1)->heaplimit != ((struct mdesc *)heap2)->heaplimit)
325     || ((((struct mdesc *)heap1)->heapsize != ((struct mdesc *)heap2)->heapsize) ))
326     return -1;
327
328   int i, j;
329
330   state->heaplimit = ((struct mdesc *)heap1)->heaplimit;
331
332   state->s_heap = (char *)mmalloc_get_current_heap() - STD_HEAP_SIZE - getpagesize();
333
334   state->heapbase1 = (char *)heap1 + BLOCKSIZE;
335   state->heapbase2 = (char *)heap2 + BLOCKSIZE;
336
337   state->heapinfo1 = (malloc_info *)((char *)heap1 + ((uintptr_t)((char *)((struct mdesc *)heap1)->heapinfo - (char *)state->s_heap)));
338   state->heapinfo2 = (malloc_info *)((char *)heap2 + ((uintptr_t)((char *)((struct mdesc *)heap2)->heapinfo - (char *)state->s_heap)));
339
340   state->heapsize1 = heap1->heapsize;
341   state->heapsize2 = heap2->heapsize;
342
343   state->to_ignore1 = i1;
344   state-> to_ignore2 = i2;
345
346   state->equals_to1 = malloc(state->heaplimit * sizeof(heap_area_t *));
347   state->types1 = malloc(state->heaplimit * sizeof(type_name *));
348   for(i=0; i<=state->heaplimit; i++){
349     state->equals_to1[i] = malloc(MAX_FRAGMENT_PER_BLOCK * sizeof(heap_area_t));
350     state->types1[i] = malloc(MAX_FRAGMENT_PER_BLOCK * sizeof(type_name));
351     for(j=0; j<MAX_FRAGMENT_PER_BLOCK; j++){
352       state->equals_to1[i][j] = NULL;
353       state->types1[i][j] = NULL;
354     }      
355   }
356
357   state->equals_to2 = malloc(state->heaplimit * sizeof(heap_area_t *));
358   state->types2 = malloc(state->heaplimit * sizeof(type_name *));
359   for(i=0; i<=state->heaplimit; i++){
360     state->equals_to2[i] = malloc(MAX_FRAGMENT_PER_BLOCK * sizeof(heap_area_t));
361     state->types2[i] = malloc(MAX_FRAGMENT_PER_BLOCK * sizeof(type_name));
362     for(j=0; j<MAX_FRAGMENT_PER_BLOCK; j++){
363       state->equals_to2[i][j] = NULL;
364       state->types2[i][j] = NULL;
365     }
366   }
367
368   if(MC_is_active()){
369     MC_ignore_global_variable("mm_diff_info");
370   }
371
372   return 0;
373
374 }
375
376 void reset_heap_information(){
377
378   struct s_mm_diff *state = mm_diff_info;
379
380   size_t i = 0, j;
381
382   for(i=0; i<=state->heaplimit; i++){
383     for(j=0; j<MAX_FRAGMENT_PER_BLOCK;j++){
384       heap_area_free(state->equals_to1[i][j]);
385       state->equals_to1[i][j] = NULL;
386       heap_area_free(state->equals_to2[i][j]);
387       state-> equals_to2[i][j] = NULL;
388       xbt_free(state->types1[i][j]);
389       state->types1[i][j] = NULL;
390       xbt_free(state->types2[i][j]);
391       state->types2[i][j] = NULL;
392     }
393     free(state->equals_to1[i]);
394     free(state->equals_to2[i]);
395     free(state->types1[i]);
396     free(state->types2[i]);
397   }
398
399   free(state->equals_to1);
400   free(state->equals_to2);
401   free(state->types1);
402   free(state->types2);
403
404   state->s_heap = NULL, state->heapbase1 = NULL, state->heapbase2 = NULL;
405   state->heapinfo1 = NULL, state->heapinfo2 = NULL;
406   state->heaplimit = 0, state->heapsize1 = 0, state->heapsize2 = 0;
407   state->to_ignore1 = NULL, state->to_ignore2 = NULL;
408   state->equals_to1 = NULL, state->equals_to2 = NULL;
409   state->types1 = NULL, state->types2 = NULL;
410
411 }
412
413 int mmalloc_compare_heap(xbt_mheap_t heap1, xbt_mheap_t heap2, mc_object_info_t info, mc_object_info_t other_info){
414
415   struct s_mm_diff *state = mm_diff_info;
416
417   if(heap1 == NULL && heap2 == NULL){
418     XBT_DEBUG("Malloc descriptors null");
419     return 0;
420   }
421
422   /* Start comparison */
423   size_t i1, i2, j1, j2, k;
424   void *addr_block1, *addr_block2, *addr_frag1, *addr_frag2;
425   int nb_diff1 = 0, nb_diff2 = 0;
426
427   xbt_dynar_t previous = xbt_dynar_new(sizeof(heap_area_pair_t), heap_area_pair_free_voidp);
428
429   int equal, res_compare = 0;
430
431   /* Check busy blocks*/
432
433   i1 = 1;
434
435   while(i1 <= state->heaplimit){
436
437     if(state->heapinfo1[i1].type == -1){ /* Free block */
438       i1++;
439       continue;
440     }
441
442     addr_block1 = ((void*) (((ADDR2UINT(i1)) - 1) * BLOCKSIZE + (char*)((xbt_mheap_t)state->s_heap)->heapbase));
443
444     if(state->heapinfo1[i1].type == 0){  /* Large block */
445       
446       if(is_stack(addr_block1)){
447         for(k=0; k < state->heapinfo1[i1].busy_block.size; k++)
448           state->equals_to1[i1+k][0] = new_heap_area(i1, -1);
449         for(k=0; k < state->heapinfo2[i1].busy_block.size; k++)
450           state->equals_to2[i1+k][0] = new_heap_area(i1, -1);
451         i1 += state->heapinfo1[i1].busy_block.size;
452         continue;
453       }
454
455       if(state->equals_to1[i1][0] != NULL){
456         i1++;
457         continue;
458       }
459     
460       i2 = 1;
461       equal = 0;
462       res_compare = 0;
463   
464       /* Try first to associate to same block in the other heap */
465       if(state->heapinfo2[i1].type == state->heapinfo1[i1].type){
466
467         if(state->equals_to2[i1][0] == NULL){
468
469           addr_block2 = ((void*) (((ADDR2UINT(i1)) - 1) * BLOCKSIZE + (char*)((xbt_mheap_t)state->s_heap)->heapbase));
470         
471           res_compare = compare_heap_area(addr_block1, addr_block2, NULL, info, other_info, NULL, 0);
472         
473           if(res_compare != 1){
474             for(k=1; k < state->heapinfo2[i1].busy_block.size; k++)
475               state->equals_to2[i1+k][0] = new_heap_area(i1, -1);
476             for(k=1; k < state->heapinfo1[i1].busy_block.size; k++)
477               state->equals_to1[i1+k][0] = new_heap_area(i1, -1);
478             equal = 1;
479             i1 += state->heapinfo1[i1].busy_block.size;
480           }
481         
482           xbt_dynar_reset(previous);
483         
484         }
485         
486       }
487
488       while(i2 <= state->heaplimit && !equal){
489
490         addr_block2 = ((void*) (((ADDR2UINT(i2)) - 1) * BLOCKSIZE + (char*)((xbt_mheap_t)state->s_heap)->heapbase));
491            
492         if(i2 == i1){
493           i2++;
494           continue;
495         }
496
497         if(state->heapinfo2[i2].type != 0){
498           i2++;
499           continue;
500         }
501     
502         if(state->equals_to2[i2][0] != NULL){
503           i2++;
504           continue;
505         }
506           
507         res_compare = compare_heap_area(addr_block1, addr_block2, NULL, info, other_info, NULL, 0);
508         
509         if(res_compare != 1 ){
510           for(k=1; k < state->heapinfo2[i2].busy_block.size; k++)
511             state->equals_to2[i2+k][0] = new_heap_area(i1, -1);
512           for(k=1; k < state->heapinfo1[i1].busy_block.size; k++)
513             state->equals_to1[i1+k][0] = new_heap_area(i2, -1);
514           equal = 1;
515           i1 += state->heapinfo1[i1].busy_block.size;
516         }
517
518         xbt_dynar_reset(previous);
519
520         i2++;
521
522       }
523
524       if(!equal){
525         XBT_DEBUG("Block %zu not found (size_used = %zu, addr = %p)", i1, state->heapinfo1[i1].busy_block.busy_size, addr_block1);
526         i1 = state->heaplimit + 1;
527         nb_diff1++;
528           //i1++;
529       }
530       
531     }else{ /* Fragmented block */
532
533       for(j1=0; j1 < (size_t) (BLOCKSIZE >> state->heapinfo1[i1].type); j1++){
534
535         if(state->heapinfo1[i1].busy_frag.frag_size[j1] == -1) /* Free fragment */
536           continue;
537
538         if(state->equals_to1[i1][j1] != NULL)
539           continue;
540
541         addr_frag1 = (void*) ((char *)addr_block1 + (j1 << state->heapinfo1[i1].type));
542
543         i2 = 1;
544         equal = 0;
545         
546         /* Try first to associate to same fragment in the other heap */
547         if(state->heapinfo2[i1].type == state->heapinfo1[i1].type){
548
549           if(state->equals_to2[i1][j1] == NULL){
550
551             addr_block2 = ((void*) (((ADDR2UINT(i1)) - 1) * BLOCKSIZE + (char*)((xbt_mheap_t)state->s_heap)->heapbase));
552             addr_frag2 = (void*) ((char *)addr_block2 + (j1 << ((xbt_mheap_t)state->s_heap)->heapinfo[i1].type));
553
554             res_compare = compare_heap_area(addr_frag1, addr_frag2, NULL, info, other_info, NULL, 0);
555
556             if(res_compare !=  1)
557               equal = 1;
558         
559             xbt_dynar_reset(previous);
560
561           }
562
563         }
564
565         while(i2 <= state->heaplimit && !equal){
566
567           if(state->heapinfo2[i2].type <= 0){
568             i2++;
569             continue;
570           }
571
572           for(j2=0; j2 < (size_t) (BLOCKSIZE >> state->heapinfo2[i2].type); j2++){
573
574             if(i2 == i1 && j2 == j1)
575               continue;
576            
577             if(state->equals_to2[i2][j2] != NULL)
578               continue;
579                           
580             addr_block2 = ((void*) (((ADDR2UINT(i2)) - 1) * BLOCKSIZE + (char*)((xbt_mheap_t)state->s_heap)->heapbase));
581             addr_frag2 = (void*) ((char *)addr_block2 + (j2 <<((xbt_mheap_t)state->s_heap)->heapinfo[i2].type));
582
583             res_compare = compare_heap_area(addr_frag1, addr_frag2, NULL, info, other_info, NULL, 0);
584             
585             if(res_compare != 1){
586               equal = 1;
587               xbt_dynar_reset(previous);
588               break;
589             }
590
591             xbt_dynar_reset(previous);
592
593           }
594
595           i2++;
596
597         }
598
599         if(!equal){
600           XBT_DEBUG("Block %zu, fragment %zu not found (size_used = %zd, address = %p)\n", i1, j1, state->heapinfo1[i1].busy_frag.frag_size[j1], addr_frag1);
601           i2 = state->heaplimit + 1;
602           i1 = state->heaplimit + 1;
603           nb_diff1++;
604           break;
605         }
606
607       }
608
609       i1++;
610       
611     }
612
613   }
614
615   /* All blocks/fragments are equal to another block/fragment ? */
616   size_t i = 1, j = 0;
617   void *real_addr_frag1 = NULL, *real_addr_block1 = NULL, *real_addr_block2 = NULL, *real_addr_frag2 = NULL;
618  
619   while(i<=state->heaplimit){
620     if(state->heapinfo1[i].type == 0){
621       if(i1 == state->heaplimit){
622         if(state->heapinfo1[i].busy_block.busy_size > 0){
623           if(state->equals_to1[i][0] == NULL){
624             if(XBT_LOG_ISENABLED(mm_diff, xbt_log_priority_debug)){
625               addr_block1 = ((void*) (((ADDR2UINT(i)) - 1) * BLOCKSIZE + (char*)state->heapbase1));
626               XBT_DEBUG("Block %zu (%p) not found (size used = %zu)", i, addr_block1, state->heapinfo1[i].busy_block.busy_size);
627               //mmalloc_backtrace_block_display((void*)heapinfo1, i);
628             }
629             nb_diff1++;
630           }
631         }
632       }
633     }
634     if(state->heapinfo1[i].type > 0){
635       addr_block1 = ((void*) (((ADDR2UINT(i)) - 1) * BLOCKSIZE + (char*)state->heapbase1));
636       real_addr_block1 =  ((void*) (((ADDR2UINT(i)) - 1) * BLOCKSIZE + (char*)((struct mdesc *)state->s_heap)->heapbase));
637       for(j=0; j < (size_t) (BLOCKSIZE >> state->heapinfo1[i].type); j++){
638         if(i1== state->heaplimit){
639           if(state->heapinfo1[i].busy_frag.frag_size[j] > 0){
640             if(state->equals_to1[i][j] == NULL){
641               if(XBT_LOG_ISENABLED(mm_diff, xbt_log_priority_debug)){
642                 addr_frag1 = (void*) ((char *)addr_block1 + (j << state->heapinfo1[i].type));
643                 real_addr_frag1 = (void*) ((char *)real_addr_block1 + (j << ((struct mdesc *)state->s_heap)->heapinfo[i].type));
644                 XBT_DEBUG("Block %zu, Fragment %zu (%p - %p) not found (size used = %zd)", i, j, addr_frag1, real_addr_frag1, state->heapinfo1[i].busy_frag.frag_size[j]);
645                 //mmalloc_backtrace_fragment_display((void*)heapinfo1, i, j);
646               }
647               nb_diff1++;
648             }
649           }
650         }
651       }
652     }
653     i++; 
654   }
655
656   if(i1 == state->heaplimit)
657     XBT_DEBUG("Number of blocks/fragments not found in heap1 : %d", nb_diff1);
658
659   i = 1;
660
661   while(i<=state->heaplimit){
662     if(state->heapinfo2[i].type == 0){
663       if(i1 == state->heaplimit){
664         if(state->heapinfo2[i].busy_block.busy_size > 0){
665           if(state->equals_to2[i][0] == NULL){
666             if(XBT_LOG_ISENABLED(mm_diff, xbt_log_priority_debug)){
667               addr_block2 = ((void*) (((ADDR2UINT(i)) - 1) * BLOCKSIZE + (char*)state->heapbase2));
668               XBT_DEBUG("Block %zu (%p) not found (size used = %zu)", i, addr_block2, state->heapinfo2[i].busy_block.busy_size);
669               //mmalloc_backtrace_block_display((void*)heapinfo2, i);
670             }
671             nb_diff2++;
672           }
673         }
674       }
675     }
676     if(state->heapinfo2[i].type > 0){
677       addr_block2 = ((void*) (((ADDR2UINT(i)) - 1) * BLOCKSIZE + (char*)state->heapbase2));
678       real_addr_block2 =  ((void*) (((ADDR2UINT(i)) - 1) * BLOCKSIZE + (char*)((struct mdesc *)state->s_heap)->heapbase));
679       for(j=0; j < (size_t) (BLOCKSIZE >> state->heapinfo2[i].type); j++){
680         if(i1 == state->heaplimit){
681           if(state->heapinfo2[i].busy_frag.frag_size[j] > 0){
682             if(state->equals_to2[i][j] == NULL){
683               if(XBT_LOG_ISENABLED(mm_diff, xbt_log_priority_debug)){
684                 addr_frag2 = (void*) ((char *)addr_block2 + (j << state->heapinfo2[i].type));
685                 real_addr_frag2 = (void*) ((char *)real_addr_block2 + (j << ((struct mdesc *)state->s_heap)->heapinfo[i].type));
686                 XBT_DEBUG( "Block %zu, Fragment %zu (%p - %p) not found (size used = %zd)", i, j, addr_frag2, real_addr_frag2, state->heapinfo2[i].busy_frag.frag_size[j]);
687                 //mmalloc_backtrace_fragment_display((void*)heapinfo2, i, j);
688               }
689               nb_diff2++;
690             }
691           }
692         }
693       }
694     }
695     i++; 
696   }
697
698   if(i1 == state->heaplimit)
699     XBT_DEBUG("Number of blocks/fragments not found in heap2 : %d", nb_diff2);
700
701   xbt_dynar_free(&previous);
702   real_addr_frag1 = NULL, real_addr_block1 = NULL, real_addr_block2 = NULL, real_addr_frag2 = NULL;
703
704   return ((nb_diff1 > 0) || (nb_diff2 > 0));
705 }
706
707 static int compare_heap_area_without_type(struct s_mm_diff *state, void *real_area1, void *real_area2, void *area1, void *area2, xbt_dynar_t previous, mc_object_info_t info, mc_object_info_t other_info, int size, int check_ignore){
708
709   int i = 0;
710   void *addr_pointed1, *addr_pointed2;
711   int pointer_align, res_compare;
712   ssize_t ignore1, ignore2;
713
714   while(i<size){
715
716     if(check_ignore > 0){
717       if((ignore1 = heap_comparison_ignore_size(state->to_ignore1, (char *)real_area1 + i)) != -1){
718         if((ignore2 = heap_comparison_ignore_size(state->to_ignore2, (char *)real_area2 + i))  == ignore1){
719           if(ignore1 == 0){
720             check_ignore--;
721             return 0;
722           }else{
723             i = i + ignore2;
724             check_ignore--;
725             continue;
726           }
727         }
728       }
729     }
730
731     if(memcmp(((char *)area1) + i, ((char *)area2) + i, 1) != 0){
732
733       pointer_align = (i / sizeof(void*)) * sizeof(void*);
734       addr_pointed1 = *((void **)((char *)area1 + pointer_align));
735       addr_pointed2 = *((void **)((char *)area2 + pointer_align));
736       
737       if(addr_pointed1 > maestro_stack_start && addr_pointed1 < maestro_stack_end && addr_pointed2 > maestro_stack_start && addr_pointed2 < maestro_stack_end){
738         i = pointer_align + sizeof(void *);
739         continue;
740       }else if((addr_pointed1 > state->s_heap) && ((char *)addr_pointed1 < (char *)state->s_heap + STD_HEAP_SIZE)
741                && (addr_pointed2 > state->s_heap) && ((char *)addr_pointed2 < (char *)state->s_heap + STD_HEAP_SIZE)){
742         res_compare = compare_heap_area(addr_pointed1, addr_pointed2, previous, info, other_info, NULL, 0);
743         if(res_compare == 1){
744           return res_compare;
745         }
746         i = pointer_align + sizeof(void *);
747         continue;
748       }else{
749         return 1;
750       }
751       
752     }
753     
754     i++;
755
756   }
757
758   return 0;
759  
760 }
761
762 // area_size is either a byte_size or an elements_count?&
763 static int compare_heap_area_with_type(struct s_mm_diff *state, void *real_area1, void *real_area2, void *area1, void *area2,
764                                        xbt_dynar_t previous, mc_object_info_t info, mc_object_info_t other_info, char *type_id,
765                                        int area_size, int check_ignore, int pointer_level){
766
767   if(is_stack(real_area1) && is_stack(real_area2))
768     return 0;
769
770   ssize_t ignore1, ignore2;
771
772   if((check_ignore > 0) && ((ignore1 = heap_comparison_ignore_size(state->to_ignore1, real_area1)) > 0) && ((ignore2 = heap_comparison_ignore_size(state->to_ignore2, real_area2))  == ignore1)){
773     return 0;
774   }
775   
776   dw_type_t type = xbt_dict_get_or_null(info->types, type_id);
777   dw_type_t subtype, subsubtype;
778   int res, elm_size, i, switch_types = 0;
779   unsigned int cursor = 0;
780   dw_type_t member;
781   void *addr_pointed1, *addr_pointed2;;
782
783   switch(type->type){
784   case DW_TAG_unspecified_type:
785     return 1;
786
787   case DW_TAG_base_type:
788     if(type->name!=NULL && strcmp(type->name, "char") == 0){ /* String, hence random (arbitrary ?) size */
789       if(real_area1 == real_area2)
790         return -1;
791       else
792         return (memcmp(area1, area2, area_size) != 0);
793     }else{
794       if(area_size != -1 && type->byte_size != area_size)
795         return -1;
796       else{
797         return  (memcmp(area1, area2, type->byte_size) != 0);
798       }
799     }
800     break;
801   case DW_TAG_enumeration_type:
802     if(area_size != -1 && type->byte_size != area_size)
803       return -1;
804     else
805       return (memcmp(area1, area2, type->byte_size) != 0);
806     break;
807   case DW_TAG_typedef:
808   case DW_TAG_const_type:
809   case DW_TAG_volatile_type:
810     return compare_heap_area_with_type(state, real_area1, real_area2, area1, area2, previous, info, other_info, type->dw_type_id, area_size, check_ignore, pointer_level);
811     break;
812   case DW_TAG_array_type:
813     subtype = xbt_dict_get_or_null(info->types, type->dw_type_id);
814     switch(subtype->type){
815     case DW_TAG_unspecified_type:
816       return 1;
817
818     case DW_TAG_base_type:
819     case DW_TAG_enumeration_type:
820     case DW_TAG_pointer_type:
821     case DW_TAG_reference_type:
822     case DW_TAG_rvalue_reference_type:
823     case DW_TAG_structure_type:
824     case DW_TAG_class_type:
825     case DW_TAG_union_type:
826       if(subtype->byte_size == 0){ /*declaration of the type, need the complete description */
827           subtype = xbt_dict_get_or_null(other_info->types_by_name, subtype->name);
828           switch_types = 1;
829       }
830       elm_size = subtype->byte_size;
831       break;
832     // TODO, just remove the type indirection?
833     case DW_TAG_const_type:
834     case DW_TAG_typedef:
835     case DW_TAG_volatile_type:
836       subsubtype = subtype->subtype;
837       if(subsubtype->byte_size == 0){ /*declaration of the type, need the complete description */
838           subsubtype = xbt_dict_get_or_null(other_info->types_by_name, subtype->name);
839           switch_types = 1;
840       }
841       elm_size = subsubtype->byte_size;
842       break;
843     default : 
844       return 0;
845       break;
846     }
847     for(i=0; i<type->element_count; i++){
848       // TODO, add support for variable stride (DW_AT_byte_stride)
849       if(switch_types)
850         res = compare_heap_area_with_type(state, (char *)real_area1 + (i*elm_size), (char *)real_area2 + (i*elm_size), (char *)area1 + (i*elm_size), (char *)area2 + (i*elm_size), previous, other_info, info, type->dw_type_id, subtype->byte_size, check_ignore, pointer_level);
851       else
852         res = compare_heap_area_with_type(state, (char *)real_area1 + (i*elm_size), (char *)real_area2 + (i*elm_size), (char *)area1 + (i*elm_size), (char *)area2 + (i*elm_size), previous, info, other_info, type->dw_type_id, subtype->byte_size, check_ignore, pointer_level);
853       if(res == 1)
854         return res;
855     }
856     break;
857   case DW_TAG_reference_type:
858   case DW_TAG_rvalue_reference_type:
859   case DW_TAG_pointer_type:
860     if(type->dw_type_id && ((dw_type_t)xbt_dict_get_or_null(info->types, type->dw_type_id))->type == DW_TAG_subroutine_type){
861       addr_pointed1 = *((void **)(area1)); 
862       addr_pointed2 = *((void **)(area2));
863       return (addr_pointed1 != addr_pointed2);;
864     }else{
865       pointer_level++;
866       if(pointer_level > 1){ /* Array of pointers */
867         for(i=0; i<(area_size/sizeof(void *)); i++){ 
868           addr_pointed1 = *((void **)((char *)area1 + (i*sizeof(void *)))); 
869           addr_pointed2 = *((void **)((char *)area2 + (i*sizeof(void *)))); 
870           if(addr_pointed1 > state->s_heap && (char *)addr_pointed1 < (char*) state->s_heap + STD_HEAP_SIZE && addr_pointed2 > state->s_heap && (char *)addr_pointed2 < (char*) state->s_heap + STD_HEAP_SIZE)
871             res =  compare_heap_area(addr_pointed1, addr_pointed2, previous, info, other_info, type->dw_type_id, pointer_level);
872           else
873             res =  (addr_pointed1 != addr_pointed2);
874           if(res == 1)
875             return res;
876         }
877       }else{
878         addr_pointed1 = *((void **)(area1)); 
879         addr_pointed2 = *((void **)(area2));
880         if(addr_pointed1 > state->s_heap && (char *)addr_pointed1 < (char*) state->s_heap + STD_HEAP_SIZE && addr_pointed2 > state->s_heap && (char *)addr_pointed2 < (char*) state->s_heap + STD_HEAP_SIZE)
881           return compare_heap_area(addr_pointed1, addr_pointed2, previous, info, other_info, type->dw_type_id, pointer_level);
882         else
883           return  (addr_pointed1 != addr_pointed2);
884       }
885     }
886     break;
887   case DW_TAG_structure_type:
888   case DW_TAG_class_type:
889     if(type->byte_size == 0){ /*declaration of the structure, need the complete description */
890       dw_type_t full_type = xbt_dict_get_or_null(info->types_by_name, type->name);
891       if(full_type){
892         type = full_type;
893       }else{
894         type = xbt_dict_get_or_null(other_info->types_by_name, type->name);
895         switch_types = 1;
896       }
897     }
898     if(area_size != -1 && type->byte_size != area_size){
899       if(area_size>type->byte_size && area_size%type->byte_size == 0){
900         for(i=0; i<(area_size/type->byte_size); i++){
901           if(switch_types)
902             res = compare_heap_area_with_type(state, (char *)real_area1 + (i*type->byte_size), (char *)real_area2 + (i*type->byte_size), (char *)area1 + (i*type->byte_size), (char *)area2 + (i*type->byte_size), previous, other_info, info, type_id, -1, check_ignore, 0);
903           else
904             res = compare_heap_area_with_type(state, (char *)real_area1 + (i*type->byte_size), (char *)real_area2 + (i*type->byte_size), (char *)area1 + (i*type->byte_size), (char *)area2 + (i*type->byte_size), previous, info, other_info, type_id, -1, check_ignore, 0);
905           if(res == 1)
906             return res;
907         }
908       }else{
909         return -1;
910       }
911     }else{
912       cursor = 0;
913       xbt_dynar_foreach(type->members, cursor, member){ 
914         if(switch_types)
915           res = compare_heap_area_with_type(state, (char *)real_area1 + member->offset, (char *)real_area2 + member->offset, (char *)area1 + member->offset, (char *)area2 + member->offset, previous, other_info, info, member->dw_type_id, -1, check_ignore, 0);
916         else
917           res = compare_heap_area_with_type(state, (char *)real_area1 + member->offset, (char *)real_area2 + member->offset, (char *)area1 + member->offset, (char *)area2 + member->offset, previous, info, other_info, member->dw_type_id, -1, check_ignore, 0);
918         if(res == 1){
919           return res;
920         }
921       }
922     }
923     break;
924   case DW_TAG_union_type:
925     return compare_heap_area_without_type(state, real_area1, real_area2, area1, area2, previous, info, other_info, type->byte_size, check_ignore);
926     break;
927   default:
928     break;
929   }
930
931   return 0;
932
933 }
934
935 static char* get_offset_type(char* type_id, int offset, mc_object_info_t info, mc_object_info_t other_info, int area_size, int *switch_type){
936   dw_type_t type = xbt_dict_get_or_null(info->types, type_id);
937   if(type == NULL){
938     type = xbt_dict_get_or_null(other_info->types, type_id);
939     *switch_type = 1;
940   }
941   switch(type->type){
942   case DW_TAG_structure_type :
943   case DW_TAG_class_type:
944     if(type->byte_size == 0){ /*declaration of the structure, need the complete description */
945       if(*switch_type == 0){
946         dw_type_t full_type = xbt_dict_get_or_null(info->types_by_name, type->name);
947         if(full_type){
948           type = full_type;
949         }else{
950           type = xbt_dict_get_or_null(other_info->types_by_name, type->name);
951           *switch_type = 1;
952         }
953       }else{
954         dw_type_t full_type = xbt_dict_get_or_null(other_info->types_by_name, type->name);
955         if(full_type){
956           type = full_type;
957         }else{
958           type = xbt_dict_get_or_null(info->types_by_name, type->name);
959           *switch_type = 0;
960         }
961       }
962     
963     }
964     if(area_size != -1 && type->byte_size != area_size){
965       if(area_size>type->byte_size && area_size%type->byte_size == 0)
966         return type_id;
967       else
968         return NULL;
969     }else{
970       unsigned int cursor = 0;
971       dw_type_t member;
972       xbt_dynar_foreach(type->members, cursor, member){ 
973         if(member->offset == offset)
974           return member->dw_type_id;
975       }
976       return NULL;
977     }
978     break;
979   default:
980     /* FIXME : other cases ? */
981     return NULL;
982     break;
983   }
984 }
985
986 int compare_heap_area(void *area1, void* area2, xbt_dynar_t previous, mc_object_info_t info, mc_object_info_t other_info, char *type_id, int pointer_level){
987
988   struct s_mm_diff* state = mm_diff_info;
989
990   int res_compare;
991   ssize_t block1, frag1, block2, frag2;
992   ssize_t size;
993   int check_ignore = 0;
994
995   void *addr_block1, *addr_block2, *addr_frag1, *addr_frag2, *real_addr_block1, *real_addr_block2,  *real_addr_frag1, *real_addr_frag2;
996   void *area1_to_compare, *area2_to_compare;
997   dw_type_t type = NULL;
998   int type_size = -1;
999   int offset1 =0, offset2 = 0;
1000   int new_size1 = -1, new_size2 = -1;
1001   char *new_type_id1 = NULL, *new_type_id2 = NULL;
1002   int switch_type = 0;
1003
1004   int match_pairs = 0;
1005
1006   if(previous == NULL){
1007     previous = xbt_dynar_new(sizeof(heap_area_pair_t), heap_area_pair_free_voidp);
1008     match_pairs = 1;
1009   }
1010
1011   block1 = ((char*)area1 - (char*)((xbt_mheap_t)state->s_heap)->heapbase) / BLOCKSIZE + 1;
1012   block2 = ((char*)area2 - (char*)((xbt_mheap_t)state->s_heap)->heapbase) / BLOCKSIZE + 1;
1013
1014   if(is_block_stack((int)block1) && is_block_stack((int)block2)){
1015     add_heap_area_pair(previous, block1, -1, block2, -1);
1016     if(match_pairs){
1017       match_equals(state, previous);
1018       xbt_dynar_free(&previous);
1019     }
1020     return 0;
1021   }
1022
1023   if(((char *)area1 < (char*)((xbt_mheap_t)state->s_heap)->heapbase)  || (block1 > state->heapsize1) || (block1 < 1) || ((char *)area2 < (char*)((xbt_mheap_t)state->s_heap)->heapbase) || (block2 > state->heapsize2) || (block2 < 1)){
1024     if(match_pairs){
1025       xbt_dynar_free(&previous);
1026     }
1027     return 1;
1028   }
1029
1030   addr_block1 = ((void*) (((ADDR2UINT(block1)) - 1) * BLOCKSIZE + (char*)state->heapbase1));
1031   addr_block2 = ((void*) (((ADDR2UINT(block2)) - 1) * BLOCKSIZE + (char*)state->heapbase2));
1032
1033   real_addr_block1 = ((void*) (((ADDR2UINT(block1)) - 1) * BLOCKSIZE + (char*)((xbt_mheap_t)state->s_heap)->heapbase));
1034   real_addr_block2 = ((void*) (((ADDR2UINT(block2)) - 1) * BLOCKSIZE + (char*)((xbt_mheap_t)state->s_heap)->heapbase));
1035
1036   if(type_id){
1037     type = xbt_dict_get_or_null(info->types, type_id);
1038     if(type->byte_size == 0){
1039       if(type->subtype == NULL){
1040         dw_type_t full_type = xbt_dict_get_or_null(info->types_by_name, type->name);
1041         if(full_type)
1042           type = full_type;
1043         else
1044           type = xbt_dict_get_or_null(other_info->types_by_name, type->name);
1045       }else{
1046         type = type->subtype;
1047       }
1048     }
1049     if((type->type == DW_TAG_pointer_type) || ((type->type == DW_TAG_base_type) && type->name!=NULL && (!strcmp(type->name, "char"))))
1050       type_size = -1;
1051     else
1052       type_size = type->byte_size;
1053   }
1054   
1055   if((state->heapinfo1[block1].type == -1) && (state->heapinfo2[block2].type == -1)){  /* Free block */
1056
1057     if(match_pairs){
1058       match_equals(state, previous);
1059       xbt_dynar_free(&previous);
1060     }
1061     return 0;
1062
1063   }else if((state->heapinfo1[block1].type == 0) && (state->heapinfo2[block2].type == 0)){ /* Complete block */
1064     
1065     if(state->equals_to1[block1][0] != NULL && state->equals_to2[block2][0] != NULL){
1066       if(equal_blocks(state, block1, block2)){
1067         if(match_pairs){
1068           match_equals(state, previous);
1069           xbt_dynar_free(&previous);
1070         }
1071         return 0;
1072       }
1073     }
1074
1075     if(type_size != -1){
1076       if(type_size != state->heapinfo1[block1].busy_block.busy_size
1077         && type_size != state->heapinfo2[block2].busy_block.busy_size
1078         && type->name!=NULL && !strcmp(type->name, "s_smx_context")){
1079         if(match_pairs){
1080           match_equals(state, previous);
1081           xbt_dynar_free(&previous);
1082         }
1083         return -1;
1084       }
1085     }
1086
1087     if(state->heapinfo1[block1].busy_block.size != state->heapinfo2[block2].busy_block.size){
1088       if(match_pairs){
1089         xbt_dynar_free(&previous);
1090       }
1091       return 1;
1092     }
1093
1094     if(state->heapinfo1[block1].busy_block.busy_size != state->heapinfo2[block2].busy_block.busy_size){
1095       if(match_pairs){
1096         xbt_dynar_free(&previous);
1097       }
1098       return 1;
1099     }
1100
1101     if(!add_heap_area_pair(previous, block1, -1, block2, -1)){
1102       if(match_pairs){
1103         match_equals(state, previous);
1104         xbt_dynar_free(&previous);
1105       }
1106       return 0;
1107     }
1108  
1109     size = state->heapinfo1[block1].busy_block.busy_size;
1110     
1111     if(type_id != NULL){
1112       xbt_free(state->types1[block1][0]);
1113       xbt_free(state->types2[block2][0]);
1114       state->types1[block1][0] = strdup(type_id);
1115       state->types2[block2][0] = strdup(type_id);
1116     }
1117
1118     if(size <= 0){
1119       if(match_pairs){
1120         match_equals(state, previous);
1121         xbt_dynar_free(&previous);
1122       }
1123       return 0;
1124     }
1125
1126     frag1 = -1;
1127     frag2 = -1;
1128
1129     area1_to_compare = addr_block1;
1130     area2_to_compare = addr_block2;
1131
1132     if((state->heapinfo1[block1].busy_block.ignore > 0) && (state->heapinfo2[block2].busy_block.ignore == state->heapinfo1[block1].busy_block.ignore))
1133       check_ignore = state->heapinfo1[block1].busy_block.ignore;
1134       
1135   }else if((state->heapinfo1[block1].type > 0) && (state->heapinfo2[block2].type > 0)){ /* Fragmented block */
1136
1137     frag1 = ((uintptr_t) (ADDR2UINT (area1) % (BLOCKSIZE))) >> state->heapinfo1[block1].type;
1138     frag2 = ((uintptr_t) (ADDR2UINT (area2) % (BLOCKSIZE))) >> state->heapinfo2[block2].type;
1139
1140     addr_frag1 = (void*) ((char *)addr_block1 + (frag1 << state->heapinfo1[block1].type));
1141     addr_frag2 = (void*) ((char *)addr_block2 + (frag2 << state->heapinfo2[block2].type));
1142
1143     real_addr_frag1 = (void*) ((char *)real_addr_block1 + (frag1 << ((xbt_mheap_t)state->s_heap)->heapinfo[block1].type));
1144     real_addr_frag2 = (void*) ((char *)real_addr_block2 + (frag2 << ((xbt_mheap_t)state->s_heap)->heapinfo[block2].type));
1145
1146     if(type_size != -1){
1147       if(state->heapinfo1[block1].busy_frag.frag_size[frag1] == -1 || state->heapinfo2[block2].busy_frag.frag_size[frag2] == -1){
1148         if(match_pairs){
1149           match_equals(state, previous);
1150           xbt_dynar_free(&previous);
1151         }
1152         return -1;
1153       }
1154       if(type_size != state->heapinfo1[block1].busy_frag.frag_size[frag1]|| type_size !=  state->heapinfo2[block2].busy_frag.frag_size[frag2]){
1155         if(match_pairs){
1156           match_equals(state, previous);
1157           xbt_dynar_free(&previous);
1158         }
1159         return -1;
1160       }
1161     }
1162
1163     if(state->equals_to1[block1][frag1] != NULL && state->equals_to2[block2][frag2] != NULL){
1164       if(equal_fragments(state, block1, frag1, block2, frag2)){
1165         if(match_pairs){
1166           match_equals(state, previous);
1167           xbt_dynar_free(&previous);
1168         }
1169         return 0;
1170       }
1171     }
1172
1173     if(state->heapinfo1[block1].busy_frag.frag_size[frag1] != state->heapinfo2[block2].busy_frag.frag_size[frag2]){
1174       if(type_size == -1){
1175          if(match_pairs){
1176           match_equals(state, previous);
1177           xbt_dynar_free(&previous);
1178         }
1179         return -1;
1180       }else{
1181         if(match_pairs){
1182           xbt_dynar_free(&previous);
1183         }
1184         return 1;
1185       }
1186     }
1187       
1188     size = state->heapinfo1[block1].busy_frag.frag_size[frag1];
1189
1190     if(type_id != NULL){
1191       xbt_free(state->types1[block1][frag1]);
1192       xbt_free(state->types2[block2][frag2]);
1193       state->types1[block1][frag1] = strdup(type_id);
1194       state->types2[block2][frag2] = strdup(type_id);
1195     }
1196
1197     if(real_addr_frag1 != area1 || real_addr_frag2 != area2){
1198       offset1 = (char *)area1 - (char *)real_addr_frag1;
1199       offset2 = (char *)area2 - (char *)real_addr_frag2;
1200       if(state->types1[block1][frag1] != NULL && state->types2[block2][frag2] != NULL){
1201         new_type_id1 = get_offset_type(state->types1[block1][frag1], offset1, info, other_info, size, &switch_type);
1202         new_type_id2 = get_offset_type(state->types2[block2][frag2], offset1, info, other_info, size, &switch_type);
1203       }else if(state->types1[block1][frag1] != NULL){
1204         new_type_id1 = get_offset_type(state->types1[block1][frag1], offset1, info, other_info, size, &switch_type);
1205         new_type_id2 = get_offset_type(state->types1[block1][frag1], offset2, info, other_info, size, &switch_type);
1206       }else if(state->types2[block2][frag2] != NULL){
1207         new_type_id1 = get_offset_type(state->types2[block2][frag2], offset1, info, other_info, size, &switch_type);
1208         new_type_id2 = get_offset_type(state->types2[block2][frag2], offset2, info, other_info, size, &switch_type);
1209       }else{
1210         if(match_pairs){
1211           match_equals(state, previous);
1212           xbt_dynar_free(&previous);
1213         }
1214         return -1;
1215       }   
1216
1217       if(new_type_id1 !=  NULL && new_type_id2 !=  NULL && !strcmp(new_type_id1, new_type_id2)){
1218         if(switch_type){
1219           type = xbt_dict_get_or_null(other_info->types, new_type_id1);
1220           while(type->byte_size == 0 && type->dw_type_id != NULL)
1221             type = xbt_dict_get_or_null(other_info->types, type->dw_type_id);
1222           new_size1 = type->byte_size;
1223           type = xbt_dict_get_or_null(other_info->types, new_type_id2);
1224           while(type->byte_size == 0 && type->dw_type_id != NULL)
1225             type = xbt_dict_get_or_null(other_info->types, type->dw_type_id);
1226           new_size2 = type->byte_size;
1227         }else{
1228           type = xbt_dict_get_or_null(info->types, new_type_id1);
1229           while(type->byte_size == 0 && type->dw_type_id != NULL)
1230             type = xbt_dict_get_or_null(info->types, type->dw_type_id);
1231           new_size1 = type->byte_size;
1232           type = xbt_dict_get_or_null(info->types, new_type_id2);
1233           while(type->byte_size == 0 && type->dw_type_id != NULL)
1234             type = xbt_dict_get_or_null(info->types, type->dw_type_id);
1235           new_size2 = type->byte_size;
1236         }
1237       }else{
1238         if(match_pairs){
1239           match_equals(state, previous);
1240           xbt_dynar_free(&previous);
1241         }
1242         return -1;
1243       }
1244     }
1245
1246     area1_to_compare = (char *)addr_frag1 + offset1;
1247     area2_to_compare = (char *)addr_frag2 + offset2;
1248     
1249     if(new_size1 > 0 && new_size1 == new_size2){
1250       type_id = new_type_id1;
1251       size = new_size1;
1252     }
1253
1254     if(offset1 == 0 && offset2 == 0){
1255       if(!add_heap_area_pair(previous, block1, frag1, block2, frag2)){
1256         if(match_pairs){
1257           match_equals(state, previous);
1258           xbt_dynar_free(&previous);
1259         }
1260         return 0;
1261       }
1262     }
1263
1264     if(size <= 0){
1265       if(match_pairs){
1266         match_equals(state, previous);
1267         xbt_dynar_free(&previous);
1268       }
1269       return 0;
1270     }
1271       
1272     if((state->heapinfo1[block1].busy_frag.ignore[frag1] > 0) && ( state->heapinfo2[block2].busy_frag.ignore[frag2] == state->heapinfo1[block1].busy_frag.ignore[frag1]))
1273       check_ignore = state->heapinfo1[block1].busy_frag.ignore[frag1];
1274     
1275   }else{
1276
1277     if(match_pairs){
1278       xbt_dynar_free(&previous);
1279     }
1280     return 1;
1281
1282   }
1283   
1284
1285   /* Start comparison*/
1286   if(type_id != NULL){
1287     if(switch_type)
1288       res_compare = compare_heap_area_with_type(state, area1, area2, area1_to_compare, area2_to_compare, previous, other_info, info, type_id, size, check_ignore, pointer_level);
1289     else
1290       res_compare = compare_heap_area_with_type(state, area1, area2, area1_to_compare, area2_to_compare, previous, info, other_info, type_id, size, check_ignore, pointer_level);
1291     if(res_compare == 1){
1292       if(match_pairs)
1293         xbt_dynar_free(&previous);
1294       return res_compare;
1295     }
1296   }else{
1297     if(switch_type)
1298       res_compare = compare_heap_area_without_type(state, area1, area2, area1_to_compare, area2_to_compare, previous, other_info, info, size, check_ignore);
1299     else
1300       res_compare = compare_heap_area_without_type(state, area1, area2, area1_to_compare, area2_to_compare, previous, info, other_info, size, check_ignore);
1301     if(res_compare == 1){
1302       if(match_pairs)
1303         xbt_dynar_free(&previous);
1304       return res_compare;
1305     }
1306   }
1307
1308   if(match_pairs){
1309     match_equals(state, previous);
1310     xbt_dynar_free(&previous);
1311   }
1312
1313   return 0;
1314 }
1315
1316 /*********************************************** Miscellaneous ***************************************************/
1317 /****************************************************************************************************************/
1318
1319 // Not used:
1320 static int get_pointed_area_size(void *area, int heap){
1321
1322   struct s_mm_diff *state = mm_diff_info;
1323
1324   int block, frag;
1325   malloc_info *heapinfo;
1326
1327   if(heap == 1)
1328     heapinfo = state->heapinfo1;
1329   else
1330     heapinfo = state->heapinfo2;
1331
1332   block = ((char*)area - (char*)((xbt_mheap_t)state->s_heap)->heapbase) / BLOCKSIZE + 1;
1333
1334   if(((char *)area < (char*)((xbt_mheap_t)state->s_heap)->heapbase)  || (block > state->heapsize1) || (block < 1))
1335     return -1;
1336
1337   if(heapinfo[block].type == -1){ /* Free block */
1338     return -1;  
1339   }else if(heapinfo[block].type == 0){ /* Complete block */
1340     return (int)heapinfo[block].busy_block.busy_size;
1341   }else{
1342     frag = ((uintptr_t) (ADDR2UINT (area) % (BLOCKSIZE))) >> heapinfo[block].type;
1343     return (int)heapinfo[block].busy_frag.frag_size[frag];
1344   }
1345
1346 }
1347
1348 // Not used:
1349 char *get_type_description(mc_object_info_t info, char *type_name){
1350
1351   xbt_dict_cursor_t dict_cursor;
1352   char *type_origin;
1353   dw_type_t type;
1354
1355   xbt_dict_foreach(info->types, dict_cursor, type_origin, type){
1356     if(type->name && (strcmp(type->name, type_name) == 0) && type->byte_size > 0){
1357       xbt_dict_cursor_free(&dict_cursor);
1358       return type_origin;
1359     }
1360   }
1361
1362   xbt_dict_cursor_free(&dict_cursor);
1363   return NULL;
1364 }
1365
1366
1367 #ifndef max
1368 #define max( a, b ) ( ((a) > (b)) ? (a) : (b) )
1369 #endif
1370
1371 // Not used:
1372 int mmalloc_linear_compare_heap(xbt_mheap_t heap1, xbt_mheap_t heap2){
1373
1374   struct s_mm_diff *state = mm_diff_info;
1375
1376   if(heap1 == NULL && heap1 == NULL){
1377     XBT_DEBUG("Malloc descriptors null");
1378     return 0;
1379   }
1380
1381   if(heap1->heaplimit != heap2->heaplimit){
1382     XBT_DEBUG("Different limit of valid info table indices");
1383     return 1;
1384   }
1385
1386   /* Heap information */
1387   state->heaplimit = ((struct mdesc *)heap1)->heaplimit;
1388
1389   state->s_heap = (char *)mmalloc_get_current_heap() - STD_HEAP_SIZE - getpagesize();
1390
1391   state->heapbase1 = (char *)heap1 + BLOCKSIZE;
1392   state->heapbase2 = (char *)heap2 + BLOCKSIZE;
1393
1394   state->heapinfo1 = (malloc_info *)((char *)heap1 + ((uintptr_t)((char *)heap1->heapinfo - (char *)state->s_heap)));
1395   state->heapinfo2 = (malloc_info *)((char *)heap2 + ((uintptr_t)((char *)heap2->heapinfo - (char *)state->s_heap)));
1396
1397   state->heapsize1 = heap1->heapsize;
1398   state->heapsize2 = heap2->heapsize;
1399
1400   /* Start comparison */
1401   size_t i, j, k;
1402   void *addr_block1, *addr_block2, *addr_frag1, *addr_frag2;
1403
1404   int distance = 0;
1405
1406   /* Check busy blocks*/
1407
1408   i = 1;
1409
1410   while(i <= state->heaplimit){
1411
1412     addr_block1 = ((void*) (((ADDR2UINT(i)) - 1) * BLOCKSIZE + (char*)state->heapbase1));
1413     addr_block2 = ((void*) (((ADDR2UINT(i)) - 1) * BLOCKSIZE + (char*)state->heapbase2));
1414
1415     if(state->heapinfo1[i].type != state->heapinfo2[i].type){
1416   
1417       distance += BLOCKSIZE;
1418       XBT_DEBUG("Different type of blocks (%zu) : %d - %d -> distance = %d", i, state->heapinfo1[i].type, state->heapinfo2[i].type, distance);
1419       i++;
1420     
1421     }else{
1422
1423       if(state->heapinfo1[i].type == -1){ /* Free block */
1424         i++;
1425         continue;
1426       }
1427
1428       if(state->heapinfo1[i].type == 0){ /* Large block */
1429        
1430         if(state->heapinfo1[i].busy_block.size != state->heapinfo2[i].busy_block.size){
1431           distance += BLOCKSIZE * max(state->heapinfo1[i].busy_block.size, state->heapinfo2[i].busy_block.size);
1432           i += max(state->heapinfo1[i].busy_block.size, state->heapinfo2[i].busy_block.size);
1433           XBT_DEBUG("Different larger of cluster at block %zu : %zu - %zu -> distance = %d", i, state->heapinfo1[i].busy_block.size, state->heapinfo2[i].busy_block.size, distance);
1434           continue;
1435         }
1436
1437         /*if(heapinfo1[i].busy_block.busy_size != heapinfo2[i].busy_block.busy_size){
1438           distance += max(heapinfo1[i].busy_block.busy_size, heapinfo2[i].busy_block.busy_size);
1439           i += max(heapinfo1[i].busy_block.size, heapinfo2[i].busy_block.size);
1440           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);
1441           continue;
1442           }*/
1443
1444         k = 0;
1445
1446         //while(k < (heapinfo1[i].busy_block.busy_size)){
1447         while(k < state->heapinfo1[i].busy_block.size * BLOCKSIZE){
1448           if(memcmp((char *)addr_block1 + k, (char *)addr_block2 + k, 1) != 0){
1449             distance ++;
1450           }
1451           k++;
1452         } 
1453
1454         i++;
1455
1456       }else { /* Fragmented block */
1457
1458         for(j=0; j < (size_t) (BLOCKSIZE >> state->heapinfo1[i].type); j++){
1459
1460           addr_frag1 = (void*) ((char *)addr_block1 + (j << state->heapinfo1[i].type));
1461           addr_frag2 = (void*) ((char *)addr_block2 + (j << state->heapinfo2[i].type));
1462
1463           if(state->heapinfo1[i].busy_frag.frag_size[j] == 0 && state->heapinfo2[i].busy_frag.frag_size[j] == 0){
1464             continue;
1465           }
1466           
1467           
1468           /*if(heapinfo1[i].busy_frag.frag_size[j] != heapinfo2[i].busy_frag.frag_size[j]){
1469             distance += max(heapinfo1[i].busy_frag.frag_size[j], heapinfo2[i].busy_frag.frag_size[j]);
1470             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); 
1471             continue;
1472             }*/
1473    
1474           k=0;
1475
1476           //while(k < max(heapinfo1[i].busy_frag.frag_size[j], heapinfo2[i].busy_frag.frag_size[j])){
1477           while(k < (BLOCKSIZE / (BLOCKSIZE >> state->heapinfo1[i].type))){
1478             if(memcmp((char *)addr_frag1 + k, (char *)addr_frag2 + k, 1) != 0){
1479               distance ++;
1480             }
1481             k++;
1482           }
1483
1484         }
1485
1486         i++;
1487
1488       }
1489       
1490     }
1491
1492   }
1493
1494   return distance;
1495   
1496 }
1497