- d_phase /= abs(d_phase); // ensure multiplication is rotation
- // FIXME. This is expensive. Maybe workaround using
- // double precision complex???
+ gr_complex z = in * d_phase; // rotate in by phase
+ d_phase *= d_phase_incr; // incr our phase (complex mult == add phases)