From terra at gnome.org  Tue Mar 23 19:10:07 2004
From: terra at gnome.org (terra@gnome.org)
Date: Tue Mar 23 19:10:14 2004
Subject: [Rd] pgeom accuracy (PR#6687)
Message-ID: <20040323181007.CA72E1043D@slim.kubism.ku.dk>

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


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

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

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


[Rd] pgeom accuracy (PR#6687)
terra@gnome.org