← Back to branch summary

~ubuntu-branches/ubuntu/precise/eglibc/precise-201308281639

« back to all changes in this revision

Viewing changes to sysdeps/ieee754/ldbl-96/w_expl.c

  • Committer: Package Import Robot
  • Author(s): Matthias Klose
  • Date: 2012-02-08 01:58:09 UTC
  • mfrom: (1.5.3) (288.1.12 precise)
  • Revision ID: package-import@ubuntu.com-20120208015809-ulscst7uteq3e22z
Tags: 2.15~pre6-0ubuntu10
Merge from Debian (r5151, 2.13-26).

Show diffs side-by-side

added added

removed removed

Lines of Context:
sysdeps/ieee754/ldbl-96/w_expl.c
1
 
/* w_expl.c -- long double version of w_exp.c.
2
 
 * Conversion to long double by Ulrich Drepper,
3
 
 * Cygnus Support, drepper@cygnus.com.
4
 
 */
5
 
 
6
 
/*
7
 
 * ====================================================
8
 
 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
9
 
 *
10
 
 * Developed at SunPro, a Sun Microsystems, Inc. business.
11
 
 * Permission to use, copy, modify, and distribute this
12
 
 * software is freely granted, provided that this notice
13
 
 * is preserved.
14
 
 * ====================================================
15
 
 */
16
 
 
17
 
#if defined(LIBM_SCCS) && !defined(lint)
18
 
static char rcsid[] = "$NetBSD: $";
19
 
#endif
20
 
 
21
 
/*
22
 
 * wrapper expl(x)
23
 
 */
24
 
 
25
 
#include "math.h"
26
 
#include "math_private.h"
27
 
 
28
 
#ifdef __STDC__
 
1
/* Copyright (C) 2011 Free Software Foundation, Inc.
 
2
   This file is part of the GNU C Library.
 
3
   Contributed by Ulrich Drepper <drepper@gmail.com>, 2011.
 
4
 
 
5
   The GNU C Library is free software; you can redistribute it and/or
 
6
   modify it under the terms of the GNU Lesser General Public
 
7
   License as published by the Free Software Foundation; either
 
8
   version 2.1 of the License, or (at your option) any later version.
 
9
 
 
10
   The GNU C Library is distributed in the hope that it will be useful,
 
11
   but WITHOUT ANY WARRANTY; without even the implied warranty of
 
12
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 
13
   Lesser General Public License for more details.
 
14
 
 
15
   You should have received a copy of the GNU Lesser General Public
 
16
   License along with the GNU C Library; if not, write to the Free
 
17
   Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
 
18
   02111-1307 USA.  */
 
19
 
 
20
#include <math.h>
 
21
#include <math_private.h>
 
22
 
29
23
static const long double
30
 
#else
31
 
static long double
32
 
#endif
33
24
o_threshold=  1.135652340629414394949193107797076489134e4,
34
25
  /* 0x400C, 0xB17217F7, 0xD1CF79AC */
35
26
u_threshold= -1.140019167866942050398521670162263001513e4;
36
27
  /* 0x400C, 0xB220C447, 0x69C201E8 */
37
28
 
38
 
#ifdef __STDC__
39
 
        long double __expl(long double x)       /* wrapper exp */
40
 
#else
41
 
        long double __expl(x)                   /* wrapper exp */
42
 
        long double x;
43
 
#endif
 
29
 
 
30
/* wrapper expl */
 
31
long double
 
32
__expl (long double x)
44
33
{
45
 
#ifdef _IEEE_LIBM
46
 
        return __ieee754_expl(x);
47
 
#else
48
 
        long double z;
49
 
        z = __ieee754_expl(x);
50
 
        if(_LIB_VERSION == _IEEE_) return z;
51
 
        if(__finitel(x)) {
52
 
            if(x>o_threshold)
53
 
                return __kernel_standard(x,x,206); /* exp overflow */
54
 
            else if(x<u_threshold)
55
 
                return __kernel_standard(x,x,207); /* exp underflow */
56
 
        }
57
 
        return z;
58
 
#endif
 
34
  if (__builtin_expect (x > o_threshold, 0))
 
35
    {
 
36
      if (_LIB_VERSION != _IEEE_)
 
37
        return __kernel_standard (x, x, 206);
 
38
    }
 
39
  else if (__builtin_expect (x < u_threshold, 0))
 
40
    {
 
41
      if (_LIB_VERSION != _IEEE_)
 
42
        return __kernel_standard (x, x, 207);
 
43
    }
 
44
 
 
45
  return __ieee754_expl (x);
59
46
}
60
47
hidden_def (__expl)
61
48
weak_alias (__expl, expl)