Things that should be checked for and moved to Rcpp: - Any trace functions should be moved from `sum(diag(x))` to `trace_mat` (~5x faster) - Any multiplications of a matrix with a diagonal matrix should be moved to `mult_diag`
Things that should be checked for and moved to Rcpp:
sum(diag(x))totrace_mat(~5x faster)mult_diag