Function principal_phase 

fn principal_phase(twice_n: i32, z: Complex) -> Complex

z^n / (-z)^n on the principal branches, for z with a signed imaginary part: (-1)^n, or for half-integer n (-1)^floor(n) times -i on and below the real axis and i above it.