  We propose a micro-macro parallel-in-time Parareal method for scalar
McKean-Vlasov stochastic differential equations (SDEs). In the algorithm, the
fine Parareal propagator is a Monte Carlo simulation of an ensemble of
particles, while an approximate ordinary differential equation (ODE)
description of the mean and the variance of the particle distribution is used
as a coarse Parareal propagator to achieve speedup. We analyse the convergence
behaviour of our method for a linear problem and provide numerical experiments
indicating the parallel weak scaling of the algorithm on a set of examples. We
show that convergence typically takes place in a low number of iterations,
depending on the quality of the ODE predictor. For bimodal SDEs, we avoid
quality deterioration of the coarse predictor (compared to unimodal SDEs)
through the usage of multiple ODEs, each describing the mean and variance of
the particle distribution in locally unimodal regions of the phase space. The
benefit of the proposed algorithm can be viewed through two lenses: (i) through
the parallel-in-time lens, speedup is obtained through the use of a very cheap
coarse integrator (an ODE moment model), and (ii) through the moment models
lens, accuracy is iteratively gained through the use of parallel machinery as a
corrector. In contrast to the isolated use of a moment model, the proposed
method (iteratively) converges to the true distribution generated by the SDE.
