Posts

The normal equations, QR, and the SVD solve the same problem with very different accuracy. The difference is the square of a condition number.

What it took to reproduce an old bootstrap analysis, and why its confidence interval moved in the third decimal place.

Recursive summation in single precision loses four digits over ten million terms. Pairwise and compensated summation keep them for little extra work.