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!

A good place to look for code examples is GitHub. For instance, I found this implementation of CP-ALS by just searching for the terms "tensor decomposition" in GitHub's search-box function (it was on page 3).

I am not familiar with Alternating Least Squares, let alone the tensor version of it, so I cannot vouch for the accuracy or correctness of the linked code. I think it's a good starting point. At the very least, this is a way for you to find some reference implementations of your algorithm.