#' \item 'jd_nvects', number of random vectors for joint-diagonalization
#' (or 0 for p=d, canonical basis by default)
#' \item 'M', moments of order 1,2,3: will be computed if not provided.
#' \item 'jd_nvects', number of random vectors for joint-diagonalization
#' (or 0 for p=d, canonical basis by default)
#' \item 'M', moments of order 1,2,3: will be computed if not provided.