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
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
|
/*
Program to test complex divide for correct results on selected values.
Checking known failure points.
*/
#include <float.h>
extern void abort (void);
extern void exit (int);
extern int ilogb (double);
int match (double _Complex, double _Complex);
#define SMALL DBL_MIN
#define MAXBIT DBL_MANT_DIG
#define ERRLIM 6
/*
Compare c (computed value) with z (expected value).
Return 0 if within allowed range. Return 1 if not.
*/
int match (double _Complex c, double _Complex z)
{
double rz, iz, rc, ic;
double rerr, ierr, rmax;
int biterr;
rz = __real__ z;
iz = __imag__ z;
rc = __real__ c;
ic = __imag__ c;
if (__builtin_fabs (rz) > SMALL)
{
rerr = __builtin_fabs (rz - rc) / __builtin_fabs (rz);
}
else if (__builtin_fabs (rz) == 0.0)
{
rerr = __builtin_fabs (rc);
}
else
{
rerr = __builtin_fabs (rz - rc) / SMALL;
}
if (__builtin_fabs (iz) > SMALL)
{
ierr = __builtin_fabs (iz - ic) / __builtin_fabs (iz);
}
else if (__builtin_fabs (iz) == 0.0)
{
ierr = __builtin_fabs (ic);
}
else
{
ierr = __builtin_fabs (iz - ic) / SMALL;
}
rmax = __builtin_fmax(rerr, ierr);
biterr = 0;
if ( rmax != 0.0)
{
biterr = ilogb (rmax) + MAXBIT + 1;
}
if (biterr >= ERRLIM)
return 0;
else
return 1;
}
int main (int argc, char** argv)
{
double _Complex a,b,c,z;
double xr[4], xi[4], yr[4], yi[4], zr[4], zi[4];
double cr, ci;
int i;
int ok = 1;
xr[0] = -0x1.16e7fad79e45ep+651;
xi[0] = -0x1.f7f75b94c6c6ap-860;
yr[0] = -0x1.2f40d8ff7e55ep+245;
yi[0] = -0x0.0000000004ebcp-1022;
zr[0] = 0x1.d6e4b0e282869p+405;
zi[0] = -0x1.e9095e311e706p-900;
xr[1] = -0x1.21ff587f953d3p-310;
xi[1] = -0x1.5a526dcc59960p+837;
yr[1] = 0x1.b88b8b552eaadp+735;
yi[1] = -0x1.873e2d6544d92p-327;
zr[1] = 0x1.65734a88b2de0p-961;
zi[1] = -0x1.927e85b8b5770p+101;
xr[2] = 0x1.4612e41aa8080p-846;
xi[2] = -0x0.0000000613e07p-1022;
yr[2] = 0x1.df9cd0d58caafp-820;
yi[2] = -0x1.e47051a9036dbp-584;
zr[2] = 0x1.9b194f3fffa32p-469;
zi[2] = 0x1.58a00ab740a6bp-263;
xr[3] = 0x1.cb27eece7c585p-355;
xi[3] = 0x0.000000223b8a8p-1022;
yr[3] = -0x1.74e7ed2b9189fp-22;
yi[3] = 0x1.3d80439e9a119p-731;
zr[3] = -0x1.3b35ed806ae5ap-333;
zi[3] = -0x0.05e01bcbfd9f6p-1022;
for (i = 0; i < 4; i++)
{
__real__ a = xr[i];
__imag__ a = xi[i];
__real__ b = yr[i];
__imag__ b = yi[i];
__real__ z = zr[i];
__imag__ z = zi[i];
c = a / b;
cr = __real__ c;
ci = __imag__ c;
if (!match (c,z)){
ok = 0;
}
}
if (!ok)
abort ();
exit (0);
}
|
