memcpy(statbuf, work + 40, 20);
}
-static unsigned int rms(const int *data, int f)
+static unsigned int rescale_rms(int rms, int energy)
+{
+ return (rms * energy) >> 10;
+}
+
+static unsigned int rms(const int *data)
{
int x;
unsigned int res = 0x10000;
res = t_sqrt(res);
res >>= (b + 10);
- res = (res * f) >> 10;
return res;
}
// coefficients
if (copynew) {
int_to_int16(decsp, ractx->lpc_coef);
- return rms(ractx->lpc_refl, energy);
+ return rescale_rms(rms(ractx->lpc_refl), energy);
} else {
int_to_int16(decsp, ractx->lpc_coef_old);
- return rms(ractx->lpc_refl_old, energy);
+ return rescale_rms(rms(ractx->lpc_refl_old), energy);
}
} else {
- return rms(work, energy);
+ return rescale_rms(rms(work), energy);
}
}
refl_rms[1] = interp(ractx, block_coefs[1], 1, energy > ractx->old_energy,
t_sqrt(energy*ractx->old_energy) >> 12);
refl_rms[2] = interp(ractx, block_coefs[2], 2, 1, energy);
- refl_rms[3] = rms(ractx->lpc_refl, energy);
+ refl_rms[3] = rescale_rms(rms(ractx->lpc_refl), energy);
int_to_int16(block_coefs[3], ractx->lpc_coef);