[SciPy-user] Numerical stability and Legendre polynomials

[Fernando Perez]

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