htswanalysis/MACS/lib/gsl/gsl-1.11/gsl-histogram.c

   1 /* histogram/gsl-histogram.c
   2  *
   3  * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2007 Brian Gough
   4  *
   5  * This program is free software; you can redistribute it and/or modify
   6  * it under the terms of the GNU General Public License as published by
   7  * the Free Software Foundation; either version 3 of the License, or (at
   8  * your option) any later version.
   9  *
  10  * This program is distributed in the hope that it will be useful, but
  11  * WITHOUT ANY WARRANTY; without even the implied warranty of
  12  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13  * General Public License for more details.
  14  *
  15  * You should have received a copy of the GNU General Public License
  16  * along with this program; if not, write to the Free Software
  17  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
  18  */
  19
  20 #include <config.h>
  21 #include <stdio.h>
  22 #include <stdlib.h>
  23 #include <gsl/gsl_histogram.h>
  24
  25 int
  26 main (int argc, char **argv)
  27 {
  28   double a = 0.0, b = 1.0;
  29   size_t n = 10;
  30
  31   if (argc != 3 && argc !=4)
  32     {
  33       printf ("Usage: gsl-histogram xmin xmax [n]\n");
  34       printf (
  35 "Computes a histogram of the data on stdin using n bins from xmin to xmax.\n"
  36 "If n is unspecified then bins of integer width are used.\n");
  37       exit (0);
  38     }
  39
  40   a = atof (argv[1]);
  41   b = atof (argv[2]);
  42
  43   if (argc == 4)
  44     {
  45       n = atoi (argv[3]);
  46     }
  47   else
  48     {
  49       n = (int)(b - a) ;
  50     }
  51
  52   {
  53     double x;
  54     gsl_histogram *h = gsl_histogram_alloc (n);
  55
  56     gsl_histogram_set_ranges_uniform (h, a, b);
  57
  58     while (fscanf(stdin, "%lg", &x) == 1)
  59       {
  60         gsl_histogram_increment(h, x);
  61       }
  62
  63 #ifdef DISPLAY_STATS
  64     {
  65       double mean = gsl_histogram_mean (h);
  66       double sigma = gsl_histogram_sigma (h);
  67       fprintf (stdout, "# mean = %g\n", mean);
  68       fprintf (stdout, "# sigma = %g\n", sigma);
  69     }
  70 #endif
  71
  72     gsl_histogram_fprintf (stdout, h, "%g", "%g");
  73
  74     gsl_histogram_free (h);
  75   }
  76
  77   return 0;
  78 }