Running 'testgsvd.R' [0s]
  Running 'testgv.R' [143s]
  Running 'testqz.R' [0s]
  Running 'tgeigen1.R' [0s]
  Comparing 'tgeigen1.Rout' to 'tgeigen1.Rout.save' ... OK
  Running 'tgeigen2.R' [0s]
  Comparing 'tgeigen2.Rout' to 'tgeigen2.Rout.save' ... OK
  Running 'tgeigen3.R' [163s]
  Running 'tgeigen4.R' [0s]
  Comparing 'tgeigen4.Rout' to 'tgeigen4.Rout.save' ... OK
  Running 'tgqz1.R' [0s]
  Comparing 'tgqz1.Rout' to 'tgqz1.Rout.save' ... OK
  Running 'tgqz2.R' [0s]
  Comparing 'tgqz2.Rout' to 'tgqz2.Rout.save' ... OK
  Running 'tgqz3.R' [0s]
  Comparing 'tgqz3.Rout' to 'tgqz3.Rout.save' ... OK
  Running 'tgqz4.R' [0s]
  Comparing 'tgqz4.Rout' to 'tgqz4.Rout.save' ... OK
  Running 'tgqz5.R' [153s]
  Running 'tgqz6.R' [156s]
  Running 'tgqz7.R' [172s]
  Running 'tgqz8.R' [0s]
  Comparing 'tgqz8.Rout' to 'tgqz8.Rout.save' ... OK
  Running 'tgqz9.R' [156s]
  Running 'tgsvd1.R' [0s]
  Comparing 'tgsvd1.Rout' to 'tgsvd1.Rout.save' ... OK
  Running 'tgsvd2.R' [165s]
  Running 'tgsvd3.R' [159s]
  Running 'tgsvd4.R' [0s]
  Comparing 'tgsvd4.Rout' to 'tgsvd4.Rout.save' ... OK
Running the tests in 'tests/testgv.R' failed.
Complete output:
  > testgv <- function(A,B,z) {
  +     V <- z$vectors         
  +     if( is.null(z$alpha) ) {
  +         ret <- A %*% V - B %*% V %*% diag(z$values)
  +     } else {
  +         ret <- A %*% V %*% diag(z$beta) - B %*% V %*% diag(z$alpha)
  +     }
  +     tol <- 100 * sqrt(.Machine$double.eps)
  +     all(abs(ret)<=tol)
  + }
  > 
  > proc.time()
     user  system elapsed 
     0.18    0.04    0.18 
Running the tests in 'tests/tgeigen3.R' failed.
Complete output:
  > 
  > library("geigen")
  > source("testgv.R")
  > 
  > # tol <-  0 triggers buglet in complex QR
  > tol <- .Machine$double.eps
  > set.seed(11)
  > 
  > n <- 5
  > A.save <- matrix(runif(n*n),nrow=n)
  > B.save <- matrix(runif(n*n),nrow=n)
  > 
  > A <- A.save
  > B <- B.save
  > B[5,] <- (B[4,]+B[3,])/2 + tol
  > z <- geigen(A,B)
  > testgv(A,B,z)
  [1] TRUE
  > 
  > A <- A.save+0i
  > B <- B.save+0i
  > B[5,] <- (B[4,]+B[3,])/2 + tol
  > z <- geigen(A,B)
  > testgv(A,B,z)
  [1] TRUE
  > 
  > A <- A.save
  > B <- B.save
  > B[1,] <- (B[2,]+B[3,])/2 + tol
  > z <- geigen(A,B)
  > testgv(A,B,z)
  [1] TRUE
  > 
  > A <- A.save+0i
  > B <- B.save+0i
  > B[1,] <- (B[2,]+B[3,])/2 + tol
  > z <- geigen(A,B)
  > testgv(A,B,z)
  [1] TRUE
  > 
  > proc.time()
     user  system elapsed 
     0.18    0.04    0.21 
