Re: [Rd] Inverting a square matrix using solve() with LAPACK=TRUE (PR#13762)

rvarad...@jhmi.edu wrote: Full_Name: Ravi Varadhan Version: 2.8.1 OS: Windows Submission from: (NULL) (162.129.251.19) Inverting a matrix with solve(), but using LAPACK=TRUE, gives erroneous results: Thanks, but there seems to be a much easier fix. Inside coef.qr, we have coef[qr$pivot, ] <

[Rd] Inverting a square matrix using solve() with LAPACK=TRUE (PR#13762)

Full_Name: Ravi Varadhan Version: 2.8.1 OS: Windows Submission from: (NULL) (162.129.251.19) Inverting a matrix with solve(), but using LAPACK=TRUE, gives erroneous results: Here is an example: hilbert <- function(n) { i <- 1:n; 1 / outer(i - 1, i, "+") } h5 <- hilbert(5) h