Higher order tensor inversion is possible for even order. This is due to the fact that a tensor group endowed with the contracted product is isomorphic to the general linear group of degree n. With these isomorphic group structures, we derive a tensor SVD which we have shown to be equivalent to well-known canonical polyadic decomposition and multilinear SVD provided that some constraints are satisfied. Moreover, within this group structure framework, multilinear systems are derived and solved for problems of high-dimensional PDEs and large discrete quantum models. We also address multilinear systems which do not fit the framework in the least-squares sense. These are cases when there is an odd number of modes or when each mode has distinct dimension. Numerically we solve multilinear systems using iterative techniques, namely, biconjugate gradient and Jacobi methods. © 2013 Society for Industrial and Applied Mathematics.