SCM

SCM Repository

[matrix] View of /pkg/R/dspMatrix.R
ViewVC logotype

View of /pkg/R/dspMatrix.R

Parent Directory Parent Directory | Revision Log Revision Log


Revision 652 - (download) (annotate)
Wed Mar 16 01:08:36 2005 UTC (14 years, 5 months ago) by bates
File size: 3181 byte(s)
Updates to multiplication and other methods
 ### Coercion and Methods for Symmetric Matrices

setAs("dspMatrix", "dsyMatrix",
      function(from) .Call("dspMatrix_as_dsyMatrix", from) )

setAs("dspMatrix", "dgeMatrix",
      function(from) as(as(from, "dsyMatrix"), "dgeMatrix"))

setAs("dspMatrix", "matrix",
      function(from) as(as(from, "dsyMatrix"), "matrix"))

setMethod("rcond", signature(x = "dspMatrix", type = "character"),
          function(x, type, ...)
          .Call("dspMatrix_rcond", x, type),
          valueClass = "numeric")

setMethod("rcond", signature(x = "dspMatrix", type = "missing"),
          function(x, type, ...)
          .Call("dspMatrix_rcond", x, "O"),
          valueClass = "numeric")

setMethod("%*%", signature(x = "dspMatrix", y = "dgeMatrix"),
          function(x, y) .Call("dspMatrix_matrix_mm", x, y, TRUE),
          valueClass = "dgeMatrix")

setMethod("%*%", signature(x = "dspMatrix", y = "matrix"),
          function(x, y) .Call("dspMatrix_matrix_mm", x, y, FALSE),
          valueClass = "dgeMatrix")

setMethod("%*%", signature(x = "dspMatrix", y = "numeric"),
          function(x, y) .Call("dspMatrix_matrix_mm", x, as.matrix(y), FALSE),
          valueClass = "dgeMatrix")

setMethod("%*%", signature(x = "dspMatrix", y = "integer"),
          function(x, y) {
              storage.mode(y) <- "double"
              .Call("dspMatrix_matrix_mm", x, as.matrix(y), FALSE)
          }, valueClass = "dgeMatrix")

setMethod("solve", signature(a = "dspMatrix", b = "missing"),
          function(a, b, ...) .Call("dspMatrix_solve", a),
          valueClass = "dspMatrix")

setMethod("solve", signature(a = "dspMatrix", b = "matrix"),
          function(a, b, ...)
          .Call("dspMatrix_matrix_solve", a, b, FALSE),
          valueClass = "dgeMatrix")

setMethod("solve", signature(a = "dspMatrix", b = "dgeMatrix"),
          function(a, b, ...)
          .Call("dspMatrix_dgeMatrix_solve", a, b, TRUE),
          valueClass = "dgeMatrix")

setMethod("solve", signature(a = "dspMatrix", b = "numeric"),
          function(a, b, ...)
          .Call("dspMatrix_dgeMatrix_solve", a, as.matrix(b), TRUE),
          valueClass = "dgeMatrix")

setMethod("solve", signature(a = "dspMatrix", b = "integer"),
          function(a, b, ...) {
              storage.mode(b) <- "double"
              .Call("dspMatrix_dgeMatrix_solve", a, as.matrix(b), TRUE)
          }, valueClass = "dgeMatrix")

setMethod("norm", signature(x = "dspMatrix", type = "character"),
          function(x, type, ...) .Call("dspMatrix_norm", x, type),
          valueClass = "numeric")

setMethod("norm", signature(x = "dspMatrix", type = "missing"),
          function(x, type, ...) .Call("dspMatrix_norm", x, "O"),
          valueClass = "numeric")

setMethod("t", signature(x = "dsyMatrix"), # need a t method for dsyMatrix
          function(x) as(t(as(x, "dsyMatrix")), "dspMatrix"),
          valueClass = "dspMatrix")

setMethod("unpack", signature(x = "dspMatrix"),
          function(x, ...) as(x, "dsyMatrix"),
          valueClass = "dsyMatrix")

setIs("dspMatrix", "dppMatrix",
      test = function(obj)
          "try-error" != class(try(.Call("dppMatrix_chol", obj), TRUE))
      )


root@r-forge.r-project.org
ViewVC Help
Powered by ViewVC 1.0.0  
Thanks to:
Vienna University of Economics and Business Powered By FusionForge