← Back to branch summary

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

« back to all changes in this revision

Viewing changes to sysdeps/ieee754/dbl-64/e_exp.c

  • Committer: Package Import Robot
  • Author(s): Matthias Klose
  • Date: 2011-10-04 17:48:26 UTC
  • mfrom: (216.1.23 oneiric)
  • Revision ID: package-import@ubuntu.com-20111004174826-2cyb9ewn3ucymlsx
Tags: 2.13-20ubuntu5
https://launchpad.net/bugs/853688
libc6-dev: Don't break the current {gnat,gcj}-4.4-base versons. LP: #853688.

Show diffs side-by-side

added added

removed removed

Lines of Context:
sysdeps/ieee754/dbl-64/e_exp.c
1
1
/*
2
2
 * IBM Accurate Mathematical Library
3
3
 * written by International Business Machines Corp.
4
 
 * Copyright (C) 2001, 2011 Free Software Foundation
 
4
 * Copyright (C) 2001 Free Software Foundation
5
5
 *
6
6
 * This program is free software; you can redistribute it and/or modify
7
7
 * it under the terms of the GNU Lesser General Public License as published by
40
40
#include "uexp.tbl"
41
41
#include "math_private.h"
42
42
 
43
 
#ifndef SECTION
44
 
# define SECTION
45
 
#endif
46
 
 
47
43
double __slowexp(double);
48
44
 
49
45
/***************************************************************************/
50
46
/* An ultimate exp routine. Given an IEEE double machine number x          */
51
47
/* it computes the correctly rounded (to nearest) value of e^x             */
52
48
/***************************************************************************/
53
 
double
54
 
SECTION
55
 
__ieee754_exp(double x) {
 
49
double __ieee754_exp(double x) {
56
50
  double bexp, t, eps, del, base, y, al, bet, res, rem, cor;
57
51
  mynumber junk1, junk2, binexp  = {{0,0}};
58
52
#if 0
151
145
    else return __slowexp(x);
152
146
  }
153
147
}
154
 
#ifndef __ieee754_exp
155
 
strong_alias (__ieee754_exp, __exp_finite)
156
 
#endif
157
148
 
158
149
/************************************************************************/
159
150
/* Compute e^(x+xx)(Double-Length number) .The routine also receive     */
162
153
/*else return   e^(x + xx)   (always positive )                         */
163
154
/************************************************************************/
164
155
 
165
 
double
166
 
SECTION
167
 
__exp1(double x, double xx, double error) {
 
156
double __exp1(double x, double xx, double error) {
168
157
  double bexp, t, eps, del, base, y, al, bet, res, rem, cor;
169
158
  mynumber junk1, junk2, binexp  = {{0,0}};
170
159
#if 0