Does someone have a code for Alternating Least Squares in C ??? I'm working on face recognition and I would like to minimize the loss function : (X-A.B)^2 where X, A and B are matrices. Knowing X, first fix A and calculate B, then fix B and calculate A and repeat until convergence. I've been searching for source codes but I haven't found this algorithm for matrices... Thanks for help