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