Stochastic variance reduced multiplicative update for nonnegative matrix factorization