Skip to content
Prev 52374 / 63424 Next
Original Raw

Numerical accuracy of matrix multiplication

Alexis Sarda 
Hello,

while testing the crossprod() function under Linux, I noticed the following:

set.seed(883)
x <- rnorm(100)
x %*% x - sum(x^2) # equal to 1.421085e-14

Is this difference normal? It seems to be rather large for double precision.

Regards,
Alexis.

Thread (5 messages)

Alexis Sarda Numerical accuracy of matrix multiplication Sep 16 Peter Dalgaard Numerical accuracy of matrix multiplication Sep 16 Martin Maechler Numerical accuracy of matrix multiplication Sep 20 Alexis Sarda Numerical accuracy of matrix multiplication Sep 20 Martin Maechler Numerical accuracy of matrix multiplication Sep 20