Sync usage with man page.
[netbsd-mini2440.git] / lib / libc / arch / m68k / gen / fpclassifyl.c
blob446622c118c37d985dfefb9f37f454957f22dcd6
1 /* $NetBSD: fpclassifyl.c,v 1.2 2007/02/17 17:50:52 mhitch Exp $ */
3 /*-
4 * Copyright (c) 2003 The NetBSD Foundation, Inc.
5 * All rights reserved.
7 * This code is derived from software contributed to The NetBSD Foundation
8 * by Klaus Klein.
10 * Redistribution and use in source and binary forms, with or without
11 * modification, are permitted provided that the following conditions
12 * are met:
13 * 1. Redistributions of source code must retain the above copyright
14 * notice, this list of conditions and the following disclaimer.
15 * 2. Redistributions in binary form must reproduce the above copyright
16 * notice, this list of conditions and the following disclaimer in the
17 * documentation and/or other materials provided with the distribution.
19 * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
20 * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
21 * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
22 * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
23 * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
24 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
25 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
26 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
27 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
28 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
29 * POSSIBILITY OF SUCH DAMAGE.
32 #include <sys/cdefs.h>
33 #if defined(LIBC_SCCS) && !defined(lint)
34 __RCSID("$NetBSD: fpclassifyl.c,v 1.2 2007/02/17 17:50:52 mhitch Exp $");
35 #endif
37 #include <machine/ieee.h>
38 #include <assert.h>
39 #include <math.h>
42 * 7.12.3.1 fpclassify - classify real floating type
43 * IEEE 754 compatible 80-bit extended-precision Motorola 68k version
45 int
46 __fpclassifyl(long double x)
48 union ieee_ext_u u;
50 u.extu_ld = x;
52 _DIAGASSERT(u.extu_ext.ext_zero == 0);
54 if (u.extu_ext.ext_exp == 0 &&
55 (u.extu_ext.ext_frach & 0x80000000) == 0) {
56 if ((u.extu_ext.ext_frach & 0x7fffffff) == 0 &&
57 u.extu_ext.ext_fracl == 0)
58 return FP_ZERO;
59 else
60 return FP_SUBNORMAL;
61 } else if (u.extu_ext.ext_exp == EXT_EXP_INFNAN) {
62 if ((u.extu_ext.ext_frach & 0x7fffffff) == 0 &&
63 u.extu_ext.ext_fracl == 0)
64 return FP_INFINITE;
65 else
66 return FP_NAN;
69 return FP_NORMAL;