if (!b)
b = -2;
+ b = 0x1000000 / b;
for (j=0; j <= i; j++)
- bp1[j] = ((bp2[j] - ((refl[i+1] * bp2[i-j]) >> 12)) * (0x1000000 / b)) >> 12;
+ bp1[j] = ((bp2[j] - ((refl[i+1] * bp2[i-j]) >> 12)) * b) >> 12;
if ((unsigned) bp1[i] + 0x1000 > 0x1fff)
return 1;