qd-0.2: double-double and quad-double number type via libqd
Numeric.QD.Bits.Raw
get_double_expn :: CDouble -> IO CIntSource