Use generalized ufunc matrix multiply to compute dot products
Implement speed up transformations (rotations) for linear state-space model