frac32.bipolar pitch_slave
frac32.bipolar pitch master
frac32buffer.bipolar out
frac32.s.map.pitch pitch slave
frac32.s.map.pitch pitch master
frac32.s.map level0
frac32.s.map level1
frac32.s.map level2
frac32.s.map level3
frac32.s.map level4
frac32.s.map level5
float x, y, z0, z_1, z_2, p, dp, _dp, pM, dpM, _dpM;
float w0, w_1, w_2;
int cpt;
float seq[6];
cpt = 0;
x = y = z_2 = z_1 = z0 = w_2 = w_1 = w0 = 0;
p = pM = 0;
int32_t idp;
MTOFEXTENDED(param_pitch_space_slave + inlet_pitch_space_slave, idp);
dp = 6.0f * (idp * (0.25f / (1 << 30)));
_dp = 1 / dp;
if (dp > 1)
dp = 1;
// _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
MTOFEXTENDED(param_pitch_space_master + inlet_pitch_space_master, idp);
dpM = (idp * (0.25f / (1 << 30)));
_dpM = 1 / dpM;
if (dpM > 1)
dpM = 1;
// _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
float dc = seq[0] = (float)param_level0;
dc += seq[1] = (float)param_level1;
dc += seq[2] = (float)param_level2;
dc += seq[3] = (float)param_level3;
dc += seq[4] = (float)param_level4;
dc += seq[5] = (float)param_level5;
dc *= 1.0f / 6;
for (int i = 0; i < 6; i++)
seq[i] -= dc;
// phase increment
pM += dpM;
// _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
// synchro ?
if (pM > 1) { // phase above 1
cpt = 0; // reset
pM -= 1; // reset phase
float alpha = pM * _dpM; // subsample time since the transition
p = alpha * dp;
//_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
// integrator evolution before the transition
float _alpha = 1 - alpha;
float alpha_x = _alpha * x;
z0 += _alpha * (y + 0.5f * alpha_x);
y += alpha_x;
//_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
// transition
x = seq[cpt];
//_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
// integrator evolution after the transition
alpha_x = alpha * x;
z0 += alpha * (y + 0.5f * alpha_x);
y += alpha_x;
} else {
// phase increment
p += dp;
// _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
// transition ?
if (p > 1) { // phase above 1
cpt++; // next step in sequence
if (cpt >= 6) // above 6 => wrap to 0
cpt = 0;
p -= 1; // reset phase
float alpha = p * _dp; // subsample time since the transition
//_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
// integrator evolution before the transition
float _alpha = 1 - alpha;
float alpha_x = _alpha * x;
z0 += _alpha * (y + 0.5f * alpha_x);
y += alpha_x;
//_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
// transition
x = seq[cpt];
//_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
// integrator evolution after the transition
alpha_x = alpha * x;
z0 += alpha * (y + 0.5f * alpha_x);
y += alpha_x;
} else {
//_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
// integrator evolution without transition
z0 += y + 0.5f * x;
y += x;
}
}
// _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
// integrator leaks
y *= 0.995f;
z0 *= 0.995f;
// _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
// 2nd order differenciation
w0 = z0 + z_2 - 2 * z_1;
// _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
// high freq compensation
outlet_out = (int32_t)(0.9f * w_1 - 0.1f * (w_2 + w0));
// _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
// previous values
z_2 = z_1;
z_1 = z0;
w_2 = w_1;
w_1 = w0;