[mirror_ubuntu-bionic-kernel.git] / arch / mips / math-emu / ieee754int.h

/*
 * IEEE754 floating point
 * common internal header file
 */
/*
 * 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.,
 *  51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA.
 */
#ifndef __IEEE754INT_H
#define __IEEE754INT_H

#include "ieee754.h"

#define CLPAIR(x, y)	((x)*6+(y))

enum maddf_flags {
	MADDF_NEGATE_PRODUCT	= 1 << 0,
};

static inline void ieee754_clearcx(void)
{
	ieee754_csr.cx = 0;
}

static inline void ieee754_setcx(const unsigned int flags)
{
	ieee754_csr.cx |= flags;
	ieee754_csr.sx |= flags;
}

static inline int ieee754_setandtestcx(const unsigned int x)
{
	ieee754_setcx(x);

	return ieee754_csr.mx & x;
}

static inline int ieee754_class_nan(int xc)
{
	return xc >= IEEE754_CLASS_SNAN;
}

#define COMPXSP \
	unsigned int xm; int xe; int xs __maybe_unused; int xc

#define COMPYSP \
	unsigned int ym; int ye; int ys; int yc

#define COMPZSP \
	unsigned int zm; int ze; int zs; int zc

#define EXPLODESP(v, vc, vs, ve, vm)					\
{									\
	vs = SPSIGN(v);							\
	ve = SPBEXP(v);							\
	vm = SPMANT(v);							\
	if (ve == SP_EMAX+1+SP_EBIAS) {					\
		if (vm == 0)						\
			vc = IEEE754_CLASS_INF;				\
		else if (ieee754_csr.nan2008 ^ !(vm & SP_MBIT(SP_FBITS - 1))) \
			vc = IEEE754_CLASS_QNAN;			\
		else							\
			vc = IEEE754_CLASS_SNAN;			\
	} else if (ve == SP_EMIN-1+SP_EBIAS) {				\
		if (vm) {						\
			ve = SP_EMIN;					\
			vc = IEEE754_CLASS_DNORM;			\
		} else							\
			vc = IEEE754_CLASS_ZERO;			\
	} else {							\
		ve -= SP_EBIAS;						\
		vm |= SP_HIDDEN_BIT;					\
		vc = IEEE754_CLASS_NORM;				\
	}								\
}
#define EXPLODEXSP EXPLODESP(x, xc, xs, xe, xm)
#define EXPLODEYSP EXPLODESP(y, yc, ys, ye, ym)
#define EXPLODEZSP EXPLODESP(z, zc, zs, ze, zm)


#define COMPXDP \
	u64 xm; int xe; int xs __maybe_unused; int xc

#define COMPYDP \
	u64 ym; int ye; int ys; int yc

#define COMPZDP \
	u64 zm; int ze; int zs; int zc

#define EXPLODEDP(v, vc, vs, ve, vm)					\
{									\
	vm = DPMANT(v);							\
	vs = DPSIGN(v);							\
	ve = DPBEXP(v);							\
	if (ve == DP_EMAX+1+DP_EBIAS) {					\
		if (vm == 0)						\
			vc = IEEE754_CLASS_INF;				\
		else if (ieee754_csr.nan2008 ^ !(vm & DP_MBIT(DP_FBITS - 1))) \
			vc = IEEE754_CLASS_QNAN;			\
		else							\
			vc = IEEE754_CLASS_SNAN;			\
	} else if (ve == DP_EMIN-1+DP_EBIAS) {				\
		if (vm) {						\
			ve = DP_EMIN;					\
			vc = IEEE754_CLASS_DNORM;			\
		} else							\
			vc = IEEE754_CLASS_ZERO;			\
	} else {							\
		ve -= DP_EBIAS;						\
		vm |= DP_HIDDEN_BIT;					\
		vc = IEEE754_CLASS_NORM;				\
	}								\
}
#define EXPLODEXDP EXPLODEDP(x, xc, xs, xe, xm)
#define EXPLODEYDP EXPLODEDP(y, yc, ys, ye, ym)
#define EXPLODEZDP EXPLODEDP(z, zc, zs, ze, zm)

#define FLUSHDP(v, vc, vs, ve, vm)					\
	if (vc==IEEE754_CLASS_DNORM) {					\
		if (ieee754_csr.nod) {					\
			ieee754_setcx(IEEE754_INEXACT);			\
			vc = IEEE754_CLASS_ZERO;			\
			ve = DP_EMIN-1+DP_EBIAS;			\
			vm = 0;						\
			v = ieee754dp_zero(vs);				\
		}							\
	}

