  A novel, optimized numerical method of modeling of an exciton-polariton
superfluid in a semiconductor microcavity was proposed. Exciton-polaritons are
spin-carrying quasiparticles formed from photons strongly coupled to excitons.
They possess unique properties, interesting from the point of view of
fundamental research as well as numerous potential applications. However, their
numerical modeling is challenging due to the structure of nonlinear
differential equations describing their evolution. In this paper, we propose to
solve the equations with a modified Runge--Kutta method of 4th order, further
optimized for efficient computations. The algorithms were implemented in form
of C++ programs fitted for parallel environments and utilizing vector
instructions. The programs form the EPCGP suite which have been used for
theoretical investigation of exciton-polaritons.
