src/npymath/npy_math.c.src File Reference

#include "npy_math_private.h"

Defines
#define	NPY_DBL_EPSILON   1.2246467991473531772E-16
#define	LOGE2   NPY_LOGE2@c@
#define	LOG2E   NPY_LOG2E@c@
#define	RAD2DEG   (180.0@c@/NPY_PI@c@)
#define	DEG2RAD   (NPY_PI@c@/180.0@c@)
Functions
double	npy_expm1 (double x)
double	npy_log1p (double x)
double	npy_atan2 (double y, double x)
double	npy_hypot (double x, double y)
double	npy_acosh (double x)
double	npy_asinh (double xx)
double	npy_atanh (double x)
double	npy_rint (double x)
double	npy_trunc (double x)
double	npy_exp2 (double x)
double	npy_log2 (double x)
type npy_ kind	c (@type @x)
type npy_ kind	c (@type @x,@type @y)
type npy_modf	c (@type @x,@type @*iptr)

---

Define Documentation

#define DEG2RAD   (NPY_PI@c@/180.0@c@)

#define LOG2E   NPY_LOG2E@c@

#define LOGE2   NPY_LOGE2@c@

end repeat*

Decorate all the math functions which are available on the current platform

begin repeat

type = npy_longdouble,double,float# c = l,,f# C = L,,F#

begin repeat1

kind = sin,cos,tan,sinh,cosh,tanh,fabs,floor,ceil,rint,trunc,sqrt,log10,

log,exp,expm1,asin,acos,atan,asinh,acosh,atanh,log1p,exp2,log2#

#KIND = SIN,COS,TAN,SINH,COSH,TANH,FABS,FLOOR,CEIL,RINT,TRUNC,SQRT,LOG10,

LOG,EXP,EXPM1,ASIN,ACOS,ATAN,ASINH,ACOSH,ATANH,LOG1P,EXP2,LOG2#

end repeat1*

begin repeat1

kind = atan2,hypot,pow,fmod,copysign# #KIND = ATAN2,HYPOT,POW,FMOD,COPYSIGN#

end repeat1*

end repeat*

Non standard functions

begin repeat

type = float, double, npy_longdouble# c = f, ,l# C = F, ,L#

#define NPY_DBL_EPSILON   1.2246467991473531772E-16

Taken from FreeBSD mlib, adapted for numpy

XXX: we could be a bit faster by reusing high/low words for inf/nan classification instead of calling npy_isinf/npy_isnan: we should have some macros for this, though, instead of doing it manually

XXX: we should have this in npy_math.h

#define RAD2DEG   (180.0@c@/NPY_PI@c@)

Referenced by c().

---

Function Documentation

type npy_exp2_m1 c

(

@type @

x

)

if C99 extensions not available then define dummy functions that use the double versions for

sin, cos, tan sinh, cosh, tanh, fabs, floor, ceil, rint, trunc sqrt, log10, log, exp, expm1 asin, acos, atan, asinh, acosh, atanh

hypot, atan2, pow, fmod, modf

We assume the above are always available in their double versions.

NOTE: some facilities may be available as macro only instead of functions. For simplicity, we define our own functions and undef the macros. We could instead test for the macro, but I am lazy to do that for now.

begin repeat

type = npy_longdouble, float# TYPE = NPY_LONGDOUBLE, FLOAT# c = l,f# C = L,F#

begin repeat1

kind = sin,cos,tan,sinh,cosh,tanh,fabs,floor,ceil,rint,trunc,sqrt,log10,

log,exp,expm1,asin,acos,atan,asinh,acosh,atanh,log1p,exp2,log2#

#KIND = SIN,COS,TAN,SINH,COSH,TANH,FABS,FLOOR,CEIL,RINT,TRUNC,SQRT,LOG10,

LOG,EXP,EXPM1,ASIN,ACOS,ATAN,ASINH,ACOSH,ATANH,LOG1P,EXP2,LOG2#

References RAD2DEG.

type npy_logaddexp2 c	(	@type @	x,
		@type @	y
	)

end repeat1*

begin repeat1

kind = atan2,hypot,pow,fmod,copysign# #KIND = ATAN2,HYPOT,POW,FMOD,COPYSIGN#

NaNs, or infinities of the same sign involved

NaNs, or infinities of the same sign involved

type npy_modf c	(	@type @	x,
		@type @*	iptr
	)

end repeat1*

double npy_acosh

(

double

x

)

double npy_asinh

(

double

xx

)

References kind(), and type().

double npy_atan2	(	double	y,
		double	x
	)

if x or y is nan, return nan

<

atan(+-0,+anything)=+-0

<

atan(+0,-anything) = pi

<

atan(-0,-anything) =-pi

<

atan(+INF,+INF)

<

atan(-INF,+INF)

<

atan(+INF,-INF)

<

atan(-INF,-INF)

<

atan(+...,+INF)

<

atan(-...,+INF)

<

atan(+...,-INF)

<

atan(-...,-INF)

compute y/x

<

|y/x| > 2**60

Docutils System Messages

System Message: ERROR/3 (<string>, line 1); backlink Undefined substitution referenced: "y/x".

<

0 > |y|/x > -2**-60

Docutils System Messages

System Message: ERROR/3 (<string>, line 1); backlink Undefined substitution referenced: "y".

<

safe to do y/x

<

atan(+,+)

<

atan(-,+)

<

atan(+,-)

<

case 3

<

atan(-,-)

Referenced by HALF_modf().

double npy_atanh

(

double

x

)

double npy_exp2

(

double

x

)

double npy_expm1

(

double

x

)

* BASIC MATH FUNCTIONS **

Original code by Konrad Hinsen.

double npy_hypot	(	double	x,
		double	y
	)

Handle the case where x or y is a NaN

double npy_log1p

(

double

x

)

double npy_log2

(

double

x

)

double npy_rint

(

double

x

)

Round to nearest even

double npy_trunc

(

double

x

)