Forgot 0.5 factor.

diff --git a/sysdeps/libm-ieee754/s_csinh.c b/sysdeps/libm-ieee754/s_csinh.c
index e377140..aab041b 100644 (file)
--- a/sysdeps/libm-ieee754/s_csinh.c
+++ b/sysdeps/libm-ieee754/s_csinh.c
@@ -37,8 +37,8 @@ __csinh (__complex__ double x)
          double exp_val = __exp (__real__ x);
          double rec_exp_val = 1.0 / exp_val;
 
-         __real__ retval = (exp_val - rec_exp_val) * __cos (__imag__ x);
-         __imag__ retval = (exp_val - rec_exp_val) * __sin (__imag__ x);
+         __real__ retval = 0.5 * (exp_val - rec_exp_val) * __cos (__imag__ x);
+         __imag__ retval = 0.5 * (exp_val - rec_exp_val) * __sin (__imag__ x);
 
          if (negate)
            __real__ retval = -__real__ retval;

diff --git a/sysdeps/libm-ieee754/s_csinhf.c b/sysdeps/libm-ieee754/s_csinhf.c
index 42d56ff..204bbfe 100644 (file)
--- a/sysdeps/libm-ieee754/s_csinhf.c
+++ b/sysdeps/libm-ieee754/s_csinhf.c
@@ -37,8 +37,10 @@ __csinhf (__complex__ float x)
          float exp_val = __expf (__real__ x);
          float rec_exp_val = 1.0 / exp_val;
 
-         __real__ retval = (exp_val - rec_exp_val) * __cosf (__imag__ x);
-         __imag__ retval = (exp_val - rec_exp_val) * __sinf (__imag__ x);
+         __real__ retval = (0.5 * (exp_val - rec_exp_val)
+                            * __cosf (__imag__ x));
+         __imag__ retval = (0.5 * (exp_val - rec_exp_val)
+                            * __sinf (__imag__ x));
 
          if (negate)
            __real__ retval = -__real__ retval;

diff --git a/sysdeps/libm-ieee754/s_csinhl.c b/sysdeps/libm-ieee754/s_csinhl.c
index ea6839b..e403dd4 100644 (file)
--- a/sysdeps/libm-ieee754/s_csinhl.c
+++ b/sysdeps/libm-ieee754/s_csinhl.c
@@ -37,8 +37,10 @@ __csinhl (__complex__ long double x)
          long double exp_val = __expl (__real__ x);
          long double rec_exp_val = 1.0 / exp_val;
 
-         __real__ retval = (exp_val - rec_exp_val) * __cosl (__imag__ x);
-         __imag__ retval = (exp_val - rec_exp_val) * __sinl (__imag__ x);
+         __real__ retval = (0.5 * (exp_val - rec_exp_val)
+                            * __cosl (__imag__ x));
+         __imag__ retval = (0.5 * (exp_val - rec_exp_val)
+                            * __sinl (__imag__ x));
 
          if (negate)
            __real__ retval = -__real__ retval;