A Flexible Optimization Framework for Regularized Matrix-Tensor Factorizations with Linear Couplings