bugGNU Scientific Library - Bugs: bug #58066, Digamma function returning...

 
 

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

bug #58066: Digamma function returning incorrect values for (most) negative integers

Submitted by:  Jackson Vanover <jacksonvanover>
Submitted on:  Sat 28 Mar 2020 05:27:05 PM UTC
 
Category:  Accuracy problem Severity:  3 - Normal
Operating System:  Ubuntu 18.04 Status:  None
Assigned to:  None Open/Closed:  Open
Release:  2.6
* Mandatory Fields

Add a New Comment (Rich Markup)
   

Comment Type & Canned Response:
       

 

Sat 28 Mar 2020 05:27:05 PM UTC, original submission:  

Using version 2.6 of GSL from ftp://ftp.gnu.org/gnu/gsl/gsl-2.6.tar.gz

OS is Ubuntu 18.04.2 LTS

Hardware is a Dell workstation with a 3.60 GHz Intel i7-4790 and 32 GB
of RAM

Compiler is gcc 7.4.0 with the -w option

Though the digamma function has poles at all negative integer values,
once the magnitude of the input is larger than 2,
GSL does not complain and returns incorrect values. For the inputs -1
and -2, the error handler is invoked with a domain error as might be
expected.

#include <gsl/gsl_sf.h>
#include <stdio.h>
#include <gsl/gsl_math.h>
#include <float.h>

int main (void){
        double out;
        out = gsl_sf_psi(-3);

        printf("%.*f\n", DBL_DIG-1, out);
        printf("%.21e\n", out);
        return 0;
}

Jackson Vanover <jacksonvanover>

 

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

Attach Files:
   
   
Comment:
   

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by jacksonvanover (Submitted the item)
  •  

     

     

     

    No changes have been made to this item

    Back to the top


    Powered by Savane 3.9