#define FLUSHSP(v, vc, vs, ve, vm)					\
	if (vc==IEEE754_CLASS_DNORM) {					\
		if (ieee754_csr.nod) {					\
			ieee754_setcx(IEEE754_INEXACT);			\
			vc = IEEE754_CLASS_ZERO;			\
			ve = SP_EMIN-1+SP_EBIAS;			\
			vm = 0;						\
			v = ieee754sp_zero(vs);				\
		}							\
	}

#define FLUSHXDP FLUSHDP(x, xc, xs, xe, xm)
#define FLUSHYDP FLUSHDP(y, yc, ys, ye, ym)
#define FLUSHZDP FLUSHDP(z, zc, zs, ze, zm)
#define FLUSHXSP FLUSHSP(x, xc, xs, xe, xm)
#define FLUSHYSP FLUSHSP(y, yc, ys, ye, ym)
#define FLUSHZSP FLUSHSP(z, zc, zs, ze, zm)

#endif /* __IEEE754INT_H  */
Commit	Line	Data
1da177e4 LT	1	/*
	2	* IEEE754 floating point
	3	* common internal header file
	4	*/
	5	/*
	6	* MIPS floating point support
	7	* Copyright (C) 1994-2000 Algorithmics Ltd.
1da177e4	8	*
1da177e4 LT	9	* This program is free software; you can distribute it and/or modify it
	10	* under the terms of the GNU General Public License (Version 2) as
	11	* published by the Free Software Foundation.
	12	*
	13	* This program is distributed in the hope it will be useful, but WITHOUT
	14	* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
	15	* FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
	16	* for more details.
	17	*
	18	* You should have received a copy of the GNU General Public License along
	19	* with this program; if not, write to the Free Software Foundation, Inc.,
3f7cac41	20	* 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA.
1da177e4	21	*/
bee16535 RB	22	#ifndef __IEEE754INT_H
bee16535 RB	23	#define __IEEE754INT_H
1da177e4 LT	24
	25	#include "ieee754.h"
	26
21a151d8	27	#define CLPAIR(x, y) ((x)*6+(y))
1da177e4	28
ae11c061 AM	29	enum maddf_flags {
	30	MADDF_NEGATE_PRODUCT = 1 << 0,
	31	};
	32
9e8bad1f RB	33	static inline void ieee754_clearcx(void)
	34	{
	35	ieee754_csr.cx = 0;
	36	}
1da177e4	37
9e8bad1f RB	38	static inline void ieee754_setcx(const unsigned int flags)
	39	{
	40	ieee754_csr.cx \|= flags;
	41	ieee754_csr.sx \|= flags;
	42	}
1da177e4	43
9e8bad1f RB	44	static inline int ieee754_setandtestcx(const unsigned int x)
	45	{
	46	ieee754_setcx(x);
1da177e4	47
9e8bad1f RB	48	return ieee754_csr.mx & x;
9e8bad1f RB	49	}
1da177e4	50
c9a10845 MR	51	static inline int ieee754_class_nan(int xc)
	52	{
	53	return xc >= IEEE754_CLASS_SNAN;
	54	}
	55
1da177e4	56	#define COMPXSP \
a58f85b5	57	unsigned int xm; int xe; int xs __maybe_unused; int xc
1da177e4 LT	58
1da177e4 LT	59	#define COMPYSP \
a58f85b5	60	unsigned int ym; int ye; int ys; int yc
47fa0c02	61
e2d11e1a	62	#define COMPZSP \
a58f85b5	63	unsigned int zm; int ze; int zs; int zc
e2d11e1a	64
47fa0c02 RB	65	#define EXPLODESP(v, vc, vs, ve, vm) \
	66	{ \
	67	vs = SPSIGN(v); \
	68	ve = SPBEXP(v); \
	69	vm = SPMANT(v); \
	70	if (ve == SP_EMAX+1+SP_EBIAS) { \
	71	if (vm == 0) \
	72	vc = IEEE754_CLASS_INF; \
90d53a91	73	else if (ieee754_csr.nan2008 ^ !(vm & SP_MBIT(SP_FBITS - 1))) \
16645743	74	vc = IEEE754_CLASS_QNAN; \
90d53a91 MR	75	else \
90d53a91 MR	76	vc = IEEE754_CLASS_SNAN; \
47fa0c02 RB	77	} else if (ve == SP_EMIN-1+SP_EBIAS) { \
	78	if (vm) { \
	79	ve = SP_EMIN; \
	80	vc = IEEE754_CLASS_DNORM; \
	81	} else \
	82	vc = IEEE754_CLASS_ZERO; \
	83	} else { \
	84	ve -= SP_EBIAS; \
	85	vm \|= SP_HIDDEN_BIT; \
	86	vc = IEEE754_CLASS_NORM; \
	87	} \
1da177e4	88	}
21a151d8 RB	89	#define EXPLODEXSP EXPLODESP(x, xc, xs, xe, xm)
21a151d8 RB	90	#define EXPLODEYSP EXPLODESP(y, yc, ys, ye, ym)
e2d11e1a	91	#define EXPLODEZSP EXPLODESP(z, zc, zs, ze, zm)
1da177e4 LT	92
	93
	94	#define COMPXDP \
