use matrix operations for SOMs wherever possible

[willb]

NB: I will squash this if it passes CI. @erikerlandson if you're interested, take a look but not over the weekend ;-)