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

Subversion Repositories openrisc

[/] [openrisc/] [trunk/] [gnu-dev/] [or1k-gcc/] [libquadmath/] [math/] [cacoshq.c] - Blame information for rev 740

Details | Compare with Previous | View Log


/* Return arc hyperbole cosine for __float128 value.
   Copyright (C) 1997, 1998, 2006 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 "quadmath-imp.h"
 
 
__complex128
cacoshq (__complex128 x)
{
  __complex128 res;
  int rcls = fpclassifyq (__real__ x);
  int icls = fpclassifyq (__imag__ x);
 
  if (rcls <= QUADFP_INFINITE || icls <= QUADFP_INFINITE)
    {
      if (icls == QUADFP_INFINITE)
        {
          __real__ res = HUGE_VALQ;
 
          if (rcls == QUADFP_NAN)
            __imag__ res = nanq ("");
          else
            __imag__ res = copysignq ((rcls == QUADFP_INFINITE
                                       ? (__real__ x < 0.0
                                          ? M_PIq - M_PI_4q : M_PI_4q)
                                       : M_PI_2q), __imag__ x);
        }
      else if (rcls == QUADFP_INFINITE)
        {
          __real__ res = HUGE_VALQ;
 
          if (icls >= QUADFP_ZERO)
            __imag__ res = copysignq (signbitq (__real__ x) ? M_PIq : 0.0,
                                      __imag__ x);
          else
            __imag__ res = nanq ("");
        }
      else
        {
          __real__ res = nanq ("");
          __imag__ res = nanq ("");
        }
    }
  else if (rcls == QUADFP_ZERO && icls == QUADFP_ZERO)
    {
      __real__ res = 0.0;
      __imag__ res = copysignq (M_PI_2q, __imag__ x);
    }
  else
    {
      __complex128 y;
 
      __real__ y = (__real__ x - __imag__ x) * (__real__ x + __imag__ x) - 1.0;
      __imag__ y = 2.0 * __real__ x * __imag__ x;
 
      y = csqrtq (y);
 
      if (__real__ x < 0.0)
        y = -y;
 
      __real__ y += __real__ x;
      __imag__ y += __imag__ x;
 
      res = clogq (y);
 
      /* We have to use the positive branch.  */
      if (__real__ res < 0.0)
        res = -res;
    }
 
  return res;
}

Line No.	Rev	Author	Line
1	740	jeremybenn	`/* Return arc hyperbole cosine for __float128 value.`
2			`Copyright (C) 1997, 1998, 2006 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 "quadmath-imp.h"`
22
23
24			`__complex128`
25			`cacoshq (__complex128 x)`
26			`{`
27			`__complex128 res;`
28			`int rcls = fpclassifyq (__real__ x);`
29			`int icls = fpclassifyq (__imag__ x);`
30
31			`if (rcls <= QUADFP_INFINITE \|\| icls <= QUADFP_INFINITE)`
32			`{`
33			`if (icls == QUADFP_INFINITE)`
34			`{`
35			`__real__ res = HUGE_VALQ;`
36
37			`if (rcls == QUADFP_NAN)`
38			`__imag__ res = nanq ("");`
39			`else`
40			`__imag__ res = copysignq ((rcls == QUADFP_INFINITE`
41			`? (__real__ x < 0.0`
42			`? M_PIq - M_PI_4q : M_PI_4q)`
43			`: M_PI_2q), __imag__ x);`
44			`}`
45			`else if (rcls == QUADFP_INFINITE)`
46			`{`
47			`__real__ res = HUGE_VALQ;`
48
49			`if (icls >= QUADFP_ZERO)`
50			`__imag__ res = copysignq (signbitq (__real__ x) ? M_PIq : 0.0,`
51			`__imag__ x);`
52			`else`
53			`__imag__ res = nanq ("");`
54			`}`
55			`else`
56			`{`
57			`__real__ res = nanq ("");`
58			`__imag__ res = nanq ("");`
59			`}`
60			`}`
61			`else if (rcls == QUADFP_ZERO && icls == QUADFP_ZERO)`
62			`{`
63			`__real__ res = 0.0;`
64			`__imag__ res = copysignq (M_PI_2q, __imag__ x);`
65			`}`
66			`else`
67			`{`
68			`__complex128 y;`
69
70			`__real__ y = (__real__ x - __imag__ x) * (__real__ x + __imag__ x) - 1.0;`
71			`__imag__ y = 2.0 * __real__ x * __imag__ x;`
72
73			`y = csqrtq (y);`
74
75			`if (__real__ x < 0.0)`
76			`y = -y;`
77
78			`__real__ y += __real__ x;`
79			`__imag__ y += __imag__ x;`
80
81			`res = clogq (y);`
82
83			`/* We have to use the positive branch. */`
84			`if (__real__ res < 0.0)`
85			`res = -res;`
86			`}`
87
88			`return res;`
89			`}`

Browse

Tools

Subversion Repositories openrisc

[/] [openrisc/] [trunk/] [gnu-dev/] [or1k-gcc/] [libquadmath/] [math/] [cacoshq.c] - Blame information for rev 740