[Top][All Lists]

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

Re: [Bug-gsl] [bug #41837] bugs in gsl_sf_hyperg_U

From: Raymond Rogers
Subject: Re: [Bug-gsl] [bug #41837] bugs in gsl_sf_hyperg_U
Date: Fri, 14 Mar 2014 18:00:27 -0400
User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:24.0) Gecko/20100101 Thunderbird/24.3.0

Hi Patrick,

Actually I didn't think so; but upon examination it seems your right about the infinite being a correction (of some sort). But calling the infinite sum with these values is causing gamma to be called with a zero; which it shouldn't be. Correcting the finite sum with an infinite one is not what I expected. Sigh.. I will look into what is going on in the infinite sum more carefully. This is not the easiest code to follow :) This time I will start altering the code to check my analysis, and write down the program flow. (:


On 03/14/2014 02:29 PM, Patrick Alken wrote:
Hi Ray,

Thanks for your report. If you have a correction you can upload it to the bug tracker or just post it to this list.

I don't know much about these functions but maybe someone else on the list can check it. Based on that code snippet, it looks like it starts off with a finite sum approximation to the true value, and then refines that approximation with an infinite sum? It seems to pass in the finite sum value to the infinite sum routine.

But you're saying under certain argument conditions, its not necessary to refine the initial approximation with the infinite sum?


On 03/14/2014 12:24 PM, Raymond Rogers wrote:
I have added an analysis to the bug report.   Somebody needs to double
check it; I usually get things a little wrong the first time through.
If it should happen that my suggestion is valid; how and where do I
submit a correction?


On 03/11/2014 06:24 PM, Patrick Alken wrote:
Additional Item Attachment, bug #41837 (project gsl):

File name: gsl_sf_hyperg_U_bug_report.cpp Size:0 KB


Reply to this item at:


    Message sent via/by Savannah

Act IV, Sc. IV
What is a man,
If his chief good and profit of his time
Be to sleep and feed.  Be a beast, no more

reply via email to

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