Efficient parallel solution of nonlinear parabolic partial differential equations by a probabilistic domain decomposition