[Top][All Lists]

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

[Bug-gsl] [bug #40116] possible error in integration routines

From: Patrick Alken
Subject: [Bug-gsl] [bug #40116] possible error in integration routines
Date: Thu, 26 Sep 2013 21:34:46 +0000
User-agent: Mozilla/5.0 (X11; Linux i686 on x86_64; rv:23.0) Gecko/20100101 Firefox/23.0


                 Summary: possible error in integration routines
                 Project: GNU Scientific Library
            Submitted by: psa
            Submitted on: Thu 26 Sep 2013 09:34:45 PM GMT
                Category: Runtime error
                Severity: 3 - Normal
        Operating System: 
                  Status: None
             Assigned to: None
             Open/Closed: Open
         Discussion Lock: Any



(from address@hidden)


I apologize if this is the wrong place. I’m currently working with the
integration utilities
of GSL and I’d like some help about this particular point.

I’m not sure the function defined in the file integration/qpsrt2.c and
invoked, for example, 
at line 418 of integration/qags.c is doing what it should, according to the
algorithm, nor according to the comment

        /* test whether the interval to be bisected next is the
           smallest interval. */

at line 415 of qags.c. 

In the QUADPACK version (http://www.netlib.org/quadpack/dqagse.f), at line 360
we can 
see the actual test being performed

        if(dabs(blist(maxerr)-alist(maxerr)).gt.small) go to 90

I think the refining loop may fail at some point, since the loop won’t cycle
when the test
evaluates to true.

I don’t have a buggy case to report but it seems to me that it isn’t

Thank you for any help,


Reply to this item at:


  Message sent via/by Savannah

reply via email to

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