47fa0c02	95	u64 xm; int xe; int xs __maybe_unused; int xc
1da177e4 LT	96
1da177e4 LT	97	#define COMPYDP \
47fa0c02 RB	98	u64 ym; int ye; int ys; int yc
47fa0c02 RB	99
e2d11e1a PB	100	#define COMPZDP \
	101	u64 zm; int ze; int zs; int zc
	102
47fa0c02 RB	103	#define EXPLODEDP(v, vc, vs, ve, vm) \
	104	{ \
	105	vm = DPMANT(v); \
	106	vs = DPSIGN(v); \
	107	ve = DPBEXP(v); \
	108	if (ve == DP_EMAX+1+DP_EBIAS) { \
	109	if (vm == 0) \
	110	vc = IEEE754_CLASS_INF; \
90d53a91	111	else if (ieee754_csr.nan2008 ^ !(vm & DP_MBIT(DP_FBITS - 1))) \
47fa0c02	112	vc = IEEE754_CLASS_QNAN; \
90d53a91 MR	113	else \
90d53a91 MR	114	vc = IEEE754_CLASS_SNAN; \
47fa0c02 RB	115	} else if (ve == DP_EMIN-1+DP_EBIAS) { \
	116	if (vm) { \
	117	ve = DP_EMIN; \
	118	vc = IEEE754_CLASS_DNORM; \
16645743 MC	119	} else \
16645743 MC	120	vc = IEEE754_CLASS_ZERO; \
47fa0c02 RB	121	} else { \
	122	ve -= DP_EBIAS; \
	123	vm \|= DP_HIDDEN_BIT; \
	124	vc = IEEE754_CLASS_NORM; \
	125	} \
1da177e4	126	}
21a151d8 RB	127	#define EXPLODEXDP EXPLODEDP(x, xc, xs, xe, xm)
21a151d8 RB	128	#define EXPLODEYDP EXPLODEDP(y, yc, ys, ye, ym)
e2d11e1a	129	#define EXPLODEZDP EXPLODEDP(z, zc, zs, ze, zm)
1da177e4	130
47fa0c02 RB	131	#define FLUSHDP(v, vc, vs, ve, vm) \
	132	if (vc==IEEE754_CLASS_DNORM) { \
	133	if (ieee754_csr.nod) { \
9e8bad1f	134	ieee754_setcx(IEEE754_INEXACT); \
47fa0c02 RB	135	vc = IEEE754_CLASS_ZERO; \
	136	ve = DP_EMIN-1+DP_EBIAS; \
	137	vm = 0; \
	138	v = ieee754dp_zero(vs); \
	139	} \
1da177e4 LT	140	}
1da177e4 LT	141
47fa0c02 RB	142	#define FLUSHSP(v, vc, vs, ve, vm) \
	143	if (vc==IEEE754_CLASS_DNORM) { \
	144	if (ieee754_csr.nod) { \
9e8bad1f	145	ieee754_setcx(IEEE754_INEXACT); \
47fa0c02 RB	146	vc = IEEE754_CLASS_ZERO; \
	147	ve = SP_EMIN-1+SP_EBIAS; \
	148	vm = 0; \
	149	v = ieee754sp_zero(vs); \
	150	} \
1da177e4 LT	151	}
1da177e4 LT	152
21a151d8 RB	153	#define FLUSHXDP FLUSHDP(x, xc, xs, xe, xm)
21a151d8 RB	154	#define FLUSHYDP FLUSHDP(y, yc, ys, ye, ym)
e2d11e1a	155	#define FLUSHZDP FLUSHDP(z, zc, zs, ze, zm)
21a151d8 RB	156	#define FLUSHXSP FLUSHSP(x, xc, xs, xe, xm)
21a151d8 RB	157	#define FLUSHYSP FLUSHSP(y, yc, ys, ye, ym)
e2d11e1a	158	#define FLUSHZSP FLUSHSP(z, zc, zs, ze, zm)
bee16535 RB	159
bee16535 RB	160	#endif /* __IEEE754INT_H */