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

   1 /* histogram/test_trap.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 <stdlib.h>
  22 #include <stdio.h>
  23 #include <gsl/gsl_histogram.h>
  24 #include <gsl/gsl_test.h>
  25 #include <gsl/gsl_errno.h>
  26 #include <gsl/gsl_ieee_utils.h>
  27
  28 #define N 397
  29
  30 static void my_error_handler (const char *reason, const char *file,
  31                               int line, int err);
  32 static int status = 0;
  33
  34 void
  35 test1d_trap (void)
  36 {
  37   gsl_histogram *h;
  38   double result, lower, upper;
  39   size_t i;
  40
  41   gsl_set_error_handler (&my_error_handler);
  42
  43   gsl_ieee_env_setup ();
  44
  45   status = 0;
  46   h = gsl_histogram_calloc (0);
  47   gsl_test (!status, "gsl_histogram_calloc traps zero-length histogram");
  48   gsl_test (h != 0,
  49             "gsl_histogram_calloc returns NULL for zero-length histogram");
  50
  51   status = 0;
  52   h = gsl_histogram_calloc_uniform (0, 0.0, 1.0);
  53   gsl_test (!status,
  54             "gsl_histogram_calloc_uniform traps zero-length histogram");
  55   gsl_test (h != 0,
  56      "gsl_histogram_calloc_uniform returns NULL for zero-length histogram");
  57
  58   status = 0;
  59   h = gsl_histogram_calloc_uniform (10, 1.0, 1.0);
  60   gsl_test (!status,
  61             "gsl_histogram_calloc_uniform traps equal endpoints");
  62   gsl_test (h != 0,
  63             "gsl_histogram_calloc_uniform returns NULL for equal endpoints");
  64
  65   status = 0;
  66   h = gsl_histogram_calloc_uniform (10, 2.0, 1.0);
  67   gsl_test (!status,
  68             "gsl_histogram_calloc_uniform traps invalid range");
  69   gsl_test (h != 0,
  70             "gsl_histogram_calloc_uniform returns NULL for invalid range");
  71
  72   h = gsl_histogram_calloc_uniform (N, 0.0, 1.0);
  73
  74   status = gsl_histogram_accumulate (h, 1.0, 10.0);
  75   gsl_test (status != GSL_EDOM, "gsl_histogram_accumulate traps x at xmax");
  76
  77   status = gsl_histogram_accumulate (h, 2.0, 100.0);
  78   gsl_test (status != GSL_EDOM, "gsl_histogram_accumulate traps x above xmax");
  79
  80   status = gsl_histogram_accumulate (h, -1.0, 1000.0);
  81   gsl_test (status != GSL_EDOM, "gsl_histogram_accumulate traps x below xmin");
  82
  83   status = gsl_histogram_increment (h, 1.0);
  84   gsl_test (status != GSL_EDOM, "gsl_histogram_increment traps x at xmax");
  85
  86   status = gsl_histogram_increment (h, 2.0);
  87   gsl_test (status != GSL_EDOM, "gsl_histogram_increment traps x above xmax");
  88
  89   status = gsl_histogram_increment (h, -1.0);
  90   gsl_test (status != GSL_EDOM, "gsl_histogram_increment traps x below xmin");
  91
  92
  93   result = gsl_histogram_get (h, N);
  94   gsl_test (result != 0, "gsl_histogram_get traps index at n");
  95
  96   result = gsl_histogram_get (h, N + 1);
  97   gsl_test (result != 0, "gsl_histogram_get traps index above n");
  98
  99   status = gsl_histogram_get_range (h, N, &lower, &upper);
 100   gsl_test (status != GSL_EDOM,
 101             "gsl_histogram_get_range traps index at n");
 102
 103   status = gsl_histogram_get_range (h, N + 1, &lower, &upper);
 104   gsl_test (status != GSL_EDOM,
 105             "gsl_histogram_get_range traps index above n");
 106
 107
 108   status = 0;
 109   gsl_histogram_find (h, -0.01, &i);
 110   gsl_test (status != GSL_EDOM, "gsl_histogram_find traps x below xmin");
 111
 112   status = 0;
 113   gsl_histogram_find (h, 1.0, &i);
 114   gsl_test (status != GSL_EDOM, "gsl_histogram_find traps x at xmax");
 115
 116   status = 0;
 117   gsl_histogram_find (h, 1.1, &i);
 118   gsl_test (status != GSL_EDOM, "gsl_histogram_find traps x above xmax");
 119
 120   gsl_histogram_free (h);
 121 }
 122
 123
 124 static void
 125 my_error_handler (const char *reason, const char *file, int line, int err)
 126 {
 127   if (0)
 128     printf ("(caught [%s:%d: %s (%d)])\n", file, line, reason, err);
 129   status = 1;
 130 }