Some internal underflows can cause substantial inaccuracy in the results of csqrt. Specifically, on x86_64, csqrt (0x1p-1074i) returns zero and csqrt (0x1p-500 + 0x1p-1074i) returns an imaginary part of zero (should be 4.469444793151709302716387622440056066334e-249).
Fixed for 2.17 by: commit cdfe2c5eb3703ed964cbfdb6906b21ace2956385 Author: Joseph Myers <joseph@codesourcery.com> Date: Thu Jul 5 11:02:13 2012 +0000 Fix csqrt underflow (bugs 14157, 14331).