einsum

Hej,

Jag har riktigt svårt att förstå hur einsum i exempelvis numpy fungerar. Finns det någon som har läst en riktigt bra förklaring eller kanske sitter på en egen riktigt bra förklaring?

Jag är med på grundläggande syntax, att om jag skriver np.einsum('ij,kl'->'ijkl', A, B) har jag döpt dimensionerna i A till i och j, och i B för k och l. Men i princip här slutar det fungera i min hjärna och det stör mig otroligt. Tex, om jag repeterar ett index som np.einsum('ij,kj'->'ikj', A, B) blir det riktig kortslutning. Jag kan såklart memorera att exempelvis np.einsum('ij,jk->ik',A,B) är matris-multiplikation, men säg att jag vill göra $B^TA$, jag kan liksom inte klura ut hur det skulle vara