Gradient-free variational inference within exponential families using least squares regression with applications on discrete distributions and likelihood-free models.