Back to index

glibc  2.9
Functions
s_sincosl.c File Reference
#include <math.h>
#include "math_private.h"

Go to the source code of this file.

Functions

void __sincosl (long double x, long double *sinx, long double *cosx)

Function Documentation

void __sincosl ( long double  x,
long double *  sinx,
long double *  cosx 
)

Definition at line 27 of file s_sincosl.c.

{
  int32_t se, i0, i1;

  /* High word of x. */
  GET_LDOUBLE_WORDS (se, i0, i1, x);

  /* |x| ~< pi/4 */
  se &= 0x7fff;
  if (se < 0x3ffe || (se == 0x3ffe && i0 <= 0xc90fdaa2))
    {
      *sinx = __kernel_sinl (x, 0.0, 0);
      *cosx = __kernel_cosl (x, 0.0);
    }
  else if (se == 0x7fff)
    {
      /* sin(Inf or NaN) is NaN */
      *sinx = *cosx = x - x;
    }
  else
    {
      /* Argument reduction needed.  */
      long double y[2];
      int n;

      n = __ieee754_rem_pio2l (x, y);
      switch (n & 3)
       {
       case 0:
         *sinx = __kernel_sinl (y[0], y[1], 1);
         *cosx = __kernel_cosl (y[0], y[1]);
         break;
       case 1:
         *sinx = __kernel_cosl (y[0], y[1]);
         *cosx = -__kernel_sinl (y[0], y[1], 1);
         break;
       case 2:
         *sinx = -__kernel_sinl (y[0], y[1], 1);
         *cosx = -__kernel_cosl (y[0], y[1]);
         break;
       default:
         *sinx = -__kernel_cosl (y[0], y[1]);
         *cosx = __kernel_sinl (y[0], y[1], 1);
         break;
       }
    }
}

Here is the call graph for this function:

Here is the caller graph for this function: