The pSSAlib provides a complete and concise implementation of the most efficient partial-propensity methods for exact stochastic chemical kinetics.