[Rd] pnorm problem (PR#7302)

[Prof Brian Ripley]

On Thu, 21 Oct 2004, Uwe Ligges wrote:

> terra at gnome.org wrote:
> > Full_Name: Morten Welinder
> > Version: 2
> > OS: Solaris(sparc)
> > Submission from: (NULL) (65.213.85.208)
> > 
> 
> Works for me with R-2.0.0, Solaris 5.7 (UltraSparc), gcc-3.2.3, 32-bit.
> So you have to specify compiler version, OS version, R version etc. much 
> more precisely.

And for me under Solaris 8, gcc-3.4.2 or Sun cc Forte 7, 32- or 64-bit, 
with options as recommended in R-admin.

I suspect this is caused by unsafe compiler options.

> > (gdb) p pnorm(-100000.1, 0.0, 1.0, 0, 1)
> > $42 = NaN(0xfffffffffffff)
> > 
> > Expected: -0
> > 
> > (gdb) p pnorm(-100000.0, 0.0, 1.0, 0, 1)
> > $43 = -0
> > 
> > Good.
> > 
> > (I know this is not a typical usage.)

-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595