libm/math: Update sinh FE_INEXACT code

Replace the code that is there to generate FE_INEXACT with __math_inexact.

Author: futog

newlib/libm/math/s_sinh.c

@@ -33,7 +33,7 @@
 
 #ifdef _NEED_FLOAT64
 
-static const __float64 one = _F_64(1.0), shuge = _F_64(1.0e307);
+static const __float64 one = _F_64(1.0);
 
 __float64
 sinh64(__float64 x)
@@ -56,8 +56,7 @@ sinh64(__float64 x)
     /* |x| in [0,22], return sign(x)*0.5*(E+E/(E+1))) */
     if (ix < 0x40360000) { /* |x|<22 */
         if (ix < 0x3e300000) /* |x|<2**-28 */
-            if (shuge + x > one)
-                return x; /* sinh(tiny) = tiny with inexact */
+            return __math_inexact64(x); /* sinh(tiny) = tiny with inexact */
         t = expm164(fabs64(x));
         if (ix < 0x3ff00000)
             return h * (_F_64(2.0) * t - t * t / (t + one));

newlib/libm/math/sf_sinh.c

@@ -15,7 +15,7 @@
 
 #include "fdlibm.h"
 
-static const float one = 1.0, shuge = 1.0e37;
+static const float one = 1.0;
 
 float
 sinhf(float x)
@@ -36,8 +36,7 @@ sinhf(float x)
     /* |x| in [0,22], return sign(x)*0.5*(E+E/(E+1))) */
     if (ix < 0x41b00000) { /* |x|<22 */
         if (ix < 0x31800000) /* |x|<2**-28 */
-            if (shuge + x > one)
-                return x; /* sinh(tiny) = tiny with inexact */
+            return __math_inexactf(x); /* sinh(tiny) = tiny with inexact */
         t = expm1f(fabsf(x));
         if (ix < 0x3f800000)
             return h * ((float)2.0 * t - t * t / (t + one));