It may be better and faster to use Ensembles for the ODE problems solved in the MPSA example. It may be nice to separate the MPSA functionality into a separate package and to modify the MPSA loop so it uses EnsembleProblems (using the full multi-core potential of Julia)