Re: [SciPy-user] Numerical stability and Legendre polynomials

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