MAINT: use the sylvester-like decomposition to compute likelihoods
MAINT: faster multiplication by diagonal matrices