Particle method and quantization-based schemes for the simulation of the McKean-Vlasov equation

In this paper, we study three numerical schemes for the McKean-Vlasov equation [\begin{cases} \;dXt=b(t, Xt, \mut) \, dt+\sigma(t, Xt, \mut) \, dBt,: \ \;\forall\, t\in[0,T],\;\mut \text{ is the probability distribution of }Xt, \end{cases}] where $X_0$ is a known random variable. Under the assumption on the Lipschitz continuity of the coefficients $b$ and $\sigma$, our first result proves the convergence rate of the particle method with respect to the Wasserstein distance, which extends a previous work [BT97] established in one-dimensional setting. In the second part, we present and analyse two quantization-based schemes, including the recursive quantization scheme (deterministic scheme) in the Vlasov setting, and the hybrid particle-quantization scheme (random scheme, inspired by the $K$-means clustering). Two examples are simulated at the end of this paper: Burger's equation and the network of FitzHugh-Nagumo neurons in dimension 3.