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
|
#include <immintrin.h>
#include "m128-check.h"
#ifndef max
#define max(a, b) (((a) > (b)) ? (a):(b))
#endif
#ifndef min
#define min(a, b) (((a) < (b)) ? (a):(b))
#endif
typedef union
{
__m256i x;
char a[32];
} union256i_b;
typedef union
{
__m256i x;
short a[16];
} union256i_w;
typedef union
{
__m256i x;
int a[8];
} union256i_d;
typedef union
{
__m256i x;
long long a[4];
} union256i_q;
typedef union
{
__m256 x;
float a[8];
} union256;
typedef union
{
__m256d x;
double a[4];
} union256d;
CHECK_EXP (union256i_b, char, "%d")
CHECK_EXP (union256i_w, short, "%d")
CHECK_EXP (union256i_d, int, "0x%x")
CHECK_EXP (union256i_q, long long, "0x%llx")
CHECK_EXP (union256, float, "%f")
CHECK_EXP (union256d, double, "%f")
#define CHECK_FP_EXP(UINON_TYPE, VALUE_TYPE, ESP, FMT) \
static int \
__attribute__((noinline, unused)) \
check_fp_##UINON_TYPE (UINON_TYPE u, const VALUE_TYPE *v) \
{ \
int i; \
int err = 0; \
\
for (i = 0; i < ARRAY_SIZE (u.a); i++) \
if (u.a[i] > (v[i] + (ESP)) || u.a[i] < (v[i] - (ESP))) \
{ \
err++; \
PRINTF ("%i: " FMT " != " FMT "\n", \
i, v[i], u.a[i]); \
} \
return err; \
}
CHECK_FP_EXP (union256, float, ESP_FLOAT, "%f")
CHECK_FP_EXP (union256d, double, ESP_DOUBLE, "%f")
|
