[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