orastron/brickworks
2
1
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases 10 Wiki Activity

made bw_allpass1 supposedly super robust

Browse Source
This commit is contained in:
Stefano D'Angelo 2022-12-13 15:37:39 +01:00
parent 7c941fa26b
commit 84c3a7d1cd
1 changed files with 9 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
include/bw_allpass1.h
View File

@ -130,7 +130,9 @@ struct _bw_allpass1_coeffs {
float t_k; float t_k;
float t; float t;
float X_k; float X_x;
float X_X_z1;
float lp_X;
// Parameters // Parameters
float cutoff; float cutoff;
@ -163,7 +165,10 @@ static inline void bw_allpass1_reset_state(const bw_allpass1_coeffs *BW_RESTRICT
static inline void bw_allpass1_update_coeffs_ctrl(bw_allpass1_coeffs *BW_RESTRICT coeffs) { static inline void bw_allpass1_update_coeffs_ctrl(bw_allpass1_coeffs *BW_RESTRICT coeffs) {
if (coeffs->cutoff != coeffs->cutoff_prev) { if (coeffs->cutoff != coeffs->cutoff_prev) {
coeffs->t = bw_tanf_3(coeffs->t_k * coeffs->cutoff); coeffs->t = bw_tanf_3(coeffs->t_k * coeffs->cutoff);
coeffs->X_k = bw_rcpf_2(1,f + coeffs->t); const float k = bw_rcpf_2(1.f + coeffs->T);
coeffs->X_x = k * coeffs->cutoff;
coeffs->X_X_z1 = k * coeffs->t;
coeffs->lp_X = bw_rcpf_2(coeffs->cutoff);
coeffs->cutoff_prev = coeffs->cutoff; coeffs->cutoff_prev = coeffs->cutoff;
} }
} }
@ -172,8 +177,8 @@ static inline void bw_allpass1_update_coeffs_audio(bw_allpass1_coeffs *BW_RESTRI
} }
static inline float bw_allpass1_process1(const bw_allpass1_coeffs *BW_RESTRICT coeffs, bw_allpass1_state *BW_RESTRICT state, float x) { static inline float bw_allpass1_process1(const bw_allpass1_coeffs *BW_RESTRICT coeffs, bw_allpass1_state *BW_RESTRICT state, float x) {
const float X = coeffs->X_k * (x - state->lp_z1 - coeffs->t * state->X_z1); const float X = coeffs->X_x * (x - state->lp_z1) - coeffs->X_X_z1 * state->X_z1;
const float lp = x - X; const float lp = x - coeffs->lp_X * X;
state->X_z1 = X; state->X_z1 = X;
state->lp_z1 = lp; state->lp_z1 = lp;
return lp + lp - x; return lp + lp - x;