1/* Complex sine function for float.
2   Copyright (C) 1997 Free Software Foundation, Inc.
3   This file is part of the GNU C Library.
4   Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
5
6   The GNU C Library is free software; you can redistribute it and/or
7   modify it under the terms of the GNU Lesser General Public
8   License as published by the Free Software Foundation; either
9   version 2.1 of the License, or (at your option) any later version.
10
11   The GNU C Library is distributed in the hope that it will be useful,
12   but WITHOUT ANY WARRANTY; without even the implied warranty of
13   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14   Lesser General Public License for more details.
15
16   You should have received a copy of the GNU Lesser General Public
17   License along with the GNU C Library; if not, write to the Free
18   Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
19   02111-1307 USA.  */
20
21#include <complex.h>
22#include <fenv.h>
23#include <math.h>
24
25#include "math_private.h"
26
27
28__complex__ float
29__csinf (__complex__ float x)
30{
31  __complex__ float retval;
32  int negate = signbit (__real__ x);
33  int rcls = fpclassify (__real__ x);
34  int icls = fpclassify (__imag__ x);
35
36  __real__ x = fabsf (__real__ x);
37
38  if (icls >= FP_ZERO)
39    {
40      /* Imaginary part is finite.  */
41      if (rcls >= FP_ZERO)
42	{
43	  /* Real part is finite.  */
44	  float sinh_val = sinhf (__imag__ x);
45	  float cosh_val = coshf (__imag__ x);
46	  float sinix, cosix;
47
48	  sincosf (__real__ x, &sinix, &cosix);
49
50	  __real__ retval = cosh_val * sinix;
51	  __imag__ retval = sinh_val * cosix;
52
53	  if (negate)
54	    __real__ retval = -__real__ retval;
55	}
56      else
57	{
58	  if (icls == FP_ZERO)
59	    {
60	      /* Imaginary part is 0.0.  */
61	      __real__ retval = nanf ("");
62	      __imag__ retval = __imag__ x;
63
64#ifdef FE_INVALID
65	      if (rcls == FP_INFINITE)
66		feraiseexcept (FE_INVALID);
67#endif
68	    }
69	  else
70	    {
71	      __real__ retval = nanf ("");
72	      __imag__ retval = nanf ("");
73
74#ifdef FE_INVALID
75	      feraiseexcept (FE_INVALID);
76#endif
77	    }
78	}
79    }
80  else if (icls == FP_INFINITE)
81    {
82      /* Imaginary part is infinite.  */
83      if (rcls == FP_ZERO)
84	{
85	  /* Real part is 0.0.  */
86	  __real__ retval = copysignf (0.0, negate ? -1.0 : 1.0);
87	  __imag__ retval = __imag__ x;
88	}
89      else if (rcls > FP_ZERO)
90	{
91	  /* Real part is finite.  */
92	  float sinix, cosix;
93
94	  sincosf (__real__ x, &sinix, &cosix);
95
96	  __real__ retval = copysignf (HUGE_VALF, sinix);
97	  __imag__ retval = copysignf (HUGE_VALF, cosix);
98
99	  if (negate)
100	    __real__ retval = -__real__ retval;
101	  if (signbit (__imag__ x))
102	    __imag__ retval = -__imag__ retval;
103	}
104      else
105	{
106	  /* The addition raises the invalid exception.  */
107	  __real__ retval = nanf ("");
108	  __imag__ retval = HUGE_VALF;
109
110#ifdef FE_INVALID
111	  if (rcls == FP_INFINITE)
112	    feraiseexcept (FE_INVALID);
113#endif
114	}
115    }
116  else
117    {
118      if (rcls == FP_ZERO)
119	__real__ retval = copysignf (0.0, negate ? -1.0 : 1.0);
120      else
121	__real__ retval = nanf ("");
122      __imag__ retval = nanf ("");
123    }
124
125  return retval;
126}
127#ifndef __csinf
128weak_alias (__csinf, csinf)
129#endif
130