qd is a port of David H. Bailey double-double and quad-double (not yet implemented) extended precision floating point arithmetic.
qd