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

Subversion Repositories openrisc

[/] [openrisc/] [tags/] [gnu-src/] [newlib-1.18.0/] [newlib-1.18.0-or32-1.0rc2/] [newlib/] [libc/] [sys/] [linux/] [cmath/] [s_cexpf.c] - Diff between revs 207 and 520

Only display areas with differences | Details | Blame | View Log

Rev 207 Rev 520
/* Return value of complex exponential function for float complex value.
/* Return value of complex exponential function for float complex value.
   Copyright (C) 1997 Free Software Foundation, Inc.
   Copyright (C) 1997 Free Software Foundation, Inc.
   This file is part of the GNU C Library.
   This file is part of the GNU C Library.
   Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
   Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
 
 
   The GNU C Library is free software; you can redistribute it and/or
   The GNU C Library is free software; you can redistribute it and/or
   modify it under the terms of the GNU Lesser General Public
   modify it under the terms of the GNU Lesser General Public
   License as published by the Free Software Foundation; either
   License as published by the Free Software Foundation; either
   version 2.1 of the License, or (at your option) any later version.
   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,
   The GNU C Library is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
   Lesser General Public License for more details.
   Lesser General Public License for more details.
 
 
   You should have received a copy of the GNU Lesser General Public
   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
   License along with the GNU C Library; if not, write to the Free
   Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
   Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
   02111-1307 USA.  */
   02111-1307 USA.  */
 
 
#include <complex.h>
#include <complex.h>
#include <fenv.h>
#include <fenv.h>
#include <math.h>
#include <math.h>
 
 
#include "math_private.h"
#include "math_private.h"
 
 
 
 
__complex__ float
__complex__ float
__cexpf (__complex__ float x)
__cexpf (__complex__ float x)
{
{
  __complex__ float retval;
  __complex__ float retval;
  int rcls = fpclassify (__real__ x);
  int rcls = fpclassify (__real__ x);
  int icls = fpclassify (__imag__ x);
  int icls = fpclassify (__imag__ x);
 
 
  if (rcls >= FP_ZERO)
  if (rcls >= FP_ZERO)
    {
    {
      /* Real part is finite.  */
      /* Real part is finite.  */
      if (icls >= FP_ZERO)
      if (icls >= FP_ZERO)
        {
        {
          /* Imaginary part is finite.  */
          /* Imaginary part is finite.  */
          float exp_val = __ieee754_expf (__real__ x);
          float exp_val = __ieee754_expf (__real__ x);
          float sinix, cosix;
          float sinix, cosix;
 
 
          __sincosf (__imag__ x, &sinix, &cosix);
          __sincosf (__imag__ x, &sinix, &cosix);
 
 
          if (isfinite (exp_val))
          if (isfinite (exp_val))
            {
            {
              __real__ retval = exp_val * cosix;
              __real__ retval = exp_val * cosix;
              __imag__ retval = exp_val * sinix;
              __imag__ retval = exp_val * sinix;
            }
            }
          else
          else
            {
            {
              __real__ retval = __copysignf (exp_val, cosix);
              __real__ retval = __copysignf (exp_val, cosix);
              __imag__ retval = __copysignf (exp_val, sinix);
              __imag__ retval = __copysignf (exp_val, sinix);
            }
            }
        }
        }
      else
      else
        {
        {
          /* If the imaginary part is +-inf or NaN and the real part
          /* If the imaginary part is +-inf or NaN and the real part
             is not +-inf the result is NaN + iNaN.  */
             is not +-inf the result is NaN + iNaN.  */
          __real__ retval = __nanf ("");
          __real__ retval = __nanf ("");
          __imag__ retval = __nanf ("");
          __imag__ retval = __nanf ("");
 
 
#ifdef FE_INVALID
#ifdef FE_INVALID
          feraiseexcept (FE_INVALID);
          feraiseexcept (FE_INVALID);
#endif
#endif
        }
        }
    }
    }
  else if (rcls == FP_INFINITE)
  else if (rcls == FP_INFINITE)
    {
    {
      /* Real part is infinite.  */
      /* Real part is infinite.  */
      if (icls >= FP_ZERO)
      if (icls >= FP_ZERO)
        {
        {
          /* Imaginary part is finite.  */
          /* Imaginary part is finite.  */
          float value = signbit (__real__ x) ? 0.0 : HUGE_VALF;
          float value = signbit (__real__ x) ? 0.0 : HUGE_VALF;
 
 
          if (icls == FP_ZERO)
          if (icls == FP_ZERO)
            {
            {
              /* Imaginary part is 0.0.  */
              /* Imaginary part is 0.0.  */
              __real__ retval = value;
              __real__ retval = value;
              __imag__ retval = __imag__ x;
              __imag__ retval = __imag__ x;
            }
            }
          else
          else
            {
            {
              float sinix, cosix;
              float sinix, cosix;
 
 
              __sincosf (__imag__ x, &sinix, &cosix);
              __sincosf (__imag__ x, &sinix, &cosix);
 
 
              __real__ retval = __copysignf (value, cosix);
              __real__ retval = __copysignf (value, cosix);
              __imag__ retval = __copysignf (value, sinix);
              __imag__ retval = __copysignf (value, sinix);
            }
            }
        }
        }
      else if (signbit (__real__ x) == 0)
      else if (signbit (__real__ x) == 0)
        {
        {
          __real__ retval = HUGE_VALF;
          __real__ retval = HUGE_VALF;
          __imag__ retval = __nanf ("");
          __imag__ retval = __nanf ("");
 
 
#ifdef FE_INVALID
#ifdef FE_INVALID
          if (icls == FP_INFINITE)
          if (icls == FP_INFINITE)
            feraiseexcept (FE_INVALID);
            feraiseexcept (FE_INVALID);
#endif
#endif
        }
        }
      else
      else
        {
        {
          __real__ retval = 0.0;
          __real__ retval = 0.0;
          __imag__ retval = __copysignf (0.0, __imag__ x);
          __imag__ retval = __copysignf (0.0, __imag__ x);
        }
        }
    }
    }
  else
  else
    {
    {
      /* If the real part is NaN the result is NaN + iNaN.  */
      /* If the real part is NaN the result is NaN + iNaN.  */
      __real__ retval = __nanf ("");
      __real__ retval = __nanf ("");
      __imag__ retval = __nanf ("");
      __imag__ retval = __nanf ("");
 
 
#ifdef FE_INVALID
#ifdef FE_INVALID
      if (rcls != FP_NAN || icls != FP_NAN)
      if (rcls != FP_NAN || icls != FP_NAN)
        feraiseexcept (FE_INVALID);
        feraiseexcept (FE_INVALID);
#endif
#endif
    }
    }
 
 
  return retval;
  return retval;
}
}
weak_alias (__cexpf, cexpf)
weak_alias (__cexpf, cexpf)
 
 

powered by: WebSVN 2.1.0

© copyright 1999-2024 OpenCores.org, equivalent to Oliscience, all rights reserved. OpenCores®, registered trademark.