log_special.c File Reference

#include <fpieee.h>
#include <excpt.h>
#include <float.h>
#include <math.h>
#include <errno.h>
#include "libm_new.h"

Include dependency graph for log_special.c:

Go to the source code of this file.

Macros
#define	LOG_X_ZERO   1
#define	LOG_X_NEG   2
#define	LOG_X_NAN   3

Functions
static float	_logf_special_common (float x, float y, U32 code, unsigned int op, char *name)
float	_logf_special (float x, float y, U32 code)
float	_log10f_special (float x, float y, U32 code)
static double	_log_special_common (double x, double y, U32 code, unsigned int op, char *name)
double	_log_special (double x, double y, U32 code)
double	_log10_special (double x, double y, U32 code)

Macro Definition Documentation

LOG_X_NAN

#define LOG_X_NAN   3

Definition at line 43 of file log_special.c.

LOG_X_NEG

#define LOG_X_NEG   2

Definition at line 42 of file log_special.c.

LOG_X_ZERO

#define LOG_X_ZERO   1

Definition at line 41 of file log_special.c.

Function Documentation

_log10_special()

double _log10_special	(	double	x,
		double	y,
		U32	code
	)

Definition at line 130 of file log_special.c.

{
    return _log_special_common(x, y, code, _FpCodeLog10, "log10");
}

_log10f_special()

float _log10f_special	(	float	x,
		float	y,
		U32	code
	)

Definition at line 91 of file log_special.c.

{
    return _logf_special_common(x, y, code, _FpCodeLog10, "log10f");
}

_log_special()

double _log_special	(	double	x,
		double	y,
		U32	code
	)

Definition at line 125 of file log_special.c.

{
    return _log_special_common(x, y, code, _FpCodeLog, "log");
}

_log_special_common()

static double _log_special_common ( double  x, double  y, U32  code, unsigned int  op, char *  name  )

static

Definition at line 96 of file log_special.c.

{
    switch(code)
    {
    case LOG_X_ZERO:
        {
            UT64 ym; ym.f64 = y;
            _handle_error(name, op, ym.u64, _SING, AMD_F_DIVBYZERO, ERANGE, x, 0.0, 1);
        }
        break;
 
    case LOG_X_NEG:
        {
            UT64 ym; ym.f64 = y;
            _handle_error(name, op, ym.u64, _DOMAIN, AMD_F_INVALID, EDOM, x, 0.0, 1);
        }
        break;
 
    case LOG_X_NAN:
        {
            UT64 ym; ym.f64 = y;
            _handle_error(name, op, ym.u64, _DOMAIN, 0, EDOM, x, 0.0, 1);
        }
        break;
    }
 
    return y;
}

Referenced by _log10_special(), and _log_special().

_logf_special()

float _logf_special	(	float	x,
		float	y,
		U32	code
	)

Definition at line 86 of file log_special.c.

{
    return _logf_special_common(x, y, code, _FpCodeLog, "logf");
}

_logf_special_common()

static float _logf_special_common ( float  x, float  y, U32  code, unsigned int  op, char *  name  )

static

Definition at line 45 of file log_special.c.

{
    switch(code)
    {
    case LOG_X_ZERO:
        {
            UT64 ym; ym.u64 = 0; ym.f32[0] = y;
            _handle_errorf(name, op, ym.u64, _SING, AMD_F_DIVBYZERO, ERANGE, x, 0.0, 1);
        }
        break;
 
    case LOG_X_NEG:
        {
            UT64 ym; ym.u64 = 0; ym.f32[0] = y;
            _handle_errorf(name, op, ym.u64, _DOMAIN, AMD_F_INVALID, EDOM, x, 0.0, 1);
        }
        break;
 
    case LOG_X_NAN:
        {
            unsigned int is_snan;
            UT32 xm; UT64 ym;
            xm.f32 = x;
            is_snan = (((xm.u32 & QNAN_MASK_32) == QNAN_SET_32) ? 0 : 1);
            ym.u64 = 0; ym.f32[0] = y;
 
            if(is_snan)
            {
                _handle_errorf(name, op, ym.u64, _DOMAIN, AMD_F_INVALID, EDOM, x, 0.0, 1);
            }
            else
            {
                _handle_errorf(name, op, ym.u64, _DOMAIN, 0, EDOM, x, 0.0, 1);
            }
        }
        break;
    }
 
    return y;
}

Referenced by _log10f_special(), and _logf_special().