integrate()

static int integrate ( unsigned  fdim, integrand_v  f, void *  fdata, unsigned  dim, const double *  xmin, const double *  xmax, unsigned  maxEval, double  reqAbsError, double  reqRelError, double *  val, double *  err, int  parallel  )

static

{
     rule *r;
     hypercube h;
     int status;
     unsigned i;
     
     if (fdim == 0) /* nothing to do */ return SUCCESS;
     if (dim == 0) { /* trivial integration */
          f(0, 1, xmin, fdata, fdim, val);
          for (i = 0; i < fdim; ++i) err[i] = 0;
          return SUCCESS;
     }
     r = dim == 1 ? make_rule15gauss(dim, fdim)
                  : make_rule75genzmalik(dim, fdim);
     if (!r) { 
          for (i = 0; i < fdim; ++i) {
               val[i] = 0;
               err[i] = HUGE_VAL; 
          }
          return FAILURE;
     }
     h = make_hypercube_range(dim, xmin, xmax);
     status = !h.data ? FAILURE
          : ruleadapt_integrate(r, fdim, f, fdata, &h,
                                maxEval, reqAbsError, reqRelError,
                                val, err, parallel);
     destroy_hypercube(&h);
     destroy_rule(r);
     return status;
}