From bug-octave-request at bevo dot che dot wisc dot edu Fri Dec 19 11:56:52 2003
Subject: Re: (-1).^.5 == NaN
From: "John W. Eaton" <jwe at bevo dot che dot wisc dot edu>
To: Paul Kienzle <pkienzle at jazz dot ncnr dot nist dot gov>
Cc: bug-octave at bevo dot che dot wisc dot edu
Date: Fri, 19 Dec 2003 11:56:50 -0600

On 19-Dec-2003, Paul Kienzle <pkienzle at jazz dot ncnr dot nist dot gov> wrote:

| > OK, how about the following change?  OTOH, I'm not sure that this
| > should be "fixed" in Octave, but maybe it can eventually be changed
| > back when the buggy libm function is no longer likely to be
| > encountered.
| 
| How about fixing it in lo-mappers instead?

I don't see any pow function in liboctave/lo-mappers.cc, so how would
you fix it there?  Do you want to add an xpow(double,double) function
that works around the bug and then use that instead of std::pow
everywhere else in Octave?

jwe



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------