[deliverable/linux.git] / arch / mips / math-emu / ieee754.h

/*
 * MIPS floating point support
 * Copyright (C) 1994-2000 Algorithmics Ltd.
 *
 *  This program is free software; you can distribute it and/or modify it
 *  under the terms of the GNU General Public License (Version 2) as
 *  published by the Free Software Foundation.
 *
 *  This program is distributed in the hope it will be useful, but WITHOUT
 *  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
 *  FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
 *  for more details.
 *
 *  You should have received a copy of the GNU General Public License along
 *  with this program; if not, write to the Free Software Foundation, Inc.,
 *  59 Temple Place - Suite 330, Boston MA 02111-1307, USA.
 *
 *  Nov 7, 2000
 *  Modification to allow integration with Linux kernel
 *
 *  Kevin D. Kissell, kevink@mips.com and Carsten Langgard, carstenl@mips.com
 *  Copyright (C) 2000 MIPS Technologies, Inc. All rights reserved.
 */
#ifndef __ARCH_MIPS_MATH_EMU_IEEE754_H
#define __ARCH_MIPS_MATH_EMU_IEEE754_H

#include <linux/compiler.h>
#include <asm/byteorder.h>
#include <linux/types.h>
#include <linux/sched.h>
#include <asm/bitfield.h>

/*
 * Not very pretty, but the Linux kernel's normal va_list definition
 * does not allow it to be used as a structure element, as it is here.
 */
#ifndef _STDARG_H
#include <stdarg.h>
#endif

struct ieee754dp_const {
	__BITFIELD_FIELD(unsigned sign:1,
	__BITFIELD_FIELD(unsigned bexp:11,
	__BITFIELD_FIELD(unsigned manthi:20,
	__BITFIELD_FIELD(unsigned mantlo:32,
	;))))
};

union ieee754dp {
	struct ieee754dp_const oparts;
	struct {
		__BITFIELD_FIELD(unsigned int sign:1,
		__BITFIELD_FIELD(unsigned int bexp:11,
		__BITFIELD_FIELD(u64 mant:52,
		;)))
	} parts;
	double d;
	u64 bits;
};

struct ieee754sp_const {
	__BITFIELD_FIELD(unsigned sign:1,
	__BITFIELD_FIELD(unsigned bexp:8,
	__BITFIELD_FIELD(unsigned mant:23,
	;)))
};

union ieee754sp {
	struct ieee754sp_const parts;
	float f;
	u32 bits;
};

/*
 * single precision (often aka float)
*/
int ieee754sp_finite(union ieee754sp x);
int ieee754sp_class(union ieee754sp x);

union ieee754sp ieee754sp_abs(union ieee754sp x);
union ieee754sp ieee754sp_neg(union ieee754sp x);
union ieee754sp ieee754sp_scalb(union ieee754sp x, int);
union ieee754sp ieee754sp_logb(union ieee754sp x);

/* x with sign of y */
union ieee754sp ieee754sp_copysign(union ieee754sp x, union ieee754sp y);

union ieee754sp ieee754sp_add(union ieee754sp x, union ieee754sp y);
union ieee754sp ieee754sp_sub(union ieee754sp x, union ieee754sp y);
union ieee754sp ieee754sp_mul(union ieee754sp x, union ieee754sp y);
union ieee754sp ieee754sp_div(union ieee754sp x, union ieee754sp y);

union ieee754sp ieee754sp_fint(int x);
union ieee754sp ieee754sp_funs(unsigned x);
union ieee754sp ieee754sp_flong(s64 x);
union ieee754sp ieee754sp_fulong(u64 x);
union ieee754sp ieee754sp_fdp(union ieee754dp x);

int ieee754sp_tint(union ieee754sp x);
unsigned int ieee754sp_tuns(union ieee754sp x);
s64 ieee754sp_tlong(union ieee754sp x);
u64 ieee754sp_tulong(union ieee754sp x);

int ieee754sp_cmp(union ieee754sp x, union ieee754sp y, int cop, int sig);
/*
 * basic sp math
 */
union ieee754sp ieee754sp_modf(union ieee754sp x, union ieee754sp * ip);
union ieee754sp ieee754sp_frexp(union ieee754sp x, int *exp);
union ieee754sp ieee754sp_ldexp(union ieee754sp x, int exp);

union ieee754sp ieee754sp_ceil(union ieee754sp x);
union ieee754sp ieee754sp_floor(union ieee754sp x);
union ieee754sp ieee754sp_trunc(union ieee754sp x);

union ieee754sp ieee754sp_sqrt(union ieee754sp x);

/*
 * double precision (often aka double)
*/
int ieee754dp_finite(union ieee754dp x);
int ieee754dp_class(union ieee754dp x);

/* x with sign of y */
union ieee754dp ieee754dp_copysign(union ieee754dp x, union ieee754dp y);

union ieee754dp ieee754dp_add(union ieee754dp x, union ieee754dp y);
union ieee754dp ieee754dp_sub(union ieee754dp x, union ieee754dp y);
union ieee754dp ieee754dp_mul(union ieee754dp x, union ieee754dp y);
union ieee754dp ieee754dp_div(union ieee754dp x, union ieee754dp y);

