Not yet reviewed by Pith; the record is open.

Modern advances in algorithms for lattice QCD calculations have steadily driven down the resources required to generate gauge field ensembles and calculate quark propagators, such that, in cases relevant to nuclear physics, performing quark contractions to assemble correlation functions from propagators has become the dominant cost. This work explores a propagator sparsening algorithm for forming correlation functions describing multi-hadron systems, such as light nuclei, with reduced computational cost. The algorithm constructs correlation functions from sparsened propagators defined on a coarsened lattice geometry, where the sparsened propagators are obtained from propagators computed on the full lattice. This algorithm is used to study the low-energy QCD ground-state spectrum using a single Wilson-clover lattice ensemble with $m_{\pi} \approx 800$ MeV. It is found that the extracted ground state masses and binding energies, as well as their statistical uncertainties, are consistent when determined from correlation functions constructed from sparsened and full propagators. In addition, while evidence of modified couplings to excited states is observed in sparsened correlation functions, it is demonstrated that these effects can be removed, if desired, with an inexpensive modification to the sparsened estimator.