R/matmult.blockmatrix.R

In ecor/blockmatrix: Tools to Solve Algebraic Systems with Partitioned Matrices

NULL

#' 
#'\code{blockmatmult} implements the implents betwwen two blockmatrix ( see \code{\link{matmult}} for \code{matrx} objects)
#' 
#' @name blockmatmult
#' @param x,y  \code{blockmatrix} objects 
#' @param ... further arguments 
#' @export
#' @rdname blockmatmult
# @method matmult blockmatrix
# @S3method matmult blockmatrix
#' @aliases blockmatmult
#'
#### @importFrom base as.matrix
#' @return The inner product between \code{x} and \code{y} as a \code{blockmatrix} object
#' 
#' @author Emanuele Cordano 
#' 
#' 



blockmatmult  <- function (x,y,...)  {
	
	ncol <- ncol(y)
	nrow <- nrow(x)
	
	xm <- as.matrix(x,...)
	ym <- as.matrix(y,...)
	
	
	
	out <- xm %*% ym 
	
	
	return(blockmatrix::as.blockmatrix(out,ncol=ncol,nrow=nrow,...))
	
}