MAINT: replace einsum, using a decomposed precision matrix

Replace einsum.  This code is slightly more obfuscated, because it stores a decomposition of the precision matrix instead of the precision matrix itself.