We present the status of our calculation of hadronic matrix elements for $D$- and $B$-meson mixing. We use a large set of the MILC collaboration's $N_f=2+1$ asqtad ensembles, which includes lattice spacings in the range $a\approx0.12$-0.045 fm, and up/down to strange quark mass ratios as low as 0.05. The asqtad action is also employed for the light valence quarks. For the heavy quarks we use the Sheikholeslami-Wohlert action with the Fermilab interpretation. Our calculation covers the complete set of five local operators needed to describe $B$-meson mixing in the Standard Model and Beyond. In the charm sector, our calculation of local mixing matrix elements may be used to constrain new physics models. We present final correlator fit results on the full data set for the $B$-meson mixing project and preliminary fit results for the $D$-meson mixing project.