← Back to branch summary

~xnox/ubuntu/trusty/gcc-arm-linux-androideabi/dima

« back to all changes in this revision

Viewing changes to android/bionic/libm/src/s_nexttowardf.c

  • Committer: Package Import Robot
  • Author(s): Dmitrijs Ledkovs
  • Date: 2013-07-05 10:12:24 UTC
  • Revision ID: package-import@ubuntu.com-20130705101224-6qo3e8jbz8p31aa1
Tags: upstream-0.20130705.1
Import upstream version 0.20130705.1

Show diffs side-by-side

added added

removed removed

Lines of Context:
android/bionic/libm/src/s_nexttowardf.c
 
1
/*
 
2
 * ====================================================
 
3
 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
 
4
 *
 
5
 * Developed at SunPro, a Sun Microsystems, Inc. business.
 
6
 * Permission to use, copy, modify, and distribute this
 
7
 * software is freely granted, provided that this notice
 
8
 * is preserved.
 
9
 * ====================================================
 
10
 */
 
11
 
 
12
#ifndef lint
 
13
static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_nexttowardf.c,v 1.1 2005/03/07 04:57:38 das Exp $";
 
14
#endif
 
15
 
 
16
#include <float.h>
 
17
 
 
18
#include "fpmath.h"
 
19
#include "math.h"
 
20
#include "math_private.h"
 
21
 
 
22
#define LDBL_INFNAN_EXP (LDBL_MAX_EXP * 2 - 1)
 
23
 
 
24
float
 
25
nexttowardf(float x, long double y)
 
26
{
 
27
        union IEEEl2bits uy;
 
28
        volatile float t;
 
29
        int32_t hx,ix;
 
30
 
 
31
        GET_FLOAT_WORD(hx,x);
 
32
        ix = hx&0x7fffffff;             /* |x| */
 
33
        uy.e = y;
 
34
 
 
35
        if((ix>0x7f800000) ||
 
36
           (uy.bits.exp == LDBL_INFNAN_EXP &&
 
37
            ((uy.bits.manh&~LDBL_NBIT)|uy.bits.manl) != 0))
 
38
           return x+y;  /* x or y is nan */
 
39
        if(x==y) return (float)y;               /* x=y, return y */
 
40
        if(ix==0) {                             /* x == 0 */
 
41
            SET_FLOAT_WORD(x,(uy.bits.sign<<31)|1);/* return +-minsubnormal */
 
42
            t = x*x;
 
43
            if(t==x) return t; else return x;   /* raise underflow flag */
 
44
        }
 
45
        if(hx>=0 ^ x < y)                       /* x -= ulp */
 
46
            hx -= 1;
 
47
        else                                    /* x += ulp */
 
48
            hx += 1;
 
49
        ix = hx&0x7f800000;
 
50
        if(ix>=0x7f800000) return x+x;  /* overflow  */
 
51
        if(ix<0x00800000) {             /* underflow */
 
52
            t = x*x;
 
53
            if(t!=x) {          /* raise underflow flag */
 
54
                SET_FLOAT_WORD(y,hx);
 
55
                return y;
 
56
            }
 
57
        }
 
58
        SET_FLOAT_WORD(x,hx);
 
59
        return x;
 
60
}