I don’t have Matlab so I’m not sure about the speed vs Matlab.
I’ve begun to look at approximating this decomposition with SGD instead of the full data on every iteration.
A convergence check is interesting, do you have any ideas how to check that? I’ll think about this too.