Variational Bayesian inference for CP tensor completion with side information