A universal source of difficulty in the numerical solution of differential equations is stiffness, which may stem from multi-scale dynamics and means that straightforward explicit discretizations are very inefficient. Consequently, choosing how to deal with stiffness (usually through the selection of the time discretization method) is often the most impactful decision made while solving these problems; it can be the difference between complete intractability and rapid solution.