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.
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
