Back to index

glibc  2.9
Defines | Functions
s_cexp.c File Reference
#include <complex.h>
#include <math.h>
#include "mathimpl.h"
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Defines

#define SUFF
#define float_type   double
#define CONCATX(a, b)   __CONCAT(a,b)
#define s(name)   CONCATX(name,SUFF)
#define m81(func)   __m81_u(s(func))
#define weak_aliasx(a, b)   weak_alias(a,b)

Functions

__complex__ float_type s() __cexp (__complex__ float_type x)

Define Documentation

#define CONCATX (   a,
  b 
)    __CONCAT(a,b)

Definition at line 32 of file s_cexp.c.

#define float_type   double

Definition at line 29 of file s_cexp.c.

#define m81 (   func)    __m81_u(s(func))

Definition at line 34 of file s_cexp.c.

#define s (   name)    CONCATX(name,SUFF)

Definition at line 33 of file s_cexp.c.

#define SUFF

Definition at line 26 of file s_cexp.c.

#define weak_aliasx (   a,
  b 
)    weak_alias(a,b)

Definition at line 116 of file s_cexp.c.


Function Documentation

__complex__ float_type s() __cexp ( __complex__ float_type  x)

Definition at line 37 of file s_cexp.c.

{
  __complex__ float_type retval;
  unsigned long ix_cond;

  ix_cond = __m81_test (__imag__ x);

  if ((ix_cond & (__M81_COND_NAN|__M81_COND_INF)) == 0)
    {
      /* Imaginary part is finite.  */
      float_type exp_val = m81(__ieee754_exp) (__real__ x);

      __real__ retval = __imag__ retval = exp_val;
      if (m81(__finite) (exp_val))
       {
         float_type sin_ix, cos_ix;
         __asm ("fsincos%.x %2,%1:%0" : "=f" (sin_ix), "=f" (cos_ix)
               : "f" (__imag__ x));
         __real__ retval *= cos_ix;
         if (ix_cond & __M81_COND_ZERO)
           __imag__ retval = __imag__ x;
         else
           __imag__ retval *= sin_ix;
       }
      else
       {
         /* Compute the sign of the result.  */
         float_type remainder, pi_2;
         int quadrant;

         __asm ("fmovecr %#0,%0\n\tfscale%.w %#-1,%0" : "=f" (pi_2));
         __asm ("fmod%.x %2,%0\n\tfmove%.l %/fpsr,%1"
               : "=f" (remainder), "=dm" (quadrant)
               : "f" (pi_2), "0" (__imag__ x));
         quadrant = (quadrant >> 16) & 0x83;
         if (quadrant & 0x80)
           quadrant ^= 0x83;
         switch (quadrant)
           {
           default:
             break;
           case 1:
             __real__ retval = -__real__ retval;
             break;
           case 2:
             __real__ retval = -__real__ retval;
           case 3:
             __imag__ retval = -__imag__ retval;
             break;
           }
         if (ix_cond & __M81_COND_ZERO && !m81(__isnan) (exp_val))
           __imag__ retval = __imag__ x;
       }
    }
  else
    {
      unsigned long rx_cond = __m81_test (__real__ x);

      if (rx_cond & __M81_COND_INF)
       {
         /* Real part is infinite.  */
         if (rx_cond & __M81_COND_NEG)
           {
             __real__ retval = __imag__ retval = 0.0;
             if (ix_cond & __M81_COND_NEG)
              __imag__ retval = -__imag__ retval;
           }
         else
           {
             __real__ retval = __real__ x;
             __imag__ retval = __imag__ x - __imag__ x;
           }
       }
      else
       __real__ retval = __imag__ retval = __imag__ x - __imag__ x;
    }

  return retval;
}

Here is the call graph for this function: