powered by
Check matrix symmetry
check_symmetric_matrix( M, tol = .Machine$double.eps^0.5, arg = as.character(substitute(M)) )