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

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

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

     4 */

     5 /*-

     6  * Copyright (c) 2004 David Schultz <das@FreeBSD.ORG>

     7  * All rights reserved.

     8  *

     9  * Redistribution and use in source and binary forms, with or without

    10  * modification, are permitted provided that the following conditions

    11  * are met:

    12  * 1. Redistributions of source code must retain the above copyright

    13  *    notice, this list of conditions and the following disclaimer.

    14  * 2. Redistributions in binary form must reproduce the above copyright

    15  *    notice, this list of conditions and the following disclaimer in the

    16  *    documentation and/or other materials provided with the distribution.

    17  *

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

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

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

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

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

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

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

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

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

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

    28  * SUCH DAMAGE.

    29  *

    30  * $FreeBSD: src/lib/libc/gen/isinf.c,v 1.2 2004/07/10 15:52:26 marcel Exp $

    31  */

    32 

    33 #include <math.h>

    34 

    35 #include "fpmath.h"

    36 

    37 /*

    38  * XXX These routines belong in libm, but they must remain in libc for

    39  *     binary compat until we can bump libm's major version number.

    40  */

    41 

    42 #ifndef __SYMBIAN32__

    43 __weak_reference(__isinf, isinf);

    44 #endif //__SYMBIAN32__

    45 

    46 EXPORT_C int __isinf(double d)

    47 {

    48 	union IEEEd2bits u;

    49 

    50 	u.d = d;

    51 	#ifdef __SYMBIAN32__

    52 	return (u.bits.sign) ? (-1)*(u.bits.exp == 2047 && u.bits.manl == 0 && u.bits.manh == 0):

    53 						   (u.bits.exp == 2047 && u.bits.manl == 0 && u.bits.manh == 0);

    54 	#else

    55 	return (u.bits.exp == 2047 && u.bits.manl == 0 && u.bits.manh == 0);

    56 	#endif //__SYMBIAN32__ 

    57 }

    58 

    59 

    60 EXPORT_C int __isinff(float f)

    61 {

    62 	union IEEEf2bits u;

    63 

    64 	u.f = f;

    65 	#ifdef __SYMBIAN32__

    66 	return (u.bits.sign) ? (-1)*(u.bits.exp == 255 && u.bits.man == 0):

    67 						   (u.bits.exp == 255 && u.bits.man == 0)	;

    68 	#else

    69 	return (u.bits.exp == 255 && u.bits.man == 0);

    70 	#endif //__SYMBIAN32__

    71 }

    72 

    73 

    74 EXPORT_C int __isinfl(long double e)

    75 {

    76 	union IEEEl2bits u;

    77 

    78 	u.e = e;

    79 	mask_nbit_l(u);

    80 #ifndef __alpha__

    81 	return (u.bits.exp == 32767 && u.bits.manl == 0 && u.bits.manh == 0);

    82 #else

    83 	return (u.bits.exp == 2047 && u.bits.manl == 0 && u.bits.manh == 0);

    84 #endif

    85 }