R/bigmemory_methods.R

In fastglm: Fast and Stable Fitting of Generalized Linear Models using 'RcppEigen'

big.crossprod <- function(x)
{
    if (!is.big.matrix(x))
    {
        stop("object must be a big.matrix object")
    }
    .Call("crossprod_big", x@address)
}

big.colSums <- function(x)
{
    if (!is.big.matrix(x))
    {
        stop("object must be a big.matrix object")
    }
    .Call("colsums_big", x@address)
}

big.colMax <- function(x)
{
    if (!is.big.matrix(x))
    {
        stop("object must be a big.matrix object")
    }
    .Call("colmax_big", x@address)
}

big.colMin <- function(x)
{
    if (!is.big.matrix(x))
    {
        stop("object must be a big.matrix object")
    }
    .Call("colmin_big", x@address)
}


#' big.matrix prod
#'
#' @param x big.matrix
#' @param y numeric vector
#' @docType methods
#' @import methods
#' @importFrom methods new
#' @importFrom bigmemory as.big.matrix is.big.matrix
#' @rdname extract-methods
setMethod("%*%",signature(x="big.matrix", y="vector"),
          function(x, y)
          {
              if(dim(x)[2] != length(y)) stop("non-conformant matrices")
              return( .Call("prod_vec_big", x@address, y) )
          },
          valueClass="vector"
)

#' big.matrix prod
#'
#' @docType methods
#' @rdname extract-methods
setMethod("%*%",signature(x="vector", y="big.matrix"),
          function(x, y)
          {
              if(dim(y)[2] != length(x)) stop("non-conformant matrices")
              return( .Call("prod_vec_big_right", x, y@address) )
          },
          valueClass="vector"
)