  High energy scattering processes involving jets generically involve matrix
elements of light- like Wilson lines, known as soft functions. These describe
the structure of soft contributions to observables and encode color and
kinematic correlations between jets. We compute the dijet soft function to
O({\alpha}_s^2) as a function of the two jet invariant masses, focusing on
terms not determined by its renormalization group evolution that have a
non-separable dependence on these masses. Our results include non-global single
and double logarithms, and analytic results for the full set of non-logarithmic
contributions as well. Using a recent result for the thrust constant, we
present the complete O({\alpha}_s^2) soft function for dijet production in both
position and momentum space.
