The R Project SVN R

Rev

Rev 74068 | Rev 82074 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 74068 Rev 75221
Line 48... Line 48...
48 
{
 48 
{
49 
    print(fixsign(E$vectors))
 49 
    print(fixsign(E$vectors))
50 
    print(zapsmall(E$values))
 50 
    print(zapsmall(E$values))
51 
    V <- E$vectors; lam <- E$values
 51 
    V <- E$vectors; lam <- E$values
52 
    stopifnot(abs(A %*% V - V %*% diag(lam)) < Eps,
 52 
    stopifnot(abs(A %*% V - V %*% diag(lam)) < Eps,
53 
              abs(lam[length(lam)]/lam[1]) < Eps || # this one not for singular A :
 53 
              abs(lam[length(lam)]/lam[1]) < Eps | # this one not for singular A :
54 
              abs(A - V %*% diag(lam) %*% t(V)) < Eps)
 54 
              abs(A - V %*% diag(lam) %*% t(V)) < Eps)
55 
}
 55 
}
56 
 
 56 
 
57 
Ceigenok <- function(A, E, Eps=1000*.Machine$double.eps)
 57 
Ceigenok <- function(A, E, Eps=1000*.Machine$double.eps)
58 
{
 58 
{