union ieee754dp ieee754dp_abs(union ieee754dp x);
union ieee754dp ieee754dp_neg(union ieee754dp x);
union ieee754dp ieee754dp_scalb(union ieee754dp x, int);

/* return exponent as integer in floating point format
 */
union ieee754dp ieee754dp_logb(union ieee754dp x);

union ieee754dp ieee754dp_fint(int x);
union ieee754dp ieee754dp_funs(unsigned x);
union ieee754dp ieee754dp_flong(s64 x);
union ieee754dp ieee754dp_fulong(u64 x);
union ieee754dp ieee754dp_fsp(union ieee754sp x);

union ieee754dp ieee754dp_ceil(union ieee754dp x);
union ieee754dp ieee754dp_floor(union ieee754dp x);
union ieee754dp ieee754dp_trunc(union ieee754dp x);

int ieee754dp_tint(union ieee754dp x);
unsigned int ieee754dp_tuns(union ieee754dp x);
s64 ieee754dp_tlong(union ieee754dp x);
u64 ieee754dp_tulong(union ieee754dp x);

int ieee754dp_cmp(union ieee754dp x, union ieee754dp y, int cop, int sig);
/*
 * basic sp math
 */
union ieee754dp ieee754dp_modf(union ieee754dp x, union ieee754dp * ip);
union ieee754dp ieee754dp_frexp(union ieee754dp x, int *exp);
union ieee754dp ieee754dp_ldexp(union ieee754dp x, int exp);

union ieee754dp ieee754dp_ceil(union ieee754dp x);
union ieee754dp ieee754dp_floor(union ieee754dp x);
union ieee754dp ieee754dp_trunc(union ieee754dp x);

union ieee754dp ieee754dp_sqrt(union ieee754dp x);


/* 5 types of floating point number
*/
enum {
	IEEE754_CLASS_NORM	= 0x00,
	IEEE754_CLASS_ZERO	= 0x01,
	IEEE754_CLASS_DNORM	= 0x02,
	IEEE754_CLASS_INF	= 0x03,
	IEEE754_CLASS_SNAN	= 0x04,
	IEEE754_CLASS_QNAN	= 0x05,
};

/* exception numbers */
#define IEEE754_INEXACT			0x01
#define IEEE754_UNDERFLOW		0x02
#define IEEE754_OVERFLOW		0x04
#define IEEE754_ZERO_DIVIDE		0x08
#define IEEE754_INVALID_OPERATION	0x10

/* cmp operators
*/
#define IEEE754_CLT	0x01
#define IEEE754_CEQ	0x02
#define IEEE754_CGT	0x04
#define IEEE754_CUN	0x08

/* rounding mode
*/
#define IEEE754_RN	0	/* round to nearest */
#define IEEE754_RZ	1	/* round toward zero  */
#define IEEE754_RD	2	/* round toward -Infinity */
#define IEEE754_RU	3	/* round toward +Infinity */

/* other naming */
#define IEEE754_RM	IEEE754_RD
#define IEEE754_RP	IEEE754_RU

/* "normal" comparisons
*/
static inline int ieee754sp_eq(union ieee754sp x, union ieee754sp y)
{
	return ieee754sp_cmp(x, y, IEEE754_CEQ, 0);
}

static inline int ieee754sp_ne(union ieee754sp x, union ieee754sp y)
{
	return ieee754sp_cmp(x, y,
			     IEEE754_CLT | IEEE754_CGT | IEEE754_CUN, 0);
}

static inline int ieee754sp_lt(union ieee754sp x, union ieee754sp y)
{
	return ieee754sp_cmp(x, y, IEEE754_CLT, 0);
}

static inline int ieee754sp_le(union ieee754sp x, union ieee754sp y)
{
	return ieee754sp_cmp(x, y, IEEE754_CLT | IEEE754_CEQ, 0);
}

static inline int ieee754sp_gt(union ieee754sp x, union ieee754sp y)
{
	return ieee754sp_cmp(x, y, IEEE754_CGT, 0);
}


static inline int ieee754sp_ge(union ieee754sp x, union ieee754sp y)
{
	return ieee754sp_cmp(x, y, IEEE754_CGT | IEEE754_CEQ, 0);
}

static inline int ieee754dp_eq(union ieee754dp x, union ieee754dp y)
{
	return ieee754dp_cmp(x, y, IEEE754_CEQ, 0);
}

static inline int ieee754dp_ne(union ieee754dp x, union ieee754dp y)
{
	return ieee754dp_cmp(x, y,
			     IEEE754_CLT | IEEE754_CGT | IEEE754_CUN, 0);
}

static inline int ieee754dp_lt(union ieee754dp x, union ieee754dp y)
{
	return ieee754dp_cmp(x, y, IEEE754_CLT, 0);
}

static inline int ieee754dp_le(union ieee754dp x, union ieee754dp y)
{
	return ieee754dp_cmp(x, y, IEEE754_CLT | IEEE754_CEQ, 0);
}

