Revision 064f27d90fcbda7122d5b74508651e2917cc3feb authored by Doug and Martin on 17 October 2011, 00:00:00 UTC, committed by Gabor Csardi on 17 October 2011, 00:00:00 UTC
1 parent a707be9
dtCMatrix.R
#### Triangular Sparse Matrices in compressed column-oriented format
setAs("dtCMatrix", "ltCMatrix",
function(from) new("ltCMatrix", i = from@i, p = from@p,
uplo = from@uplo, diag = from@diag,
x = as.logical(from@x),
## FIXME?: use from@factors smartly
Dim = from@Dim, Dimnames = from@Dimnames))
setAs("dtCMatrix", "ntCMatrix", # just drop 'x' slot:
function(from) new("ntCMatrix", i = from@i, p = from@p,
uplo = from@uplo, diag = from@diag,
## FIXME?: use from@factors smartly
Dim = from@Dim, Dimnames = from@Dimnames))
setAs("matrix", "dtCMatrix",
function(from) as(as(from, "dtTMatrix"), "dtCMatrix"))
setAs("dtCMatrix", "dgCMatrix",
function(from) {
if (from@diag == "U")
from <- .Call(Csparse_diagU2N, from)
new("dgCMatrix",
i = from@i, p = from@p, x = from@x,
Dim = from@Dim, Dimnames = from@Dimnames)
})
setAs("dtCMatrix", "dsCMatrix", function(from) as(from, "symmetricMatrix"))
setAs("dtCMatrix", "dgTMatrix",
function(from) {
if (from@diag == "U") from <- .Call(Csparse_diagU2N, from)
## ignore triangularity in conversion to TsparseMatrix
.Call(Csparse_to_Tsparse, from, FALSE)
})
## FIXME: make more efficient
## ----- and as(., "triangularMatrix") is even worse via as_Sp()
setAs("dgCMatrix", "dtCMatrix", # to triangular, needed for triu,..
function(from) as(.Call(Csparse_to_Tsparse, from, FALSE), "dtCMatrix"))
setAs("dtCMatrix", "dgeMatrix",
function(from) as(as(from, "dgTMatrix"), "dgeMatrix"))
## These are all needed because cholmod doesn't support triangular:
## (see end of ./Csparse.R ), e.g. for triu()
setAs("dtCMatrix", "dtTMatrix",
function(from) .Call(Csparse_to_Tsparse, from, TRUE))
## {# and this is not elegant:
## x <- as(from, "dgTMatrix")
## if (from@diag == "U") { ## drop diagonal entries '1':
## i <- x@i; j <- x@j
## nonD <- i != j
## xx <- x@x[nonD] ; i <- i[nonD] ; j <- j[nonD]
## } else {
## xx <- x@x; i <- x@i; j <- x@j
## }
## new("dtTMatrix", x = xx, i = i, j = j, Dim = x@Dim,
## Dimnames = x@Dimnames, uplo = from@uplo, diag = from@diag)
## })
## Now that we support triangular matrices use the inherited method.
## setAs("dtCMatrix", "TsparseMatrix", function(from) as(from, "dtTMatrix"))
setAs("dtCMatrix", "dtrMatrix",
function(from) as(as(from, "dtTMatrix"), "dtrMatrix"))
setMethod("determinant", signature(x = "dtCMatrix", logarithm = "logical"),
function(x, logarithm = TRUE, ...) {
if(x@diag == "N")
mkDet(diag(x), logarithm)
else
structure(list(modulus = structure(if (logarithm) 0 else 1,
"logarithm" = logarithm),
sign = 1L),
class = "det")
})
setMethod("solve", signature(a = "dtCMatrix", b = "missing"),
function(a, b, ...) {
stopifnot((n <- nrow(a)) == ncol(a))
as(.Call(dtCMatrix_sparse_solve, a, .trDiagonal(n)),
"dtCMatrix")
}, valueClass = "dtCMatrix")
setMethod("solve", signature(a = "dtCMatrix", b = "dgeMatrix"),
function(a, b, ...) .Call(dtCMatrix_matrix_solve, a, b, TRUE),
valueClass = "dgeMatrix")
setMethod("solve", signature(a = "dtCMatrix", b = "CsparseMatrix"),
function(a, b, ...) .Call(dtCMatrix_sparse_solve, a, b),
valueClass = "dgCMatrix")
setMethod("solve", signature(a = "dtCMatrix", b = "matrix"),
function(a, b, ...) {
storage.mode(b) <- "double"
.Call(dtCMatrix_matrix_solve, a, b, FALSE)
}, valueClass = "dgeMatrix")
## Isn't this case handled by the method for (a = "Matrix', b =
## "numeric") in ./Matrix.R? Or is this method defined here for
## the as.double coercion?
setMethod("solve", signature(a = "dtCMatrix", b = "numeric"),
function(a, b, ...) .Call(dtCMatrix_matrix_solve, a,
as.matrix(as.double(b)), FALSE),
valueClass = "dgeMatrix")
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...