00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023 #ifndef O2SCL_GSL_INTE_QAGI_H
00024 #define O2SCL_GSL_INTE_QAGI_H
00025
00026 #include <o2scl/inte.h>
00027 #include <o2scl/gsl_inte_qags.h>
00028
00029 #ifndef DOXYGENP
00030 namespace o2scl {
00031 #endif
00032
00033
00034
00035
00036 template<class param_t, class func_t> class gsl_inte_qagi :
00037 public gsl_inte_transform<param_t,func_t> {
00038
00039 public:
00040
00041
00042
00043
00044
00045
00046 virtual double integ(func_t &func, double a, double b, param_t &pa) {
00047 double res, err;
00048 integ_err(func,a,b,pa,res,err);
00049 this->interror=err;
00050 return res;
00051 }
00052
00053
00054
00055
00056
00057
00058 virtual int integ_err(func_t &func, double a, double b,
00059 param_t &pa, double &res, double &err2) {
00060
00061 double fv1[8], fv2[8];
00062
00063 int status=qags(func,8,o2scl_inte_qag_coeffs::qk15_xgk,
00064 o2scl_inte_qag_coeffs::qk15_wg,
00065 o2scl_inte_qag_coeffs::qk15_wgk,
00066 fv1,fv2,0.0,1.0,this->tolx,this->tolf,
00067 this->wkspace,&res,&err2,pa);
00068 return status;
00069
00070 }
00071
00072 #ifndef DOXYGEN_INTERNAL
00073
00074 protected:
00075
00076
00077 virtual double transform(func_t &func, double t, param_t &pa) {
00078 double x=(1-t)/t, y1, y2;
00079 func(x,y1,pa);
00080 func(-x,y2,pa);
00081 return (y1+y2)/t/t;
00082 }
00083
00084 #endif
00085
00086 };
00087
00088 #ifndef DOXYGENP
00089 }
00090 #endif
00091
00092 #endif