From maechler at stat.math.ethz.ch  Fri Jun 27 16:45:18 2008
From: maechler at stat.math.ethz.ch (maechler at stat.math.ethz.ch)
Date: Fri, 27 Jun 2008 16:45:18 +0200 (CEST)
Subject: [Rd] phyper hang (PR#11813)
Message-ID: <20080627144518.BD20B2834337@mail.pubhealth.ku.dk>

>>>>> "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



[Rd] phyper hang (PR#11813)
maechler at stat.math.ethz.ch