4 #ifndef THIRDPARTY_QUADPACK_H_ 
    5 #define THIRDPARTY_QUADPACK_H_ 
   35         double (*
m_f)(
const double);
 
   45         const double *, 
const double *, 
const double *, 
const double *, 
 
   46         const double *, 
const long *);
 
   48                                                                         const sys_float *, 
const sys_float *, 
 
   49                                                                         const sys_float *, 
const long *);
 
   51 void dqage_(
const D_fp& f, 
const double *a, 
const double *b, 
const double *epsabs, 
 
   52                                 const double *epsrel, 
const long *key, 
const long *limit, 
 
   53                                 double *result, 
double *abserr, 
long *neval, 
long *ier, 
double *
 
   54                                 alist__, 
double *blist, 
double *rlist, 
double *elist, 
 
   55                                 long *iord, 
long *last);
 
   56 void dqag_(
const D_fp& f, 
const double *a, 
const double *b, 
 
   57                           const double *epsabs, 
const double *epsrel, 
const long *key, 
 
   59                           double *abserr, 
long *neval, 
long *ier, 
long *limit, 
 
   60                           const long *lenw, 
long *last, 
long *iwork, 
double *work);
 
   61 void dqagie_(
const D_fp& f, 
const double *bound, 
const long *inf, 
 
   62                                  const double *epsabs, 
const double *epsrel, 
const long *limit, 
 
   63                                  double *result, 
double *abserr, 
long *neval, 
long *ier, 
double *
 
   64                                  alist__, 
double *blist, 
double *rlist, 
double *elist, 
 
   65                                  long *iord, 
long *last);
 
   66 void dqagi_(
const D_fp& f, 
const double *bound, 
const long *inf, 
 
   67                                 const double *epsabs, 
const double *epsrel, 
double *result, 
 
   68                                 double *abserr, 
long *neval, 
long *ier, 
long *limit, 
 
   69                                 const long *lenw, 
long *last, 
long *iwork, 
double *work);
 
   70 void dqagpe_(
const D_fp& f, 
const double *a, 
const double *b, 
 
   71                                  const long *npts2, 
const double *polongs, 
const double *epsabs, 
 
   73                                  const long *limit, 
double *result, 
double *abserr, 
long *
 
   74                                  neval, 
long *ier, 
double *alist__, 
double *blist, 
 
   75                                  double *rlist, 
double *elist, 
double *pts, 
long *iord, 
 
   76                                  long *level, 
long *ndin, 
long *last);
 
   77 void dqagp_(
const D_fp& f, 
const double *a, 
const double *b, 
const long *npts2, 
 
   78                                 const double *polongs, 
const double *epsabs, 
const double *epsrel, 
 
   79                                 double *result, 
double *abserr, 
long *neval, 
long *ier, 
 
   80                                 const long *leniw, 
const long *lenw, 
long *last, 
long *iwork, 
 
   82 void dqagse_(
const D_fp& f, 
const double *a, 
const double *b, 
const double *epsabs, 
 
   83                                  const double *epsrel, 
const long *limit, 
double *result, 
 
   84                                  double *abserr, 
long *neval, 
long *ier, 
double *alist__,
 
   85                                  double *blist, 
double *rlist, 
double *elist, 
long *
 
   87 void dqags_(
const D_fp& f, 
const double *a, 
const double *b, 
const double *epsabs, 
 
   88                                 const double *epsrel, 
double *result, 
double *abserr, 
 
   89                                 long *neval, 
long *ier, 
const long *limit, 
const long *lenw, 
long *last, 
 
   90                                 long *iwork, 
double *work);
 
   91 void dqawce_(
const D_fp& f, 
const double *a, 
const double *b, 
const double *c__, 
 
   92                                  const double *epsabs, 
const double *epsrel, 
const long *limit, 
 
   93                                  double *result, 
double *abserr, 
long *neval, 
long *ier, 
 
   94                                  double *alist__, 
double *blist, 
double *rlist, 
double  
   95                                  *elist, 
long *iord, 
long *last);
 
   96 void dqawc_(
const D_fp& f, 
const double *a, 
const double *b, 
const double *c__, 
 
   97                                 const double *epsabs, 
const double *epsrel, 
double *result, 
 
   98                                 double *abserr, 
long *neval, 
long *ier, 
long *limit, 
 
   99                                 const long *lenw, 
long *last, 
long *iwork, 
double *work);
 
  100 void dqawfe_(
const D_fp& f, 
const double *a, 
const double *omega, 
 
  101                                  const long *integr, 
const double *epsabs, 
const long *limlst, 
 
  102                                  const long *limit, 
const long *maxp1, 
 
  103                                  double *result, 
double *abserr, 
long *
 
  104                                  neval, 
long *ier, 
double *rslst, 
double *erlst, 
long *
 
  105                                  ierlst, 
long *lst, 
double *alist__, 
double *blist, 
 
  106                                  double *rlist, 
double *elist, 
long *iord, 
long *nnlog, 
 
  108 void dqawf_(
const D_fp& f, 
const double *a, 
const double *omega, 
const long *integr, 
 
  109                                 const double *epsabs, 
double *result, 
double *abserr, 
 
  110                                 long *neval, 
long *ier, 
long *limlst, 
long *lst, 
 
  111                                 const long *leniw, 
const long *maxp1, 
const long *lenw, 
 
  112                                 long *iwork, 
double *work);
 
  113 void dqawoe_(
const D_fp& f, 
const double *a, 
const double *b, 
const double *omega, 
 
  114                                  const long *integr, 
const double *epsabs, 
const double *epsrel, 
 
  115                                  const long *limit, 
const long *icall, 
const long *maxp1, 
 
  117                                  double *abserr, 
long *neval, 
long *ier, 
long *last, 
 
  118                                  double *alist__, 
double *blist, 
double *rlist, 
double  
  119                                  *elist, 
long *iord, 
long *nnlog, 
long *momcom, 
double *
 
  121 void dqawo_(
const D_fp& f, 
const double *a, 
const double *b, 
const double *omega, 
 
  122                                 const long *integr, 
const double *epsabs, 
const double *epsrel, 
 
  123                                 double *result, 
double *abserr, 
long *neval, 
long *ier, 
 
  124                                 const long *leniw, 
long *maxp1, 
const long *lenw, 
long *last, 
long  
  125                                 *iwork, 
double *work);
 
  126 void dqawse_(
const D_fp& f, 
const double *a, 
const double *b, 
const double *alfa, 
 
  127                                  const double *beta, 
const long *integr, 
const double *epsabs, 
 
  128                                  const double *epsrel, 
const long *limit, 
double *result, 
 
  129                                  double *abserr, 
long *neval, 
long *ier, 
double *alist__, 
 
  130                                  double *blist, 
double *rlist, 
double *elist, 
long *iord, 
 
  132 void dqaws_(
const D_fp& f, 
const double *a, 
const double *b, 
const double *alfa, 
 
  133                                 const double *beta, 
const long *integr, 
const double *epsabs, 
 
  134                                 const double *epsrel, 
double *result, 
double *abserr, 
 
  135                                 long *neval, 
long *ier, 
long *limit, 
const long *lenw, 
long *last, 
 
  136                                 long *iwork, 
double *work);
 
  137 void dqc25c_(
const D_fp& f, 
const double *a, 
const double *b, 
const double *c__, 
 
  138                                  double *result, 
double *abserr, 
long *krul, 
long *neval);
 
  139 void dqc25f_(
const D_fp& f, 
const double *a, 
const double *b, 
const double *omega, 
 
  140                                  const long *integr, 
const long *nrmom, 
const long *maxp1, 
 
  141                                  const long *ksave, 
double *result, 
double *abserr, 
long *neval, 
 
  142                                  double *resabs, 
double *resasc, 
long *momcom, 
double *
 
  144 void dqc25s_(
const D_fp& f, 
const double *a, 
const double *b, 
const double *bl, 
 
  145                                  const double *br, 
const double *alfa, 
const double *beta, 
 
  146                                  const double *ri, 
const double *rj, 
const double *rg, 
 
  148                                  double *result, 
double *abserr, 
double *resasc, 
 
  149                                  const long *integr, 
long *nev);
 
  150 void dqcheb_(
const double *x, 
double *fval, 
double *cheb12, 
double *cheb24);
 
  151 void dqelg_(
long *n, 
double *epstab, 
double *result, 
 
  152                                 double *abserr, 
double *res3la, 
long *nres);
 
  153 void dqk15_(
const D_fp& f, 
const double *a, 
const double *b, 
double *
 
  154                                 result, 
double *abserr, 
double *resabs, 
double *resasc);
 
  155 void dqk15i_(
const D_fp& f, 
const double *boun, 
const long *inf, 
 
  156                                  const double *a, 
const double *b, 
double *result, 
double *abserr, 
 
  157                                  double *resabs, 
double *resasc);
 
  159                                  const double *p3, 
const double *p4, 
const long *kp, 
 
  161                                  const double *b, 
double *result, 
double *abserr, 
double *
 
  162                                  resabs, 
double *resasc);
 
  163 void dqk21_(
const D_fp& f, 
const double *a, 
const double *b, 
double *
 
  164                                 result, 
double *abserr, 
double *resabs, 
double *resasc);
 
  165 void dqk31_(
const D_fp& f, 
const double *a, 
const double *b, 
double *
 
  166                                 result, 
double *abserr, 
double *resabs, 
double *resasc);
 
  167 void dqk41_(
const D_fp& f, 
const double *a, 
const double *b, 
double *
 
  168                                 result, 
double *abserr, 
double *resabs, 
double *resasc);
 
  169 void dqk51_(
const D_fp& f, 
const double *a, 
const double *b, 
double *
 
  170                                 result, 
double *abserr, 
double *resabs, 
double *resasc);
 
  171 void dqk61_(
const D_fp& f, 
const double *a, 
const double *b, 
double *
 
  172                                 result, 
double *abserr, 
double *resabs, 
double *resasc);
 
  173 void dqmomo_(
const double *alfa, 
const double *beta, 
double *
 
  174                                  ri, 
double *rj, 
double *rg, 
double *rh, 
const long *integr);
 
  175 void dqng_(
const D_fp& f, 
const double *a, 
const double *b, 
const double *epsabs, 
 
  176                           const double *epsrel, 
double *result, 
double *abserr, 
 
  177                           long *neval, 
long *ier);
 
  178 void dqpsrt_(
const long *limit, 
long *last, 
long *maxerr, 
 
  179                                  double *ermax, 
double *elist, 
long *iord, 
long *nrmax);
 
  180 void qage_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
 
  181                           const sys_float *epsabs, 
const sys_float *epsrel, 
 
  182                           const long *key, 
const long *limit, 
 
  183                           sys_float *result, sys_float *abserr, 
 
  184                           long *neval, 
long *ier, sys_float *alist__, sys_float *blist, sys_float *rlist,
 
  185                           sys_float *elist, 
long *iord, 
long *last);
 
  186 void qag_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
 
  187                          const sys_float *epsabs, 
const sys_float *epsrel, 
const long *key, 
 
  188                          sys_float *result, sys_float *abserr, 
long *neval, 
 
  189                          long *ier, 
long *limit, 
const long *lenw, 
long *last, 
long *
 
  190                          iwork, sys_float *work);
 
  191 void qagie_(
const E_fp& f, 
const sys_float *bound, 
const long *inf, 
 
  192                                 const sys_float *epsabs, 
const sys_float *epsrel, 
 
  193                                 const long *limit, sys_float *result, sys_float *abserr, 
 
  194                                 long *neval, 
long *ier, sys_float *alist__, sys_float *blist, 
 
  195                                 sys_float *rlist, sys_float *elist, 
long *iord, 
long *last);
 
  196 void qagi_(
const E_fp& f, 
const sys_float *bound, 
const long *inf, 
 
  197                           const sys_float *epsabs, 
const sys_float *epsrel, 
 
  198                           sys_float *result, sys_float *abserr, 
 
  199                           long *neval, 
long *ier, 
const long *limit, 
const long *lenw, 
 
  200                           long *last, 
long *iwork, sys_float *work);
 
  201 void qagpe_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
const long *npts2, 
 
  202                                 const sys_float *polongs, 
const sys_float *epsabs, 
 
  203                                 const sys_float *epsrel, 
const long *limit, sys_float *result, 
 
  204                                 sys_float *abserr, 
long *neval, 
long *ier, sys_float *alist__, 
 
  205                                 sys_float *blist, sys_float *rlist, sys_float *elist, 
 
  206                                 sys_float *pts, 
long *iord, 
long *level, 
long *ndin, 
long *last);
 
  207 void qagp_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
const long *npts2,
 
  208                           const sys_float *polongs, 
const sys_float *epsabs, 
 
  209                           const sys_float *epsrel, sys_float *result, sys_float *abserr, 
 
  210                           long *neval, 
long *ier, 
const long *leniw, 
const long *lenw, 
 
  211                           long *last, 
long *iwork, sys_float *work);
 
  212 void qagse_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
 
  213                                 const sys_float *epsabs, 
const sys_float *epsrel, 
 
  214                                 const long *limit, sys_float *result, sys_float *abserr, 
 
  216                                 long *ier, sys_float *alist__, sys_float *blist, sys_float *rlist, sys_float *elist, 
 
  217                                 long *iord, 
long *last);
 
  218 void qags_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
 
  219                           const sys_float *epsabs, 
