sparc64/gen/fpclassifyl.c

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