Although this can be computed in CPUs within several minutes, more memories will be required with the growth of the amount of training data.
we present an approximate solution, where W is optimized with gradient decent (GD) and is not constrained to be orthogonal.