Running the tests in 'tests/tgqz5.R' failed.
Complete output:
  > 
  > library("geigen")
  > source("testqz.R")
  > 
  > # example from NAG (F08XNF)
  > 
  > A <- matrix(c( -21.1-22.5i , 53.5-50.5i, -34.5+127.5i ,   7.5+.5i,
  +                  -.46-7.78i, -3.5-37.5i, -15.5+ 58.5i , -10.5-1.5i,
  +                  4.3-5.5i  , 39.7-17.1i, -68.5+12.5i  ,  -7.5-3.5i,
  +                  5.5+4.4i  , 14.4+43.3i, -32.5-46.0i  , -19.0-32.5i), nrow=4, byrow=TRUE)
  > 
  > B <- matrix(c( 1.0-5i,  1.6+1.2i , -3.0+0i,       -1i,
  +                .8-.6i,  3-5i     , -4+3i  , -2.4-3.2i,
  +                1.0+0i,  2.4+1.8i , -4-5i  ,       -3i,
  +                    1i, -1.8+2.4i ,   -4i  ,      4-5i ), nrow=4, byrow=TRUE)
  > 
  > z <- gqz(A, B,"R")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > 
  > z <- gqz(A, B)
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > 
  > zn <- gqz(A, B, "N")
  > testqz(A,B,zn)
  [1] TRUE TRUE TRUE TRUE
  > 
  > identical(z,zn)
  [1] TRUE
  > 
  > proc.time()
     user  system elapsed 
     0.21    0.03    0.25 
Running the tests in 'tests/tgqz6.R' failed.
Complete output:
  > 
  > library("geigen")
  > source("testqz.R")
  > 
  > set.seed(11)
  > 
  > n <- 5
  > A <- matrix(runif(n*n),nrow=n)+0i
  > B <- matrix(runif(n*n),nrow=n)+0i
  > 
  > # Test interface to zgges (QZ method)
  > 
  > z <- gqz(A, B,"-")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(Re(gev)<0) == seq_len(z$sdim))
  [1] TRUE
  > 
  > z <- gqz(A, B,"+")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(Re(gev)>0) == seq_len(z$sdim))
  [1] TRUE
  > 
  > z <- gqz(A, B,"S")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(abs(gev)<1) == seq_len(z$sdim))
  [1] TRUE
  > 
  > z <- gqz(A, B,"B")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(abs(gev)>1) == seq_len(z$sdim))
  [1] TRUE
  > 
  > z <- gqz(A, B,"R")
  > testqz(A,B,z) 
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(abs(Im(gev)) <= 100*.Machine$double.eps) == seq_len(z$sdim))
  [1] TRUE
  > 
  > proc.time()
     user  system elapsed 
     0.26    0.06    0.31 
Running the tests in 'tests/tgqz7.R' failed.
Complete output:
  > 
  > library("geigen")
  > source("testqz.R")
  > 
  > set.seed(11)
  > 
  > n <- 5
  > A <- matrix(runif(n*n),nrow=n)
  > B <- matrix(runif(n*n),nrow=n)
  > 
  > B[5,] <- (B[4,]+B[3,])/2
  > 
  > # Test interface to dgges (QZ method)
  > 
  > z <- gqz(A, B,"-")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(ifelse(is.finite(gev),Re(gev)<0,FALSE)) == seq_len(z$sdim))
  [1] TRUE
  > 
  > z <- gqz(A, B,"+")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(ifelse(is.finite(gev),Re(gev)>0,FALSE)) == seq_len(z$sdim))
  [1] TRUE
  > 
  > z <- gqz(A, B,"S")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(ifelse(is.finite(gev),abs(gev)<1,FALSE)) == seq_len(z$sdim))
  [1] TRUE
  > 
  > z <- gqz(A, B,"B")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(ifelse(is.finite(gev),abs(gev)>1,FALSE)) == seq_len(z$sdim))
  [1] TRUE
  > 
  > z <- gqz(A, B,"R")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(ifelse(is.finite(gev),abs(Im(gev)) <= 100*.Machine$double.eps,FALSE)) == seq_len(z$sdim))
  [1] TRUE
  > 
  > proc.time()
     user  system elapsed 
     0.23    0.06    0.28 