const sys_float *epsrel, 
 
  220                           sys_float *result, sys_float *abserr, 
long *neval, 
long *ier, 
 
  221                           const long *limit, 
const long *lenw, 
long *last, 
long *iwork,
 
  223 void qawce_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
 
  224                                 const sys_float *c__, 
const sys_float *epsabs,
 
  225                                 const sys_float *epsrel, 
const long *limit, 
 
  226                                 sys_float *result, sys_float *abserr, 
long *
 
  227                                 neval, 
long *ier, sys_float *alist__, sys_float *blist, 
 
  228                                 sys_float *rlist, sys_float *elist, 
long *iord, 
long *last);
 
  229 void qawc_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
 
  230                           const sys_float *c__, 
const sys_float *epsabs, 
 
  231                           const sys_float *epsrel, sys_float *result, sys_float *abserr, 
 
  232                           long *neval, 
long *ier, 
long *limit, 
const long *lenw, 
 
  233                           long *last, 
long *iwork, 
 
  235 void qawfe_(
const E_fp& f, 
const sys_float *a, 
const sys_float *omega, 
 
  236                                 const long *integr, 
const sys_float *epsabs, 
const long *limlst, 
 
  237                                 const long *limit, 
const long *maxp1, 
 
  238                                 sys_float *result, sys_float *abserr, 
long *neval, 
 
  239                                 long *ier, sys_float *rslst, sys_float *erlst, 
 
  240                                 long *ierlst, 
long *lst, sys_float *alist__, sys_float *blist, 
 
  241                                 sys_float *rlist, sys_float *elist, 
long *iord, 
 
  242                                 long *nnlog, sys_float *chebmo);
 
  243 void qawf_(
const E_fp& f, 
const sys_float *a, 
const sys_float *omega, 
 
  245                           const sys_float *epsabs, sys_float *result, 
 
  246                           sys_float *abserr, 
long *neval, 
long *ier, 
 
  247                           const long *limlst, 
long *lst, 
const long *leniw, 
const long *maxp1, 
 
  248                           const long *lenw, 
long *iwork, sys_float *work);
 
  249 void qawoe_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
 
  250                                 const sys_float *omega, 
