A Contour Stochastic Gradient Langevin Dynamics Algorithm for Simulations of Multi-modal Distributions