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
56
57
58
59
60
61
62
63
64
65
66
67
68
|
/* { dg-do run } */
/* { dg-skip-if "-mpowerpc-gpopt not supported" { powerpc*-*-darwin* } } */
/* { dg-options "-O2 -ffast-math -fno-inline -fno-unroll-loops -lm -mpowerpc-gpopt" } */
#include <math.h>
extern void abort (void);
#define NVALS 6
static double
convert_it_1 (double x)
{
return pow (x, 10.0 / 3.0);
}
static double
convert_it_2 (double x)
{
return pow (x, 11.0 / 3.0);
}
static double
convert_it_3 (double x)
{
return pow (x, -7.0 / 3.0);
}
static double
convert_it_4 (double x)
{
return pow (x, -8.0 / 3.0);
}
int
main (int argc, char *argv[])
{
double values[NVALS] = { 3.0, 1.95, 2.227, 4.0, 256.0, .0008797 };
double PREC = .999999;
unsigned i;
for (i = 0; i < NVALS; i++)
{
volatile double x, y;
x = __builtin_powi (values[i], 3);
y = __builtin_powi (cbrt (values[i]), 1);
if (fabs (convert_it_1 (values[i]) / (x * y)) < PREC)
abort ();
x = __builtin_powi (values[i], 3);
y = __builtin_powi (cbrt (values[i]), 2);
if (fabs (convert_it_2 (values[i]) / (x * y)) < PREC)
abort ();
x = __builtin_powi (values[i], -3);
y = __builtin_powi (cbrt (values[i]), 2);
if (fabs (convert_it_3 (values[i]) / (x * y)) < PREC)
abort ();
x = __builtin_powi (values[i], -3);
y = __builtin_powi (cbrt (values[i]), 1);
if (fabs (convert_it_4 (values[i]) / (x * y)) < PREC)
abort ();
}
return 0;
}
|
