1 files changed, 110 insertions, 0 deletions
diff --git a/libstdc++-v3/libmath/csinh.c b/libstdc++-v3/libmath/csinh.c
new file mode 100644
index 00000000000..9da924b1c66
--- /dev/null
+++ b/libstdc++-v3/libmath/csinh.c
@@ -0,0 +1,110 @@
+/* Complex sine hyperbole function for double. */
+
+/* Copyright (C) 1997-1999 Free Software Foundation, Inc.
+
+   This file is part of the GNU ISO C++ Library.  This library is free
+   software; you can redistribute it and/or modify it under the
+   terms of the GNU General Public License as published by the
+   Free Software Foundation; either version 2, or (at your option)
+   any later version.
+
+   This library is distributed in the hope that it will be useful,
+   but WITHOUT ANY WARRANTY; without even the implied warranty of
+   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+   GNU General Public License for more details.
+
+   You should have received a copy of the GNU General Public License along
+   with this library; see the file COPYING.  If not, write to the Free
+   Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307,
+   USA.
+
+   As a special exception, you may use this file as part of a free software
+   library without restriction.  Specifically, if other files instantiate
+   templates or use macros or inline functions from this file, or you compile
+   this file and link it with other files to produce an executable, this
+   file does not by itself cause the resulting executable to be covered by
+   the GNU General Public License.  This exception does not however
+   invalidate any other reasons why the executable file might be covered by
+   the GNU General Public License.  */
+
+
+#include <math.h>
+#include "mathconf.h"
+
+
+__complex__ double
+csinh (__complex__ double x)
+{
+  __complex__ double retval;
+  int negate = signbit (__real__ x);
+
+  __real__ x = fabs (__real__ x);
+
+  if (FINITE_P (__real__ x))
+    {
+      /* Real part is finite.  */
+      if (FINITE_P (__imag__ x))
+	{
+	  /* Imaginary part is finite.  */
+	  double sinh_val = sinh (__real__ x);
+	  double cosh_val = cosh (__real__ x);
+	  double sinix = sin (__imag__ x);
+	  double cosix = cos (__imag__ x);
+
+	  __real__ retval = sinh_val * cosix;
+	  __imag__ retval = cosh_val * sinix;
+
+	  if (negate)
+	    __real__ retval = -__real__ retval;
+	}
+      else
+	{
+	  if (__real__ x == 0.0)
+	    {
+	      /* Real part is 0.0.  */
+	      __real__ retval = copysign (0.0, negate ? -1.0 : 1.0);
+	      __imag__ retval = NAN + NAN;
+	    }
+	  else
+	    {
+	      __real__ retval = NAN;
+	      __imag__ retval = NAN;
+	    }
+	}
+    }
+  else if (INFINITE_P (__real__ x))
+    {
+      /* Real part is infinite.  */
+      if (__imag__ x == 0.0)
+	{
+	  /* Imaginary part is 0.0.  */
+	  __real__ retval = negate ? -HUGE_VAL : HUGE_VAL;
+	  __imag__ retval = __imag__ x;
+	}
+      else if (FINITE_P (__imag__ x))
+	{
+	  /* Imaginary part is finite.  */
+	  double sinix = sin (__imag__ x);
+	  double cosix = cos (__imag__ x);
+
+	  __real__ retval = copysign (HUGE_VAL, cosix);
+	  __imag__ retval = copysign (HUGE_VAL, sinix);
+
+	  if (negate)
+	    __real__ retval = -__real__ retval;
+	}
+      else
+	{
+	  /* The addition raises the invalid exception.  */
+	  __real__ retval = HUGE_VAL;
+	  __imag__ retval = NAN + NAN;
+	}
+    }
+  else
+    {
+      __real__ retval = NAN;
+      __imag__ retval = __imag__ x == 0.0 ? __imag__ x : NAN;
+    }
+
+  return retval;
+}