00001 /* 00002 ------------------------------------------------------------------- 00003 00004 Copyright (C) 2006, 2007, 2008, 2009, Andrew W. Steiner 00005 00006 This file is part of O2scl. 00007 00008 O2scl is free software; you can redistribute it and/or modify 00009 it under the terms of the GNU General Public License as published by 00010 the Free Software Foundation; either version 3 of the License, or 00011 (at your option) any later version. 00012 00013 O2scl is distributed in the hope that it will be useful, 00014 but WITHOUT ANY WARRANTY; without even the implied warranty of 00015 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00016 GNU General Public License for more details. 00017 00018 You should have received a copy of the GNU General Public License 00019 along with O2scl. If not, see <http://www.gnu.org/licenses/>. 00020 00021 ------------------------------------------------------------------- 00022 */ 00023 #ifndef O2SCL_GSL_INTE_QAGIL_H 00024 #define O2SCL_GSL_INTE_QAGIL_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 /** \brief Integrate a function from \f$ -\infty \f$ to 00034 \f$ \mathrm{b} \f$ (GSL) 00035 00036 The number of subdivisions of the original interval which 00037 this class is allowed to make is dictated by the workspace 00038 size for the integration class, which can be set using 00039 \ref gsl_inte_table::set_wkspace() . 00040 00041 */ 00042 template<class param_t, class func_t=funct<param_t> > class gsl_inte_qagil : 00043 public gsl_inte_transform<param_t,func_t> { 00044 00045 #ifndef DOXYGEN_INTERNAL 00046 00047 protected: 00048 00049 /// Store the upper limit 00050 double lb; 00051 00052 #endif 00053 00054 public: 00055 00056 /** 00057 \brief Integrate function \c func from \f$ -\infty \f$ to \c b. 00058 00059 The value given in \c a is ignored. 00060 */ 00061 virtual double integ(func_t &func, double a, double b, param_t &pa) { 00062 double res, err; 00063 integ_err(func,a,b,pa,res,err); 00064 this->interror=err; 00065 return res; 00066 } 00067 00068 /** 00069 \brief Integrate function \c func from \f$ -\infty \f$ to \c b 00070 and place the result in \c res and the error in \c err2 00071 00072 The value given in \c a is ignored. 00073 */ 00074 virtual int integ_err(func_t &func, double a, double b, 00075 param_t &pa, double &res, double &err2) { 00076 00077 lb=b; 00078 00079 double fv1[8], fv2[8]; 00080 00081 int status=qags(func,8,o2scl_inte_qag_coeffs::qk15_xgk, 00082 o2scl_inte_qag_coeffs::qk15_wg, 00083 o2scl_inte_qag_coeffs::qk15_wgk, 00084 fv1,fv2,0.0,1.0,this->tolx,this->tolf, 00085 this->wkspace,&res,&err2,pa); 00086 return status; 00087 00088 } 00089 00090 protected: 00091 00092 /// Transform to \f$ t \in (0,1] \f$ 00093 virtual double transform(func_t &func, double t, param_t &pa) { 00094 double x=lb-(1-t)/t, y; 00095 func(x,y,pa); 00096 return y/t/t; 00097 } 00098 00099 }; 00100 00101 #ifndef DOXYGENP 00102 } 00103 #endif 00104 00105 #endif
Documentation generated with Doxygen and provided under the GNU Free Documentation License. See License Information for details.
Project hosting provided by
,
O2scl Sourceforge Project Page