[SciPy-user] Numerical stability and Legendre polynomials

Fernando Perez Fernando.Perez at colorado.edu
Fri Sep 9 13:58:59 CDT 2005


Yaroslav Bulatov wrote:

> I'm using 64 bit arithmetic (floats), scipy version 0.3.2, why do you
> get higher accuracy?

I'm not really.  Relative differences O(1e-15) in floating point (these are 
absoulute, but we're comparing against zero) are essentially numerical noise. 
  Compiler differences, libc differences, whatever.  For all intents and 
purposes, the results are quite similar and I would not worry about the 
discrepancy.

Cheers,

f



More information about the SciPy-user mailing list