swh:1:snp:887de40c312a4259a22e92c95293ef0f709b045b
Tip revision: a7001ff1805634d18a10fa371b38dbe3e48f8c9e authored by HwB on 30 October 2011, 00:00:00 UTC
version 0.8.1
version 0.8.1
Tip revision: a7001ff
mtrace.R
##
## m t r a c e . R Matrix trace
##
mtrace <- function(a) {
if (length(a) <= 1) return(a)
if ((!is.numeric(a) && !is.complex(a)) || !is.matrix(a))
stop("Argument 'a' must be a real or complex matrix.")
if (nrow(a) != ncol(a))
stop("Matrix 'a' must be square.")
return(sum(diag(a)))
}