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?