Complex numbers

Hey,

Here a first implementation of complex numbers. I tested them with CUDA and they work there, at least.

I am not sure what other functions might be needed on complex numbers (probably arg and norm at least), but I need to see how to implement them so that they work with both backends.

Another thing is speed, I don't if we should maybe mark functions as inline (or __forceinline__ for CUDA?).