1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
|
/* Copyright (C) 2003 Free Software Foundation.
Check that constant folding of built-in math functions doesn't
break anything and produces the expected results.
Written by Roger Sayle, 2nd April 2003. */
/* { dg-do link } */
/* { dg-options "-O2 -ffast-math" } */
extern void link_error(void);
extern double exp(double);
extern double log(double);
extern double sqrt(double);
extern double pow(double,double);
extern double fabs(double);
void test(double x)
{
if (sqrt(pow(x,4.0)) != x*x)
link_error ();
if (pow(sqrt(x),4.0) != x*x)
link_error ();
if (pow(pow(x,4.0),0.25) != x)
/* XFAIL. PR41098. */;
}
void test2(double x, double y, double z)
{
if (sqrt(pow(x,y)) != pow(fabs(x),y*0.5))
link_error ();
if (log(pow(x,y)) != y*log(x))
link_error ();
if (pow(exp(x),y) != exp(x*y))
link_error ();
if (pow(sqrt(x),y) != pow(x,y*0.5))
link_error ();
if (pow(pow(fabs(x),y),z) != pow(fabs(x),y*z))
link_error ();
}
int main()
{
test (2.0);
test2 (2.0, 3.0, 4.0);
return 0;
}
|