const long *integr, 
 
  251                                 const sys_float *epsabs, 
const sys_float *epsrel, 
 
  252                                 const long *limit, 
const long *icall, 
const long *maxp1, 
 
  253                                 sys_float *result, sys_float *abserr, 
long *neval, 
long *
 
  254                                 ier, 
long *last, sys_float *alist__, sys_float *blist, sys_float *rlist, sys_float *
 
  255                                 elist, 
long *iord, 
long *nnlog, 
long *momcom, sys_float *chebmo);
 
  256 void qawo_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
 
  257                           const sys_float *omega, 
const long *integr, 
 
  258                           const sys_float *epsabs, 
const sys_float *epsrel, 
 
  259                           sys_float *result, sys_float *abserr, 
 
  260                           long *neval, 
long *ier, 
const long *leniw, 
const long *maxp1, 
 
  261                           const long *lenw, 
long *last, 
long *iwork, sys_float *work);
 
  262 void qawse_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
 
  263                                 const sys_float *alfa, 
const sys_float *beta, 
 
  264                                 const long *integr, 
const sys_float *epsabs, 
 
  265                                 const sys_float *epsrel, 
const long *limit, 
 
  266                                 sys_float *result, sys_float *abserr, 
long *neval, 
 
  267                                 long *ier, sys_float *alist__, 
 
  268                                 sys_float *blist, sys_float *rlist, sys_float *elist, 
