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))