http://gcc.gnu.org/bugzilla/show_bug.cgi?id=52413
--- Comment #7 from Francois-Xavier Coudert <fxcoudert at gcc dot gnu.org> --- I forgot in the last comment to say: handling of sign for non-zero cases, in old MPFR versions, is done by this line which was missing in the existing code: + mpfr_copysign (result->value.real, result->value.real, + x->value.real, GFC_RND_MODE);