1  /* Compute sine and cosine of argument.
       2     Copyright (C) 1997-2023 Free Software Foundation, Inc.
       3     This file is part of the GNU C Library.
       4  
       5     The GNU C Library is free software; you can redistribute it and/or
       6     modify it under the terms of the GNU Lesser General Public
       7     License as published by the Free Software Foundation; either
       8     version 2.1 of the License, or (at your option) any later version.
       9  
      10     The GNU C Library is distributed in the hope that it will be useful,
      11     but WITHOUT ANY WARRANTY; without even the implied warranty of
      12     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
      13     Lesser General Public License for more details.
      14  
      15     You should have received a copy of the GNU Lesser General Public
      16     License along with the GNU C Library; if not, see
      17     <https://www.gnu.org/licenses/>.  */
      18  
      19  #include <errno.h>
      20  #include <math.h>
      21  
      22  #include <math_private.h>
      23  #include <libm-alias-ldouble.h>
      24  
      25  void
      26  __sincosl (_Float128 x, _Float128 *sinx, _Float128 *cosx)
      27  {
      28    int64_t ix;
      29  
      30    /* High word of x. */
      31    GET_LDOUBLE_MSW64 (ix, x);
      32  
      33    /* |x| ~< pi/4 */
      34    ix &= 0x7fffffffffffffffLL;
      35    if (ix <= 0x3ffe921fb54442d1LL)
      36      __kernel_sincosl (x, 0, sinx, cosx, 0);
      37    else if (ix >= 0x7fff000000000000LL)
      38      {
      39        /* sin(Inf or NaN) is NaN */
      40        *sinx = *cosx = x - x;
      41        if (isinf (x))
      42  	__set_errno (EDOM);
      43      }
      44    else
      45      {
      46        /* Argument reduction needed.  */
      47        _Float128 y[2];
      48        int n;
      49  
      50        n = __ieee754_rem_pio2l (x, y);
      51        switch (n & 3)
      52  	{
      53  	case 0:
      54  	  __kernel_sincosl (y[0], y[1], sinx, cosx, 1);
      55  	  break;
      56  	case 1:
      57  	  __kernel_sincosl (y[0], y[1], cosx, sinx, 1);
      58  	  *cosx = -*cosx;
      59  	  break;
      60  	case 2:
      61  	  __kernel_sincosl (y[0], y[1], sinx, cosx, 1);
      62  	  *sinx = -*sinx;
      63  	  *cosx = -*cosx;
      64  	  break;
      65  	default:
      66  	  __kernel_sincosl (y[0], y[1], cosx, sinx, 1);
      67  	  *sinx = -*sinx;
      68  	  break;
      69  	}
      70      }
      71  }
      72  libm_alias_ldouble (__sincos, sincos)