histotest.cpp

统计软件包

/** \file

   Test the histogram code

 */
#include <stdio.h>

#include "histogram.h"

double data[] = {68.23,
		 69.40,
		 69.47,
		 69.75,
		 70.00,
		 70.00,
		 70.00,
		 70.59,
		 70.70,
		 70.71,
		 70.99,
		 71.00,
		 71.25,
		 71.39,
		 71.40,
		 71.45,
		 71.49,
		 71.60,
		 71.60,
		 71.74,
		 72.15,
		 72.19,
		 72.25,
		 72.41,
		 72.70,
		 72.70,
		 73.48,
		 73.62,
		 73.90,
		 74.09,
		 74.20,
		 75.20,
		 75.95,
		 76.77,
		 76.90,
		 77.40,
		 77.50,
		 77.75,
		 78.25,
		 78.80,
		 78.85,
		 79.65,
		 80.50,
		 80.71,
		 80.91,
		 80.95,
		 81.87,
		 82.00,
		 82.05,
		 82.25,
		 82.29,
		 82.80,
		 83.00,
		 83.11,
		 83.50,
		 83.75,
		 84.39,
		 84.65,
		 84.80,
		 85.00,
		 85.05,
		 85.46,
		 85.48,
		 86.00,
		 86.40,
		 86.49 };

int
main()
{
  const size_t N = sizeof( data ) / sizeof( double );
  const size_t num_bins = 20;
  histogram::bin_vec binz( num_bins );
  histogram histo;

  histo.calculate( data, N, binz );

  for (size_t i = 0; i < num_bins; i++) {
    printf("%7.4f %2d\n", binz.start(i), binz[i] );
  }
  return 0;
}