search for: eigenok

Hi all, I've only now noticed that on our Opteron boxes running Suse Enterprise Server 8 that R-2.0.1 did not pass make chek-all. It bombed in tests/lapack.R at: > eigenok <- function(A, E, Eps = 1000 * .Machine$double.eps) { V <- E$vect lam <- E$values stopifnot(abs(A %*% V - V %*% diag(lam)) < Eps, .... [TRUNCATED] > Ceigenok <- function(A, E, Eps = 1000 * .Machine$double.eps) { V <- E$vect lam <- E$values stopifno...

...t some 64bit combinations failing. Well, my configuration fails. Is this really something to worry about? for R-patched: fi1->cat SVN-REVISION Revision: 39576 Last Changed Date: 2006-10-03 fi1-> > ## failed for some 64bit-Lapack-gcc combinations: > sm <- cbind(1, 3:1, 1:3) > eigenok(sm, eigen(sm)) Error: abs(A %*% V - V %*% diag(lam)) < Eps is not all TRUE > eigenok(sm, eigen(sm, sym=FALSE)) Error: abs(A %*% V - V %*% diag(lam)) < Eps is not all TRUE > I'm not using a nonstandard blas or lapack, just the ones included in R. Here is my system information:...

Hallo, I've just installed R on our Opteron System (AMD Opteron 246). Running make check-all results in a failure during the lapack test: > ## ------- tests of random real and complex matrices ------ > > ## 100 may cause failures here. > eigenok <- function(A, E, Eps=1000*.Machine$double.eps) + { + V <- E$vect; lam <- E$values + stopifnot(abs(A %*% V - V %*% diag(lam)) < Eps, + abs(lam[length(lam)]/lam[1]) < Eps || # this one not for singular A : + abs(A - V %*% diag(lam) %*% t(V)) < Ep...

...tar.gz with the following: RHEL3 on x86_64 gcc (GCC) 3.2.3 20030502 (Red Hat Linux 3.2.3-47) GNU Fortran (GCC 3.2.3 20030502 (Red Hat Linux 3.2.3-47)) 3.2.3 20030502 (Red Hat Linux 3.2.3-47) code snip: > ## failed for some 64bit-Lapack-gcc combinations: > sm <- cbind(1, 3:1, 1:3) > eigenok(sm, eigen(sm)) Error: abs(A %*% V - V %*% diag(lam)) < Eps is not all TRUE Execution halted Thanks in advance, Whit This e-mail message is intended only for the named recipient(s) above. It may contain confidential information. If you are not the intended recipient you are hereby notified...