[SciPy-Dev] PR: Implement higher-order SVD in scipy.linalg