static inline int ieee754dp_gt(union ieee754dp x, union ieee754dp y)
{
	return ieee754dp_cmp(x, y, IEEE754_CGT, 0);
}

static inline int ieee754dp_ge(union ieee754dp x, union ieee754dp y)
{
	return ieee754dp_cmp(x, y, IEEE754_CGT | IEEE754_CEQ, 0);
}


/*
 * Like strtod
 */
union ieee754dp ieee754dp_fstr(const char *s, char **endp);
char *ieee754dp_tstr(union ieee754dp x, int prec, int fmt, int af);


/*
 * The control status register
 */
struct _ieee754_csr {
	__BITFIELD_FIELD(unsigned pad0:7,
	__BITFIELD_FIELD(unsigned nod:1,	/* set 1 for no denormalised numbers */
	__BITFIELD_FIELD(unsigned c:1,		/* condition */
	__BITFIELD_FIELD(unsigned pad1:5,
	__BITFIELD_FIELD(unsigned cx:6,		/* exceptions this operation */
	__BITFIELD_FIELD(unsigned mx:5,		/* exception enable  mask */
	__BITFIELD_FIELD(unsigned sx:5,		/* exceptions total */
	__BITFIELD_FIELD(unsigned rm:2,		/* current rounding mode */
	;))))))))
};
#define ieee754_csr (*(struct _ieee754_csr *)(&current->thread.fpu.fcr31))

static inline unsigned ieee754_getrm(void)
{
	return (ieee754_csr.rm);
}
static inline unsigned ieee754_setrm(unsigned rm)
{
	return (ieee754_csr.rm = rm);
}

/*
 * get current exceptions
 */
static inline unsigned ieee754_getcx(void)
{
	return (ieee754_csr.cx);
}

/* test for current exception condition
 */
static inline int ieee754_cxtest(unsigned n)
{
	return (ieee754_csr.cx & n);
}

/*
 * get sticky exceptions
 */
static inline unsigned ieee754_getsx(void)
{
	return (ieee754_csr.sx);
}

/* clear sticky conditions
*/
static inline unsigned ieee754_clrsx(void)
{
	return (ieee754_csr.sx = 0);
}

/* test for sticky exception condition
 */
static inline int ieee754_sxtest(unsigned n)
{
	return (ieee754_csr.sx & n);
}

/* debugging */
union ieee754sp ieee754sp_dump(char *s, union ieee754sp x);
union ieee754dp ieee754dp_dump(char *s, union ieee754dp x);

#define IEEE754_SPCVAL_PZERO	0
#define IEEE754_SPCVAL_NZERO	1
#define IEEE754_SPCVAL_PONE	2
#define IEEE754_SPCVAL_NONE	3
#define IEEE754_SPCVAL_PTEN	4
#define IEEE754_SPCVAL_NTEN	5
#define IEEE754_SPCVAL_PINFINITY	6
#define IEEE754_SPCVAL_NINFINITY	7
#define IEEE754_SPCVAL_INDEF	8
#define IEEE754_SPCVAL_PMAX	9	/* +max norm */
#define IEEE754_SPCVAL_NMAX	10	/* -max norm */
#define IEEE754_SPCVAL_PMIN	11	/* +min norm */
#define IEEE754_SPCVAL_NMIN	12	/* +min norm */
#define IEEE754_SPCVAL_PMIND	13	/* +min denorm */
#define IEEE754_SPCVAL_NMIND	14	/* +min denorm */
#define IEEE754_SPCVAL_P1E31	15	/* + 1.0e31 */
#define IEEE754_SPCVAL_P1E63	16	/* + 1.0e63 */

extern const struct ieee754dp_const __ieee754dp_spcvals[];
extern const struct ieee754sp_const __ieee754sp_spcvals[];
#define ieee754dp_spcvals ((const union ieee754dp *)__ieee754dp_spcvals)
#define ieee754sp_spcvals ((const union ieee754sp *)__ieee754sp_spcvals)

/*
 * Return infinity with given sign
 */
#define ieee754dp_inf(sn)     (ieee754dp_spcvals[IEEE754_SPCVAL_PINFINITY+(sn)])
#define ieee754dp_zero(sn)	(ieee754dp_spcvals[IEEE754_SPCVAL_PZERO+(sn)])
#define ieee754dp_one(sn)	(ieee754dp_spcvals[IEEE754_SPCVAL_PONE+(sn)])
#define ieee754dp_ten(sn)	(ieee754dp_spcvals[IEEE754_SPCVAL_PTEN+(sn)])
#define ieee754dp_indef()	(ieee754dp_spcvals[IEEE754_SPCVAL_INDEF])
#define ieee754dp_max(sn)	(ieee754dp_spcvals[IEEE754_SPCVAL_PMAX+(sn)])
#define ieee754dp_min(sn)	(ieee754dp_spcvals[IEEE754_SPCVAL_PMIN+(sn)])
#define ieee754dp_mind(sn)	(ieee754dp_spcvals[IEEE754_SPCVAL_PMIND+(sn)])
#define ieee754dp_1e31()	(ieee754dp_spcvals[IEEE754_SPCVAL_P1E31])
#define ieee754dp_1e63()	(ieee754dp_spcvals[IEEE754_SPCVAL_P1E63])

