3
<title>QAGP adaptive integration with known singular points - GNU Scientific Library -- Reference Manual</title>
4
<meta http-equiv="Content-Type" content="text/html">
5
<meta name="description" content="GNU Scientific Library -- Reference Manual">
6
<meta name="generator" content="makeinfo 4.8">
7
<link title="Top" rel="start" href="index.html#Top">
8
<link rel="up" href="Numerical-Integration.html#Numerical-Integration" title="Numerical Integration">
9
<link rel="prev" href="QAGS-adaptive-integration-with-singularities.html#QAGS-adaptive-integration-with-singularities" title="QAGS adaptive integration with singularities">
10
<link rel="next" href="QAGI-adaptive-integration-on-infinite-intervals.html#QAGI-adaptive-integration-on-infinite-intervals" title="QAGI adaptive integration on infinite intervals">
11
<link href="http://www.gnu.org/software/texinfo/" rel="generator-home" title="Texinfo Homepage">
13
Copyright (C) 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006 The GSL Team.
15
Permission is granted to copy, distribute and/or modify this document
16
under the terms of the GNU Free Documentation License, Version 1.2 or
17
any later version published by the Free Software Foundation; with the
18
Invariant Sections being ``GNU General Public License'' and ``Free Software
19
Needs Free Documentation'', the Front-Cover text being ``A GNU Manual'',
20
and with the Back-Cover Text being (a) (see below). A copy of the
21
license is included in the section entitled ``GNU Free Documentation
24
(a) The Back-Cover Text is: ``You have freedom to copy and modify this
25
GNU Manual, like GNU software.''-->
26
<meta http-equiv="Content-Style-Type" content="text/css">
27
<style type="text/css"><!--
28
pre.display { font-family:inherit }
29
pre.format { font-family:inherit }
30
pre.smalldisplay { font-family:inherit; font-size:smaller }
31
pre.smallformat { font-family:inherit; font-size:smaller }
32
pre.smallexample { font-size:smaller }
33
pre.smalllisp { font-size:smaller }
34
span.sc { font-variant:small-caps }
35
span.roman { font-family:serif; font-weight:normal; }
36
span.sansserif { font-family:sans-serif; font-weight:normal; }
42
<a name="QAGP-adaptive-integration-with-known-singular-points"></a>
43
Next: <a rel="next" accesskey="n" href="QAGI-adaptive-integration-on-infinite-intervals.html#QAGI-adaptive-integration-on-infinite-intervals">QAGI adaptive integration on infinite intervals</a>,
44
Previous: <a rel="previous" accesskey="p" href="QAGS-adaptive-integration-with-singularities.html#QAGS-adaptive-integration-with-singularities">QAGS adaptive integration with singularities</a>,
45
Up: <a rel="up" accesskey="u" href="Numerical-Integration.html#Numerical-Integration">Numerical Integration</a>
49
<h3 class="section">16.5 QAGP adaptive integration with known singular points</h3>
51
<p><a name="index-singular-points_002c-specifying-positions-in-quadrature-1352"></a>
54
— Function: int <b>gsl_integration_qagp</b> (<var>const gsl_function * f, double * pts, size_t npts, double epsabs, double epsrel, size_t limit, gsl_integration_workspace * workspace, double * result, double * abserr</var>)<var><a name="index-gsl_005fintegration_005fqagp-1353"></a></var><br>
56
<p>This function applies the adaptive integration algorithm QAGS taking
57
account of the user-supplied locations of singular points. The array
58
<var>pts</var> of length <var>npts</var> should contain the endpoints of the
59
integration ranges defined by the integration region and locations of
60
the singularities. For example, to integrate over the region
61
(a,b) with break-points at x_1, x_2, x_3 (where
62
a < x_1 < x_2 < x_3 < b) the following <var>pts</var> array should be used
64
<pre class="example"> pts[0] = a
70
<p class="noindent">with <var>npts</var> = 5.
72
<p class="noindent">If you know the locations of the singular points in the integration
73
region then this routine will be faster than <code>QAGS</code>.