function real,x
type = size(x,/type)
format = intarr(20)
format[[5,9]] = 1
return, format[type] ? double(x) : float(x)
end
 


function abs2,x
return, real(x)^2 + imaginary(x)^2
end