re PR libfortran/59313 (gfortran.dg/erf_3.F90 FAILs on Solaris/SPARC)

	PR libfortran/59313
	* intrinsics/erfc_scaled.c (erfc_scaled_r16): Also provide for
	quadruple precision long double variant.

From-SVN: r205574

diff --git a/libgfortran/ChangeLog b/libgfortran/ChangeLog
index 5bded877e0f6f7ae84cdde84576a0a1cc39b1257..a428205e5626cb5f94f70e67fe7ab5f26a43f5b0 100644 (file)
--- a/libgfortran/ChangeLog
+++ b/libgfortran/ChangeLog
@@ -1,3 +1,9 @@
+2013-12-01  Uros Bizjak  <ubizjak@gmail.com>
+
+       PR libfortran/59313
+       * intrinsics/erfc_scaled.c (erfc_scaled_r16): Also provide for
+       quadruple precision long double variant.
+
 2013-11-20  Francois-Xavier Coudert  <fxcoudert@gcc.gnu.org>
 
        * intrinsics/erfc_scaled.c (erfc_scaled_r16): Don't define if

diff --git a/libgfortran/intrinsics/erfc_scaled.c b/libgfortran/intrinsics/erfc_scaled.c
index f69e1fc7357d7dd58a847d14d98dbd851b9f9b7b..0498508a5bbcf7f8e909ed026e0d00c88be3b2ca 100644 (file)
--- a/libgfortran/intrinsics/erfc_scaled.c
+++ b/libgfortran/intrinsics/erfc_scaled.c
@@ -45,17 +45,34 @@ see the files COPYING3 and COPYING.RUNTIME respectively.  If not, see
 #include "erfc_scaled_inc.c"
 #endif
 
-#if defined(HAVE_GFC_REAL_16) && defined(GFC_REAL_16_IS_LONG_DOUBLE)
-#undef KIND
-#define KIND 16
-#include "erfc_scaled_inc.c"
-#endif
+#ifdef HAVE_GFC_REAL_16
 
+/* For quadruple-precision, netlib's implementation is
+   not accurate enough.  We provide another one.  */
 
 #ifdef GFC_REAL_16_IS_FLOAT128
 
-/* For quadruple-precision (__float128), netlib's implementation is
-   not accurate enough.  We provide another one.  */
+# define _THRESH -106.566990228185312813205074546585730Q
+# define _M_2_SQRTPI M_2_SQRTPIq
+# define _INF __builtin_infq()
+# define _ERFC(x) erfcq(x)
+# define _EXP(x) expq(x)
+
+#else
+
+# define _THRESH -106.566990228185312813205074546585730L
+# define _M_2_SQRTPI M_2_SQRTPIl
+# define _INF __builtin_infl()
+# ifdef HAVE_ERFCL
+#  define _ERFC(x) erfcl(x)
+# endif
+# ifdef HAVE_EXPL
+#  define _EXP(x) expl(x)
+# endif
+
+#endif
+
+#if defined(_ERFC) && defined(_EXP)
 
 extern GFC_REAL_16 erfc_scaled_r16 (GFC_REAL_16);
 export_proto(erfc_scaled_r16);
@@ -63,15 +80,15 @@ export_proto(erfc_scaled_r16);
 GFC_REAL_16
 erfc_scaled_r16 (GFC_REAL_16 x)
 {
-  if (x < -106.566990228185312813205074546585730Q)
+  if (x < _THRESH)
     {
-      return __builtin_infq();
+      return _INF;
     }
   if (x < 12)
     {
       /* Compute directly as ERFC_SCALED(x) = ERFC(x) * EXP(X**2).
         This is not perfect, but much better than netlib.  */
-      return erfcq(x) * expq(x * x);
+      return _ERFC(x) * _EXP(x * x);
     }
   else
     {
@@ -97,9 +114,10 @@ erfc_scaled_r16 (GFC_REAL_16 x)
          n++;
        }
 
-      return (1 + sum) / x * (M_2_SQRTPIq / 2);
+      return (1 + sum) / x * (_M_2_SQRTPI / 2);
     }
 }
 
 #endif
 
+#endif

diff --git a/libgfortran/intrinsics/erfc_scaled_inc.c b/libgfortran/intrinsics/erfc_scaled_inc.c
index 107d91a6c9db2c2c03cfcab637901ba349bcf3ae..b643edf43a4802ab0a39d62c4ac21ce1919f602b 100644 (file)
--- a/libgfortran/intrinsics/erfc_scaled_inc.c
+++ b/libgfortran/intrinsics/erfc_scaled_inc.c
@@ -39,7 +39,7 @@ see the files COPYING3 and COPYING.RUNTIME respectively.  If not, see
 # define EXP(x) exp(x)
 # define TRUNC(x) trunc(x)
 
-#elif (KIND == 10) || (KIND == 16 && defined(GFC_REAL_16_IS_LONG_DOUBLE))
+#elif (KIND == 10)
 
 # ifdef HAVE_EXPL
 #  define EXP(x) expl(x)