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
|
#include <smmintrin.h>
#include <fenv.h>
#include "sse4_1-check.h"
#define DIM(a) (sizeof (a) / sizeof (a)[0])
static int roundings[] =
{
_MM_FROUND_TO_NEAREST_INT,
_MM_FROUND_TO_NEG_INF,
_MM_FROUND_TO_POS_INF,
_MM_FROUND_TO_ZERO,
_MM_FROUND_CUR_DIRECTION
};
static int modes[] =
{
FE_TONEAREST,
FE_UPWARD,
FE_DOWNWARD,
FE_TOWARDZERO
};
static void
TEST (void)
{
int i, j, ri, mi, round_save;
round_save = fegetround ();
for (mi = 0; mi < DIM (modes); mi++) {
fesetround (modes[mi]);
for (i = 0; i < DIM (data); i++) {
for (ri = 0; ri < DIM (roundings); ri++) {
union value guess;
union value *current_answers = answers[ri];
switch ( roundings[ri] ) {
case _MM_FROUND_TO_NEAREST_INT:
guess.x = ROUND_INTRIN (data[i].value1.x, data[i].value2.x,
_MM_FROUND_TO_NEAREST_INT);
break;
case _MM_FROUND_TO_NEG_INF:
guess.x = ROUND_INTRIN (data[i].value1.x, data[i].value2.x,
_MM_FROUND_TO_NEG_INF);
break;
case _MM_FROUND_TO_POS_INF:
guess.x = ROUND_INTRIN (data[i].value1.x, data[i].value2.x,
_MM_FROUND_TO_POS_INF);
break;
case _MM_FROUND_TO_ZERO:
guess.x = ROUND_INTRIN (data[i].value1.x, data[i].value2.x,
_MM_FROUND_TO_ZERO);
break;
case _MM_FROUND_CUR_DIRECTION:
guess.x = ROUND_INTRIN (data[i].value1.x, data[i].value2.x,
_MM_FROUND_CUR_DIRECTION);
switch ( modes[mi] ) {
case FE_TONEAREST:
current_answers = answers_NEAREST_INT;
break;
case FE_UPWARD:
current_answers = answers_POS_INF;
break;
case FE_DOWNWARD:
current_answers = answers_NEG_INF;
break;
case FE_TOWARDZERO:
current_answers = answers_ZERO;
break;
}
break;
default:
abort ();
}
for (j = 0; j < DIM (guess.f); j++)
if (guess.f[j] != current_answers[i].f[j])
abort ();
}
}
}
fesetround (round_save);
}
|
