  We present two cubature on Wiener space algorithms for the numerical solution
of McKean-Vlasov SDEs with smooth scalar interaction. The analysis hinges on
sharp gradient to time-inhomogeneous parabolic PDEs bounds. These bounds may be
of independent interest. They extend the classical results of Kusuoka \&
Stroock. Both algorithms are tested through two numerical examples.
