From maechler at stat.math.ethz.ch  Wed Mar 24 09:06:20 2004
From: maechler at stat.math.ethz.ch (Martin Maechler)
Date: Wed Mar 24 09:06:22 2004
Subject: cvs.r-project.org offline [Re: [Rd] pgeom accuracy (PR#6687)]
In-Reply-To: <20040324073614.ABE8E1043D@slim.kubism.ku.dk>
References: <20040324073614.ABE8E1043D@slim.kubism.ku.dk>
Message-ID: <16481.16764.340338.672269@gargle.gargle.HOWL>

>>>>> "MM" == Martin Maechler <maechler@stat.math.ethz.ch>
>>>>>     on Wed, 24 Mar 2004 08:36:14 +0100 (CET) writes:

>>>>> "Morten" == Morten Welinder <terra@gnome.org>
>>>>>     on Tue, 23 Mar 2004 19:10:07 +0100 (CET) writes:

    Morten> Full_Name: M. Welinder Version: 1.6.2 and 1.8.1 OS:
    Morten> Solaris Submission from: (NULL) (65.213.85.223)


    Morten> The common lower_tail==TRUE && log_p==FALSE case can
    Morten> be improved as

    Morten> if (lower_tail && !log_p) return -expm1 ((x + 1) *
    Morten> log1p (-p));

    Morten> This improves precision when p is very small and x
    Morten> not too big by avoiding two problematic "1-sth"
    Morten> calculations.

    MM> That's entirely correct, Thank you, Morten!

    MM> Improvement committed for R 1.9.0 beta.  Martin

eehm;  this meant I had done the committment command, but
it timed out :

cvs.r-project.org is currently not available from here.
Martin


cvs.r-project.org offline [Re: [Rd] pgeom accuracy (PR#6687)]
Martin Maechler