/*****************************
 * heap_stats.c
 *
 * Measure various heap performance statistics using 
 * heaps of uniformly distributed randoms.
 *
 *   $ ./heap_stats
 *   -- analyze heap creation --
 *        size     steps 
 *         100        92 
 *         100        90 
 *         100        84 
 *         300       277 
 *         300       276 
 *         300       280 
 *        1000       938 
 *        1000       932 
 *        1000       944 
 *        3000      2811 
 *        3000      2787 
 *        3000      2807 
 *       10000      9389 
 *       10000      9380 
 *       10000      9419 
 *       30000     28257 
 *       30000     28194 
 *       30000     28102 
 *   ... which is pretty clearly O(n) 
 * 
 * $Id: heap_stats.c 12567 2007-04-09 16:12:05Z mahoney $
 ***************/

#include <stdlib.h>
#include <stdio.h>
#include <time.h>
#include "utility.h"
#include "heap.h"

// Make random() start in a different place each time this program runs.
void randomize(){
  // I'm not confident of this recipe, but it seems to work;
  // using time(NULL) alone seemed to give nearly the same starting rand().
  srand(314159*time(NULL));
}

node* make_random_nodelist(int size){
  int i, r;
  node* nodes = (node*) malloc(size*sizeof(node));
  heap h;
  for (i=0; i<size; i++){
    r = random();
    nodes[i] = new_node(r, int2string(r));
  }
  return nodes;
}

int main(){
  int n, i, trial;
  bool descending = TRUE;
  int n_trials = 3;
  int n_sizes = 6;
  int sizes[6] = {100, 300, 1000, 3000, 10000, 30000};
  heap h;
  node* nodes;
  randomize(); 

  printf("-- analyze heap creation --\n");
  printf(" %8s  %8s \n", "size", "steps");
  for (i=0; i<n_sizes; i++){
    for (trial=0; trial<n_trials; trial++){
      nodes = make_random_nodelist(sizes[i]);
      reset_step_counter();
      h = new_heap(sizes[i], nodes, descending);
      printf(" %8i  %8i \n", sizes[i], get_step_counter());
    }
  }

  // // see what rand() is doing ...
  // printf("RAND_MAX = %i \n", RAND_MAX);
  // for (i=1; i<n_times; i++){
  //  random_int = rand();
  //  random_fraction = 1.0*random_int/RAND_MAX;
  //  printf("random_int = %i; fraction=%f \n", random_int, random_fraction);
  // }

}