bugGNU Scientific Library - Bugs: bug #26706, gsl_sf_hyperg_U(-2.0, 0.5, 1.13)...

 
 

You are not allowed to post comments on this tracker with your current authentication level.

bug #26706: gsl_sf_hyperg_U(-2.0, 0.5, 1.13) gives incorrect result

Submitter:  -Deleted Account- <bjg>
Submitted:  Thu 28 May 2009 03:20:24 PM UTC
   
 
Category:  Runtime error Severity:  3 - Normal
Operating System:  Status:  Fixed
Assigned to:  None Open/Closed:  Closed
Release:  1.12

Fri 17 Jul 2009 07:07:07 PM UTC, comment #2: 

fixed by 48989e142b764a631d674f5c0ef56b25f8c73a6b

-Deleted Account- <bjg>
Tue 02 Jun 2009 10:11:39 AM UTC, comment #1: 

Added reference

"Computation of Special Functions", Shanjie Zhang and    
    Jianming Jin, John Wiley and Sons, 1996, chapter 12.  

-----------------------------------------------------------
   I simply used the fortran implementation from the attached file. It should be
relatively easy to port that implementation to C.

Kovid.
----------------------------------------

(file #18217)

-Deleted Account- <bjg>
Thu 28 May 2009 03:20:24 PM UTC, original submission:  

From: Kovid Goyal <kovid@kovidgoyal.net>
To: -email is unavailable-
Subject: [Bug-gsl] Bug in gsl_sf_hyperg_U
Date: Wed, 27 May 2009 18:38:37 -0700

[1  <multipart/signed (7bit)>]
[1.1  <text/plain; us-ascii (quoted-printable)>]
Hi,

The Hyper geometric U function gsl_sf_hyperg_U gives incorrect results on a
linux x86 machine (gsl version 1.12).

For example:

gsl_sf_hyperg_U(-2.0, 0.5, 1.13) gives 0.75 which is the wrong value (it is the
expected value for x = 0.

Running gsl_sf_hyperg_U(-2.0, 0.5, 0) gives a domain error.

I'm not subscribed to the list, so please CC replies to -email is unavailable-

Thanks,

Kovid.
--
___________________________________

Dr. Kovid Goyal
http://www.kovidgoyal.net
http://calibre.kovidgoyal.net

-Deleted Account- <bjg>

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attached Files
file #18217:  hyper.f added by bjg (19KiB - application/octet-stream)

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by bjg
  •  

    Follow 4 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2009-07-17 bjg StatusConfirmed Fixed
        Open/ClosedOpen Closed
    2009-06-02 bjg Attached File- Added hyper.f, #18217
    2009-05-28 bjg Carbon-Copy- Added -email is unavailable-

    Back to the top

    Powered by Savane 3.13-3230.
    Corresponding source code