------- Comment #5 from tobias dot burnus at physik dot fu-berlin dot de 2006-09-10 22:43 ------- > Can you try this patch? It required MPFR 2.2.0. Tested on openSUSE/AMD64 with current gfortran trunk + patch and mpfr-2.2.0-9.
It works ok with kind = 4 and kind = 8. However, with kind = 10 I get a NaN for print *, scale (fraction (a), exponent (a)) / a My test case: ---------------- integer, parameter :: k = 10 real(k), parameter :: one = 1.0_k real(k) :: a a = one print *, "This ratio should be 1:" print *, scale (fraction (a), exponent (a)) / a end -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=28276