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]] ______________________________________________ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel