Parallel methods in the numerical treatment of population dynamic models