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 

    31 #include <sys/cdefs.h>

    32 #ifndef __SYMBIAN32__

    33 __FBSDID("$FreeBSD: src/lib/msun/src/s_fmaxf.c,v 1.1 2004/06/30 07:04:01 das Exp $");

    34 #endif //__SYMBIAN32__

    35 

    36 #include <math.h>

    37 

    38 #include "fpmath.h"

    39 

    40 EXPORT_C float

    41 fmaxf(float x, float y)

    42 {

    43 	union IEEEf2bits u[2];

    44 

    45 	u[0].f = x;

    46 	u[1].f = y;

    47 

    48 	/* Check for NaNs to avoid raising spurious exceptions. */

    49 	if (u[0].bits.exp == 255 && u[0].bits.man != 0)

    50 		return (y);

    51 	if (u[1].bits.exp == 255 && u[1].bits.man != 0)

    52 		return (x);

    53 

    54 	/* Handle comparisons of signed zeroes. */

    55 	if (u[0].bits.sign != u[1].bits.sign)

    56 		return (u[u[0].bits.sign].f);

    57 

    58 	return (x > y ? x : y);

    59 }