URL https://opencores.org/ocsvn/scarts/scarts/trunk

Subversion Repositories scarts

[/] [scarts/] [trunk/] [toolchain/] [scarts-newlib/] [newlib-1.17.0/] [newlib/] [libc/] [sys/] [linux/] [cmath/] [s_csinhl.c] - Blame information for rev 9

Details | Compare with Previous | View Log


/* Complex sine hyperbole function for long double.
   Copyright (C) 1997 Free Software Foundation, Inc.
   This file is part of the GNU C Library.
   Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
 
   The GNU C Library is free software; you can redistribute it and/or
   modify it under the terms of the GNU Lesser General Public
   License as published by the Free Software Foundation; either
   version 2.1 of the License, or (at your option) any later version.
 
   The GNU C 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
   Lesser General Public License for more details.
 
   You should have received a copy of the GNU Lesser General Public
   License along with the GNU C Library; if not, write to the Free
   Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
   02111-1307 USA.  */
 
#include <complex.h>
#include <fenv.h>
#include <math.h>
 
#include "math_private.h"
 
 
__complex__ long double
__csinhl (__complex__ long double x)
{
  __complex__ long double retval;
  int negate = signbit (__real__ x);
  int rcls = fpclassify (__real__ x);
  int icls = fpclassify (__imag__ x);
 
  __real__ x = fabsl (__real__ x);
 
  if (rcls >= FP_ZERO)
    {
      /* Real part is finite.  */
      if (icls >= FP_ZERO)
        {
          /* Imaginary part is finite.  */
          long double sinh_val = __ieee754_sinhl (__real__ x);
          long double cosh_val = __ieee754_coshl (__real__ x);
          long double sinix, cosix;
 
          __sincosl (__imag__ x, &sinix, &cosix);
 
          __real__ retval = sinh_val * cosix;
          __imag__ retval = cosh_val * sinix;
 
          if (negate)
            __real__ retval = -__real__ retval;
        }
      else
        {
          if (rcls == FP_ZERO)
            {
              /* Real part is 0.0.  */
              __real__ retval = __copysignl (0.0, negate ? -1.0 : 1.0);
              __imag__ retval = __nanl ("") + __nanl ("");
 
#ifdef FE_INVALID
              if (icls == FP_INFINITE)
                feraiseexcept (FE_INVALID);
#endif
            }
          else
            {
              __real__ retval = __nanl ("");
              __imag__ retval = __nanl ("");
 
#ifdef FE_INVALID
              feraiseexcept (FE_INVALID);
#endif
            }
        }
    }
  else if (rcls == FP_INFINITE)
    {
      /* Real part is infinite.  */
      if (icls == FP_ZERO)
        {
          /* Imaginary part is 0.0.  */
          __real__ retval = negate ? -HUGE_VALL : HUGE_VALL;
          __imag__ retval = __imag__ x;
        }
      else if (icls > FP_ZERO)
        {
          /* Imaginary part is finite.  */
          long double sinix, cosix;
 
          __sincosl (__imag__ x, &sinix, &cosix);
 
          __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 = __nanl ("") + __nanl ("");
 
#ifdef FE_INVALID
          if (icls == FP_INFINITE)
            feraiseexcept (FE_INVALID);
#endif
        }
    }
  else
    {
      __real__ retval = __nanl ("");
      __imag__ retval = __imag__ x == 0.0 ? __imag__ x : __nanl ("");
    }
 
  return retval;
}
weak_alias (__csinhl, csinhl)

Line No.	Rev	Author	Line
1	9	jlechner	`/* Complex sine hyperbole function for long double.`
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__ long double`
29			`__csinhl (__complex__ long double x)`
30			`{`
31			`__complex__ long double retval;`
32			`int negate = signbit (__real__ x);`
33			`int rcls = fpclassify (__real__ x);`
34			`int icls = fpclassify (__imag__ x);`
35
36			`__real__ x = fabsl (__real__ x);`
37
38			`if (rcls >= FP_ZERO)`
39			`{`
40			`/* Real part is finite. */`
41			`if (icls >= FP_ZERO)`
42			`{`
43			`/* Imaginary part is finite. */`
44			`long double sinh_val = __ieee754_sinhl (__real__ x);`
45			`long double cosh_val = __ieee754_coshl (__real__ x);`
46			`long double sinix, cosix;`
47
48			`__sincosl (__imag__ x, &sinix, &cosix);`
49
50			`__real__ retval = sinh_val * cosix;`
51			`__imag__ retval = cosh_val * sinix;`
52
53			`if (negate)`
54			`__real__ retval = -__real__ retval;`
55			`}`
56			`else`
57			`{`
58			`if (rcls == FP_ZERO)`
59			`{`
60			`/* Real part is 0.0. */`
61			`__real__ retval = __copysignl (0.0, negate ? -1.0 : 1.0);`
62			`__imag__ retval = __nanl ("") + __nanl ("");`
63
64			`#ifdef FE_INVALID`
65			`if (icls == FP_INFINITE)`
66			`feraiseexcept (FE_INVALID);`
67			`#endif`
68			`}`
69			`else`
70			`{`
71			`__real__ retval = __nanl ("");`
72			`__imag__ retval = __nanl ("");`
73
74			`#ifdef FE_INVALID`
75			`feraiseexcept (FE_INVALID);`
76			`#endif`
77			`}`
78			`}`
79			`}`
80			`else if (rcls == FP_INFINITE)`
81			`{`
82			`/* Real part is infinite. */`
83			`if (icls == FP_ZERO)`
84			`{`
85			`/* Imaginary part is 0.0. */`
86			`__real__ retval = negate ? -HUGE_VALL : HUGE_VALL;`
87			`__imag__ retval = __imag__ x;`
88			`}`
89			`else if (icls > FP_ZERO)`
90			`{`
91			`/* Imaginary part is finite. */`
92			`long double sinix, cosix;`
93
94			`__sincosl (__imag__ x, &sinix, &cosix);`
95
96			`__real__ retval = __copysignl (HUGE_VALL, cosix);`
97			`__imag__ retval = __copysignl (HUGE_VALL, sinix);`
98
99			`if (negate)`
100			`__real__ retval = -__real__ retval;`
101			`}`
102			`else`
103			`{`
104			`/* The addition raises the invalid exception. */`
105			`__real__ retval = HUGE_VALL;`
106			`__imag__ retval = __nanl ("") + __nanl ("");`
107
108			`#ifdef FE_INVALID`
109			`if (icls == FP_INFINITE)`
110			`feraiseexcept (FE_INVALID);`
111			`#endif`
112			`}`
113			`}`
114			`else`
115			`{`
116			`__real__ retval = __nanl ("");`
117			`__imag__ retval = __imag__ x == 0.0 ? __imag__ x : __nanl ("");`
118			`}`
119
120			`return retval;`
121			`}`
122			`weak_alias (__csinhl, csinhl)`

Browse

Tools

Subversion Repositories scarts

[/] [scarts/] [trunk/] [toolchain/] [scarts-newlib/] [newlib-1.17.0/] [newlib/] [libc/] [sys/] [linux/] [cmath/] [s_csinhl.c] - Blame information for rev 9