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_ccoshl.c] - Blame information for rev 9

Details | Compare with Previous | View Log


/* Complex cosine 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
__ccoshl (__complex__ long double x)
{
  __complex__ long double retval;
  int rcls = fpclassify (__real__ x);
  int icls = fpclassify (__imag__ 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 = cosh_val * cosix;
          __imag__ retval = sinh_val * sinix;
        }
      else
        {
          __imag__ retval = __real__ x == 0.0 ? 0.0 : __nanl ("");
          __real__ retval = __nanl ("") + __nanl ("");
 
#ifdef FE_INVALID
          if (icls == FP_INFINITE)
            feraiseexcept (FE_INVALID);
#endif
        }
    }
  else if (rcls == FP_INFINITE)
    {
      /* Real part is infinite.  */
      if (icls == FP_ZERO)
        {
          /* Imaginary part is 0.0.  */
          __real__ retval = HUGE_VALL;
          __imag__ retval = __imag__ x * __copysignl (1.0, __real__ 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)
                             * __copysignl (1.0, __real__ x));
        }
      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 (__ccoshl, ccoshl)

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

Browse

Tools

Subversion Repositories scarts

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