Symaptic: os/ossrv/genericopenlibs/openenvcore/libm/src/fpclassify.c@260cb5ec6c19

     1 /*--------------------------------------------------------------------

     2  *© Portions copyright (c) 2006 Nokia Corporation.  All rights reserved.

     3  *--------------------------------------------------------------------

     4 */

     5 /*-

     6  * Copyright (c) 2003 Mike Barcroft <mike@FreeBSD.org>

     7  * Copyright (c) 2002, 2003 David Schultz <das@FreeBSD.ORG>

     8  * All rights reserved.

     9  *

    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.

    18  *

    19  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND

    20  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE

    21  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE

    22  * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE

    23  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL

    24  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS

    25  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)

    26  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT

    27  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY

    28  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF

    29  * SUCH DAMAGE.

    30  *

    31  * $FreeBSD: src/lib/libc/gen/fpclassify.c,v 1.2 2005/02/06 03:23:31 das Exp $

    32  */

    34 #include <sys/endian.h>

    36 #include <math.h>

    37 #include <stdint.h>

    39 #include "fpmath.h"

    42 EXPORT_C int __fpclassifyf(float f)

    43 {

    44 	union IEEEf2bits u;

    46 	u.f = f;

    47 	if (u.bits.exp == 0) {

    48 		if (u.bits.man == 0)

    49 			return (FP_ZERO);

    50 		return (FP_SUBNORMAL);

    51 	}

    52 	if (u.bits.exp == 255) {

    53 		if (u.bits.man == 0)

    54 			return (FP_INFINITE);

    55 		return (FP_NAN);

    56 	}

    57 	return (FP_NORMAL);

    58 }

    60 EXPORT_C int __fpclassifyd(double d)

    61 {

    62 	union IEEEd2bits u;

    64 	u.d = d;

    65 	if (u.bits.exp == 0) {

    66 		if ((u.bits.manl | u.bits.manh) == 0)

    67 			return (FP_ZERO);

    68 		return (FP_SUBNORMAL);

    69 	}

    70 	if (u.bits.exp == 2047) {

    71 		if ((u.bits.manl | u.bits.manh) == 0)

    72 			return (FP_INFINITE);

    73 		return (FP_NAN);

    74 	}

    75 	return (FP_NORMAL);

    76 }

    79 EXPORT_C int __fpclassifyl(long double e)

    80 {

    81 	union IEEEl2bits u;

    83 	u.e = e;

    84 	if (u.bits.exp == 0) {

    85 		if ((u.bits.manl | u.bits.manh) == 0)

    86 			return (FP_ZERO);

    87 		return (FP_SUBNORMAL);

    88 	}

    89 	mask_nbit_l(u);		/* Mask normalization bit if applicable. */

    90 	if (u.bits.exp == 32767) {

    91 		if ((u.bits.manl | u.bits.manh) == 0)

    92 			return (FP_INFINITE);

    93 		return (FP_NAN);

    94 	}

    95 	return (FP_NORMAL);

    96 }

author	sl
	Tue, 10 Jun 2014 14:32:02 +0200
changeset 1	260cb5ec6c19
permissions	-rw-r--r--