#define ieee754sp_inf(sn)     (ieee754sp_spcvals[IEEE754_SPCVAL_PINFINITY+(sn)])
#define ieee754sp_zero(sn)	(ieee754sp_spcvals[IEEE754_SPCVAL_PZERO+(sn)])
#define ieee754sp_one(sn)	(ieee754sp_spcvals[IEEE754_SPCVAL_PONE+(sn)])
#define ieee754sp_ten(sn)	(ieee754sp_spcvals[IEEE754_SPCVAL_PTEN+(sn)])
#define ieee754sp_indef()	(ieee754sp_spcvals[IEEE754_SPCVAL_INDEF])
#define ieee754sp_max(sn)	(ieee754sp_spcvals[IEEE754_SPCVAL_PMAX+(sn)])
#define ieee754sp_min(sn)	(ieee754sp_spcvals[IEEE754_SPCVAL_PMIN+(sn)])
#define ieee754sp_mind(sn)	(ieee754sp_spcvals[IEEE754_SPCVAL_PMIND+(sn)])
#define ieee754sp_1e31()	(ieee754sp_spcvals[IEEE754_SPCVAL_P1E31])
#define ieee754sp_1e63()	(ieee754sp_spcvals[IEEE754_SPCVAL_P1E63])

/*
 * Indefinite integer value
 */
#define ieee754si_indef()	INT_MAX
#ifdef LONG_LONG_MAX
#define ieee754di_indef()	LONG_LONG_MAX
#else
#define ieee754di_indef()	((s64)(~0ULL>>1))
#endif

/* IEEE exception context, passed to handler */
struct ieee754xctx {
	const char *op;		/* operation name */
	int rt;			/* result type */
	union {
		union ieee754sp sp;	/* single precision */
		union ieee754dp dp;	/* double precision */
		int si;		/* standard signed integer (32bits) */
		s64 di;		/* extended signed integer (64bits) */
	} rv;			/* default result format implied by op */
	va_list ap;
};

/* result types for xctx.rt */
#define IEEE754_RT_SP	0
#define IEEE754_RT_DP	1
#define IEEE754_RT_XP	2
#define IEEE754_RT_SI	3
#define IEEE754_RT_DI	4

extern void __cold ieee754_xcpt(struct ieee754xctx *xcp);

/* compat */
#define ieee754dp_fix(x)	ieee754dp_tint(x)
#define ieee754sp_fix(x)	ieee754sp_tint(x)

#endif /* __ARCH_MIPS_MATH_EMU_IEEE754_H */
Commit	Line	Data
1da177e4 LT	1	/*
	2	* MIPS floating point support
	3	* Copyright (C) 1994-2000 Algorithmics Ltd.
1da177e4	4	*
1da177e4 LT	5	* This program is free software; you can distribute it and/or modify it
	6	* under the terms of the GNU General Public License (Version 2) as
	7	* published by the Free Software Foundation.
	8	*
	9	* This program is distributed in the hope it will be useful, but WITHOUT
	10	* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
	11	* FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
	12	* for more details.
	13	*
	14	* You should have received a copy of the GNU General Public License along
	15	* with this program; if not, write to the Free Software Foundation, Inc.,
	16	* 59 Temple Place - Suite 330, Boston MA 02111-1307, USA.
	17	*
1da177e4 LT	18	* Nov 7, 2000
	19	* Modification to allow integration with Linux kernel
	20	*
	21	* Kevin D. Kissell, kevink@mips.com and Carsten Langgard, carstenl@mips.com
	22	* Copyright (C) 2000 MIPS Technologies, Inc. All rights reserved.
cd21dfcf	23	*/
a663bf90 RB	24	#ifndef __ARCH_MIPS_MATH_EMU_IEEE754_H
a663bf90 RB	25	#define __ARCH_MIPS_MATH_EMU_IEEE754_H
1da177e4	26
cae55066	27	#include <linux/compiler.h>
cd21dfcf	28	#include <asm/byteorder.h>
1da177e4	29	#include <linux/types.h>
cd21dfcf	30	#include <linux/sched.h>
f80cc08d	31	#include <asm/bitfield.h>
1da177e4 LT	32
	33	/*
	34	* Not very pretty, but the Linux kernel's normal va_list definition
	35	* does not allow it to be used as a structure element, as it is here.
	36	*/
	37	#ifndef _STDARG_H
	38	#include <stdarg.h>
	39	#endif
	40
