R devel - Jun 2008 - phyper hang (PR#11813)

Full_Name: Morten Welinder
Version: 
OS: Linux
Submission from: (NULL) (69.113.253.99)


It appears that phyper hangs for insanely large x, such as the call

    phyper(1e67,0,0,0)

Adding something like this after the tail swap seems to cure that.

    if (x >= NR)
	return R_DT_1;

>>>>> "MW" == Morten Welinder <terra at
gnome.org>
>>>>>     on Thu, 26 Jun 2008 03:55:07 +0200 (CEST) writes:
    MW> Full_Name: Morten Welinder
    MW> Version: 
    MW> OS: Linux
    MW> Submission from: (NULL) (69.113.253.99)


    MW> It appears that phyper hangs for insanely large x, such as the call

    MW> phyper(1e67,0,0,0)

Yes.  I think it only happens for  (x, 0,0,0)  though.

    MW> Adding something like this after the tail swap seems to cure that.

    MW> if (x >= NR)
    MW>    return R_DT_1;

yes, and that also cures  phyper(1,0,0,0)  which gave NaN
instead of 1.
Is going to be fixed for R-devel and R-patched.

Thank you, Matt!

Martin Maechler, ETH Zurich