re PR libfortran/34699 (FAIL: gfortran.fortran-torture/execute/intrinsic_set_exponent...

	PR libfortran/34699
	Backport:
	2007-09-21  Francois-Xavier Coudert  <fxcoudert@gcc.gnu.org>

	PR libfortran/26253
	* intrinsics/c99_functions.c (scalbn): Use ldexp if appopriate.

From-SVN: r131699

diff --git a/libgfortran/ChangeLog b/libgfortran/ChangeLog
index a71d215b8ce7a4abcedd65c35e2f03c97c53a336..b60bfe33484c514a5f116038c00f5bc28f0b17a1 100644 (file)
--- a/libgfortran/ChangeLog
+++ b/libgfortran/ChangeLog
@@ -1,3 +1,12 @@
+2008-01-21  John David Anglin  <dave.anglin@nrc-cnrc.gc.ca>
+
+       PR libfortran/34699
+       Backport:
+       2007-09-21  Francois-Xavier Coudert  <fxcoudert@gcc.gnu.org>
+
+       PR libfortran/26253
+       * intrinsics/c99_functions.c (scalbn): Use ldexp if appopriate.
+
 2007-05-23  Jerry DeLisle  <jvdelisle@gcc.gnu.org>
 
        PR libfortran/31964

diff --git a/libgfortran/intrinsics/c99_functions.c b/libgfortran/intrinsics/c99_functions.c
index 3eabd0dc5a7ebdcb450df1cd709a9f687b2c6801..40e0ba763465b0b0249cebe0a15cc90e7bb0faf9 100644 (file)
--- a/libgfortran/intrinsics/c99_functions.c
+++ b/libgfortran/intrinsics/c99_functions.c
@@ -226,7 +226,11 @@ log10f(float x)
 double
 scalbn(double x, int y)
 {
+#if (FLT_RADIX == 2) && defined(HAVE_LDEXP)
+  return ldexp (x, y);
+#else
   return x * pow(FLT_RADIX, y);
+#endif
 }
 #endif