m68k/gen/fpclassifyl.c

1.3.8.2  martin /*	$NetBSD: fpclassifyl.c,v 1.3.8.2 2008/04/28 20:22:57 martin Exp $	*/
1.3.8.2  martin
1.3.8.2  martin /*-
1.3.8.2  martin  * Copyright (c) 2003 The NetBSD Foundation, Inc.
1.3.8.2  martin  * All rights reserved.
1.3.8.2  martin  *
1.3.8.2  martin  * This code is derived from software contributed to The NetBSD Foundation
1.3.8.2  martin  * by Klaus Klein.
1.3.8.2  martin  *
1.3.8.2  martin  * Redistribution and use in source and binary forms, with or without
1.3.8.2  martin  * modification, are permitted provided that the following conditions
1.3.8.2  martin  * are met:
1.3.8.2  martin  * 1. Redistributions of source code must retain the above copyright
1.3.8.2  martin  *    notice, this list of conditions and the following disclaimer.
1.3.8.2  martin  * 2. Redistributions in binary form must reproduce the above copyright
1.3.8.2  martin  *    notice, this list of conditions and the following disclaimer in the
1.3.8.2  martin  *    documentation and/or other materials provided with the distribution.
1.3.8.2  martin  *
1.3.8.2  martin  * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
1.3.8.2  martin  * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
1.3.8.2  martin  * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
1.3.8.2  martin  * PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
1.3.8.2  martin  * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
1.3.8.2  martin  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
1.3.8.2  martin  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
1.3.8.2  martin  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
1.3.8.2  martin  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
1.3.8.2  martin  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
1.3.8.2  martin  * POSSIBILITY OF SUCH DAMAGE.
1.3.8.2  martin  */
1.3.8.2  martin
1.3.8.2  martin #include <sys/cdefs.h>
1.3.8.2  martin #if defined(LIBC_SCCS) && !defined(lint)
1.3.8.2  martin __RCSID("$NetBSD: fpclassifyl.c,v 1.3.8.2 2008/04/28 20:22:57 martin Exp $");
1.3.8.2  martin #endif
1.3.8.2  martin
1.3.8.2  martin #include <machine/ieee.h>
1.3.8.2  martin #include <assert.h>
1.3.8.2  martin #include <math.h>
1.3.8.2  martin
1.3.8.2  martin /*
1.3.8.2  martin  * 7.12.3.1 fpclassify - classify real floating type
1.3.8.2  martin  *          IEEE 754 compatible 80-bit extended-precision Motorola 68k version
1.3.8.2  martin  */
1.3.8.2  martin int
1.3.8.2  martin __fpclassifyl(long double x)
1.3.8.2  martin {
1.3.8.2  martin 	union ieee_ext_u u;
1.3.8.2  martin
1.3.8.2  martin 	u.extu_ld = x;
1.3.8.2  martin
1.3.8.2  martin 	_DIAGASSERT(u.extu_ext.ext_zero == 0);
1.3.8.2  martin
1.3.8.2  martin 	if (u.extu_ext.ext_exp == 0 &&
1.3.8.2  martin 	    (u.extu_ext.ext_frach & 0x80000000) == 0) {
1.3.8.2  martin 		if ((u.extu_ext.ext_frach & 0x7fffffff) == 0 &&
1.3.8.2  martin 		    u.extu_ext.ext_fracl == 0)
1.3.8.2  martin 			return FP_ZERO;
1.3.8.2  martin 		else
1.3.8.2  martin 			return FP_SUBNORMAL;
1.3.8.2  martin 	} else if (u.extu_ext.ext_exp == EXT_EXP_INFNAN) {
1.3.8.2  martin 		if ((u.extu_ext.ext_frach & 0x7fffffff) == 0 &&
1.3.8.2  martin 		    u.extu_ext.ext_fracl == 0)
1.3.8.2  martin 			return FP_INFINITE;
1.3.8.2  martin 		else
1.3.8.2  martin 			return FP_NAN;
1.3.8.2  martin 	}
1.3.8.2  martin
1.3.8.2  martin 	return FP_NORMAL;
1.3.8.2  martin }