From terra at gnome.org  Mon Apr 19 17:12:20 2004
From: terra at gnome.org (terra@gnome.org)
Date: Mon Apr 19 17:13:14 2004
Subject: [Rd] pgeom accuracy (PR#6792)
Message-ID: <20040419151220.33B8D104B1@slim.kubism.ku.dk>

Full_Name: Morten Welinder
Version: snapshot
OS: 
Submission from: (NULL) (65.213.85.218)


This should fix the remaining two 1-p cancellation issues.

  double l_rt = log1p (-p) * (x + 1);
  if (log_p)
    return R_DT_Clog (l_rt);
  else
    return lower_tail ? -expm1 (l_rt) : exp (l_rt);


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