inst/ex_verify/dexA2_dgees.r

library(QZ, quiet = TRUE)

### https://www.nag.com/lapack-ex/node89.html
A <- exA2$A
ret <- qz.dgees(A)

# Verify 1
A.new <- ret$Q %*% ret$T %*% solve(ret$Q)
round(A - A.new)

# verify 2
round(ret$Q %*% solve(ret$Q))
snoweye/QZ documentation built on Sept. 12, 2023, 4:59 a.m.