void
test_quantile(const double p, const double data[], const size_t n,
const double expected, const double tol, const char *desc)
{
gsl_rstat_quantile_workspace *w = gsl_rstat_quantile_alloc(p);
double result;
size_t i;
for (i = 0; i < n; ++i)
gsl_rstat_quantile_add(data[i], w);
result = gsl_rstat_quantile_get(w);
if (fabs(expected) < 1.0e-4)
gsl_test_abs(result, expected, tol, "%s p=%g", desc, p);
else
gsl_test_rel(result, expected, tol, "%s p=%g", desc, p);
gsl_rstat_quantile_free(w);
}
gsl_rstat_workspace *
gsl_rstat_alloc(void)
{
gsl_rstat_workspace *w;
w = calloc(1, sizeof(gsl_rstat_workspace));
if (w == 0)
{
GSL_ERROR_NULL ("failed to allocate space for workspace", GSL_ENOMEM);
}
w->median_workspace_p = gsl_rstat_quantile_alloc(0.5);
if (w->median_workspace_p == 0)
{
GSL_ERROR_NULL ("failed to allocate space for median workspace",
GSL_ENOMEM);
}
gsl_rstat_reset(w);
return w;
} /* gsl_rstat_alloc() */
