Nonnegative CP Tensor Decomposition by Hierarchical Alternating Least Squares (HALS)ΒΆ