MAINT: work in progress -- use the structure of the rate matrices to speed up the likelihood calculations