Parallel solvers for simulating morphogenesis