23708818	41	struct ieee754dp_const {
f80cc08d RB	42	__BITFIELD_FIELD(unsigned sign:1,
	43	__BITFIELD_FIELD(unsigned bexp:11,
	44	__BITFIELD_FIELD(unsigned manthi:20,
	45	__BITFIELD_FIELD(unsigned mantlo:32,
	46	;))))
1da177e4 LT	47	};
1da177e4 LT	48
2209bcb1	49	union ieee754dp {
23708818	50	struct ieee754dp_const oparts;
1da177e4	51	struct {
f80cc08d RB	52	__BITFIELD_FIELD(unsigned int sign:1,
	53	__BITFIELD_FIELD(unsigned int bexp:11,
	54	__BITFIELD_FIELD(u64 mant:52,
	55	;)))
1da177e4 LT	56	} parts;
	57	double d;
	58	u64 bits;
2209bcb1	59	};
1da177e4	60
23708818	61	struct ieee754sp_const {
f80cc08d RB	62	__BITFIELD_FIELD(unsigned sign:1,
	63	__BITFIELD_FIELD(unsigned bexp:8,
	64	__BITFIELD_FIELD(unsigned mant:23,
	65	;)))
1da177e4 LT	66	};
1da177e4 LT	67
2209bcb1	68	union ieee754sp {
23708818	69	struct ieee754sp_const parts;
1da177e4 LT	70	float f;
1da177e4 LT	71	u32 bits;
2209bcb1	72	};
1da177e4 LT	73
	74	/*
	75	* single precision (often aka float)
	76	*/
2209bcb1 RB	77	int ieee754sp_finite(union ieee754sp x);
2209bcb1 RB	78	int ieee754sp_class(union ieee754sp x);
1da177e4	79
2209bcb1 RB	80	union ieee754sp ieee754sp_abs(union ieee754sp x);
	81	union ieee754sp ieee754sp_neg(union ieee754sp x);
	82	union ieee754sp ieee754sp_scalb(union ieee754sp x, int);
	83	union ieee754sp ieee754sp_logb(union ieee754sp x);
1da177e4 LT	84
1da177e4 LT	85	/* x with sign of y */
2209bcb1	86	union ieee754sp ieee754sp_copysign(union ieee754sp x, union ieee754sp y);
1da177e4	87
2209bcb1 RB	88	union ieee754sp ieee754sp_add(union ieee754sp x, union ieee754sp y);
	89	union ieee754sp ieee754sp_sub(union ieee754sp x, union ieee754sp y);
	90	union ieee754sp ieee754sp_mul(union ieee754sp x, union ieee754sp y);
	91	union ieee754sp ieee754sp_div(union ieee754sp x, union ieee754sp y);
1da177e4	92
2209bcb1 RB	93	union ieee754sp ieee754sp_fint(int x);
	94	union ieee754sp ieee754sp_funs(unsigned x);
	95	union ieee754sp ieee754sp_flong(s64 x);
	96	union ieee754sp ieee754sp_fulong(u64 x);
	97	union ieee754sp ieee754sp_fdp(union ieee754dp x);
1da177e4	98
2209bcb1 RB	99	int ieee754sp_tint(union ieee754sp x);
	100	unsigned int ieee754sp_tuns(union ieee754sp x);
	101	s64 ieee754sp_tlong(union ieee754sp x);
	102	u64 ieee754sp_tulong(union ieee754sp x);
1da177e4	103
2209bcb1	104	int ieee754sp_cmp(union ieee754sp x, union ieee754sp y, int cop, int sig);
1da177e4 LT	105	/*
	106	* basic sp math
	107	*/
2209bcb1 RB	108	union ieee754sp ieee754sp_modf(union ieee754sp x, union ieee754sp * ip);
	109	union ieee754sp ieee754sp_frexp(union ieee754sp x, int *exp);
	110	union ieee754sp ieee754sp_ldexp(union ieee754sp x, int exp);
1da177e4	111
2209bcb1 RB	112	union ieee754sp ieee754sp_ceil(union ieee754sp x);
	113	union ieee754sp ieee754sp_floor(union ieee754sp x);
	114	union ieee754sp ieee754sp_trunc(union ieee754sp x);
1da177e4	115
2209bcb1	116	union ieee754sp ieee754sp_sqrt(union ieee754sp x);
1da177e4 LT	117
	118	/*
	119	* double precision (often aka double)
	120	*/
2209bcb1 RB	121	int ieee754dp_finite(union ieee754dp x);
2209bcb1 RB	122	int ieee754dp_class(union ieee754dp x);
1da177e4 LT	123
1da177e4 LT	124	/* x with sign of y */
2209bcb1	125	union ieee754dp ieee754dp_copysign(union ieee754dp x, union ieee754dp y);
1da177e4	126
2209bcb1 RB	127	union ieee754dp ieee754dp_add(union ieee754dp x, union ieee754dp y);
	128	union ieee754dp ieee754dp_sub(union ieee754dp x, union ieee754dp y);
	129	union ieee754dp ieee754dp_mul(union ieee754dp x, union ieee754dp y);
	130	union ieee754dp ieee754dp_div(union ieee754dp x, union ieee754dp y);
1da177e4	131
2209bcb1 RB	132	union ieee754dp ieee754dp_abs(union ieee754dp x);
	133	union ieee754dp ieee754dp_neg(union ieee754dp x);
	134	union ieee754dp ieee754dp_scalb(union ieee754dp x, int);
