Hyperspectral Unmixing Under Endmember Variability: A Variational Inference Framework