Follow-up Comment #1, bug #58066 (group gsl): he implementation in the function psi_x() in the file psi.c checks only 0, -1 and -2 for poles and fails to do this for other negative integers. So instead of
if(x == 0.0 || x == -1.0 || x == -2.0) {
DOMAIN_ERROR(result);
}
one would have to write
if(x <= 0.0 && round(x) == x) {
DOMAIN_ERROR(result);
}
_______________________________________________________
Reply to this item at:
<https://savannah.gnu.org/bugs/?58066>
_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/
signature.asc
Description: PGP signature
