Dear R-helpers, I have a question regarding LU-decomposition with function lu in package Matrix. The following simple example confuses me: Why is as.matrix(elu$U) not an upper triangular matrix?
u3 <- matrix(c(1,1,1,1,1,1,-1,1,0,0,0,0,0,-1,1,0,0,0,-1,0,1,0,0,0,0,0,-1,1,0,0),5,6,byrow=T) elu <- expand(lu(Matrix(u3,sparse=F))) as.matrix(elu$U) I only have very limited experience with the package and its different types of matrices, and I am lost where to start looking for a reason. Regards, Ulrike Grömping -- View this message in context: http://www.nabble.com/Question-regarding-lu-in-package-Matrix-tp18366291p18366291.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
