sqrtf.c File Reference

#include "libm.h"
#include "libm_util.h"
#include "libm_inlines.h"
#include "libm_errno.h"

Include dependency graph for sqrtf.c:

Go to the source code of this file.

Macros
#define	USE_NANF_WITH_FLAGS
#define	USE_HANDLE_ERRORF

Functions
float	sqrtf (float x)

Macro Definition Documentation

USE_HANDLE_ERRORF

#define USE_HANDLE_ERRORF

Definition at line 34 of file sqrtf.c.

USE_NANF_WITH_FLAGS

#define USE_NANF_WITH_FLAGS

Definition at line 33 of file sqrtf.c.

Function Documentation

sqrtf()

float sqrtf

(

float

x

)

Definition at line 50 of file sqrtf.c.

{
#if USE_SOFTWARE_SQRT
  return sqrtf_amd_inline(x);
#else
  float r;
  unsigned int ux;
  GET_BITS_SP32(x, ux);
  /* Check for special cases for Microsoft error handling */
  if ((ux & PINFBITPATT_SP32) == PINFBITPATT_SP32)
    {
      /* x is infinity, or NaN */
      if (ux & MANTBITS_SP32)
        {
          /* NaN  of some sort */
          /* If it's a signaling NaN, convert to QNaN */
            return _handle_errorf("sqrtf", OP_SQRT, ux|0x00400000, _DOMAIN, 0,
                               EDOM, x, 0.0F, 1);
        }
      else
        {
          /* +/-infinity  */
          if (ux & SIGNBIT_SP32)
            {
              /* - infinity */
                return _handle_errorf("sqrtf", OP_SQRT, INDEFBITPATT_SP32, 
                            _DOMAIN, AMD_F_INVALID, EDOM, x, 0.0F, 1);
            }
          /* positive infinite is not a problem */
        }
    }
  if ((ux & SIGNBIT_SP32)&&(ux & ~SIGNBIT_SP32))  /* if x < zero */
    {
        return _handle_errorf("sqrtf", OP_SQRT, INDEFBITPATT_SP32, 
                    _DOMAIN, AMD_F_INVALID, EDOM, x, 0.0F, 1);
    }
 
      /* VC++ intrinsic call */
      _mm_store_ss(&r, _mm_sqrt_ss(_mm_load_ss(&x)));
  return r;
#endif
}