# CP-ALS Tensor Decomposition

I'm trying to implement CP-ALS (alternating least squares algorithm for canonical polyadic decomposition) tensor rank decomposition, but I cannot find any references for good guesses for the matrix initialization. For all intents and purposes, third-order-tensors should suffice. At this point, I'm almost thinking of just using unit matrices.

I have searched Google and Youtube and Wikipedia in vain, and have been unable to find either pseudocode, reference implementations, or just walk-through examples with real numbers. Am I looking in the wrong places?

Thanks!

If you have an N-dimensional tensor, the factor matrices can be initialized as, $$A^{(n)} = R \ leading\ left\ singular\ vectors\ of \ X_{(n)}$$ as outlined on page 472 in this excellent review paper.