1 files changed, 110 insertions, 0 deletions
diff --git a/libstdc++-v3/libmath/csinhl.c b/libstdc++-v3/libmath/csinhl.c
new file mode 100644
index 00000000000..efd45a24323
--- /dev/null
+++ b/libstdc++-v3/libmath/csinhl.c
@@ -0,0 +1,110 @@
+/* Complex sine hyperbole function for long 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__ long double
+csinhl (__complex__ long double x)
+{
+  __complex__ long double retval;
+  int negate = signbit (__real__ x);
+
+  __real__ x = fabsl (__real__ x);
+
+  if (FINITEL_P (__real__ x))
+    {
+      /* Real part is finite.  */
+      if (FINITEL_P (__imag__ x))
+	{
+	  /* Imaginary part is finite.  */
+	  long double sinh_val = sinhl (__real__ x);
+	  long double cosh_val = coshl (__real__ x);
+	  long double sinix = sinl (__imag__ x);
+	  long double cosix = cosl (__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 = copysignl (0.0, negate ? -1.0 : 1.0);
+	      __imag__ retval = NAN + NAN;
+	    }
+	  else
+	    {
+	      __real__ retval = NAN;
+	      __imag__ retval = NAN;
+	    }
+	}
+    }
+  else if (INFINITEL_P (__real__ x))
+    {
+      /* Real part is infinite.  */
+      if (__imag__ x == 0.0)
+	{
+	  /* Imaginary part is 0.0.  */
+	  __real__ retval = negate ? -HUGE_VALL : HUGE_VALL;
+	  __imag__ retval = __imag__ x;
+	}
+      else if (FINITEL_P (__imag__ x))
+	{
+	  /* Imaginary part is finite.  */
+	  long double sinix = sinl (__imag__ x);
+	  long double cosix = cosl (__imag__ x);
+
+	  __real__ retval = copysignl (HUGE_VALL, cosix);
+	  __imag__ retval = copysignl (HUGE_VALL, sinix);
+
+	  if (negate)
+	    __real__ retval = -__real__ retval;
+	}
+      else
+	{
+	  /* The addition raises the invalid exception.  */
+	  __real__ retval = HUGE_VALL;
+	  __imag__ retval = NAN + NAN;
+	}
+    }
+  else
+    {
+      __real__ retval = NAN;
+      __imag__ retval = __imag__ x == 0.0 ? __imag__ x : NAN;
+    }
+
+  return retval;
+}