1da177e4 LT	135
	136	/* return exponent as integer in floating point format
	137	*/
2209bcb1	138	union ieee754dp ieee754dp_logb(union ieee754dp x);
1da177e4	139
2209bcb1 RB	140	union ieee754dp ieee754dp_fint(int x);
	141	union ieee754dp ieee754dp_funs(unsigned x);
	142	union ieee754dp ieee754dp_flong(s64 x);
	143	union ieee754dp ieee754dp_fulong(u64 x);
	144	union ieee754dp ieee754dp_fsp(union ieee754sp x);
1da177e4	145
2209bcb1 RB	146	union ieee754dp ieee754dp_ceil(union ieee754dp x);
	147	union ieee754dp ieee754dp_floor(union ieee754dp x);
	148	union ieee754dp ieee754dp_trunc(union ieee754dp x);
1da177e4	149
2209bcb1 RB	150	int ieee754dp_tint(union ieee754dp x);
	151	unsigned int ieee754dp_tuns(union ieee754dp x);
	152	s64 ieee754dp_tlong(union ieee754dp x);
	153	u64 ieee754dp_tulong(union ieee754dp x);
1da177e4	154
2209bcb1	155	int ieee754dp_cmp(union ieee754dp x, union ieee754dp y, int cop, int sig);
1da177e4 LT	156	/*
	157	* basic sp math
	158	*/
2209bcb1 RB	159	union ieee754dp ieee754dp_modf(union ieee754dp x, union ieee754dp * ip);
	160	union ieee754dp ieee754dp_frexp(union ieee754dp x, int *exp);
	161	union ieee754dp ieee754dp_ldexp(union ieee754dp x, int exp);
1da177e4	162
2209bcb1 RB	163	union ieee754dp ieee754dp_ceil(union ieee754dp x);
	164	union ieee754dp ieee754dp_floor(union ieee754dp x);
	165	union ieee754dp ieee754dp_trunc(union ieee754dp x);
1da177e4	166
2209bcb1	167	union ieee754dp ieee754dp_sqrt(union ieee754dp x);
1da177e4 LT	168
	169
	170
	171	/* 5 types of floating point number
	172	*/
9e8bad1f RB	173	enum {
	174	IEEE754_CLASS_NORM = 0x00,
	175	IEEE754_CLASS_ZERO = 0x01,
	176	IEEE754_CLASS_DNORM = 0x02,
	177	IEEE754_CLASS_INF = 0x03,
	178	IEEE754_CLASS_SNAN = 0x04,
	179	IEEE754_CLASS_QNAN = 0x05,
	180	};
1da177e4 LT	181
	182	/* exception numbers */
	183	#define IEEE754_INEXACT 0x01
	184	#define IEEE754_UNDERFLOW 0x02
	185	#define IEEE754_OVERFLOW 0x04
	186	#define IEEE754_ZERO_DIVIDE 0x08
	187	#define IEEE754_INVALID_OPERATION 0x10
	188
	189	/* cmp operators
	190	*/
	191	#define IEEE754_CLT 0x01
	192	#define IEEE754_CEQ 0x02
	193	#define IEEE754_CGT 0x04
	194	#define IEEE754_CUN 0x08
	195
	196	/* rounding mode
	197	*/
	198	#define IEEE754_RN 0 /* round to nearest */
	199	#define IEEE754_RZ 1 /* round toward zero */
	200	#define IEEE754_RD 2 /* round toward -Infinity */
	201	#define IEEE754_RU 3 /* round toward +Infinity */
	202
	203	/* other naming */
	204	#define IEEE754_RM IEEE754_RD
	205	#define IEEE754_RP IEEE754_RU
	206
	207	/* "normal" comparisons
	208	*/
2209bcb1	209	static inline int ieee754sp_eq(union ieee754sp x, union ieee754sp y)
1da177e4 LT	210	{
	211	return ieee754sp_cmp(x, y, IEEE754_CEQ, 0);
	212	}
	213
2209bcb1	214	static inline int ieee754sp_ne(union ieee754sp x, union ieee754sp y)
1da177e4 LT	215	{
	216	return ieee754sp_cmp(x, y,
	217	IEEE754_CLT \| IEEE754_CGT \| IEEE754_CUN, 0);
	218	}
	219
2209bcb1	220	static inline int ieee754sp_lt(union ieee754sp x, union ieee754sp y)
1da177e4 LT	221	{
	222	return ieee754sp_cmp(x, y, IEEE754_CLT, 0);
	223	}
	224
2209bcb1	225	static inline int ieee754sp_le(union ieee754sp x, union ieee754sp y)
1da177e4 LT	226	{
	227	return ieee754sp_cmp(x, y, IEEE754_CLT \| IEEE754_CEQ, 0);
	228	}
	229
2209bcb1	230	static inline int ieee754sp_gt(union ieee754sp x, union ieee754sp y)
1da177e4 LT	231	{
	232	return ieee754sp_cmp(x, y, IEEE754_CGT, 0);
	233	}
	234
	235
