pub fn peak_error_db(wdf: &[f64], reference: &[f64]) -> f64
Compute peak absolute error in dB relative to reference peak.
error_db = 20 * log10(max|wdf - ref| / max|ref|)