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

Subversion Repositories openrisc

[/] [openrisc/] [trunk/] [gnu-old/] [newlib-1.17.0/] [newlib/] [libc/] [sys/] [linux/] [cmath/] [s_csinf.c] - Blame information for rev 158

Go to most recent revision | Details | Compare with Previous | View Log


/* Complex sine function for float.
   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__ float
__csinf (__complex__ float x)
{
  __complex__ float retval;
  int negate = signbit (__real__ x);
  int rcls = fpclassify (__real__ x);
  int icls = fpclassify (__imag__ x);
 
  __real__ x = fabsf (__real__ x);
 
  if (icls >= FP_ZERO)
    {
      /* Imaginary part is finite.  */
      if (rcls >= FP_ZERO)
        {
          /* Real part is finite.  */
          float sinh_val = __ieee754_sinhf (__imag__ x);
          float cosh_val = __ieee754_coshf (__imag__ x);
          float sinix, cosix;
 
          __sincosf (__real__ x, &sinix, &cosix);
 
          __real__ retval = cosh_val * sinix;
          __imag__ retval = sinh_val * cosix;
 
          if (negate)
            __real__ retval = -__real__ retval;
        }
      else
        {
          if (icls == FP_ZERO)
            {
              /* Imaginary part is 0.0.  */
              __real__ retval = __nanf ("");
              __imag__ retval = __imag__ x;
 
#ifdef FE_INVALID
              if (rcls == FP_INFINITE)
                feraiseexcept (FE_INVALID);
#endif
            }
          else
            {
              __real__ retval = __nanf ("");
              __imag__ retval = __nanf ("");
 
#ifdef FE_INVALID
              feraiseexcept (FE_INVALID);
#endif
            }
        }
    }
  else if (icls == FP_INFINITE)
    {
      /* Imaginary part is infinite.  */
      if (rcls == FP_ZERO)
        {
          /* Real part is 0.0.  */
          __real__ retval = __copysignf (0.0, negate ? -1.0 : 1.0);
          __imag__ retval = __imag__ x;
        }
      else if (rcls > FP_ZERO)
        {
          /* Real part is finite.  */
          float sinix, cosix;
 
          __sincosf (__real__ x, &sinix, &cosix);
 
          __real__ retval = __copysignf (HUGE_VALF, sinix);
          __imag__ retval = __copysignf (HUGE_VALF, cosix);
 
          if (negate)
            __real__ retval = -__real__ retval;
          if (signbit (__imag__ x))
            __imag__ retval = -__imag__ retval;
        }
      else
        {
          /* The addition raises the invalid exception.  */
          __real__ retval = __nanf ("");
          __imag__ retval = HUGE_VALF;
 
#ifdef FE_INVALID
          if (rcls == FP_INFINITE)
            feraiseexcept (FE_INVALID);
#endif
        }
    }
  else
    {
      if (rcls == FP_ZERO)
        __real__ retval = __copysignf (0.0, negate ? -1.0 : 1.0);
      else
        __real__ retval = __nanf ("");
      __imag__ retval = __nanf ("");
    }
 
  return retval;
}
weak_alias (__csinf, csinf)

Line No.	Rev	Author	Line
1	148	jeremybenn	`/* 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 = __ieee754_sinhf (__imag__ x);`
45			`float cosh_val = __ieee754_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			`weak_alias (__csinf, csinf)`

Browse

Tools

Subversion Repositories openrisc

[/] [openrisc/] [trunk/] [gnu-old/] [newlib-1.17.0/] [newlib/] [libc/] [sys/] [linux/] [cmath/] [s_csinf.c] - Blame information for rev 158