2209bcb1	236	static inline int ieee754sp_ge(union ieee754sp x, union ieee754sp y)
1da177e4 LT	237	{
	238	return ieee754sp_cmp(x, y, IEEE754_CGT \| IEEE754_CEQ, 0);
	239	}
	240
2209bcb1	241	static inline int ieee754dp_eq(union ieee754dp x, union ieee754dp y)
1da177e4 LT	242	{
	243	return ieee754dp_cmp(x, y, IEEE754_CEQ, 0);
	244	}
	245
2209bcb1	246	static inline int ieee754dp_ne(union ieee754dp x, union ieee754dp y)
1da177e4 LT	247	{
	248	return ieee754dp_cmp(x, y,
	249	IEEE754_CLT \| IEEE754_CGT \| IEEE754_CUN, 0);
	250	}
	251
2209bcb1	252	static inline int ieee754dp_lt(union ieee754dp x, union ieee754dp y)
1da177e4 LT	253	{
	254	return ieee754dp_cmp(x, y, IEEE754_CLT, 0);
	255	}
	256
2209bcb1	257	static inline int ieee754dp_le(union ieee754dp x, union ieee754dp y)
1da177e4 LT	258	{
	259	return ieee754dp_cmp(x, y, IEEE754_CLT \| IEEE754_CEQ, 0);
	260	}
	261
2209bcb1	262	static inline int ieee754dp_gt(union ieee754dp x, union ieee754dp y)
1da177e4 LT	263	{
	264	return ieee754dp_cmp(x, y, IEEE754_CGT, 0);
	265	}
	266
2209bcb1	267	static inline int ieee754dp_ge(union ieee754dp x, union ieee754dp y)
1da177e4 LT	268	{
	269	return ieee754dp_cmp(x, y, IEEE754_CGT \| IEEE754_CEQ, 0);
	270	}
	271
	272
cd21dfcf RB	273	/*
	274	* Like strtod
	275	*/
2209bcb1 RB	276	union ieee754dp ieee754dp_fstr(const char s, char *endp);
2209bcb1 RB	277	char *ieee754dp_tstr(union ieee754dp x, int prec, int fmt, int af);
1da177e4 LT	278
1da177e4 LT	279
cd21dfcf RB	280	/*
	281	* The control status register
	282	*/
	283	struct _ieee754_csr {
f80cc08d RB	284	__BITFIELD_FIELD(unsigned pad0:7,
	285	__BITFIELD_FIELD(unsigned nod:1, /* set 1 for no denormalised numbers */
	286	__BITFIELD_FIELD(unsigned c:1, /* condition */
	287	__BITFIELD_FIELD(unsigned pad1:5,
	288	__BITFIELD_FIELD(unsigned cx:6, /* exceptions this operation */
	289	__BITFIELD_FIELD(unsigned mx:5, /* exception enable mask */
	290	__BITFIELD_FIELD(unsigned sx:5, /* exceptions total */
	291	__BITFIELD_FIELD(unsigned rm:2, /* current rounding mode */
	292	;))))))))
1da177e4	293	};
eae89076	294	#define ieee754_csr ((struct _ieee754_csr )(&current->thread.fpu.fcr31))
1da177e4	295
cd21dfcf	296	static inline unsigned ieee754_getrm(void)
1da177e4 LT	297	{
	298	return (ieee754_csr.rm);
	299	}
cd21dfcf	300	static inline unsigned ieee754_setrm(unsigned rm)
1da177e4 LT	301	{
	302	return (ieee754_csr.rm = rm);
	303	}
	304
	305	/*
	306	* get current exceptions
	307	*/
cd21dfcf	308	static inline unsigned ieee754_getcx(void)
1da177e4 LT	309	{
	310	return (ieee754_csr.cx);
	311	}
	312
	313	/* test for current exception condition
	314	*/
cd21dfcf	315	static inline int ieee754_cxtest(unsigned n)
1da177e4 LT	316	{
	317	return (ieee754_csr.cx & n);
	318	}
	319
	320	/*
	321	* get sticky exceptions
	322	*/
cd21dfcf	323	static inline unsigned ieee754_getsx(void)
1da177e4 LT	324	{
	325	return (ieee754_csr.sx);
	326	}
	327
	328	/* clear sticky conditions
	329	*/
cd21dfcf	330	static inline unsigned ieee754_clrsx(void)
1da177e4 LT	331	{
	332	return (ieee754_csr.sx = 0);
	333	}
	334
	335	/* test for sticky exception condition
	336	*/
cd21dfcf	337	static inline int ieee754_sxtest(unsigned n)
1da177e4 LT	338	{
	339	return (ieee754_csr.sx & n);
	340	}
	341
	342	/* debugging */
