gsl_histogram *
gsl_histogram_calloc_uniform (const size_t n, const double xmin,
const double xmax)
{
gsl_histogram *h;
if (xmin >= xmax)
{
GSL_ERROR_VAL ("xmin must be less than xmax", GSL_EINVAL, 0);
}
h = gsl_histogram_calloc (n);
if (h == 0)
{
return h;
}
make_uniform (h->range, n, xmin, xmax);
return h;
}
void
test1d_trap (void)
{
gsl_histogram *h;
double result, lower, upper;
size_t i;
gsl_set_error_handler (&my_error_handler);
gsl_ieee_env_setup ();
status = 0;
h = gsl_histogram_calloc (0);
gsl_test (!status, "gsl_histogram_calloc traps zero-length histogram");
gsl_test (h != 0,
"gsl_histogram_calloc returns NULL for zero-length histogram");
status = 0;
h = gsl_histogram_calloc_uniform (0, 0.0, 1.0);
gsl_test (!status,
"gsl_histogram_calloc_uniform traps zero-length histogram");
gsl_test (h != 0,
"gsl_histogram_calloc_uniform returns NULL for zero-length histogram");
status = 0;
h = gsl_histogram_calloc_uniform (10, 1.0, 1.0);
gsl_test (!status,
"gsl_histogram_calloc_uniform traps equal endpoints");
gsl_test (h != 0,
"gsl_histogram_calloc_uniform returns NULL for equal endpoints");
status = 0;
h = gsl_histogram_calloc_uniform (10, 2.0, 1.0);
gsl_test (!status,
"gsl_histogram_calloc_uniform traps invalid range");
gsl_test (h != 0,
"gsl_histogram_calloc_uniform returns NULL for invalid range");
h = gsl_histogram_calloc_uniform (N, 0.0, 1.0);
status = gsl_histogram_accumulate (h, 1.0, 10.0);
gsl_test (status != GSL_EDOM, "gsl_histogram_accumulate traps x at xmax");
status = gsl_histogram_accumulate (h, 2.0, 100.0);
gsl_test (status != GSL_EDOM, "gsl_histogram_accumulate traps x above xmax");
status = gsl_histogram_accumulate (h, -1.0, 1000.0);
gsl_test (status != GSL_EDOM, "gsl_histogram_accumulate traps x below xmin");
status = gsl_histogram_increment (h, 1.0);
gsl_test (status != GSL_EDOM, "gsl_histogram_increment traps x at xmax");
status = gsl_histogram_increment (h, 2.0);
gsl_test (status != GSL_EDOM, "gsl_histogram_increment traps x above xmax");
status = gsl_histogram_increment (h, -1.0);
gsl_test (status != GSL_EDOM, "gsl_histogram_increment traps x below xmin");
result = gsl_histogram_get (h, N);
gsl_test (result != 0, "gsl_histogram_get traps index at n");
result = gsl_histogram_get (h, N + 1);
gsl_test (result != 0, "gsl_histogram_get traps index above n");
status = gsl_histogram_get_range (h, N, &lower, &upper);
gsl_test (status != GSL_EDOM,
"gsl_histogram_get_range traps index at n");
status = gsl_histogram_get_range (h, N + 1, &lower, &upper);
gsl_test (status != GSL_EDOM,
"gsl_histogram_get_range traps index above n");
status = 0;
gsl_histogram_find (h, -0.01, &i);
gsl_test (status != GSL_EDOM, "gsl_histogram_find traps x below xmin");
status = 0;
gsl_histogram_find (h, 1.0, &i);
gsl_test (status != GSL_EDOM, "gsl_histogram_find traps x at xmax");
status = 0;
gsl_histogram_find (h, 1.1, &i);
gsl_test (status != GSL_EDOM, "gsl_histogram_find traps x above xmax");
gsl_histogram_free (h);
}
