| | 1 | /*- |
| | 2 | * Copyright (c) 2003 The NetBSD Foundation, Inc. |
| | 3 | * All rights reserved. |
| | 4 | * |
| | 5 | * This code is derived from software contributed to The NetBSD Foundation |
| | 6 | * by Klaus Klein. |
| | 7 | * |
| | 8 | * Redistribution and use in source and binary forms, with or without |
| | 9 | * modification, are permitted provided that the following conditions |
| | 10 | * are met: |
| | 11 | * 1. Redistributions of source code must retain the above copyright |
| | 12 | * notice, this list of conditions and the following disclaimer. |
| | 13 | * 2. Redistributions in binary form must reproduce the above copyright |
| | 14 | * notice, this list of conditions and the following disclaimer in the |
| | 15 | * documentation and/or other materials provided with the distribution. |
| | 16 | * 3. All advertising materials mentioning features or use of this software |
| | 17 | * must display the following acknowledgement: |
| | 18 | * This product includes software developed by the NetBSD |
| | 19 | * Foundation, Inc. and its contributors. |
| | 20 | * 4. Neither the name of The NetBSD Foundation nor the names of its |
| | 21 | * contributors may be used to endorse or promote products derived |
| | 22 | * from this software without specific prior written permission. |
| | 23 | * |
| | 24 | * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS |
| | 25 | * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED |
| | 26 | * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR |
| | 27 | * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS |
| | 28 | * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR |
| | 29 | * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF |
| | 30 | * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS |
| | 31 | * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN |
| | 32 | * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) |
| | 33 | * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE |
| | 34 | * POSSIBILITY OF SUCH DAMAGE. |
| | 35 | * |
| | 36 | * $NetBSD: fpclassifyl.c,v 1.2 2004/01/18 19:33:01 matt Exp $ |
| | 37 | */ |
| | 38 | |
| | 39 | #include <machine/ieee.h> |
| | 40 | #include <assert.h> |
| | 41 | #include <math.h> |
| | 42 | |
| | 43 | /* |
| | 44 | * 7.12.3.1 fpclassify - classify real floating type |
| | 45 | * IEEE 754 compatible 80-bit extended-precision Intel 386 version |
| | 46 | */ |
| | 47 | int |
| | 48 | __fpclassifyl(long double x) |
| | 49 | { |
| | 50 | union ieee_ext_u u; |
| | 51 | |
| | 52 | u.extu_ld = x; |
| | 53 | |
| | 54 | _DIAGASSERT(u.extu_ext.ext_exp == 0 || u.extu_ext.ext_int == 1); |
| | 55 | |
| | 56 | if (u.extu_ext.ext_exp == 0) { |
| | 57 | if (u.extu_ext.ext_frach == 0 && 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 == 0 && u.extu_ext.ext_fracl == 0) |
| | 63 | return FP_INFINITE; |
| | 64 | else |
| | 65 | return FP_NAN; |
| | 66 | } |
| | 67 | |
| | 68 | return FP_NORMAL; |
| | 69 | } |
projects / dragonfly.git / blame_incremental