Multivariate (and univariate) algorithms for log-likelihood estimation of arbitrary covariance matrix/trees