Home | History | Annotate | Line # | Download | only in gen
fpclassifyd.c revision 1.3.8.2 
      1  1.3.8.2  martin /*	$NetBSD: fpclassifyd.c,v 1.3.8.2 2008/04/28 20:22:58 martin Exp $	*/
      2  1.3.8.2  martin 
      3  1.3.8.2  martin /*-
      4  1.3.8.2  martin  * Copyright (c) 2004 The NetBSD Foundation, Inc.
      5  1.3.8.2  martin  * All rights reserved.
      6  1.3.8.2  martin  *
      7  1.3.8.2  martin  * This code is derived from software contributed to The NetBSD Foundation
      8  1.3.8.2  martin  * by Klaus Klein.
      9  1.3.8.2  martin  *
     10  1.3.8.2  martin  * Redistribution and use in source and binary forms, with or without
     11  1.3.8.2  martin  * modification, are permitted provided that the following conditions
     12  1.3.8.2  martin  * are met:
     13  1.3.8.2  martin  * 1. Redistributions of source code must retain the above copyright
     14  1.3.8.2  martin  *    notice, this list of conditions and the following disclaimer.
     15  1.3.8.2  martin  * 2. Redistributions in binary form must reproduce the above copyright
     16  1.3.8.2  martin  *    notice, this list of conditions and the following disclaimer in the
     17  1.3.8.2  martin  *    documentation and/or other materials provided with the distribution.
     18  1.3.8.2  martin  *
     19  1.3.8.2  martin  * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
     20  1.3.8.2  martin  * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
     21  1.3.8.2  martin  * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
     22  1.3.8.2  martin  * PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
     23  1.3.8.2  martin  * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
     24  1.3.8.2  martin  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
     25  1.3.8.2  martin  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
     26  1.3.8.2  martin  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
     27  1.3.8.2  martin  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
     28  1.3.8.2  martin  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
     29  1.3.8.2  martin  * POSSIBILITY OF SUCH DAMAGE.
     30  1.3.8.2  martin  */
     31  1.3.8.2  martin 
     32  1.3.8.2  martin #include <sys/cdefs.h>
     33  1.3.8.2  martin #if defined(LIBC_SCCS) && !defined(lint)
     34  1.3.8.2  martin __RCSID("$NetBSD: fpclassifyd.c,v 1.3.8.2 2008/04/28 20:22:58 martin Exp $");
     35  1.3.8.2  martin #endif
     36  1.3.8.2  martin 
     37  1.3.8.2  martin #include <machine/vaxfp.h>
     38  1.3.8.2  martin #include <math.h>
     39  1.3.8.2  martin 
     40  1.3.8.2  martin /*
     41  1.3.8.2  martin  * 7.12.3.1 fpclassify - classify real floating type
     42  1.3.8.2  martin  *          VAX D_floating version
     43  1.3.8.2  martin  *
     44  1.3.8.2  martin  *          Implementation notes:
     45  1.3.8.2  martin  *          Reserved operand -> FP_ROP
     46  1.3.8.2  martin  *          True zero        -> FP_ZERO
     47  1.3.8.2  martin  *          Dirty zero       -> FP_DIRTY_ZERO
     48  1.3.8.2  martin  *          Finite           -> FP_NORMAL
     49  1.3.8.2  martin  */
     50  1.3.8.2  martin int
     51  1.3.8.2  martin __fpclassifyd(double x)
     52  1.3.8.2  martin {
     53  1.3.8.2  martin 	union vax_dfloating_u u;
     54  1.3.8.2  martin 
     55  1.3.8.2  martin 	u.dfltu_d = x;
     56  1.3.8.2  martin 
     57  1.3.8.2  martin 	if (u.dfltu_dflt.dflt_exp == 0) {
     58  1.3.8.2  martin 		if (u.dfltu_dflt.dflt_sign != 0)
     59  1.3.8.2  martin 			return FP_ROP;
     60  1.3.8.2  martin 		else if (u.dfltu_dflt.dflt_frach != 0 ||
     61  1.3.8.2  martin 			 u.dfltu_dflt.dflt_fracm != 0 ||
     62  1.3.8.2  martin 			 u.dfltu_dflt.dflt_fracl != 0)
     63  1.3.8.2  martin 			return FP_DIRTYZERO;
     64  1.3.8.2  martin 		else
     65  1.3.8.2  martin 			return FP_ZERO;
     66  1.3.8.2  martin 	} else {
     67  1.3.8.2  martin 		return FP_NORMAL;
     68  1.3.8.2  martin 	}
     69  1.3.8.2  martin }
     70