← Back to branch summary

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

« back to all changes in this revision

Viewing changes to math/test-powl.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/test-powl.c
 
1
/* Test for powl
 
2
   Copyright (C) 2011 Free Software Foundation, Inc.
 
3
   This file is part of the GNU C Library.
 
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 <stdio.h>
 
21
#include <math.h>
 
22
#include <float.h>
 
23
#include <ieee754.h>
 
24
 
 
25
int
 
26
main (void)
 
27
{
 
28
  int result = 0;
 
29
 
 
30
#ifndef NO_LONG_DOUBLE
 
31
# if LDBL_MANT_DIG == 64
 
32
    {
 
33
      long double x = 1e-20;
 
34
      union ieee854_long_double u;
 
35
      u.ieee.mantissa0 = 1;
 
36
      u.ieee.mantissa1 = 1;
 
37
      u.ieee.exponent = 0;
 
38
      u.ieee.negative = 0;
 
39
      (void) powl (0.2, u.d);
 
40
      x = powl (x, 1.5);
 
41
      if (fabsl (x - 1e-30) > 1e-10)
 
42
        {
 
43
          printf ("powl (1e-20, 1.5): wrong result: %Lg\n", x);
 
44
          result = 1;
 
45
        }
 
46
    }
 
47
# endif
 
48
#endif
 
49
 
 
50
  return result;
 
51
}