Only down side in the moment is that complex<double> and complex<float> must never be mixed in a kernel (real scalars of the other type are mostly ok due to manually implemented templates).
Should work on CPU and GPU.
Another thing that this PR changes is that also the headers attribute of SymPy Expression is checked to determine necessary headers.