From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 30228 invoked by alias); 10 Sep 2006 22:43:28 -0000 Received: (qmail 30201 invoked by uid 48); 10 Sep 2006 22:43:21 -0000 Date: Sun, 10 Sep 2006 22:43:00 -0000 Message-ID: <20060910224321.30200.qmail@sourceware.org> X-Bugzilla-Reason: CC References: Subject: [Bug fortran/28276] EXPONENT() broken for real constants In-Reply-To: Reply-To: gcc-bugzilla@gcc.gnu.org To: gcc-bugs@gcc.gnu.org From: "tobias dot burnus at physik dot fu-berlin dot de" Mailing-List: contact gcc-bugs-help@gcc.gnu.org; run by ezmlm Precedence: bulk List-Archive: List-Post: List-Help: Sender: gcc-bugs-owner@gcc.gnu.org X-SW-Source: 2006-09/txt/msg00849.txt.bz2 List-Id: ------- 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