#define RE(x)   (__real__ (x)) // real part
#define IM(x)   (__imag__ (x)) // image part

#define CC(x,y) ((x) + (y)*1i) // complex constant