Raw Message

Message-ID: <20040419151220.33B8D104B1@slim.kubism.ku.dk>
Date: 2004-04-19T17:13:14Z
From: terra@diku.dk
Subject: pgeom accuracy (PR#6792)

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);