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
|
/* Copyright (C) 2003,2007 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, 5th April 2003. */
/* { dg-do link } */
/* { dg-options "-O2 -ffast-math" } */
extern void link_error(void);
extern double exp(double);
extern double sqrt(double);
extern double cbrt(double);
extern double pow(double,double);
void test(double x, double y, double z)
{
if (sqrt(x)*sqrt(x) != x)
link_error ();
if (sqrt(x)*sqrt(y) != sqrt(x*y))
link_error ();
if (exp(x)*exp(y) != exp(x+y))
link_error ();
if (pow(x,y)*pow(z,y) != pow(z*x,y))
link_error ();
if (pow(x,y)*pow(x,z) != pow(x,y+z))
link_error ();
if (x/exp(y) != x*exp(-y))
link_error ();
if (x/pow(y,z) != x*pow(y,-z))
link_error ();
if (x/sqrt(y/z) != x*sqrt(z/y))
link_error ();
if (x/cbrt(y/z) != x*cbrt(z/y))
link_error ();
}
int main()
{
test (2.0, 3.0, 4.0);
return 0;
}
|
