pub fn normalized_rms_error_db(wdf: &[f64], reference: &[f64]) -> f64
Compute normalized RMS error in dB.
error_db = 20 * log10(rms(wdf - ref) / rms(ref))
Lower is better. -60 dB means error is 0.1% of signal.