  We present large $N_f$ QCD/QED as a test bed for improved pressure
calculations, show how to apply the hints obtained on optimized renormalization
scales at large $N_f$ to finite $N_f=2$, and compare the results to recent
lattice data.
