#### Logical Symmetric Sparse Matrices in Compressed column-oriented format ### contains = "lsparseMatrix" setAs("lsCMatrix", "matrix", function(from) as(as(from, "lgCMatrix"), "matrix")) setAs("lsCMatrix", "lgCMatrix", function(from) .Call(Csparse_symmetric_to_general, from)) ## Specific conversions, should they be necessary. Better to convert as ## as(x, "TsparseMatrix") or as(x, "denseMatrix") setAs("lsCMatrix", "lsTMatrix", function(from) .Call(Csparse_to_Tsparse, from, FALSE)) setAs("lsCMatrix", "dsCMatrix", function(from) new("dsCMatrix", i = from@i, p = from@p, x = rep(1, length(from@i)), uplo = from@uplo, Dim = from@Dim, Dimnames = from@Dimnames)) setAs("lsCMatrix", "dgTMatrix", function(from) as(as(x, "dsCMatrix"), "dgTMatrix")) ## have rather tril() and triu() methods than ## setAs("lsCMatrix", "ltCMatrix", ....) setMethod("tril", "lsCMatrix", function(x, k = 0, ...) { if(x@uplo == "L" && k == 0) ## same internal structure (speedup potential !?) new("ltCMatrix", uplo = x@uplo, i = x@i, p = x@p, Dim = x@Dim, Dimnames = x@Dimnames) else tril(as(x, "lgCMatrix"), k = k, ...) }) setMethod("triu", "lsCMatrix", function(x, k = 0, ...) { if(x@uplo == "U" && k == 0) ## same internal structure (speedup potential !?) new("ltCMatrix", uplo = x@uplo, i = x@i, p = x@p, Dim = x@Dim, Dimnames = x@Dimnames) else triu(as(x, "lgCMatrix"), k = k, ...) }) ## FIXME: generalize to "lsparseMatrix" or (class union) "symmetric sparse" setMethod("image", "lsCMatrix", function(x, ...) { x <- as(as(x, "dsCMatrix"), "dgTMatrix") callGeneric() }) setMethod("chol", signature(x = "lsCMatrix", pivot = "missing"), function(x, pivot, LINPACK) .Call(lsCMatrix_chol, x, TRUE)) setMethod("chol", signature(x = "lsCMatrix", pivot = "logical"), function(x, pivot, LINPACK) .Call(lsCMatrix_chol, x, pivot)) setMethod("t", signature(x = "lsCMatrix"), function(x) .Call(lsCMatrix_trans, x), valueClass = "lsCMatrix")