The free spectra are KDEs of the PSD in a given PTA dataset without any imposed spectral shape. We need to compare a discrete distribution of gwb spectra to these. I see two methods for this
Compute KDEs for holodeck realizations (slow) NOT SLOW. Going with this.
Bin the free spectra and holodeck realizations similarly so that we can extract features with some kind of embedding network