long *iord, 
long *last);
 
  269 void qaws_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
 
  270                           const sys_float *alfa, 
const sys_float *beta, 
 
  271                           const long *integr, 
const sys_float *epsabs, 
 
  272                           const sys_float *epsrel, sys_float *result, sys_float *abserr, 
 
  273                           long *neval, 
long *ier, 
const long *limit, 
const long *lenw, 
 
  274                           long *last, 
long *iwork, sys_float *work);
 
  275 void qc25c_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
 
  276                                 const sys_float *c__, sys_float *result,
 
  277                                 sys_float *abserr, 
long *krul, 
long *neval);
 
  278 void qc25f_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
 
  279                                 const sys_float *omega, 
const long *integr, 
const long *nrmom, 
 
  280                                 const long *maxp1, 
const long *ksave, sys_float *result, 
 
  281                                 sys_float *abserr, 
long *neval, sys_float *resabs, sys_float *resasc, 
long *
 
  282                                 momcom, sys_float *chebmo);
 
  283 void qc25s_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, 
 
  284                                 const sys_float *bl, 
const sys_float *br, 
 
  285                                 const sys_float *alfa, 
const sys_float *beta, 
 
  286                                 sys_float *ri, sys_float *rj, sys_float *rg, sys_float *rh, 
 
  287                                 sys_float *result, sys_float *abserr, sys_float *resasc, 
 
  288                                 const long *integr, 
