[Bug-gsl]Bug in gsl_integration_qawc and correction

Brian Gough brian.gough@network-theory.co.uk
Sun Jan 12 19:07:00 GMT 2003


Axel Hutt writes:
 > > Thanks for the bug report -- can you send a small example program
 > > which demonstrates the problem.
 > here comes a tar.gz with two pieces of source code: a main routine
 > for the call and a commented qawc.c . In a README, I have added
 > some comments.

Thank you. Now fixed in CVS.  

I used the following change (I think the first <= is sufficient):

Brian

Index: integration/qawc.c
===================================================================
RCS file: /cvs/gsl/gsl/integration/qawc.c,v
retrieving revision 1.9
diff -r1.9 qawc.c
133c133
<       if (c > a1 && c < b1) 
---
>       if (c > a1 && c <= b1) 



More information about the Gsl-discuss mailing list