lib/libc/amd64/gen/fpclassifyl.c

   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  * $DragonFly: src/lib/libc/amd64/gen/fpclassifyl.c,v 1.1 2006/07/27 00:46:57 corecode Exp $
  38  */
  39
  40 #include <sys/cdefs.h>
  41
  42 #include <machine/ieee.h>
  43 #include <assert.h>
  44 #include <math.h>
  45
  46 /*
  47  * 7.12.3.1 fpclassify - classify real floating type
  48  *          IEEE 754 compatible 80-bit extended-precision Intel 386 version
  49  */
  50 int
  51 __fpclassifyl(long double x)
  52 {
  53         union ieee_ext_u u;
  54
  55         u.extu_ld = x;
  56
  57         _DIAGASSERT(u.extu_ext.ext_exp == 0 || u.extu_ext.ext_int == 1);
  58
  59         if (u.extu_ext.ext_exp == 0) {
  60                 if (u.extu_ext.ext_frach == 0 && u.extu_ext.ext_fracl == 0)
  61                         return FP_ZERO;
  62                 else
  63                         return FP_SUBNORMAL;
  64         } else if (u.extu_ext.ext_exp == EXT_EXP_INFNAN) {
  65                 if (u.extu_ext.ext_frach == 0 && u.extu_ext.ext_fracl == 0)
  66                         return FP_INFINITE;
  67                 else
  68                         return FP_NAN;
  69         }
  70
  71         return FP_NORMAL;
  72 }