pub fn residue(x: Cx, y: Cx) -> Cx
x | y: y reduced modulo x, with the complex floor doing the rounding.
x | y