frac32.bipolar density
frac32buffer.bipolar impulses
frac32.u.map density
uint32_t r;
int32_t y;
int32_t density;
y = 1 << 27;
r = 137 + GenerateRandomNumber();
density = param_density + inlet_density;
density = ___SMMLA(density << 2, density << 3, density >> 4); // x^2 + x/16
density = ___SMMLA(density << 2, density << 3,
density >> 2); // (x^2 + x/16)^2 + (x^2 + x/16)/4
density >>= 6; // q31
// a
r *= 69069UL;
if ((r >> 5) < density)
y = -y;
outlet_out = y;