Running the tests in 'tests/tgqz9.R' failed.
Complete output:
  > 
  > library("geigen")
  > source("testqz.R")
  > 
  > set.seed(11)
  > 
  > n <- 5
  > A <- matrix(runif(n*n),nrow=n)+0i
  > B <- matrix(runif(n*n),nrow=n)+0i
  > 
  > B[1,] <- (B[2,]+B[3,])/2
  > 
  > # Test interface to zgges (QZ method)
  > 
  > z <- gqz(A, B,"N")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > z$sdim == 0
  [1] TRUE
  > 
  > z <- gqz(A, B,"-")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(ifelse(is.finite(gev),Re(gev)<0,FALSE)) == seq_len(z$sdim))
  [1] TRUE
  > 
  > z <- gqz(A, B,"+")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(ifelse(is.finite(gev),Re(gev)>0,FALSE)) == seq_len(z$sdim))
  [1] TRUE
  > 
  > z <- gqz(A, B,"S")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(ifelse(is.finite(gev),abs(gev)<1,FALSE)) == seq_len(z$sdim))
  [1] TRUE
  > 
  > z <- gqz(A, B,"B")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(ifelse(is.finite(gev),abs(gev)>1,FALSE)) == seq_len(z$sdim))
  [1] TRUE
  > 
  > z <- gqz(A, B,"R")
  > testqz(A,B,z)
  [1] TRUE TRUE TRUE TRUE
  > gev <- gevalues(z)
  > all(which(ifelse(is.finite(gev),abs(Im(gev)) <= 100*.Machine$double.eps,FALSE)) == seq_len(z$sdim))
  [1] TRUE
  > 
  > proc.time()
     user  system elapsed 
     0.26    0.07    0.32 
Running the tests in 'tests/tgsvd2.R' failed.
Complete output:
  > 
  > library("geigen")
  > source("testgsvd.R")
  > 
  > # Example from NAG F08VAF
  > 
  > A <- matrix(c(1,2,3,3,2,1,4,5,6,7,8,8), nrow=4, ncol=3, byrow=TRUE)
  > 
  > B <- matrix(c(-2,-3,3,4,6,5), nrow=2, ncol=3, byrow=TRUE)
  > 
  > A
       [,1] [,2] [,3]
  [1,]    1    2    3
  [2,]    3    2    1
  [3,]    4    5    6
  [4,]    7    8    8
  > B
       [,1] [,2] [,3]
  [1,]   -2   -3    3
  [2,]    4    6    5
  > 
  > z <- gsvd(A,B)
  > testgsvd(z,A,B)
  [1] TRUE TRUE TRUE TRUE
  > 
  > # B modified
  > 
  > B <- matrix(1:18,byrow=TRUE, ncol=3)
  > B
       [,1] [,2] [,3]
  [1,]    1    2    3
  [2,]    4    5    6
  [3,]    7    8    9
  [4,]   10   11   12
  [5,]   13   14   15
  [6,]   16   17   18
  > 
  > z <- gsvd(A,B)
  > testgsvd(z,A,B)
  [1] TRUE TRUE TRUE TRUE
  > 
  > proc.time()
     user  system elapsed 
     0.21    0.03    0.25 
Running the tests in 'tests/tgsvd3.R' failed.
Complete output:
  > 
  > library("geigen")
  > source("testgsvd.R")
  > 
  > # double
  > set.seed(11)
  > A <- matrix(round(runif(18),2), nrow=3, byrow=TRUE)
  > B <- matrix(round(runif(30),2), ncol=6, byrow=TRUE)
  > z <- gsvd(A,B)
  > testgsvd(z,A,B)
  [1] TRUE TRUE TRUE TRUE
  > 
  > set.seed(11)
  > A <- matrix(round(runif(18),2), nrow=3, byrow=TRUE)
  > B <- matrix(round(runif(24),2), ncol=6, byrow=TRUE)
  > z <- gsvd(A,B)
  > testgsvd(z,A,B)
  [1] TRUE TRUE TRUE TRUE
  > 
  > set.seed(11)
  > A <- matrix(round(runif(36),2), ncol=6, byrow=TRUE)
  > B <- matrix(round(runif(36),2), ncol=6, byrow=TRUE)
  > z <- gsvd(A,B)
  > testgsvd(z,A,B)
  [1] TRUE TRUE TRUE TRUE
  > 
  > set.seed(11)
  > A <- matrix(round(runif(36),2), ncol=6, byrow=TRUE)
  > B <- matrix(round(runif(24),2), ncol=6, byrow=TRUE)
  > z <- gsvd(A,B)
  > testgsvd(z,A,B)
  [1] TRUE TRUE TRUE TRUE
  > 
  > set.seed(11)
  > A <- matrix(round(runif(12),2), nrow=2, byrow=TRUE)
  > B <- matrix(round(runif(24),2), ncol=6, byrow=TRUE)
  > z <- gsvd(A,B)
  > testgsvd(z,A,B)
  [1] TRUE TRUE TRUE TRUE
  > 
  > proc.time()
     user  system elapsed 
     0.18    0.12    0.26