We present a (2+1)-dimensional partial differential equation model for spatial-lateral dynamics of edge-emitting broad-area semiconductor devices and several extensions of this model describing different physical effects. MPI-based parallelization of the resulting middlesize numerical problem is implemented and tested on the blade cluster and separate multi-core computers at the Weierstrass Institute in Berlin. It was found, that an application of 25-30 parallel processes on all considered platforms was guaranteeing a nearly optimal performance of the algorithm with the speedup around 20-25 and the efficiency of 0.7-0.8. It was also shown, that a simultaneous usage of several in-house available multi-core computers allows a further increase of the speedup without a significant loss of the efficiency. Finally, an importance of the considered problem and the efficient numerical simulations of this problem were illustrated by a few examples occurring in real world applications
