search for: ceigenok

Displaying 3 results from an estimated 3 matches for "ceigenok".

Did you mean: eigenok
2004 Dec 14
3
make check-all failed on SLES8 for x86_64
...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 stopifnot(Mod(A %*% V - V %*% diag(lam)) < Eps, .... [TRUNCATED] > sm <- cbind(1, 3:1, 1:3) > eigenok(sm, eigen(sm)) Error: abs(A %*% V - V %*% diag(lam)) < Eps is not TRUE This...
2004 Oct 19
1
lapack test failed (PR#7294)
..., 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)) < Eps) + } > > Ceigenok <- function(A, E, Eps=1000*.Machine$double.eps) + { + V <- E$vect; lam <- E$values + stopifnot(Mod(A %*% V - V %*% diag(lam)) < Eps, + Mod(A - V %*% diag(lam) %*% Conj(t(V))) < Eps) + } > > ## failed for some 64bit-Lapack-gcc combinations: > sm <- cb...
2007 Nov 23
1
MacOSX 10.4.11 update breaks tests/lapack.R (R 2.6.0)? (PR#10454)
...digits=4) curr.fun: symbol options .. after 'expression(options(digits=4))' <SNIP> ##> sm[] <- as.complex(sm) curr.fun: symbol <- .. after 'expression(sm[] <- as.complex(sm))' >>>> eval(expression_nr. 35 ) ================= ##> Ceigenok(sm, eigen(sm)) *** caught bus error *** address 0x1, cause 'invalid alignment' Traceback: 1: .Call("La_rs_cmplx", x, only.values, PACKAGE = "base") 2: eigen(sm) 3: apply(A, 2, function(x) x * sign(Re(x[1]))) 4: fixsign(E$vectors) 5: print(fixsign(E$vectors)) 6:...