  Numerical studies of shock waves in large scale systems via kinetic
simulations with millions of particles are too computationally demanding to be
processed in serial. In this work we focus on optimizing the parallel
performance of a kinetic Monte Carlo code for astrophysical simulations such as
core-collapse supernovae. Our goal is to attain a flexible program that scales
well with the architecture of modern supercomputers. This approach requires a
hybrid model of programming that combines a message passing interface (MPI)
with a multithreading model (OpenMP) in C++. We report on our approach to
implement the hybrid design into the kinetic code and show first results which
demonstrate a significant gain in performance when many processors are applied.
