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
88
89
90
91
92
93
94
95
96
|
/* Test to see if floating point state is being properly maintained
for each thread. Different threads doing floating point operations
simultaneously should not interfere with one another. This
includes operations that might change some FPU flags, such as
rounding modes, at least implicitly. */
#include <pthread.h>
#include <math.h>
#include <stdio.h>
int limit = 2;
int float_passed = 0;
int float_failed = 1;
void *log_loop (void *x) {
int i;
double d, d1, d2;
/* sleep (1); */
for (i = 0; i < limit; i++) {
d = 42.0;
d = log (exp (d));
d = (d + 39.0) / d;
if (i == 0)
d1 = d;
else {
d2 = d;
d = sin(d);
/* if (d2 != d1) { */
if (memcmp (&d2, &d1, 8)) {
pthread_exit(&float_failed);
}
}
}
pthread_exit(&float_passed);
}
void *trig_loop (void *x) {
int i;
double d, d1, d2;
/* sleep (1); */
for (i = 0; i < limit; i++) {
d = 35.0;
d *= M_PI;
d /= M_LN2;
d = sin (d);
d = cos (1 / d);
if (i == 0)
d1 = d;
else {
d2 = d;
d = sin(d);
/* if (d2 != d1) { */
if (memcmp (&d2, &d1, 8)) {
pthread_exit(&float_failed);
}
}
}
pthread_exit(&float_passed);
}
#define N 10
int main () {
pthread_t thread[2];
pthread_attr_t attr;
int *x, *y;
pthread_attr_init(&attr);
/* pthread_attr_setfloatstate(&attr, PTHREAD_NOFLOAT); */
while(limit < 100000) {
pthread_create (&thread[0], &attr, trig_loop, 0);
pthread_create (&thread[1], &attr, log_loop, 0);
pthread_join(thread[0], (void **) &x);
pthread_join(thread[1], (void **) &y);
if ((*x == float_failed) || (*y == float_failed)) {
limit *= 4;
break;
}
limit *= 4;
}
if ((*x == float_passed) && (*y == float_passed)) {
printf("test_preemption_float INDETERMINATE\n");
return(0);
}
pthread_create (&thread[0], NULL, trig_loop, 0);
pthread_create (&thread[1], NULL, log_loop, 0);
pthread_join(thread[0], (void **) &x);
pthread_join(thread[1], (void **) &y);
if ((*x == float_failed) || (*y == float_failed)) {
printf("test_preemption_float FAILED\n");
return(1);
}
printf("test_preemption_float PASSED\n");
return(0);
}
|
