[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Bug-gsl]Bug in gsl_integration_qawc and correction

From: Axel Hutt
Subject: [Bug-gsl]Bug in gsl_integration_qawc and correction
Date: Wed, 08 Jan 2003 18:20:41 +0100

Dear all,

there is a bug in the quadrature routine for Cauchy principal
values in gsl_integration_qawc. 

The problem is, that for a special relation of integration interval 
(a,b) and singularity, the bisection routine in 
$(GSL-ROOT-DIR)/integration/qawc.c might lead to integration borders 
a1=c or b1=c, which yields NaN from qc25c(..) for the integral.
This is easily resolved by changing the check in which subinterval
c lies (line 133) from

if (c > a1 && c < b1)


if (c >= a1 && c <= b1)

(or equivalently for a2 and b2 in line 138), just add equation signs 
in qawc.c in line 133.

Unfortunately, I do not know how to handle patch files, sorry for


reply via email to

[Prev in Thread] Current Thread [Next in Thread]