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
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
|
/* mpcheck-longdouble -- compare mpc functions against "long double complex"
from the GNU libc implementation
Copyright (C) 2020 INRIA
This file is part of GNU MPC.
GNU MPC is free software; you can redistribute it and/or modify it under
the terms of the GNU Lesser General Public License as published by the
Free Software Foundation; either version 3 of the License, or (at your
option) any later version.
GNU MPC is distributed in the hope that it will be useful, but WITHOUT ANY
WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
more details.
You should have received a copy of the GNU Lesser General Public License
along with this program. If not, see http://www.gnu.org/licenses/ .
*/
/* the GNU libc provides the following functions (as of 2.31),
with 'f' suffix for the float/binary32 version, with no suffix
for the double/binary64 version, with 'l' suffix for the long double
version, and with 'f128' suffix for the __float128 version:
cabs casinh cexp csinh
cacos catan clog csqrt
cacosh catanh clog10 ctan
carg ccos cpow ctanh
casin ccosh csin
*/
#define _GNU_SOURCE /* for clog10 */
#include <stdlib.h>
#include <string.h>
#include <complex.h>
#include <sys/types.h>
#include <unistd.h>
#include <assert.h>
#include "mpc-tests.h"
#ifdef __GNUC__
#include <gnu/libc-version.h>
#endif
#define PRECISION 64
#define EMAX 16384
#define TYPE long double
#define SUFFIX l
#define mpc_get_type mpc_get_ldc
#define mpc_set_type mpc_set_ldc
#define mpfr_set_type mpfr_set_ld
gmp_randstate_t state;
mpz_t expz; /* global variable used in mpcheck_random */
unsigned long seed = 0;
int verbose = 0;
mpfr_exp_t emin, emax;
#include "mpcheck-common.c"
#define FOO add
#define CFOO(x,y) (x+y)
#include "mpcheck-template3.c"
#define FOO sub
#define CFOO(x,y) (x-y)
#include "mpcheck-template3.c"
#define FOO mul
#define CFOO(x,y) (x*y)
#include "mpcheck-template3.c"
#define FOO div
#define CFOO(x,y) (x/y)
#include "mpcheck-template3.c"
#define FOO pow
#include "mpcheck-template3.c"
#define FOO abs
#include "mpcheck-template2.c"
#define FOO arg
#include "mpcheck-template2.c"
#define FOO sqrt
#include "mpcheck-template1.c"
#define FOO acos
#include "mpcheck-template1.c"
#define FOO acosh
#include "mpcheck-template1.c"
#define FOO asin
#include "mpcheck-template1.c"
#define FOO asinh
#include "mpcheck-template1.c"
#define FOO atan
#include "mpcheck-template1.c"
#define FOO atanh
#include "mpcheck-template1.c"
#define FOO cos
#include "mpcheck-template1.c"
#define FOO cosh
#include "mpcheck-template1.c"
#define FOO exp
#include "mpcheck-template1.c"
#define FOO log
#include "mpcheck-template1.c"
#define FOO log10
#include "mpcheck-template1.c"
#define FOO sin
#include "mpcheck-template1.c"
#define FOO sinh
#include "mpcheck-template1.c"
#define FOO tan
#include "mpcheck-template1.c"
#define FOO tanh
#include "mpcheck-template1.c"
int
main (int argc, char *argv[])
{
mpfr_prec_t p = PRECISION; /* precision of 'double' */
unsigned long n = 1000000; /* default number of random tests per function */
while (argc >= 2 && argv[1][0] == '-')
{
if (argc >= 3 && strcmp (argv[1], "-p") == 0)
{
p = atoi (argv[2]);
argc -= 2;
argv += 2;
}
else if (argc >= 3 && strcmp (argv[1], "-seed") == 0)
{
seed = atoi (argv[2]);
argc -= 2;
argv += 2;
}
else if (argc >= 3 && strcmp (argv[1], "-num") == 0)
{
n = atoi (argv[2]);
argc -= 2;
argv += 2;
}
else if (strcmp (argv[1], "-v") == 0)
{
verbose ++;
argc --;
argv ++;
}
else
{
fprintf (stderr, "Unknown option %s\n", argv[1]);
exit (1);
}
}
/* set exponent range */
emin = -EMAX - PRECISION + 4; /* should be -16444 */
emax = EMAX;
mpfr_set_emin (emin);
mpfr_set_emax (emax);
gmp_randinit_default (state);
mpz_init (expz);
#ifdef __GNUC__
printf ("GNU libc version: %s\n", gnu_get_libc_version ());
printf ("GNU libc release: %s\n", gnu_get_libc_release ());
#endif
if (seed == 0)
seed = getpid ();
printf ("Using random seed %lu\n", seed);
/* (complex,complex) -> complex */
test_add (p, n);
test_sub (p, n);
test_mul (p, n);
test_div (p, n);
test_pow (p, n);
/* complex -> real */
test_abs (p, n);
test_arg (p, n);
/* complex -> complex */
test_sqrt (p, n);
test_acos (p, n);
test_acosh (p, n);
test_asin (p, n);
test_asinh (p, n);
test_atan (p, n);
test_atanh (p, n);
test_cos (p, n);
test_cosh (p, n);
test_exp (p, n);
test_log (p, n);
test_log10 (p, n);
test_sin (p, n);
test_sinh (p, n);
test_tan (p, n);
test_tanh (p, n);
gmp_randclear (state);
mpz_clear (expz);
return 0;
}
|