long *nev);
 
  289 void qcheb_(sys_float *x, sys_float *fval, sys_float *cheb12, sys_float *cheb24);
 
  290 void qelg_(
long *n, sys_float *epstab, sys_float *result, sys_float *
 
  291                           abserr, sys_float *res3la, 
long *nres);
 
  292 void qk15_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, sys_float *result, 
 
  293                           sys_float *abserr, sys_float *resabs, sys_float *resasc);
 
  294 void qk15i_(
const E_fp& f, 
const sys_float *boun, 
const long *inf, 
const sys_float *a,
 
  295                                 const sys_float *b, sys_float *result, sys_float *abserr, 
 
  296                                 sys_float *resabs, sys_float *resasc);
 
  297 void qk15w_(
const E_fp& f, 
E_fp1 w, 
const sys_float *p1, 
const sys_float *p2, 
 
  298                                 const sys_float *p3, 
const sys_float *p4, 
const long *kp, 
 
  299                                 const sys_float *a, 
const sys_float *b, sys_float *result, 
 
  300                                 sys_float *abserr, sys_float *resabs, sys_float *resasc);
 
  301 void qk21_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, sys_float *result, 
 
  302                           sys_float *abserr, sys_float *resabs, sys_float *resasc);
 
  303 void qk31_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, sys_float *result, 
 
  304                           sys_float *abserr, sys_float *resabs, sys_float *resasc);
 
  305 void qk41_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, sys_float *result, 
 
  306                           sys_float *abserr, sys_float *resabs, sys_float *resasc);
 
  307 void qk51_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, sys_float *result, 
 
  308                           sys_float *abserr, sys_float *resabs, sys_float *resasc);
 
  309 void qk61_(
const E_fp& f, 
const sys_float *a, 
const sys_float *b, sys_float *result, 
 
  310                           sys_float *abserr, sys_float *resabs, sys_float *resasc);
 
  311 void qmomo_(
const sys_float *alfa, 
const sys_float *beta, sys_float *ri, 
 
  312                                 sys_float *rj, sys_float *rg, sys_float *rh, 
const long *integr);
 
  313 void qng_(
const E_fp& f, sys_float *a, sys_float *b, sys_float *epsabs, sys_float *
 
  314                          epsrel, sys_float *result, sys_float *abserr, 
long *neval, 
long *ier);
 
  315 void qpsrt_(
const long *limit, 
long *last, 
long *maxerr, 
 
  316                                 sys_float *ermax, sys_float *elist, 
long *iord, 
long *nrmax);
 
  317 double dqwgtc_(
const double *x, 
const double *c__, 
const double *p2, 
 
  318                                         const double *p3, 
const double *p4, 
const long *kp);
 
  319 double dqwgtf_(
const double *x, 
const double *omega, 
const double *p2, 
 
  320                                         const double *p3, 
const double *p4, 
const long *integr);
 
  321 double dqwgts_(
const double *x, 
const double *a, 
const double *b, 
 
  322                                         const double *alfa, 
const double *beta, 
const long *integr);
 
  323 sys_float 
qwgtc_(
const sys_float *x, 
const sys_float *c__, 
const sys_float *p2,
 
  324                                           const sys_float *p3, 
const sys_float *p4, 
const long *kp);
 
  325 sys_float 
