[net-next-2.6.git] / arch / mips / math-emu / sp_simple.c

/* IEEE754 floating point arithmetic
 * single precision
 */
/*
 * MIPS floating point support
 * Copyright (C) 1994-2000 Algorithmics Ltd.
 * http://www.algor.co.uk
 *
 * ########################################################################
 *
 *  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.
 *
 * ########################################################################
 */


#include "ieee754sp.h"

int ieee754sp_finite(ieee754sp x)
{
	return SPBEXP(x) != SP_EMAX + 1 + SP_EBIAS;
}

ieee754sp ieee754sp_copysign(ieee754sp x, ieee754sp y)
{
	CLEARCX;
	SPSIGN(x) = SPSIGN(y);
	return x;
}


ieee754sp ieee754sp_neg(ieee754sp x)
{
	COMPXSP;

	EXPLODEXSP;
	CLEARCX;
	FLUSHXSP;

	/*
	 * Invert the sign ALWAYS to prevent an endless recursion on
	 * pow() in libc.
	 */
	/* quick fix up */
	SPSIGN(x) ^= 1;

	if (xc == IEEE754_CLASS_SNAN) {
		ieee754sp y = ieee754sp_indef();
		SETCX(IEEE754_INVALID_OPERATION);
		SPSIGN(y) = SPSIGN(x);
		return ieee754sp_nanxcpt(y, "neg");
	}

	return x;
}


ieee754sp ieee754sp_abs(ieee754sp x)
{
	COMPXSP;

	EXPLODEXSP;
	CLEARCX;
	FLUSHXSP;

	/* Clear sign ALWAYS, irrespective of NaN */
	SPSIGN(x) = 0;

	if (xc == IEEE754_CLASS_SNAN) {
		SETCX(IEEE754_INVALID_OPERATION);
		return ieee754sp_nanxcpt(ieee754sp_indef(), "abs");
	}

	return x;
}
Commit	Line	Data
1da177e4 LT	1	/* IEEE754 floating point arithmetic
	2	* single precision
	3	*/
	4	/*
	5	* MIPS floating point support
	6	* Copyright (C) 1994-2000 Algorithmics Ltd.
	7	* http://www.algor.co.uk
	8	*
	9	* ########################################################################
	10	*
	11	* This program is free software; you can distribute it and/or modify it
	12	* under the terms of the GNU General Public License (Version 2) as
	13	* published by the Free Software Foundation.
	14	*
	15	* This program is distributed in the hope it will be useful, but WITHOUT
	16	* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
	17	* FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
	18	* for more details.
	19	*
	20	* You should have received a copy of the GNU General Public License along
	21	* with this program; if not, write to the Free Software Foundation, Inc.,
	22	* 59 Temple Place - Suite 330, Boston MA 02111-1307, USA.
	23	*
	24	* ########################################################################
	25	*/
	26
	27
	28	#include "ieee754sp.h"
	29
	30	int ieee754sp_finite(ieee754sp x)
	31	{
	32	return SPBEXP(x) != SP_EMAX + 1 + SP_EBIAS;
	33	}
	34
	35	ieee754sp ieee754sp_copysign(ieee754sp x, ieee754sp y)
	36	{
	37	CLEARCX;
	38	SPSIGN(x) = SPSIGN(y);
	39	return x;
	40	}
	41
	42
	43	ieee754sp ieee754sp_neg(ieee754sp x)
	44	{
	45	COMPXSP;
	46
	47	EXPLODEXSP;
	48	CLEARCX;
	49	FLUSHXSP;
	50
c5033d78 AN	51	/*
	52	* Invert the sign ALWAYS to prevent an endless recursion on
	53	* pow() in libc.
	54	*/
	55	/* quick fix up */
	56	SPSIGN(x) ^= 1;
	57
1da177e4	58	if (xc == IEEE754_CLASS_SNAN) {
c5033d78	59	ieee754sp y = ieee754sp_indef();
1da177e4	60	SETCX(IEEE754_INVALID_OPERATION);
c5033d78 AN	61	SPSIGN(y) = SPSIGN(x);
c5033d78 AN	62	return ieee754sp_nanxcpt(y, "neg");
1da177e4 LT	63	}
1da177e4 LT	64
1da177e4 LT	65	return x;
	66	}
	67
	68
	69	ieee754sp ieee754sp_abs(ieee754sp x)
	70	{
	71	COMPXSP;
	72
	73	EXPLODEXSP;
	74	CLEARCX;
	75	FLUSHXSP;
	76
cea2be44 NS	77	/* Clear sign ALWAYS, irrespective of NaN */
	78	SPSIGN(x) = 0;
	79
1da177e4	80	if (xc == IEEE754_CLASS_SNAN) {
e48682dd	81	SETCX(IEEE754_INVALID_OPERATION);
1da177e4 LT	82	return ieee754sp_nanxcpt(ieee754sp_indef(), "abs");
	83	}
	84
1da177e4 LT	85	return x;
1da177e4 LT	86	}