Non-negative Tensor Factorization Based on Alternating Large-scale Non-negativity-constrained Least Squares