qwgtf_(
const sys_float *x, 
const sys_float *omega, 
 
  326                                           const sys_float *p2, 
const sys_float *p3, 
const sys_float *p4,
 
  328 sys_float 
qwgts_(
const sys_float *x, 
const sys_float *a, 
const sys_float *b, 
 
  329                                           const sys_float *alfa, 
const sys_float *beta, 
 
void qk15i_(const E_fp &f, const sys_float *boun, const long *inf, const sys_float *a, const sys_float *b, sys_float *result, sys_float *abserr, sys_float *resabs, sys_float *resasc)
sys_float(* E_fp1)(const sys_float *, const sys_float *, const sys_float *, const sys_float *, const sys_float *, const long *)
sys_float(* m_f)(sys_float)
void qc25f_(const E_fp &f, const sys_float *a, const sys_float *b, const sys_float *omega, const long *integr, const long *nrmom, const long *maxp1, const long *ksave, sys_float *result, sys_float *abserr, long *neval, sys_float *resabs, sys_float *resasc, long *momcom, sys_float *chebmo)
double dqwgtc_(const double *x, const double *c__, const double *, const double *, const double *, const long *)
void dqk61_(const D_fp &f, const double *a, const double *b, double *result, double *abserr, double *resabs, double *resasc)
virtual double operator()(double) const =0
void dqk15w_(const D_fp &f, D_fp1 w, const double *p1, const double *p2, const double *p3, const double *p4, const long *kp, const double *a, const double *b, double *result, double *abserr, double *resabs, double *resasc)
void qk51_(const E_fp &f, const sys_float *a, const sys_float *b, sys_float *result, sys_float *abserr, sys_float *resabs, sys_float *resasc)
void qelg_(long *n, sys_float *epstab, sys_float *result, sys_float *abserr, sys_float *res3la, long *nres)
void qk15w_(const E_fp &f, E_fp1 w, const sys_float *p1, const sys_float *p2, const sys_float *p3, const sys_float *p4, const long *kp, const sys_float *a, const sys_float *b, sys_float *result, sys_float *abserr, sys_float *resabs, sys_float *resasc)
void qagse_(const E_fp &f, const sys_float *a, const sys_float *b, const sys_float *epsabs, const sys_float *epsrel, const long *limit, sys_float *result, sys_float *abserr, long *neval, long *ier, sys_float *alist__, sys_float *blist, sys_float *rlist, sys_float *elist, long *iord, long *last)
sys_float operator()(sys_float x) const 
void qk21_(const E_fp &f, const sys_float *a, const sys_float *b, sys_float *result, sys_float *abserr, sys_float *resabs, sys_float *resasc)
void dqawoe_(const D_fp &f, const double *a, const double *b, const double *omega, const long *integr, const double *epsabs, const double *epsrel, const long *limit, const long *icall, const long *maxp1, double *result, double *abserr, long *neval, long *ier, long *last, double *alist__, double *blist, double *rlist, double *elist, long *iord, long *nnlog, long *momcom, double *chebmo)
double(* D_fp1)(const double *, const double *, const double *, const double *, const double *, const long *)
void dqagse_(const D_fp &f, const double *a, const double *b, const double *epsabs, const double *epsrel, const long *limit, double *result, double *abserr, long *neval, long *ier, double *alist__, double *blist, double *rlist, double *elist, long *iord, long *last)
void dqag_(const D_fp &f, const double *a, const double *b, const double *epsabs, const double *epsrel, const long *key, double *result, double *abserr, long *neval, long *ier, long *limit, const long *lenw, long *last, long *iwork, double *work)
void dqagie_(const D_fp &f, const double *bound, const long *inf, const double *epsabs, const double *epsrel, const long *limit, double *result, double *abserr, long *neval, long *ier, double *alist__, double *blist, double *rlist, double *elist, long *iord, long *last)
void qage_(const E_fp &f, const sys_float *a, const sys_float *b, const sys_float *epsabs, const sys_float *epsrel, const long *key, const long *limit, sys_float *result, sys_float *abserr, long *neval, long *ier, sys_float *alist__, sys_float *blist, sys_float *rlist, sys_float *elist, long *iord, long *last)
void dqk15_(const D_fp &f, const double *a, const double *b, double *result, double *abserr, double *resabs, double *resasc)
void dqmomo_(const double *alfa, const double *beta, double *ri, double *rj, double *rg, double *rh, const long *integr)
void qk41_(const E_fp &f, const sys_float *a, const sys_float *b, sys_float *result, sys_float *abserr, sys_float *resabs, sys_float *resasc)
double dqwgtf_(const double *x, const double *omega, const double *, const double *, const double *, const long *integr)
void qawse_(const E_fp &f, const sys_float *a, const sys_float *b, const sys_float *alfa, const sys_float *beta, const long *integr, const sys_float *epsabs, const sys_float *epsrel, const long *limit, sys_float *result, sys_float *abserr, long *neval, long *ier, sys_float *alist__, sys_float *blist, sys_float *rlist, sys_float *elist, long *iord, long *last)
void qc25c_(const E_fp &f, const sys_float *a, const sys_float *b, const sys_float *c__, sys_float *result, sys_float *abserr, long *krul, long *neval)
void dqagp_(const D_fp &f, const double *a, const double *b, const long *npts2, const double *points, const double *epsabs, const double *epsrel, double *result, double *abserr, long *neval, long *ier, const long *leniw, const long *lenw, long *last, long *iwork, double *work)
sys_float qwgtf_(const sys_float *x, const sys_float *omega, const sys_float *, const sys_float *, const sys_float *, const long *integr)
void dqawc_(const D_fp &f, const double *a, const double *b, const double *c__, const double *epsabs, const double *epsrel, double *result, double *abserr, long *neval, long *ier, long *limit, const long *lenw, long *last, long *iwork, double *work)
void qagie_(const E_fp &f, const sys_float *bound, const long *inf, const sys_float *epsabs, const sys_float *epsrel, const long *limit, sys_float *result, sys_float *abserr, long *neval, long *ier, sys_float *alist__, sys_float *blist, sys_float *rlist, sys_float *elist, long *iord, long *last)
D_fp_fp(double(*f)(double))
void qk15_(const E_fp &f, const sys_float *a, const sys_float *b, sys_float *result, sys_float *abserr, sys_float *resabs, sys_float *resasc)
void dqage_(const D_fp &f, const double *a, const double *b, const double *epsabs, const double *epsrel, const long *key, const long *limit, double *result, double *abserr, long *neval, long *ier, double *alist__, double *blist, double *rlist, double *elist, long *iord, long *last)
sys_float qwgtc_(const sys_float *x, const sys_float *c__, const sys_float *, const sys_float *, const sys_float *, const long *)
void qaws_(const E_fp &f, const sys_float *a, const sys_float *b, const sys_float *alfa, const sys_float *beta, const long *integr, const sys_float *epsabs, const sys_float *epsrel, sys_float *result, sys_float *abserr, long *neval, long *ier, const long *limit, const long *lenw, long *last, long *iwork, sys_float *work)
double(* m_f)(const double)
void dqawo_(const D_fp &f, const double *a, const double *b, const double *omega, const long *integr, const double *epsabs, const double *epsrel, double *result, double *abserr, long *neval, long *ier, const long *leniw, long *maxp1, const long *lenw, long *last, long *iwork, double *work)
void dqawfe_(const D_fp &f, const double *a, const double *omega, const long *integr, const double *epsabs, const long *limlst, const long *limit, const long *maxp1, double *result, double *abserr, long *neval, long *ier, double *rslst, double *erlst, long *ierlst, long *lst, double *alist__, double *blist, double *rlist, double *elist, long *iord, long *nnlog, double *chebmo)
void qk61_(const E_fp &f, const sys_float *a, const sys_float *b, sys_float *result, sys_float *abserr, sys_float *resabs, sys_float *resasc)
void dqc25s_(const D_fp &f, const double *a, const double *b, const double *bl, const double *br, const double *alfa, const double *beta, const double *ri, const double *rj, const double *rg, const double *rh, double *result, double *abserr, double *resasc, const long *integr, long *nev)
void dqags_(const D_fp &f, const double *a, const double *b, const double *epsabs, const double *epsrel, double *result, double *abserr, long *neval, long *ier, const long *limit, const long *lenw, long *last, long *iwork, double *work)
void dqawf_(const D_fp &f, const double *a, const double *omega, const long *integr, const double *epsabs, double *result, double *abserr, long *neval, long *ier, long *limlst, long *lst, const long *leniw, const long *maxp1, const long *lenw, long *iwork, double *work)
void qc25s_(const E_fp &f, const sys_float *a, const sys_float *b, const sys_float *bl, const sys_float *br, const sys_float *alfa, const sys_float *beta, sys_float *ri, sys_float *rj, sys_float *rg, sys_float *rh, sys_float *result, sys_float *abserr, sys_float *resasc, const long *integr, long *nev)
void qawo_(const E_fp &f, const sys_float *a, const sys_float *b, const sys_float *omega, const long *integr, const sys_float *epsabs, const sys_float *epsrel, sys_float *result, sys_float *abserr, long *neval, long *ier, const long *leniw, const long *maxp1, const long *lenw, long *last, long *iwork, sys_float *work)
sys_float qwgts_(const sys_float *x, const sys_float *a, const sys_float *b, const sys_float *alfa, const sys_float *beta, const long *integr)
void dqaws_(const D_fp &f, const double *a, const double *b, const double *alfa, const double *beta, const long *integr, const double *epsabs, const double *epsrel, double *result, double *abserr, long *neval, long *ier, long *limit, const long *lenw, long *last, long *iwork, double *work)
void dqawce_(const D_fp &f, const double *a, const double *b, const double *c__, const double *epsabs, const double *epsrel, const long *limit, double *result, double *abserr, long *neval, long *ier, double *alist__, double *blist, double *rlist, double *elist, long *iord, long *last)
void qcheb_(sys_float *x, sys_float *fval, sys_float *cheb12, sys_float *cheb24)
void qawoe_(const E_fp &f, const sys_float *a, const sys_float *b, const sys_float *omega, const long *integr, const sys_float *epsabs, const sys_float *epsrel, const long *limit, const long *icall, const long *maxp1, sys_float *result, sys_float *abserr, long *neval, long *ier, long *last, sys_float *alist__, sys_float *blist, sys_float *rlist, sys_float *elist, long *iord, long *nnlog, long *momcom, sys_float *chebmo)
void dqng_(const D_fp &f, const double *a, const double *b, const double *epsabs, const double *epsrel, double *result, double *abserr, long *neval, long *ier)
void qawc_(const E_fp &f, const sys_float *a, const sys_float *b, const sys_float *c__, const sys_float *epsabs, const sys_float *epsrel, sys_float *result, sys_float *abserr, long *neval, long *ier, long *limit, const long *lenw, long *last, long *iwork, sys_float *work)
void qmomo_(const sys_float *alfa, const sys_float *beta, sys_float *ri, sys_float *rj, sys_float *rg, sys_float *rh, const long *integr)
void qagpe_(const E_fp &f, const sys_float *a, const sys_float *b, const long *npts2, const sys_float *points, const sys_float *epsabs, const sys_float *epsrel, const long *limit, sys_float *result, sys_float *abserr, long *neval, long *ier, sys_float *alist__, sys_float *blist, sys_float *rlist, sys_float *elist, sys_float *pts, long *iord, long *level, long *ndin, long *last)
void dqk51_(const D_fp &f, const double *a, const double *b, double *result, double *abserr, double *resabs, double *resasc)
void dqpsrt_(const long *limit, long *last, long *maxerr, double *ermax, double *elist, long *iord, long *nrmax)
void dqc25f_(const D_fp &f, const double *a, const double *b, const double *omega, const long *integr, const long *nrmom, const long *maxp1, const long *ksave, double *result, double *abserr, long *neval, double *resabs, double *resasc, long *momcom, double *chebmo)
virtual sys_float operator()(sys_float) const =0
void dqk15i_(const D_fp &f, const double *boun, const long *inf, const double *a, const double *b, double *result, double *abserr, double *resabs, double *resasc)
void dqelg_(long *n, double *epstab, double *result, double *abserr, double *res3la, long *nres)
double dqwgts_(const double *x, const double *a, const double *b, const double *alfa, const double *beta, const long *integr)
void qpsrt_(const long *limit, long *last, long *maxerr, sys_float *ermax, sys_float *elist, long *iord, long *nrmax)
void qags_(const E_fp &f, const sys_float *a, const sys_float *b, const sys_float *epsabs, const sys_float *epsrel, sys_float *result, sys_float *abserr, long *neval, long *ier, const long *limit, const long *lenw, long *last, long *iwork, sys_float *work)
double operator()(double x) const 
void qagp_(const E_fp &f, const sys_float *a, const sys_float *b, const long *npts2, const sys_float *points, const sys_float *epsabs, const sys_float *epsrel, sys_float *result, sys_float *abserr, long *neval, long *ier, const long *leniw, const long *lenw, long *last, long *iwork, sys_float *work)
void dqk21_(const D_fp &f, const double *a, const double *b, double *result, double *abserr, double *resabs, double *resasc)
void qawf_(const E_fp &f, const sys_float *a, const sys_float *omega, const long *integr, const sys_float *epsabs, sys_float *result, sys_float *abserr, long *neval, long *ier, const long *limlst, long *lst, const long *leniw, const long *maxp1, const long *lenw, long *iwork, sys_float *work)
void dqagi_(const D_fp &f, const double *bound, const long *inf, const double *epsabs, const double *epsrel, double *result, double *abserr, long *neval, long *ier, long *limit, const long *lenw, long *last, long *iwork, double *work)
void qag_(const E_fp &f, const sys_float *a, const sys_float *b, const sys_float *epsabs, const sys_float *epsrel, const long *key, sys_float *result, sys_float *abserr, long *neval, long *ier, long *limit, const long *lenw, long *last, long *iwork, sys_float *work)
void dqk41_(const D_fp &f, const double *a, const double *b, double *result, double *abserr, double *resabs, double *resasc)
void dqk31_(const D_fp &f, const double *a, const double *b, double *result, double *abserr, double *resabs, double *resasc)
void dqcheb_(const double *x, double *fval, double *cheb12, double *cheb24)
void qagi_(const E_fp &f, const sys_float *bound, const long *inf, const sys_float *epsabs, const sys_float *epsrel, sys_float *result, sys_float *abserr, long *neval, long *ier, const long *limit, const long *lenw, long *last, long *iwork, sys_float *work)
void qawfe_(const E_fp &f, const sys_float *a, const sys_float *omega, const long *integr, const sys_float *epsabs, const long *limlst, const long *limit, const long *maxp1, sys_float *result, sys_float *abserr, long *neval, long *ier, sys_float *rslst, sys_float *erlst, long *ierlst, long *lst, sys_float *alist__, sys_float *blist, sys_float *rlist, sys_float *elist, long *iord, long *nnlog, sys_float *chebmo)
void qawce_(const E_fp &f, const sys_float *a, const sys_float *b, const sys_float *c__, const sys_float *epsabs, const sys_float *epsrel, const long *limit, sys_float *result, sys_float *abserr, long *neval, long *ier, sys_float *alist__, sys_float *blist, sys_float *rlist, sys_float *elist, long *iord, long *last)
void dqc25c_(const D_fp &f, const double *a, const double *b, const double *c__, double *result, double *abserr, long *krul, long *neval)
void qng_(const E_fp &f, sys_float *a, sys_float *b, sys_float *epsabs, sys_float *epsrel, sys_float *result, sys_float *abserr, long *neval, long *ier)
void qk31_(const E_fp &f, const sys_float *a, const sys_float *b, sys_float *result, sys_float *abserr, sys_float *resabs, sys_float *resasc)
void dqagpe_(const D_fp &f, const double *a, const double *b, const long *npts2, const double *points, const double *epsabs, const double *epsrel, const long *limit, double *result, double *abserr, long *neval, long *ier, double *alist__, double *blist, double *rlist, double *elist, double *pts, long *iord, long *level, long *ndin, long *last)
void dqawse_(const D_fp &f, const double *a, const double *b, const double *alfa, const double *beta, const long *integr, const double *epsabs, const double *epsrel, const long *limit, double *result, double *abserr, long *neval, long *ier, double *alist__, double *blist, double *rlist, double *elist, long *iord, long *last)
E_fp_fp(sys_float(*f)(sys_float))