From b2c75e6881c709bad359f71fd32bcbde591ae239 Mon Sep 17 00:00:00 2001 From: Benjamin Stahl Date: Sat, 28 Mar 2020 20:13:43 +0100 Subject: [PATCH] fixed error in calculation of B_f by adding parenthesis --- pydsm/iso226.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/pydsm/iso226.py b/pydsm/iso226.py index 15a8bf1..e4a3cf8 100644 --- a/pydsm/iso226.py +++ b/pydsm/iso226.py @@ -269,8 +269,8 @@ def tabled_B_f(L_p, hfe=False): reference pressure level P0 (close to the hearing threshold at 1 kHz and set to 20 uPa RMS). """ - B_f = ((0.4*10**(L_p+tbl_L_U)/10.-9.)**tbl_alpha_f - - (0.4*10**(tbl_T_f+tbl_L_U)/10.-9.)**tbl_alpha_f + + B_f = ((0.4*10**((L_p+tbl_L_U)/10.-9.))**tbl_alpha_f - + (0.4*10**((tbl_T_f+tbl_L_U)/10.-9.))**tbl_alpha_f + 0.005135) return np.append(B_f, B_f[0]) if hfe else B_f