In the usual approach to the dark matter thermal relic abundance calculation an assumption oflocal thermal equilibrium is made. This proceedings presents a numerical precision tool DRAKE (drake.hepforge.org) that goes beyond this assumption by tracing not only the DM relic density, but also its velocity dispersion and full phase space distribution function. We review the generalmotivation for this approach and, for illustration, highlight two concrete classes of models wherekinetic and chemical decoupling are intertwined in a way that can impact the value of the relicdensity by as much as an order of magnitude: dark matter annihilation via a narrow resonance and‘forbidden’ annihilation to final states that are kinematically inaccessible at threshold.