R/tr.R

In vlulla/vlutils: Common utilities that I find useful

Defines functions tr

Documented in tr

tr <-
function(x) {
  stopifnot(is.matrix(x), dim(x)[[1L]] == dim(x)[[2L]])  # Ensure x is square matrix
  sum(diag(x))
}