dbl_reduce
double dbl_reduce(double init, const double *s, const double *e, double (*merge)(double, double));
double dbl_reduce(double, const double *, const double *, double (*)(double, double));