HOSVD higher order singular value decomposition implemented from chapter 8 of Matrix Methods in Data Mining and Pattern Recognition - Lars Eldén