Monte Carlo simulation routines for high-performance parallelization of differential equation solvers and scientific machine learning