Re: [boost] Numerical Quadrature Library

From: Neal Becker (ndbecker2_at_[hidden])
Date: 2006-12-07 10:22:52

Starting with basics.  Here is my idea for a family of 1-d integrators:

template<typename func_t>
struct arg {
  typedef typename func_t::arg_t type;
};

template<typename func_t>
struct ret {
  typedef typename func_t::ret_t type;
};

template<typename func_t, typename arg_t=typename arg<func_t>::type, typename ret_t=typename ret<func_t>::type>
struct trap_1d {
  ret_t operator() (arg_t lower, arg_t upper);
};

I thought about having the integrator type (romberg, trap, etc) be specified as the
1st template parameter, but offhand I don't see what value that adds.  Any thoughts?