Sorry to have only just noticed this; I see you're using PDL::Complex, which is now officially deprecated. If you delete that line, then with a recent (2.047+) PDL "i" will return a "cdouble" (a "native complex" type) with the obvious value and all the type stuff will do the right thing.