for each sample:
    t0 = x * sample;
    real = t0 + pole * (real - t0);

    t0 = -y * sample;
    imag = t0 + pole * (imag - t0);

    t0 = x * c1 - y * c2;
    y = x * c2 + y * c1;
    x = t0;