all.equal-methods {Matrix} | R Documentation |
Methods for function all.equal()
(from R package
base) are defined for all Matrix
classes.
\
\
these three methods are
simply using all.equal.numeric
directly and work via
as.vector()
.
There are more methods, notably also for
"sparseVector"
's, see showMethods("all.equal")
.
showMethods("all.equal") (A <- spMatrix(3,3, i= c(1:3,2:1), j=c(3:1,1:2), x = 1:5)) ex <- expand(lu. <- lu(A)) stopifnot( all.equal(as(A[lu.@p + 1L, lu.@q + 1L], "CsparseMatrix"), lu.@L %*% lu.@U), with(ex, all.equal(as(P %*% A %*% Q, "CsparseMatrix"), L %*% U)), with(ex, all.equal(as(A, "CsparseMatrix"), t(P) %*% L %*% U %*% t(Q))))