Function thd_error_db 

pub fn thd_error_db(
    wdf: &[f64],
    reference: &[f64],
    fundamental_hz: f64,
    sample_rate: f64,
) -> f64

Compute THD error between WDF and reference in dB.