← Back to branch summary

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

« back to all changes in this revision

Viewing changes to math/w_j0l.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:
math/w_j0l.c
1
 
/* w_j0l.c -- long double version of w_j0.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 j0l(long double x), y0l(long double x)
23
 
 */
24
 
 
 
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 <fenv.h>
25
21
#include <math.h>
26
 
#include "math_private.h"
27
 
 
28
 
#ifdef __STDC__
29
 
        long double __j0l(long double x)                /* wrapper j0l */
30
 
#else
31
 
        long double __j0l(x)                            /* wrapper j0 */
32
 
        long double x;
33
 
#endif
 
22
#include <math_private.h>
 
23
 
 
24
 
 
25
/* wrapper j0l */
 
26
long double
 
27
__j0l (long double x)
34
28
{
35
 
#ifdef _IEEE_LIBM
36
 
        return __ieee754_j0l(x);
37
 
#else
38
 
        long double z = __ieee754_j0l(x);
39
 
        if(_LIB_VERSION == _IEEE_ || __isnanl(x)) return z;
40
 
        if(fabsl(x)>X_TLOSS) {
41
 
                return __kernel_standard(x,x,234); /* j0(|x|>X_TLOSS) */
42
 
        } else
43
 
            return z;
44
 
#endif
 
29
  if (__builtin_expect (fabsl (x) > X_TLOSS, 0) && _LIB_VERSION != _IEEE_)
 
30
    /* j0(|x|>X_TLOSS) */
 
31
    return __kernel_standard (x, x, 234);
 
32
 
 
33
  return __ieee754_j0l (x);
45
34
}
46
35
weak_alias (__j0l, j0l)
47
36
 
48
 
#ifdef __STDC__
49
 
        long double __y0l(long double x)                /* wrapper y0l */
50
 
#else
51
 
        long double __y0l(x)                            /* wrapper y0 */
52
 
        long double x;
53
 
#endif
 
37
 
 
38
/* wrapper y0l */
 
39
long double
 
40
__y0l (long double x)
54
41
{
55
 
#ifdef _IEEE_LIBM
56
 
        return __ieee754_y0l(x);
57
 
#else
58
 
        long double z;
59
 
        z = __ieee754_y0l(x);
60
 
        if(_LIB_VERSION == _IEEE_ || __isnanl(x) ) return z;
61
 
        if(x <= 0.0){
62
 
                if(x==0.0)
63
 
                    /* d= -one/(x-x); */
64
 
                    return __kernel_standard(x,x,208);
65
 
                else
66
 
                    /* d = zero/(x-x); */
67
 
                    return __kernel_standard(x,x,209);
68
 
        }
69
 
        if(x>X_TLOSS) {
70
 
                return __kernel_standard(x,x,235); /* y0(x>X_TLOSS) */
71
 
        } else
72
 
            return z;
73
 
#endif
 
42
  if (__builtin_expect (x <= 0.0L || x > X_TLOSS, 0) && _LIB_VERSION != _IEEE_)
 
43
    {
 
44
      if (x < 0.0L)
 
45
        {
 
46
          /* d = zero/(x-x) */
 
47
          feraiseexcept (FE_INVALID);
 
48
          return __kernel_standard (x, x, 209);
 
49
        }
 
50
      else if (x == 0.0L)
 
51
        /* d = -one/(x-x) */
 
52
        return __kernel_standard (x, x, 208);
 
53
      else
 
54
        /* y0(x>X_TLOSS) */
 
55
        return __kernel_standard (x, x, 235);
 
56
    }
 
57
 
 
58
  return __ieee754_y0l (x);
74
59
}
75
 
 
76
60
weak_alias (__y0l, y0l)