xref: /netbsd-src/external/gpl3/gcc.old/dist/libquadmath/math/casinq.c (revision 627f7eb200a4419d89b531d55fccd2ee3ffdcde0)
1*627f7eb2Smrg /* Return arc sine of a complex float type.
2*627f7eb2Smrg    Copyright (C) 1997-2018 Free Software Foundation, Inc.
3*627f7eb2Smrg    This file is part of the GNU C Library.
4*627f7eb2Smrg    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
5*627f7eb2Smrg 
6*627f7eb2Smrg    The GNU C Library is free software; you can redistribute it and/or
7*627f7eb2Smrg    modify it under the terms of the GNU Lesser General Public
8*627f7eb2Smrg    License as published by the Free Software Foundation; either
9*627f7eb2Smrg    version 2.1 of the License, or (at your option) any later version.
10*627f7eb2Smrg 
11*627f7eb2Smrg    The GNU C Library is distributed in the hope that it will be useful,
12*627f7eb2Smrg    but WITHOUT ANY WARRANTY; without even the implied warranty of
13*627f7eb2Smrg    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14*627f7eb2Smrg    Lesser General Public License for more details.
15*627f7eb2Smrg 
16*627f7eb2Smrg    You should have received a copy of the GNU Lesser General Public
17*627f7eb2Smrg    License along with the GNU C Library; if not, see
18*627f7eb2Smrg    <http://www.gnu.org/licenses/>.  */
19*627f7eb2Smrg 
20*627f7eb2Smrg #include "quadmath-imp.h"
21*627f7eb2Smrg 
22*627f7eb2Smrg __complex128
casinq(__complex128 x)23*627f7eb2Smrg casinq (__complex128 x)
24*627f7eb2Smrg {
25*627f7eb2Smrg   __complex128 res;
26*627f7eb2Smrg 
27*627f7eb2Smrg   if (isnanq (__real__ x) || isnanq (__imag__ x))
28*627f7eb2Smrg     {
29*627f7eb2Smrg       if (__real__ x == 0)
30*627f7eb2Smrg 	{
31*627f7eb2Smrg 	  res = x;
32*627f7eb2Smrg 	}
33*627f7eb2Smrg       else if (isinfq (__real__ x) || isinfq (__imag__ x))
34*627f7eb2Smrg 	{
35*627f7eb2Smrg 	  __real__ res = nanq ("");
36*627f7eb2Smrg 	  __imag__ res = copysignq (HUGE_VALQ, __imag__ x);
37*627f7eb2Smrg 	}
38*627f7eb2Smrg       else
39*627f7eb2Smrg 	{
40*627f7eb2Smrg 	  __real__ res = nanq ("");
41*627f7eb2Smrg 	  __imag__ res = nanq ("");
42*627f7eb2Smrg 	}
43*627f7eb2Smrg     }
44*627f7eb2Smrg   else
45*627f7eb2Smrg     {
46*627f7eb2Smrg       __complex128 y;
47*627f7eb2Smrg 
48*627f7eb2Smrg       __real__ y = -__imag__ x;
49*627f7eb2Smrg       __imag__ y = __real__ x;
50*627f7eb2Smrg 
51*627f7eb2Smrg       y = casinhq (y);
52*627f7eb2Smrg 
53*627f7eb2Smrg       __real__ res = __imag__ y;
54*627f7eb2Smrg       __imag__ res = -__real__ y;
55*627f7eb2Smrg     }
56*627f7eb2Smrg 
57*627f7eb2Smrg   return res;
58*627f7eb2Smrg }
59