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
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
|
/* PR c/64868 */
/* { dg-do run } */
/* { dg-options "-O2 -fopenmp" } */
float f = 2.0f;
double d = 4.0;
long double ld = 8.0L;
void
foo ()
{
#pragma omp atomic
f = 1.0f / f;
#pragma omp atomic
f = 1 / f;
#pragma omp atomic
f = f / 2.0f;
#pragma omp atomic
f = f / 2;
#pragma omp atomic
f /= 2.0f;
#pragma omp atomic
f /= 2;
#pragma omp atomic
d = 1.0 / d;
#pragma omp atomic
d = 1 / d;
#pragma omp atomic
d = d / 2.0;
#pragma omp atomic
d = d / 2;
#pragma omp atomic
d /= 2.0;
#pragma omp atomic
d /= 2;
#pragma omp atomic
ld = 1.0L / ld;
#pragma omp atomic
ld = 1 / ld;
#pragma omp atomic
ld = ld / 2.0L;
#pragma omp atomic
ld = ld / 2;
#pragma omp atomic
ld /= 2.0L;
#pragma omp atomic
ld /= 2;
if (f != 0.125f || d != 0.25 || ld != 0.5L)
__builtin_abort ();
}
#ifdef __cplusplus
template <typename T, int N1, int N2>
void
bar ()
{
T v = ::d;
#pragma omp atomic
v *= 16;
#pragma omp atomic
v = 1.0 / v;
#pragma omp atomic
v = N1 / v;
#pragma omp atomic
v = v / 2.0;
#pragma omp atomic
v = v / N2;
#pragma omp atomic
v /= 2.0;
#pragma omp atomic
v /= N2;
if (v != 0.25)
__builtin_abort ();
}
#endif
int
main ()
{
foo ();
#ifdef __cplusplus
bar<float, 1, 2> ();
bar<double, 1, 2> ();
bar<long double, 1, 2> ();
#endif
return 0;
}
|