2209bcb1 RB	343	union ieee754sp ieee754sp_dump(char *s, union ieee754sp x);
2209bcb1 RB	344	union ieee754dp ieee754dp_dump(char *s, union ieee754dp x);
1da177e4 LT	345
	346	#define IEEE754_SPCVAL_PZERO 0
	347	#define IEEE754_SPCVAL_NZERO 1
	348	#define IEEE754_SPCVAL_PONE 2
	349	#define IEEE754_SPCVAL_NONE 3
	350	#define IEEE754_SPCVAL_PTEN 4
	351	#define IEEE754_SPCVAL_NTEN 5
	352	#define IEEE754_SPCVAL_PINFINITY 6
	353	#define IEEE754_SPCVAL_NINFINITY 7
	354	#define IEEE754_SPCVAL_INDEF 8
	355	#define IEEE754_SPCVAL_PMAX 9 /* +max norm */
	356	#define IEEE754_SPCVAL_NMAX 10 /* -max norm */
	357	#define IEEE754_SPCVAL_PMIN 11 /* +min norm */
	358	#define IEEE754_SPCVAL_NMIN 12 /* +min norm */
	359	#define IEEE754_SPCVAL_PMIND 13 /* +min denorm */
	360	#define IEEE754_SPCVAL_NMIND 14 /* +min denorm */
	361	#define IEEE754_SPCVAL_P1E31 15 /* + 1.0e31 */
	362	#define IEEE754_SPCVAL_P1E63 16 /* + 1.0e63 */
	363
23708818 RB	364	extern const struct ieee754dp_const __ieee754dp_spcvals[];
23708818 RB	365	extern const struct ieee754sp_const __ieee754sp_spcvals[];
2209bcb1 RB	366	#define ieee754dp_spcvals ((const union ieee754dp *)__ieee754dp_spcvals)
2209bcb1 RB	367	#define ieee754sp_spcvals ((const union ieee754sp *)__ieee754sp_spcvals)
1da177e4	368
cd21dfcf RB	369	/*
	370	* Return infinity with given sign
	371	*/
	372	#define ieee754dp_inf(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PINFINITY+(sn)])
	373	#define ieee754dp_zero(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PZERO+(sn)])
	374	#define ieee754dp_one(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PONE+(sn)])
	375	#define ieee754dp_ten(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PTEN+(sn)])
	376	#define ieee754dp_indef() (ieee754dp_spcvals[IEEE754_SPCVAL_INDEF])
	377	#define ieee754dp_max(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PMAX+(sn)])
	378	#define ieee754dp_min(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PMIN+(sn)])
	379	#define ieee754dp_mind(sn) (ieee754dp_spcvals[IEEE754_SPCVAL_PMIND+(sn)])
	380	#define ieee754dp_1e31() (ieee754dp_spcvals[IEEE754_SPCVAL_P1E31])
	381	#define ieee754dp_1e63() (ieee754dp_spcvals[IEEE754_SPCVAL_P1E63])
	382
	383	#define ieee754sp_inf(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PINFINITY+(sn)])
	384	#define ieee754sp_zero(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PZERO+(sn)])
	385	#define ieee754sp_one(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PONE+(sn)])
	386	#define ieee754sp_ten(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PTEN+(sn)])
	387	#define ieee754sp_indef() (ieee754sp_spcvals[IEEE754_SPCVAL_INDEF])
	388	#define ieee754sp_max(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PMAX+(sn)])
	389	#define ieee754sp_min(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PMIN+(sn)])
	390	#define ieee754sp_mind(sn) (ieee754sp_spcvals[IEEE754_SPCVAL_PMIND+(sn)])
	391	#define ieee754sp_1e31() (ieee754sp_spcvals[IEEE754_SPCVAL_P1E31])
	392	#define ieee754sp_1e63() (ieee754sp_spcvals[IEEE754_SPCVAL_P1E63])
	393
	394	/*
	395	* Indefinite integer value
	396	*/
1da177e4 LT	397	#define ieee754si_indef() INT_MAX
	398	#ifdef LONG_LONG_MAX
	399	#define ieee754di_indef() LONG_LONG_MAX
	400	#else
	401	#define ieee754di_indef() ((s64)(~0ULL>>1))
	402	#endif
	403
	404	/* IEEE exception context, passed to handler */
	405	struct ieee754xctx {
	406	const char op; / operation name */
	407	int rt; /* result type */
	408	union {
2209bcb1 RB	409	union ieee754sp sp; /* single precision */
2209bcb1 RB	410	union ieee754dp dp; /* double precision */
1da177e4 LT	411	int si; /* standard signed integer (32bits) */
	412	s64 di; /* extended signed integer (64bits) */
	413	} rv; /* default result format implied by op */
	414	va_list ap;
	415	};
	416
	417	/* result types for xctx.rt */
	418	#define IEEE754_RT_SP 0
	419	#define IEEE754_RT_DP 1
	420	#define IEEE754_RT_XP 2
	421	#define IEEE754_RT_SI 3
	422	#define IEEE754_RT_DI 4
	423
cae55066	424	extern void __cold ieee754_xcpt(struct ieee754xctx *xcp);
1da177e4 LT	425
	426	/* compat */
	427	#define ieee754dp_fix(x) ieee754dp_tint(x)
	428	#define ieee754sp_fix(x) ieee754sp_tint(x)
a663bf90 RB	429
a663bf90 RB	430	#endif /* __ARCH_MIPS_MATH_EMU_IEEE754_H */