OpenCores
URL https://opencores.org/ocsvn/openrisc/openrisc/trunk

Subversion Repositories openrisc

[/] [openrisc/] [trunk/] [gnu-dev/] [or1k-gcc/] [libgo/] [go/] [math/] [all_test.go] - Blame information for rev 858

Go to most recent revision | Details | Compare with Previous | View Log

Line No. Rev Author Line
1 747 jeremybenn 
// Copyright 2009 The Go Authors. All rights reserved.
2 
// Use of this source code is governed by a BSD-style
3 
// license that can be found in the LICENSE file.
4 
 
5 
package math_test
6 
 
7 
import (
8 
        "fmt"
9 
        . "math"
10 
        "testing"
11 
)
12 
 
13 
var vf = []float64{
14 
        4.9790119248836735e+00,
15 
        7.7388724745781045e+00,
16 
        -2.7688005719200159e-01,
17 
        -5.0106036182710749e+00,
18 
        9.6362937071984173e+00,
19 
        2.9263772392439646e+00,
20 
        5.2290834314593066e+00,
21 
        2.7279399104360102e+00,
22 
        1.8253080916808550e+00,
23 
        -8.6859247685756013e+00,
24 
}
25 
 
26 
// The expected results below were computed by the high precision calculators
27 
// at http://keisan.casio.com/.  More exact input values (array vf[], above)
28 
// were obtained by printing them with "%.26f".  The answers were calculated
29 
// to 26 digits (by using the "Digit number" drop-down control of each
30 
// calculator).
31 
var acos = []float64{
32 
        1.0496193546107222142571536e+00,
33 
        6.8584012813664425171660692e-01,
34 
        1.5984878714577160325521819e+00,
35 
        2.0956199361475859327461799e+00,
36 
        2.7053008467824138592616927e-01,
37 
        1.2738121680361776018155625e+00,
38 
        1.0205369421140629186287407e+00,
39 
        1.2945003481781246062157835e+00,
40 
        1.3872364345374451433846657e+00,
41 
        2.6231510803970463967294145e+00,
42 
}
43 
var acosh = []float64{
44 
        2.4743347004159012494457618e+00,
45 
        2.8576385344292769649802701e+00,
46 
        7.2796961502981066190593175e-01,
47 
        2.4796794418831451156471977e+00,
48 
        3.0552020742306061857212962e+00,
49 
        2.044238592688586588942468e+00,
50 
        2.5158701513104513595766636e+00,
51 
        1.99050839282411638174299e+00,
52 
        1.6988625798424034227205445e+00,
53 
        2.9611454842470387925531875e+00,
54 
}
55 
var asin = []float64{
56 
        5.2117697218417440497416805e-01,
57 
        8.8495619865825236751471477e-01,
58 
        -02.769154466281941332086016e-02,
59 
        -5.2482360935268931351485822e-01,
60 
        1.3002662421166552333051524e+00,
61 
        2.9698415875871901741575922e-01,
62 
        5.5025938468083370060258102e-01,
63 
        2.7629597861677201301553823e-01,
64 
        1.83559892257451475846656e-01,
65 
        -1.0523547536021497774980928e+00,
66 
}
67 
var asinh = []float64{
68 
        2.3083139124923523427628243e+00,
69 
        2.743551594301593620039021e+00,
70 
        -2.7345908534880091229413487e-01,
71 
        -2.3145157644718338650499085e+00,
72 
        2.9613652154015058521951083e+00,
73 
        1.7949041616585821933067568e+00,
74 
        2.3564032905983506405561554e+00,
75 
        1.7287118790768438878045346e+00,
76 
        1.3626658083714826013073193e+00,
77 
        -2.8581483626513914445234004e+00,
78 
}
79 
var atan = []float64{
80 
        1.372590262129621651920085e+00,
81 
        1.442290609645298083020664e+00,
82 
        -2.7011324359471758245192595e-01,
83 
        -1.3738077684543379452781531e+00,
84 
        1.4673921193587666049154681e+00,
85 
        1.2415173565870168649117764e+00,
86 
        1.3818396865615168979966498e+00,
87 
        1.2194305844639670701091426e+00,
88 
        1.0696031952318783760193244e+00,
89 
        -1.4561721938838084990898679e+00,
90 
}
91 
var atanh = []float64{
92 
        5.4651163712251938116878204e-01,
93 
        1.0299474112843111224914709e+00,
94 
        -2.7695084420740135145234906e-02,
95 
        -5.5072096119207195480202529e-01,
96 
        1.9943940993171843235906642e+00,
97 
        3.01448604578089708203017e-01,
98 
        5.8033427206942188834370595e-01,
99 
        2.7987997499441511013958297e-01,
100 
        1.8459947964298794318714228e-01,
101 
        -1.3273186910532645867272502e+00,
102 
}
103 
var atan2 = []float64{
104 
        1.1088291730037004444527075e+00,
105 
        9.1218183188715804018797795e-01,
106 
        1.5984772603216203736068915e+00,
107 
        2.0352918654092086637227327e+00,
108 
        8.0391819139044720267356014e-01,
109 
        1.2861075249894661588866752e+00,
110 
        1.0889904479131695712182587e+00,
111 
        1.3044821793397925293797357e+00,
112 
        1.3902530903455392306872261e+00,
113 
        2.2859857424479142655411058e+00,
114 
}
115 
var cbrt = []float64{
116 
        1.7075799841925094446722675e+00,
117 
        1.9779982212970353936691498e+00,
118 
        -6.5177429017779910853339447e-01,
119 
        -1.7111838886544019873338113e+00,
120 
        2.1279920909827937423960472e+00,
121 
        1.4303536770460741452312367e+00,
122 
        1.7357021059106154902341052e+00,
123 
        1.3972633462554328350552916e+00,
124 
        1.2221149580905388454977636e+00,
125 
        -2.0556003730500069110343596e+00,
126 
}
127 
var ceil = []float64{
128 
        5.0000000000000000e+00,
129 
        8.0000000000000000e+00,
130 
        0.0000000000000000e+00,
131 
        -5.0000000000000000e+00,
132 
        1.0000000000000000e+01,
133 
        3.0000000000000000e+00,
134 
        6.0000000000000000e+00,
135 
        3.0000000000000000e+00,
136 
        2.0000000000000000e+00,
137 
        -8.0000000000000000e+00,
138 
}
139 
var copysign = []float64{
140 
        -4.9790119248836735e+00,
141 
        -7.7388724745781045e+00,
142 
        -2.7688005719200159e-01,
143 
        -5.0106036182710749e+00,
144 
        -9.6362937071984173e+00,
145 
        -2.9263772392439646e+00,
146 
        -5.2290834314593066e+00,
147 
        -2.7279399104360102e+00,
148 
        -1.8253080916808550e+00,
149 
        -8.6859247685756013e+00,
150 
}
151 
var cos = []float64{
152 
        2.634752140995199110787593e-01,
153 
        1.148551260848219865642039e-01,
154 
        9.6191297325640768154550453e-01,
155 
        2.938141150061714816890637e-01,
156 
        -9.777138189897924126294461e-01,
157 
        -9.7693041344303219127199518e-01,
158 
        4.940088096948647263961162e-01,
159 
        -9.1565869021018925545016502e-01,
160 
        -2.517729313893103197176091e-01,
161 
        -7.39241351595676573201918e-01,
162 
}
163 
 
164 
// Results for 100000 * Pi + vf[i]
165 
var cosLarge = []float64{
166 
        2.634752141185559426744e-01,
167 
        1.14855126055543100712e-01,
168 
        9.61912973266488928113e-01,
169 
        2.9381411499556122552e-01,
170 
        -9.777138189880161924641e-01,
171 
        -9.76930413445147608049e-01,
172 
        4.940088097314976789841e-01,
173 
        -9.15658690217517835002e-01,
174 
        -2.51772931436786954751e-01,
175 
        -7.3924135157173099849e-01,
176 
}
177 
var cosh = []float64{
178 
        7.2668796942212842775517446e+01,
179 
        1.1479413465659254502011135e+03,
180 
        1.0385767908766418550935495e+00,
181 
        7.5000957789658051428857788e+01,
182 
        7.655246669605357888468613e+03,
183 
        9.3567491758321272072888257e+00,
184 
        9.331351599270605471131735e+01,
185 
        7.6833430994624643209296404e+00,
186 
        3.1829371625150718153881164e+00,
187 
        2.9595059261916188501640911e+03,
188 
}
189 
var erf = []float64{
190 
        5.1865354817738701906913566e-01,
191 
        7.2623875834137295116929844e-01,
192 
        -3.123458688281309990629839e-02,
193 
        -5.2143121110253302920437013e-01,
194 
        8.2704742671312902508629582e-01,
195 
        3.2101767558376376743993945e-01,
196 
        5.403990312223245516066252e-01,
197 
        3.0034702916738588551174831e-01,
198 
        2.0369924417882241241559589e-01,
199 
        -7.8069386968009226729944677e-01,
200 
}
201 
var erfc = []float64{
202 
        4.8134645182261298093086434e-01,
203 
        2.7376124165862704883070156e-01,
204 
        1.0312345868828130999062984e+00,
205 
        1.5214312111025330292043701e+00,
206 
        1.7295257328687097491370418e-01,
207 
        6.7898232441623623256006055e-01,
208 
        4.596009687776754483933748e-01,
209 
        6.9965297083261411448825169e-01,
210 
        7.9630075582117758758440411e-01,
211 
        1.7806938696800922672994468e+00,
212 
}
213 
var exp = []float64{
214 
        1.4533071302642137507696589e+02,
215 
        2.2958822575694449002537581e+03,
216 
        7.5814542574851666582042306e-01,
217 
        6.6668778421791005061482264e-03,
218 
        1.5310493273896033740861206e+04,
219 
        1.8659907517999328638667732e+01,
220 
        1.8662167355098714543942057e+02,
221 
        1.5301332413189378961665788e+01,
222 
        6.2047063430646876349125085e+00,
223 
        1.6894712385826521111610438e-04,
224 
}
225 
var expm1 = []float64{
226 
        5.105047796122957327384770212e-02,
227 
        8.046199708567344080562675439e-02,
228 
        -2.764970978891639815187418703e-03,
229 
        -4.8871434888875355394330300273e-02,
230 
        1.0115864277221467777117227494e-01,
231 
        2.969616407795910726014621657e-02,
232 
        5.368214487944892300914037972e-02,
233 
        2.765488851131274068067445335e-02,
234 
        1.842068661871398836913874273e-02,
235 
        -8.3193870863553801814961137573e-02,
236 
}
237 
var exp2 = []float64{
238 
        3.1537839463286288034313104e+01,
239 
        2.1361549283756232296144849e+02,
240 
        8.2537402562185562902577219e-01,
241 
        3.1021158628740294833424229e-02,
242 
        7.9581744110252191462569661e+02,
243 
        7.6019905892596359262696423e+00,
244 
        3.7506882048388096973183084e+01,
245 
        6.6250893439173561733216375e+00,
246 
        3.5438267900243941544605339e+00,
247 
        2.4281533133513300984289196e-03,
248 
}
249 
var fabs = []float64{
250 
        4.9790119248836735e+00,
251 
        7.7388724745781045e+00,
252 
        2.7688005719200159e-01,
253 
        5.0106036182710749e+00,
254 
        9.6362937071984173e+00,
255 
        2.9263772392439646e+00,
256 
        5.2290834314593066e+00,
257 
        2.7279399104360102e+00,
258 
        1.8253080916808550e+00,
259 
        8.6859247685756013e+00,
260 
}
261 
var fdim = []float64{
262 
        4.9790119248836735e+00,
263 
        7.7388724745781045e+00,
264 
        0.0000000000000000e+00,
265 
        0.0000000000000000e+00,
266 
        9.6362937071984173e+00,
267 
        2.9263772392439646e+00,
268 
        5.2290834314593066e+00,
269 
        2.7279399104360102e+00,
270 
        1.8253080916808550e+00,
271 
        0.0000000000000000e+00,
272 
}
273 
var floor = []float64{
274 
        4.0000000000000000e+00,
275 
        7.0000000000000000e+00,
276 
        -1.0000000000000000e+00,
277 
        -6.0000000000000000e+00,
278 
        9.0000000000000000e+00,
279 
        2.0000000000000000e+00,
280 
        5.0000000000000000e+00,
281 
        2.0000000000000000e+00,
282 
        1.0000000000000000e+00,
283 
        -9.0000000000000000e+00,
284 
}
285 
var fmod = []float64{
286 
        4.197615023265299782906368e-02,
287 
        2.261127525421895434476482e+00,
288 
        3.231794108794261433104108e-02,
289 
        4.989396381728925078391512e+00,
290 
        3.637062928015826201999516e-01,
291 
        1.220868282268106064236690e+00,
292 
        4.770916568540693347699744e+00,
293 
        1.816180268691969246219742e+00,
294 
        8.734595415957246977711748e-01,
295 
        1.314075231424398637614104e+00,
296 
}
297 
 
298 
type fi struct {
299 
        f float64
300 
        i int
301 
}
302 
 
303 
var frexp = []fi{
304 
        {6.2237649061045918750e-01, 3},
305 
        {9.6735905932226306250e-01, 3},
306 
        {-5.5376011438400318000e-01, -1},
307 
        {-6.2632545228388436250e-01, 3},
308 
        {6.02268356699901081250e-01, 4},
309 
        {7.3159430981099115000e-01, 2},
310 
        {6.5363542893241332500e-01, 3},
311 
        {6.8198497760900255000e-01, 2},
312 
        {9.1265404584042750000e-01, 1},
313 
        {-5.4287029803597508250e-01, 4},
314 
}
315 
var gamma = []float64{
316 
        2.3254348370739963835386613898e+01,
317 
        2.991153837155317076427529816e+03,
318 
        -4.561154336726758060575129109e+00,
319 
        7.719403468842639065959210984e-01,
320 
        1.6111876618855418534325755566e+05,
321 
        1.8706575145216421164173224946e+00,
322 
        3.4082787447257502836734201635e+01,
323 
        1.579733951448952054898583387e+00,
324 
        9.3834586598354592860187267089e-01,
325 
        -2.093995902923148389186189429e-05,
326 
}
327 
var j0 = []float64{
328 
        -1.8444682230601672018219338e-01,
329 
        2.27353668906331975435892e-01,
330 
        9.809259936157051116270273e-01,
331 
        -1.741170131426226587841181e-01,
332 
        -2.1389448451144143352039069e-01,
333 
        -2.340905848928038763337414e-01,
334 
        -1.0029099691890912094586326e-01,
335 
        -1.5466726714884328135358907e-01,
336 
        3.252650187653420388714693e-01,
337 
        -8.72218484409407250005360235e-03,
338 
}
339 
var j1 = []float64{
340 
        -3.251526395295203422162967e-01,
341 
        1.893581711430515718062564e-01,
342 
        -1.3711761352467242914491514e-01,
343 
        3.287486536269617297529617e-01,
344 
        1.3133899188830978473849215e-01,
345 
        3.660243417832986825301766e-01,
346 
        -3.4436769271848174665420672e-01,
347 
        4.329481396640773768835036e-01,
348 
        5.8181350531954794639333955e-01,
349 
        -2.7030574577733036112996607e-01,
350 
}
351 
var j2 = []float64{
352 
        5.3837518920137802565192769e-02,
353 
        -1.7841678003393207281244667e-01,
354 
        9.521746934916464142495821e-03,
355 
        4.28958355470987397983072e-02,
356 
        2.4115371837854494725492872e-01,
357 
        4.842458532394520316844449e-01,
358 
        -3.142145220618633390125946e-02,
359 
        4.720849184745124761189957e-01,
360 
        3.122312022520957042957497e-01,
361 
        7.096213118930231185707277e-02,
362 
}
363 
var jM3 = []float64{
364 
        -3.684042080996403091021151e-01,
365 
        2.8157665936340887268092661e-01,
366 
        4.401005480841948348343589e-04,
367 
        3.629926999056814081597135e-01,
368 
        3.123672198825455192489266e-02,
369 
        -2.958805510589623607540455e-01,
370 
        -3.2033177696533233403289416e-01,
371 
        -2.592737332129663376736604e-01,
372 
        -1.0241334641061485092351251e-01,
373 
        -2.3762660886100206491674503e-01,
374 
}
375 
var lgamma = []fi{
376 
        {3.146492141244545774319734e+00, 1},
377 
        {8.003414490659126375852113e+00, 1},
378 
        {1.517575735509779707488106e+00, -1},
379 
        {-2.588480028182145853558748e-01, 1},
380 
        {1.1989897050205555002007985e+01, 1},
381 
        {6.262899811091257519386906e-01, 1},
382 
        {3.5287924899091566764846037e+00, 1},
383 
        {4.5725644770161182299423372e-01, 1},
384 
        {-6.363667087767961257654854e-02, 1},
385 
        {-1.077385130910300066425564e+01, -1},
386 
}
387 
var log = []float64{
388 
        1.605231462693062999102599e+00,
389 
        2.0462560018708770653153909e+00,
390 
        -1.2841708730962657801275038e+00,
391 
        1.6115563905281545116286206e+00,
392 
        2.2655365644872016636317461e+00,
393 
        1.0737652208918379856272735e+00,
394 
        1.6542360106073546632707956e+00,
395 
        1.0035467127723465801264487e+00,
396 
        6.0174879014578057187016475e-01,
397 
        2.161703872847352815363655e+00,
398 
}
399 
var logb = []float64{
400 
        2.0000000000000000e+00,
401 
        2.0000000000000000e+00,
402 
        -2.0000000000000000e+00,
403 
        2.0000000000000000e+00,
404 
        3.0000000000000000e+00,
405 
        1.0000000000000000e+00,
406 
        2.0000000000000000e+00,
407 
        1.0000000000000000e+00,
408 
        0.0000000000000000e+00,
409 
        3.0000000000000000e+00,
410 
}
411 
var log10 = []float64{
412 
        6.9714316642508290997617083e-01,
413 
        8.886776901739320576279124e-01,
414 
        -5.5770832400658929815908236e-01,
415 
        6.998900476822994346229723e-01,
416 
        9.8391002850684232013281033e-01,
417 
        4.6633031029295153334285302e-01,
418 
        7.1842557117242328821552533e-01,
419 
        4.3583479968917773161304553e-01,
420 
        2.6133617905227038228626834e-01,
421 
        9.3881606348649405716214241e-01,
422 
}
423 
var log1p = []float64{
424 
        4.8590257759797794104158205e-02,
425 
        7.4540265965225865330849141e-02,
426 
        -2.7726407903942672823234024e-03,
427 
        -5.1404917651627649094953380e-02,
428 
        9.1998280672258624681335010e-02,
429 
        2.8843762576593352865894824e-02,
430 
        5.0969534581863707268992645e-02,
431 
        2.6913947602193238458458594e-02,
432 
        1.8088493239630770262045333e-02,
433 
        -9.0865245631588989681559268e-02,
434 
}
435 
var log2 = []float64{
436 
        2.3158594707062190618898251e+00,
437 
        2.9521233862883917703341018e+00,
438 
        -1.8526669502700329984917062e+00,
439 
        2.3249844127278861543568029e+00,
440 
        3.268478366538305087466309e+00,
441 
        1.5491157592596970278166492e+00,
442 
        2.3865580889631732407886495e+00,
443 
        1.447811865817085365540347e+00,
444 
        8.6813999540425116282815557e-01,
445 
        3.118679457227342224364709e+00,
446 
}
447 
var modf = [][2]float64{
448 
        {4.0000000000000000e+00, 9.7901192488367350108546816e-01},
449 
        {7.0000000000000000e+00, 7.3887247457810456552351752e-01},
450 
        {0.0000000000000000e+00, -2.7688005719200159404635997e-01},
451 
        {-5.0000000000000000e+00, -1.060361827107492160848778e-02},
452 
        {9.0000000000000000e+00, 6.3629370719841737980004837e-01},
453 
        {2.0000000000000000e+00, 9.2637723924396464525443662e-01},
454 
        {5.0000000000000000e+00, 2.2908343145930665230025625e-01},
455 
        {2.0000000000000000e+00, 7.2793991043601025126008608e-01},
456 
        {1.0000000000000000e+00, 8.2530809168085506044576505e-01},
457 
        {-8.0000000000000000e+00, -6.8592476857560136238589621e-01},
458 
}
459 
var nextafter = []float64{
460 
        4.97901192488367438926388786e+00,
461 
        7.73887247457810545370193722e+00,
462 
        -2.7688005719200153853520874e-01,
463 
        -5.01060361827107403343006808e+00,
464 
        9.63629370719841915615688777e+00,
465 
        2.92637723924396508934364647e+00,
466 
        5.22908343145930754047867595e+00,
467 
        2.72793991043601069534929593e+00,
468 
        1.82530809168085528249036997e+00,
469 
        -8.68592476857559958602905681e+00,
470 
}
471 
var pow = []float64{
472 
        9.5282232631648411840742957e+04,
473 
        5.4811599352999901232411871e+07,
474 
        5.2859121715894396531132279e-01,
475 
        9.7587991957286474464259698e-06,
476 
        4.328064329346044846740467e+09,
477 
        8.4406761805034547437659092e+02,
478 
        1.6946633276191194947742146e+05,
479 
        5.3449040147551939075312879e+02,
480 
        6.688182138451414936380374e+01,
481 
        2.0609869004248742886827439e-09,
482 
}
483 
var remainder = []float64{
484 
        4.197615023265299782906368e-02,
485 
        2.261127525421895434476482e+00,
486 
        3.231794108794261433104108e-02,
487 
        -2.120723654214984321697556e-02,
488 
        3.637062928015826201999516e-01,
489 
        1.220868282268106064236690e+00,
490 
        -4.581668629186133046005125e-01,
491 
        -9.117596417440410050403443e-01,
492 
        8.734595415957246977711748e-01,
493 
        1.314075231424398637614104e+00,
494 
}
495 
var signbit = []bool{
496 
        false,
497 
        false,
498 
        true,
499 
        true,
500 
        false,
501 
        false,
502 
        false,
503 
        false,
504 
        false,
505 
        true,
506 
}
507 
var sin = []float64{
508 
        -9.6466616586009283766724726e-01,
509 
        9.9338225271646545763467022e-01,
510 
        -2.7335587039794393342449301e-01,
511 
        9.5586257685042792878173752e-01,
512 
        -2.099421066779969164496634e-01,
513 
        2.135578780799860532750616e-01,
514 
        -8.694568971167362743327708e-01,
515 
        4.019566681155577786649878e-01,
516 
        9.6778633541687993721617774e-01,
517 
        -6.734405869050344734943028e-01,
518 
}
519 
 
520 
// Results for 100000 * Pi + vf[i]
521 
var sinLarge = []float64{
522 
        -9.646661658548936063912e-01,
523 
        9.933822527198506903752e-01,
524 
        -2.7335587036246899796e-01,
525 
        9.55862576853689321268e-01,
526 
        -2.099421066862688873691e-01,
527 
        2.13557878070308981163e-01,
528 
        -8.694568970959221300497e-01,
529 
        4.01956668098863248917e-01,
530 
        9.67786335404528727927e-01,
531 
        -6.7344058693131973066e-01,
532 
}
533 
var sinh = []float64{
534 
        7.2661916084208532301448439e+01,
535 
        1.1479409110035194500526446e+03,
536 
        -2.8043136512812518927312641e-01,
537 
        -7.499429091181587232835164e+01,
538 
        7.6552466042906758523925934e+03,
539 
        9.3031583421672014313789064e+00,
540 
        9.330815755828109072810322e+01,
541 
        7.6179893137269146407361477e+00,
542 
        3.021769180549615819524392e+00,
543 
        -2.95950575724449499189888e+03,
544 
}
545 
var sqrt = []float64{
546 
        2.2313699659365484748756904e+00,
547 
        2.7818829009464263511285458e+00,
548 
        5.2619393496314796848143251e-01,
549 
        2.2384377628763938724244104e+00,
550 
        3.1042380236055381099288487e+00,
551 
        1.7106657298385224403917771e+00,
552 
        2.286718922705479046148059e+00,
553 
        1.6516476350711159636222979e+00,
554 
        1.3510396336454586262419247e+00,
555 
        2.9471892997524949215723329e+00,
556 
}
557 
var tan = []float64{
558 
        -3.661316565040227801781974e+00,
559 
        8.64900232648597589369854e+00,
560 
        -2.8417941955033612725238097e-01,
561 
        3.253290185974728640827156e+00,
562 
        2.147275640380293804770778e-01,
563 
        -2.18600910711067004921551e-01,
564 
        -1.760002817872367935518928e+00,
565 
        -4.389808914752818126249079e-01,
566 
        -3.843885560201130679995041e+00,
567 
        9.10988793377685105753416e-01,
568 
}
569 
 
570 
// Results for 100000 * Pi + vf[i]
571 
var tanLarge = []float64{
572 
        -3.66131656475596512705e+00,
573 
        8.6490023287202547927e+00,
574 
        -2.841794195104782406e-01,
575 
        3.2532901861033120983e+00,
576 
        2.14727564046880001365e-01,
577 
        -2.18600910700688062874e-01,
578 
        -1.760002817699722747043e+00,
579 
        -4.38980891453536115952e-01,
580 
        -3.84388555942723509071e+00,
581 
        9.1098879344275101051e-01,
582 
}
583 
var tanh = []float64{
584 
        9.9990531206936338549262119e-01,
585 
        9.9999962057085294197613294e-01,
586 
        -2.7001505097318677233756845e-01,
587 
        -9.9991110943061718603541401e-01,
588 
        9.9999999146798465745022007e-01,
589 
        9.9427249436125236705001048e-01,
590 
        9.9994257600983138572705076e-01,
591 
        9.9149409509772875982054701e-01,
592 
        9.4936501296239685514466577e-01,
593 
        -9.9999994291374030946055701e-01,
594 
}
595 
var trunc = []float64{
596 
        4.0000000000000000e+00,
597 
        7.0000000000000000e+00,
598 
        -0.0000000000000000e+00,
599 
        -5.0000000000000000e+00,
600 
        9.0000000000000000e+00,
601 
        2.0000000000000000e+00,
602 
        5.0000000000000000e+00,
603 
        2.0000000000000000e+00,
604 
        1.0000000000000000e+00,
605 
        -8.0000000000000000e+00,
606 
}
607 
var y0 = []float64{
608 
        -3.053399153780788357534855e-01,
609 
        1.7437227649515231515503649e-01,
610 
        -8.6221781263678836910392572e-01,
611 
        -3.100664880987498407872839e-01,
612 
        1.422200649300982280645377e-01,
613 
        4.000004067997901144239363e-01,
614 
        -3.3340749753099352392332536e-01,
615 
        4.5399790746668954555205502e-01,
616 
        4.8290004112497761007536522e-01,
617 
        2.7036697826604756229601611e-01,
618 
}
619 
var y1 = []float64{
620 
        0.15494213737457922210218611,
621 
        -0.2165955142081145245075746,
622 
        -2.4644949631241895201032829,
623 
        0.1442740489541836405154505,
624 
        0.2215379960518984777080163,
625 
        0.3038800915160754150565448,
626 
        0.0691107642452362383808547,
627 
        0.2380116417809914424860165,
628 
        -0.20849492979459761009678934,
629 
        0.0242503179793232308250804,
630 
}
631 
var y2 = []float64{
632 
        0.3675780219390303613394936,
633 
        -0.23034826393250119879267257,
634 
        -16.939677983817727205631397,
635 
        0.367653980523052152867791,
636 
        -0.0962401471767804440353136,
637 
        -0.1923169356184851105200523,
638 
        0.35984072054267882391843766,
639 
        -0.2794987252299739821654982,
640 
        -0.7113490692587462579757954,
641 
        -0.2647831587821263302087457,
642 
}
643 
var yM3 = []float64{
644 
        -0.14035984421094849100895341,
645 
        -0.097535139617792072703973,
646 
        242.25775994555580176377379,
647 
        -0.1492267014802818619511046,
648 
        0.26148702629155918694500469,
649 
        0.56675383593895176530394248,
650 
        -0.206150264009006981070575,
651 
        0.64784284687568332737963658,
652 
        1.3503631555901938037008443,
653 
        0.1461869756579956803341844,
654 
}
655 
 
656 
// arguments and expected results for special cases
657 
var vfacosSC = []float64{
658 
        -Pi,
659 
        1,
660 
        Pi,
661 
        NaN(),
662 
}
663 
var acosSC = []float64{
664 
        NaN(),
665 
        0,
666 
        NaN(),
667 
        NaN(),
668 
}
669 
 
670 
var vfacoshSC = []float64{
671 
        Inf(-1),
672 
        0.5,
673 
        1,
674 
        Inf(1),
675 
        NaN(),
676 
}
677 
var acoshSC = []float64{
678 
        NaN(),
679 
        NaN(),
680 
        0,
681 
        Inf(1),
682 
        NaN(),
683 
}
684 
 
685 
var vfasinSC = []float64{
686 
        -Pi,
687 
        Copysign(0, -1),
688 
        0,
689 
        Pi,
690 
        NaN(),
691 
}
692 
var asinSC = []float64{
693 
        NaN(),
694 
        Copysign(0, -1),
695 
        0,
696 
        NaN(),
697 
        NaN(),
698 
}
699 
 
700 
var vfasinhSC = []float64{
701 
        Inf(-1),
702 
        Copysign(0, -1),
703 
        0,
704 
        Inf(1),
705 
        NaN(),
706 
}
707 
var asinhSC = []float64{
708 
        Inf(-1),
709 
        Copysign(0, -1),
710 
        0,
711 
        Inf(1),
712 
        NaN(),
713 
}
714 
 
715 
var vfatanSC = []float64{
716 
        Inf(-1),
717 
        Copysign(0, -1),
718 
        0,
719 
        Inf(1),
720 
        NaN(),
721 
}
722 
var atanSC = []float64{
723 
        -Pi / 2,
724 
        Copysign(0, -1),
725 
        0,
726 
        Pi / 2,
727 
        NaN(),
728 
}
729 
 
730 
var vfatanhSC = []float64{
731 
        Inf(-1),
732 
        -Pi,
733 
        -1,
734 
        Copysign(0, -1),
735 
        0,
736 
        1,
737 
        Pi,
738 
        Inf(1),
739 
        NaN(),
740 
}
741 
var atanhSC = []float64{
742 
        NaN(),
743 
        NaN(),
744 
        Inf(-1),
745 
        Copysign(0, -1),
746 
        0,
747 
        Inf(1),
748 
        NaN(),
749 
        NaN(),
750 
        NaN(),
751 
}
752 
var vfatan2SC = [][2]float64{
753 
        {Inf(-1), Inf(-1)},
754 
        {Inf(-1), -Pi},
755 
        {Inf(-1), 0},
756 
        {Inf(-1), +Pi},
757 
        {Inf(-1), Inf(1)},
758 
        {Inf(-1), NaN()},
759 
        {-Pi, Inf(-1)},
760 
        {-Pi, 0},
761 
        {-Pi, Inf(1)},
762 
        {-Pi, NaN()},
763 
        {Copysign(0, -1), Inf(-1)},
764 
        {Copysign(0, -1), -Pi},
765 
        {Copysign(0, -1), Copysign(0, -1)},
766 
        {Copysign(0, -1), 0},
767 
        {Copysign(0, -1), +Pi},
768 
        {Copysign(0, -1), Inf(1)},
769 
        {Copysign(0, -1), NaN()},
770 
        {0, Inf(-1)},
771 
        {0, -Pi},
772 
        {0, Copysign(0, -1)},
773 
        {0, 0},
774 
        {0, +Pi},
775 
        {0, Inf(1)},
776 
        {0, NaN()},
777 
        {+Pi, Inf(-1)},
778 
        {+Pi, 0},
779 
        {+Pi, Inf(1)},
780 
        {+Pi, NaN()},
781 
        {Inf(1), Inf(-1)},
782 
        {Inf(1), -Pi},
783 
        {Inf(1), 0},
784 
        {Inf(1), +Pi},
785 
        {Inf(1), Inf(1)},
786 
        {Inf(1), NaN()},
787 
        {NaN(), NaN()},
788 
}
789 
var atan2SC = []float64{
790 
        -3 * Pi / 4,     // atan2(-Inf, -Inf)
791 
        -Pi / 2,         // atan2(-Inf, -Pi)
792 
        -Pi / 2,         // atan2(-Inf, +0)
793 
        -Pi / 2,         // atan2(-Inf, +Pi)
794 
        -Pi / 4,         // atan2(-Inf, +Inf)
795 
        NaN(),           // atan2(-Inf, NaN)
796 
        -Pi,             // atan2(-Pi, -Inf)
797 
        -Pi / 2,         // atan2(-Pi, +0)
798 
        Copysign(0, -1), // atan2(-Pi, Inf)
799 
        NaN(),           // atan2(-Pi, NaN)
800 
        -Pi,             // atan2(-0, -Inf)
801 
        -Pi,             // atan2(-0, -Pi)
802 
        -Pi,             // atan2(-0, -0)
803 
        Copysign(0, -1), // atan2(-0, +0)
804 
        Copysign(0, -1), // atan2(-0, +Pi)
805 
        Copysign(0, -1), // atan2(-0, +Inf)
806 
        NaN(),           // atan2(-0, NaN)
807 
        Pi,              // atan2(+0, -Inf)
808 
        Pi,              // atan2(+0, -Pi)
809 
        Pi,              // atan2(+0, -0)
810 
        0,               // atan2(+0, +0)
811 
        0,               // atan2(+0, +Pi)
812 
        0,               // atan2(+0, +Inf)
813 
        NaN(),           // atan2(+0, NaN)
814 
        Pi,              // atan2(+Pi, -Inf)
815 
        Pi / 2,          // atan2(+Pi, +0)
816 
        0,               // atan2(+Pi, +Inf)
817 
        NaN(),           // atan2(+Pi, NaN)
818 
        3 * Pi / 4,      // atan2(+Inf, -Inf)
819 
        Pi / 2,          // atan2(+Inf, -Pi)
820 
        Pi / 2,          // atan2(+Inf, +0)
821 
        Pi / 2,          // atan2(+Inf, +Pi)
822 
        Pi / 4,          // atan2(+Inf, +Inf)
823 
        NaN(),           // atan2(+Inf, NaN)
824 
        NaN(),           // atan2(NaN, NaN)
825 
}
826 
 
827 
var vfcbrtSC = []float64{
828 
        Inf(-1),
829 
        Copysign(0, -1),
830 
        0,
831 
        Inf(1),
832 
        NaN(),
833 
}
834 
var cbrtSC = []float64{
835 
        Inf(-1),
836 
        Copysign(0, -1),
837 
        0,
838 
        Inf(1),
839 
        NaN(),
840 
}
841 
 
842 
var vfceilSC = []float64{
843 
        Inf(-1),
844 
        Copysign(0, -1),
845 
        0,
846 
        Inf(1),
847 
        NaN(),
848 
}
849 
var ceilSC = []float64{
850 
        Inf(-1),
851 
        Copysign(0, -1),
852 
        0,
853 
        Inf(1),
854 
        NaN(),
855 
}
856 
 
857 
var vfcopysignSC = []float64{
858 
        Inf(-1),
859 
        Inf(1),
860 
        NaN(),
861 
}
862 
var copysignSC = []float64{
863 
        Inf(-1),
864 
        Inf(-1),
865 
        NaN(),
866 
}
867 
 
868 
var vfcosSC = []float64{
869 
        Inf(-1),
870 
        Inf(1),
871 
        NaN(),
872 
}
873 
var cosSC = []float64{
874 
        NaN(),
875 
        NaN(),
876 
        NaN(),
877 
}
878 
 
879 
var vfcoshSC = []float64{
880 
        Inf(-1),
881 
        Copysign(0, -1),
882 
        0,
883 
        Inf(1),
884 
        NaN(),
885 
}
886 
var coshSC = []float64{
887 
        Inf(1),
888 
        1,
889 
        1,
890 
        Inf(1),
891 
        NaN(),
892 
}
893 
 
894 
var vferfSC = []float64{
895 
        Inf(-1),
896 
        Copysign(0, -1),
897 
        0,
898 
        Inf(1),
899 
        NaN(),
900 
}
901 
var erfSC = []float64{
902 
        -1,
903 
        Copysign(0, -1),
904 
        0,
905 
        1,
906 
        NaN(),
907 
}
908 
 
909 
var vferfcSC = []float64{
910 
        Inf(-1),
911 
        Inf(1),
912 
        NaN(),
913 
}
914 
var erfcSC = []float64{
915 
        2,
916 
        0,
917 
        NaN(),
918 
}
919 
 
920 
var vfexpSC = []float64{
921 
        Inf(-1),
922 
        -2000,
923 
        2000,
924 
        Inf(1),
925 
        NaN(),
926 
}
927 
var expSC = []float64{
928 
        0,
929 
        0,
930 
        Inf(1),
931 
        Inf(1),
932 
        NaN(),
933 
}
934 
 
935 
var vfexpm1SC = []float64{
936 
        Inf(-1),
937 
        Copysign(0, -1),
938 
        0,
939 
        Inf(1),
940 
        NaN(),
941 
}
942 
var expm1SC = []float64{
943 
        -1,
944 
        Copysign(0, -1),
945 
        0,
946 
        Inf(1),
947 
        NaN(),
948 
}
949 
 
950 
var vffabsSC = []float64{
951 
        Inf(-1),
952 
        Copysign(0, -1),
953 
        0,
954 
        Inf(1),
955 
        NaN(),
956 
}
957 
var fabsSC = []float64{
958 
        Inf(1),
959 
        0,
960 
        0,
961 
        Inf(1),
962 
        NaN(),
963 
}
964 
 
965 
var vffdimSC = [][2]float64{
966 
        {Inf(-1), Inf(-1)},
967 
        {Inf(-1), Inf(1)},
968 
        {Inf(-1), NaN()},
969 
        {Copysign(0, -1), Copysign(0, -1)},
970 
        {Copysign(0, -1), 0},
971 
        {0, Copysign(0, -1)},
972 
        {0, 0},
973 
        {Inf(1), Inf(-1)},
974 
        {Inf(1), Inf(1)},
975 
        {Inf(1), NaN()},
976 
        {NaN(), Inf(-1)},
977 
        {NaN(), Copysign(0, -1)},
978 
        {NaN(), 0},
979 
        {NaN(), Inf(1)},
980 
        {NaN(), NaN()},
981 
}
982 
var fdimSC = []float64{
983 
        NaN(),
984 
        0,
985 
        NaN(),
986 
        0,
987 
        0,
988 
        0,
989 
        0,
990 
        Inf(1),
991 
        NaN(),
992 
        NaN(),
993 
        NaN(),
994 
        NaN(),
995 
        NaN(),
996 
        NaN(),
997 
        NaN(),
998 
}
999 
var fmaxSC = []float64{
1000 
        Inf(-1),
1001 
        Inf(1),
1002 
        NaN(),
1003 
        Copysign(0, -1),
1004 
        0,
1005 
        0,
1006 
        0,
1007 
        Inf(1),
1008 
        Inf(1),
1009 
        Inf(1),
1010 
        NaN(),
1011 
        NaN(),
1012 
        NaN(),
1013 
        Inf(1),
1014 
        NaN(),
1015 
}
1016 
var fminSC = []float64{
1017 
        Inf(-1),
1018 
        Inf(-1),
1019 
        Inf(-1),
1020 
        Copysign(0, -1),
1021 
        Copysign(0, -1),
1022 
        Copysign(0, -1),
1023 
        0,
1024 
        Inf(-1),
1025 
        Inf(1),
1026 
        NaN(),
1027 
        Inf(-1),
1028 
        NaN(),
1029 
        NaN(),
1030 
        NaN(),
1031 
        NaN(),
1032 
}
1033 
 
1034 
var vffmodSC = [][2]float64{
1035 
        {Inf(-1), Inf(-1)},
1036 
        {Inf(-1), -Pi},
1037 
        {Inf(-1), 0},
1038 
        {Inf(-1), Pi},
1039 
        {Inf(-1), Inf(1)},
1040 
        {Inf(-1), NaN()},
1041 
        {-Pi, Inf(-1)},
1042 
        {-Pi, 0},
1043 
        {-Pi, Inf(1)},
1044 
        {-Pi, NaN()},
1045 
        {Copysign(0, -1), Inf(-1)},
1046 
        {Copysign(0, -1), 0},
1047 
        {Copysign(0, -1), Inf(1)},
1048 
        {Copysign(0, -1), NaN()},
1049 
        {0, Inf(-1)},
1050 
        {0, 0},
1051 
        {0, Inf(1)},
1052 
        {0, NaN()},
1053 
        {Pi, Inf(-1)},
1054 
        {Pi, 0},
1055 
        {Pi, Inf(1)},
1056 
        {Pi, NaN()},
1057 
        {Inf(1), Inf(-1)},
1058 
        {Inf(1), -Pi},
1059 
        {Inf(1), 0},
1060 
        {Inf(1), Pi},
1061 
        {Inf(1), Inf(1)},
1062 
        {Inf(1), NaN()},
1063 
        {NaN(), Inf(-1)},
1064 
        {NaN(), -Pi},
1065 
        {NaN(), 0},
1066 
        {NaN(), Pi},
1067 
        {NaN(), Inf(1)},
1068 
        {NaN(), NaN()},
1069 
}
1070 
var fmodSC = []float64{
1071 
        NaN(),           // fmod(-Inf, -Inf)
1072 
        NaN(),           // fmod(-Inf, -Pi)
1073 
        NaN(),           // fmod(-Inf, 0)
1074 
        NaN(),           // fmod(-Inf, Pi)
1075 
        NaN(),           // fmod(-Inf, +Inf)
1076 
        NaN(),           // fmod(-Inf, NaN)
1077 
        -Pi,             // fmod(-Pi, -Inf)
1078 
        NaN(),           // fmod(-Pi, 0)
1079 
        -Pi,             // fmod(-Pi, +Inf)
1080 
        NaN(),           // fmod(-Pi, NaN)
1081 
        Copysign(0, -1), // fmod(-0, -Inf)
1082 
        NaN(),           // fmod(-0, 0)
1083 
        Copysign(0, -1), // fmod(-0, Inf)
1084 
        NaN(),           // fmod(-0, NaN)
1085 
        0,               // fmod(0, -Inf)
1086 
        NaN(),           // fmod(0, 0)
1087 
        0,               // fmod(0, +Inf)
1088 
        NaN(),           // fmod(0, NaN)
1089 
        Pi,              // fmod(Pi, -Inf)
1090 
        NaN(),           // fmod(Pi, 0)
1091 
        Pi,              // fmod(Pi, +Inf)
1092 
        NaN(),           // fmod(Pi, NaN)
1093 
        NaN(),           // fmod(+Inf, -Inf)
1094 
        NaN(),           // fmod(+Inf, -Pi)
1095 
        NaN(),           // fmod(+Inf, 0)
1096 
        NaN(),           // fmod(+Inf, Pi)
1097 
        NaN(),           // fmod(+Inf, +Inf)
1098 
        NaN(),           // fmod(+Inf, NaN)
1099 
        NaN(),           // fmod(NaN, -Inf)
1100 
        NaN(),           // fmod(NaN, -Pi)
1101 
        NaN(),           // fmod(NaN, 0)
1102 
        NaN(),           // fmod(NaN, Pi)
1103 
        NaN(),           // fmod(NaN, +Inf)
1104 
        NaN(),           // fmod(NaN, NaN)
1105 
}
1106 
 
1107 
var vffrexpSC = []float64{
1108 
        Inf(-1),
1109 
        Copysign(0, -1),
1110 
        0,
1111 
        Inf(1),
1112 
        NaN(),
1113 
}
1114 
var frexpSC = []fi{
1115 
        {Inf(-1), 0},
1116 
        {Copysign(0, -1), 0},
1117 
        {0, 0},
1118 
        {Inf(1), 0},
1119 
        {NaN(), 0},
1120 
}
1121 
 
1122 
var vfgammaSC = []float64{
1123 
        Inf(-1),
1124 
        -3,
1125 
        Copysign(0, -1),
1126 
        0,
1127 
        Inf(1),
1128 
        NaN(),
1129 
}
1130 
var gammaSC = []float64{
1131 
        Inf(-1),
1132 
        Inf(1),
1133 
        Inf(1),
1134 
        Inf(1),
1135 
        Inf(1),
1136 
        NaN(),
1137 
}
1138 
 
1139 
var vfhypotSC = [][2]float64{
1140 
        {Inf(-1), Inf(-1)},
1141 
        {Inf(-1), 0},
1142 
        {Inf(-1), Inf(1)},
1143 
        {Inf(-1), NaN()},
1144 
        {Copysign(0, -1), Copysign(0, -1)},
1145 
        {Copysign(0, -1), 0},
1146 
        {0, Copysign(0, -1)},
1147 
        {0, 0}, // +0, +0
1148 
        {0, Inf(-1)},
1149 
        {0, Inf(1)},
1150 
        {0, NaN()},
1151 
        {Inf(1), Inf(-1)},
1152 
        {Inf(1), 0},
1153 
        {Inf(1), Inf(1)},
1154 
        {Inf(1), NaN()},
1155 
        {NaN(), Inf(-1)},
1156 
        {NaN(), 0},
1157 
        {NaN(), Inf(1)},
1158 
        {NaN(), NaN()},
1159 
}
1160 
var hypotSC = []float64{
1161 
        Inf(1),
1162 
        Inf(1),
1163 
        Inf(1),
1164 
        Inf(1),
1165 
        0,
1166 
        0,
1167 
        0,
1168 
        0,
1169 
        Inf(1),
1170 
        Inf(1),
1171 
        NaN(),
1172 
        Inf(1),
1173 
        Inf(1),
1174 
        Inf(1),
1175 
        Inf(1),
1176 
        Inf(1),
1177 
        NaN(),
1178 
        Inf(1),
1179 
        NaN(),
1180 
}
1181 
 
1182 
var vfilogbSC = []float64{
1183 
        Inf(-1),
1184 
        0,
1185 
        Inf(1),
1186 
        NaN(),
1187 
}
1188 
var ilogbSC = []int{
1189 
        MaxInt32,
1190 
        MinInt32,
1191 
        MaxInt32,
1192 
        MaxInt32,
1193 
}
1194 
 
1195 
var vfj0SC = []float64{
1196 
        Inf(-1),
1197 
        0,
1198 
        Inf(1),
1199 
        NaN(),
1200 
}
1201 
var j0SC = []float64{
1202 
        0,
1203 
        1,
1204 
        0,
1205 
        NaN(),
1206 
}
1207 
var j1SC = []float64{
1208 
        0,
1209 
        0,
1210 
        0,
1211 
        NaN(),
1212 
}
1213 
var j2SC = []float64{
1214 
        0,
1215 
        0,
1216 
        0,
1217 
        NaN(),
1218 
}
1219 
var jM3SC = []float64{
1220 
        0,
1221 
        0,
1222 
        0,
1223 
        NaN(),
1224 
}
1225 
 
1226 
var vfldexpSC = []fi{
1227 
        {0, 0},
1228 
        {0, -1075},
1229 
        {0, 1024},
1230 
        {Copysign(0, -1), 0},
1231 
        {Copysign(0, -1), -1075},
1232 
        {Copysign(0, -1), 1024},
1233 
        {Inf(1), 0},
1234 
        {Inf(1), -1024},
1235 
        {Inf(-1), 0},
1236 
        {Inf(-1), -1024},
1237 
        {NaN(), -1024},
1238 
}
1239 
var ldexpSC = []float64{
1240 
        0,
1241 
        0,
1242 
        0,
1243 
        Copysign(0, -1),
1244 
        Copysign(0, -1),
1245 
        Copysign(0, -1),
1246 
        Inf(1),
1247 
        Inf(1),
1248 
        Inf(-1),
1249 
        Inf(-1),
1250 
        NaN(),
1251 
}
1252 
 
1253 
var vflgammaSC = []float64{
1254 
        Inf(-1),
1255 
        -3,
1256 
        0,
1257 
        1,
1258 
        2,
1259 
        Inf(1),
1260 
        NaN(),
1261 
}
1262 
var lgammaSC = []fi{
1263 
        {Inf(-1), 1},
1264 
        {Inf(1), 1},
1265 
        {Inf(1), 1},
1266 
        {0, 1},
1267 
        {0, 1},
1268 
        {Inf(1), 1},
1269 
        {NaN(), 1},
1270 
}
1271 
 
1272 
var vflogSC = []float64{
1273 
        Inf(-1),
1274 
        -Pi,
1275 
        Copysign(0, -1),
1276 
        0,
1277 
        1,
1278 
        Inf(1),
1279 
        NaN(),
1280 
}
1281 
var logSC = []float64{
1282 
        NaN(),
1283 
        NaN(),
1284 
        Inf(-1),
1285 
        Inf(-1),
1286 
        0,
1287 
        Inf(1),
1288 
        NaN(),
1289 
}
1290 
 
1291 
var vflogbSC = []float64{
1292 
        Inf(-1),
1293 
        0,
1294 
        Inf(1),
1295 
        NaN(),
1296 
}
1297 
var logbSC = []float64{
1298 
        Inf(1),
1299 
        Inf(-1),
1300 
        Inf(1),
1301 
        NaN(),
1302 
}
1303 
 
1304 
var vflog1pSC = []float64{
1305 
        Inf(-1),
1306 
        -Pi,
1307 
        -1,
1308 
        Copysign(0, -1),
1309 
        0,
1310 
        Inf(1),
1311 
        NaN(),
1312 
}
1313 
var log1pSC = []float64{
1314 
        NaN(),
1315 
        NaN(),
1316 
        Inf(-1),
1317 
        Copysign(0, -1),
1318 
        0,
1319 
        Inf(1),
1320 
        NaN(),
1321 
}
1322 
 
1323 
var vfmodfSC = []float64{
1324 
        Inf(-1),
1325 
        Inf(1),
1326 
        NaN(),
1327 
}
1328 
var modfSC = [][2]float64{
1329 
        {Inf(-1), NaN()}, // [2]float64{Copysign(0, -1), Inf(-1)},
1330 
        {Inf(1), NaN()},  // [2]float64{0, Inf(1)},
1331 
        {NaN(), NaN()},
1332 
}
1333 
 
1334 
var vfnextafterSC = [][2]float64{
1335 
        {0, 0},
1336 
        {0, Copysign(0, -1)},
1337 
        {0, -1},
1338 
        {0, NaN()},
1339 
        {Copysign(0, -1), 1},
1340 
        {Copysign(0, -1), 0},
1341 
        {Copysign(0, -1), Copysign(0, -1)},
1342 
        {Copysign(0, -1), -1},
1343 
        {NaN(), 0},
1344 
        {NaN(), NaN()},
1345 
}
1346 
var nextafterSC = []float64{
1347 
        0,
1348 
        0,
1349 
        -4.9406564584124654418e-324, // Float64frombits(0x8000000000000001)
1350 
        NaN(),
1351 
        4.9406564584124654418e-324, // Float64frombits(0x0000000000000001)
1352 
        Copysign(0, -1),
1353 
        Copysign(0, -1),
1354 
        -4.9406564584124654418e-324, // Float64frombits(0x8000000000000001)
1355 
        NaN(),
1356 
        NaN(),
1357 
}
1358 
 
1359 
var vfpowSC = [][2]float64{
1360 
        {Inf(-1), -Pi},
1361 
        {Inf(-1), -3},
1362 
        {Inf(-1), Copysign(0, -1)},
1363 
        {Inf(-1), 0},
1364 
        {Inf(-1), 1},
1365 
        {Inf(-1), 3},
1366 
        {Inf(-1), Pi},
1367 
        {Inf(-1), NaN()},
1368 
 
1369 
        {-Pi, Inf(-1)},
1370 
        {-Pi, -Pi},
1371 
        {-Pi, Copysign(0, -1)},
1372 
        {-Pi, 0},
1373 
        {-Pi, 1},
1374 
        {-Pi, Pi},
1375 
        {-Pi, Inf(1)},
1376 
        {-Pi, NaN()},
1377 
 
1378 
        {-1, Inf(-1)},
1379 
        {-1, Inf(1)},
1380 
        {-1, NaN()},
1381 
        {-1 / 2, Inf(-1)},
1382 
        {-1 / 2, Inf(1)},
1383 
        {Copysign(0, -1), Inf(-1)},
1384 
        {Copysign(0, -1), -Pi},
1385 
        {Copysign(0, -1), -3},
1386 
        {Copysign(0, -1), 3},
1387 
        {Copysign(0, -1), Pi},
1388 
        {Copysign(0, -1), Inf(1)},
1389 
 
1390 
        {0, Inf(-1)},
1391 
        {0, -Pi},
1392 
        {0, -3},
1393 
        {0, Copysign(0, -1)},
1394 
        {0, 0},
1395 
        {0, 3},
1396 
        {0, Pi},
1397 
        {0, Inf(1)},
1398 
        {0, NaN()},
1399 
 
1400 
        {1 / 2, Inf(-1)},
1401 
        {1 / 2, Inf(1)},
1402 
        {1, Inf(-1)},
1403 
        {1, Inf(1)},
1404 
        {1, NaN()},
1405 
 
1406 
        {Pi, Inf(-1)},
1407 
        {Pi, Copysign(0, -1)},
1408 
        {Pi, 0},
1409 
        {Pi, 1},
1410 
        {Pi, Inf(1)},
1411 
        {Pi, NaN()},
1412 
        {Inf(1), -Pi},
1413 
        {Inf(1), Copysign(0, -1)},
1414 
        {Inf(1), 0},
1415 
        {Inf(1), 1},
1416 
        {Inf(1), Pi},
1417 
        {Inf(1), NaN()},
1418 
        {NaN(), -Pi},
1419 
        {NaN(), Copysign(0, -1)},
1420 
        {NaN(), 0},
1421 
        {NaN(), 1},
1422 
        {NaN(), Pi},
1423 
        {NaN(), NaN()},
1424 
}
1425 
var powSC = []float64{
1426 
        0,               // pow(-Inf, -Pi)
1427 
        Copysign(0, -1), // pow(-Inf, -3)
1428 
        1,               // pow(-Inf, -0)
1429 
        1,               // pow(-Inf, +0)
1430 
        Inf(-1),         // pow(-Inf, 1)
1431 
        Inf(-1),         // pow(-Inf, 3)
1432 
        Inf(1),          // pow(-Inf, Pi)
1433 
        NaN(),           // pow(-Inf, NaN)
1434 
        0,               // pow(-Pi, -Inf)
1435 
        NaN(),           // pow(-Pi, -Pi)
1436 
        1,               // pow(-Pi, -0)
1437 
        1,               // pow(-Pi, +0)
1438 
        -Pi,             // pow(-Pi, 1)
1439 
        NaN(),           // pow(-Pi, Pi)
1440 
        Inf(1),          // pow(-Pi, +Inf)
1441 
        NaN(),           // pow(-Pi, NaN)
1442 
        1,               // pow(-1, -Inf) IEEE 754-2008
1443 
        1,               // pow(-1, +Inf) IEEE 754-2008
1444 
        NaN(),           // pow(-1, NaN)
1445 
        Inf(1),          // pow(-1/2, -Inf)
1446 
        0,               // pow(-1/2, +Inf)
1447 
        Inf(1),          // pow(-0, -Inf)
1448 
        Inf(1),          // pow(-0, -Pi)
1449 
        Inf(-1),         // pow(-0, -3) IEEE 754-2008
1450 
        Copysign(0, -1), // pow(-0, 3) IEEE 754-2008
1451 
        0,               // pow(-0, +Pi)
1452 
        0,               // pow(-0, +Inf)
1453 
        Inf(1),          // pow(+0, -Inf)
1454 
        Inf(1),          // pow(+0, -Pi)
1455 
        Inf(1),          // pow(+0, -3)
1456 
        1,               // pow(+0, -0)
1457 
        1,               // pow(+0, +0)
1458 
        0,               // pow(+0, 3)
1459 
        0,               // pow(+0, +Pi)
1460 
        0,               // pow(+0, +Inf)
1461 
        NaN(),           // pow(+0, NaN)
1462 
        Inf(1),          // pow(1/2, -Inf)
1463 
        0,               // pow(1/2, +Inf)
1464 
        1,               // pow(1, -Inf) IEEE 754-2008
1465 
        1,               // pow(1, +Inf) IEEE 754-2008
1466 
        1,               // pow(1, NaN) IEEE 754-2008
1467 
        0,               // pow(+Pi, -Inf)
1468 
        1,               // pow(+Pi, -0)
1469 
        1,               // pow(+Pi, +0)
1470 
        Pi,              // pow(+Pi, 1)
1471 
        Inf(1),          // pow(+Pi, +Inf)
1472 
        NaN(),           // pow(+Pi, NaN)
1473 
        0,               // pow(+Inf, -Pi)
1474 
        1,               // pow(+Inf, -0)
1475 
        1,               // pow(+Inf, +0)
1476 
        Inf(1),          // pow(+Inf, 1)
1477 
        Inf(1),          // pow(+Inf, Pi)
1478 
        NaN(),           // pow(+Inf, NaN)
1479 
        NaN(),           // pow(NaN, -Pi)
1480 
        1,               // pow(NaN, -0)
1481 
        1,               // pow(NaN, +0)
1482 
        NaN(),           // pow(NaN, 1)
1483 
        NaN(),           // pow(NaN, +Pi)
1484 
        NaN(),           // pow(NaN, NaN)
1485 
}
1486 
 
1487 
var vfpow10SC = []int{
1488 
        MinInt32,
1489 
        MaxInt32,
1490 
        -325,
1491 
        309,
1492 
}
1493 
 
1494 
var pow10SC = []float64{
1495 
        0,      // pow10(MinInt32)
1496 
        Inf(1), // pow10(MaxInt32)
1497 
        0,      // pow10(-325)
1498 
        Inf(1), // pow10(309)
1499 
}
1500 
 
1501 
var vfsignbitSC = []float64{
1502 
        Inf(-1),
1503 
        Copysign(0, -1),
1504 
        0,
1505 
        Inf(1),
1506 
        NaN(),
1507 
}
1508 
var signbitSC = []bool{
1509 
        true,
1510 
        true,
1511 
        false,
1512 
        false,
1513 
        false,
1514 
}
1515 
 
1516 
var vfsinSC = []float64{
1517 
        Inf(-1),
1518 
        Copysign(0, -1),
1519 
        0,
1520 
        Inf(1),
1521 
        NaN(),
1522 
}
1523 
var sinSC = []float64{
1524 
        NaN(),
1525 
        Copysign(0, -1),
1526 
        0,
1527 
        NaN(),
1528 
        NaN(),
1529 
}
1530 
 
1531 
var vfsinhSC = []float64{
1532 
        Inf(-1),
1533 
        Copysign(0, -1),
1534 
        0,
1535 
        Inf(1),
1536 
        NaN(),
1537 
}
1538 
var sinhSC = []float64{
1539 
        Inf(-1),
1540 
        Copysign(0, -1),
1541 
        0,
1542 
        Inf(1),
1543 
        NaN(),
1544 
}
1545 
 
1546 
var vfsqrtSC = []float64{
1547 
        Inf(-1),
1548 
        -Pi,
1549 
        Copysign(0, -1),
1550 
        0,
1551 
        Inf(1),
1552 
        NaN(),
1553 
}
1554 
var sqrtSC = []float64{
1555 
        NaN(),
1556 
        NaN(),
1557 
        Copysign(0, -1),
1558 
        0,
1559 
        Inf(1),
1560 
        NaN(),
1561 
}
1562 
 
1563 
var vftanhSC = []float64{
1564 
        Inf(-1),
1565 
        Copysign(0, -1),
1566 
        0,
1567 
        Inf(1),
1568 
        NaN(),
1569 
}
1570 
var tanhSC = []float64{
1571 
        -1,
1572 
        Copysign(0, -1),
1573 
        0,
1574 
        1,
1575 
        NaN(),
1576 
}
1577 
 
1578 
var vfy0SC = []float64{
1579 
        Inf(-1),
1580 
        0,
1581 
        Inf(1),
1582 
        NaN(),
1583 
}
1584 
var y0SC = []float64{
1585 
        NaN(),
1586 
        Inf(-1),
1587 
        0,
1588 
        NaN(),
1589 
}
1590 
var y1SC = []float64{
1591 
        NaN(),
1592 
        Inf(-1),
1593 
        0,
1594 
        NaN(),
1595 
}
1596 
var y2SC = []float64{
1597 
        NaN(),
1598 
        Inf(-1),
1599 
        0,
1600 
        NaN(),
1601 
}
1602 
var yM3SC = []float64{
1603 
        NaN(),
1604 
        Inf(1),
1605 
        0,
1606 
        NaN(),
1607 
}
1608 
 
1609 
// arguments and expected results for boundary cases
1610 
const (
1611 
        SmallestNormalFloat64   = 2.2250738585072014e-308 // 2**-1022
1612 
        LargestSubnormalFloat64 = SmallestNormalFloat64 - SmallestNonzeroFloat64
1613 
)
1614 
 
1615 
var vffrexpBC = []float64{
1616 
        SmallestNormalFloat64,
1617 
        LargestSubnormalFloat64,
1618 
        SmallestNonzeroFloat64,
1619 
        MaxFloat64,
1620 
        -SmallestNormalFloat64,
1621 
        -LargestSubnormalFloat64,
1622 
        -SmallestNonzeroFloat64,
1623 
        -MaxFloat64,
1624 
}
1625 
var frexpBC = []fi{
1626 
        {0.5, -1021},
1627 
        {0.99999999999999978, -1022},
1628 
        {0.5, -1073},
1629 
        {0.99999999999999989, 1024},
1630 
        {-0.5, -1021},
1631 
        {-0.99999999999999978, -1022},
1632 
        {-0.5, -1073},
1633 
        {-0.99999999999999989, 1024},
1634 
}
1635 
 
1636 
var vfldexpBC = []fi{
1637 
        {SmallestNormalFloat64, -52},
1638 
        {LargestSubnormalFloat64, -51},
1639 
        {SmallestNonzeroFloat64, 1074},
1640 
        {MaxFloat64, -(1023 + 1074)},
1641 
        {1, -1075},
1642 
        {-1, -1075},
1643 
        {1, 1024},
1644 
        {-1, 1024},
1645 
}
1646 
var ldexpBC = []float64{
1647 
        SmallestNonzeroFloat64,
1648 
        1e-323, // 2**-1073
1649 
        1,
1650 
        1e-323, // 2**-1073
1651 
        0,
1652 
        Copysign(0, -1),
1653 
        Inf(1),
1654 
        Inf(-1),
1655 
}
1656 
 
1657 
var logbBC = []float64{
1658 
        -1022,
1659 
        -1023,
1660 
        -1074,
1661 
        1023,
1662 
        -1022,
1663 
        -1023,
1664 
        -1074,
1665 
        1023,
1666 
}
1667 
 
1668 
func tolerance(a, b, e float64) bool {
1669 
        d := a - b
1670 
        if d < 0 {
1671 
                d = -d
1672 
        }
1673 
 
1674 
        if a != 0 {
1675 
                e = e * a
1676 
                if e < 0 {
1677 
                        e = -e
1678 
                }
1679 
        }
1680 
        return d < e
1681 
}
1682 
func kindaclose(a, b float64) bool { return tolerance(a, b, 1e-8) }
1683 
func close(a, b float64) bool      { return tolerance(a, b, 1e-14) }
1684 
func veryclose(a, b float64) bool  { return tolerance(a, b, 4e-16) }
1685 
func soclose(a, b, e float64) bool { return tolerance(a, b, e) }
1686 
func alike(a, b float64) bool {
1687 
        switch {
1688 
        case IsNaN(a) && IsNaN(b):
1689 
                return true
1690 
        case a == b:
1691 
                return Signbit(a) == Signbit(b)
1692 
        }
1693 
        return false
1694 
}
1695 
 
1696 
func TestAcos(t *testing.T) {
1697 
        for i := 0; i < len(vf); i++ {
1698 
                a := vf[i] / 10
1699 
                if f := Acos(a); !close(acos[i], f) {
1700 
                        t.Errorf("Acos(%g) = %g, want %g", a, f, acos[i])
1701 
                }
1702 
        }
1703 
        for i := 0; i < len(vfacosSC); i++ {
1704 
                if f := Acos(vfacosSC[i]); !alike(acosSC[i], f) {
1705 
                        t.Errorf("Acos(%g) = %g, want %g", vfacosSC[i], f, acosSC[i])
1706 
                }
1707 
        }
1708 
}
1709 
 
1710 
func TestAcosh(t *testing.T) {
1711 
        for i := 0; i < len(vf); i++ {
1712 
                a := 1 + Abs(vf[i])
1713 
                if f := Acosh(a); !veryclose(acosh[i], f) {
1714 
                        t.Errorf("Acosh(%g) = %g, want %g", a, f, acosh[i])
1715 
                }
1716 
        }
1717 
        for i := 0; i < len(vfacoshSC); i++ {
1718 
                if f := Acosh(vfacoshSC[i]); !alike(acoshSC[i], f) {
1719 
                        t.Errorf("Acosh(%g) = %g, want %g", vfacoshSC[i], f, acoshSC[i])
1720 
                }
1721 
        }
1722 
}
1723 
 
1724 
func TestAsin(t *testing.T) {
1725 
        for i := 0; i < len(vf); i++ {
1726 
                a := vf[i] / 10
1727 
                if f := Asin(a); !veryclose(asin[i], f) {
1728 
                        t.Errorf("Asin(%g) = %g, want %g", a, f, asin[i])
1729 
                }
1730 
        }
1731 
        for i := 0; i < len(vfasinSC); i++ {
1732 
                if f := Asin(vfasinSC[i]); !alike(asinSC[i], f) {
1733 
                        t.Errorf("Asin(%g) = %g, want %g", vfasinSC[i], f, asinSC[i])
1734 
                }
1735 
        }
1736 
}
1737 
 
1738 
func TestAsinh(t *testing.T) {
1739 
        for i := 0; i < len(vf); i++ {
1740 
                if f := Asinh(vf[i]); !veryclose(asinh[i], f) {
1741 
                        t.Errorf("Asinh(%g) = %g, want %g", vf[i], f, asinh[i])
1742 
                }
1743 
        }
1744 
        for i := 0; i < len(vfasinhSC); i++ {
1745 
                if f := Asinh(vfasinhSC[i]); !alike(asinhSC[i], f) {
1746 
                        t.Errorf("Asinh(%g) = %g, want %g", vfasinhSC[i], f, asinhSC[i])
1747 
                }
1748 
        }
1749 
}
1750 
 
1751 
func TestAtan(t *testing.T) {
1752 
        for i := 0; i < len(vf); i++ {
1753 
                if f := Atan(vf[i]); !veryclose(atan[i], f) {
1754 
                        t.Errorf("Atan(%g) = %g, want %g", vf[i], f, atan[i])
1755 
                }
1756 
        }
1757 
        for i := 0; i < len(vfatanSC); i++ {
1758 
                if f := Atan(vfatanSC[i]); !alike(atanSC[i], f) {
1759 
                        t.Errorf("Atan(%g) = %g, want %g", vfatanSC[i], f, atanSC[i])
1760 
                }
1761 
        }
1762 
}
1763 
 
1764 
func TestAtanh(t *testing.T) {
1765 
        for i := 0; i < len(vf); i++ {
1766 
                a := vf[i] / 10
1767 
                if f := Atanh(a); !veryclose(atanh[i], f) {
1768 
                        t.Errorf("Atanh(%g) = %g, want %g", a, f, atanh[i])
1769 
                }
1770 
        }
1771 
        for i := 0; i < len(vfatanhSC); i++ {
1772 
                if f := Atanh(vfatanhSC[i]); !alike(atanhSC[i], f) {
1773 
                        t.Errorf("Atanh(%g) = %g, want %g", vfatanhSC[i], f, atanhSC[i])
1774 
                }
1775 
        }
1776 
}
1777 
 
1778 
func TestAtan2(t *testing.T) {
1779 
        for i := 0; i < len(vf); i++ {
1780 
                if f := Atan2(10, vf[i]); !veryclose(atan2[i], f) {
1781 
                        t.Errorf("Atan2(10, %g) = %g, want %g", vf[i], f, atan2[i])
1782 
                }
1783 
        }
1784 
        for i := 0; i < len(vfatan2SC); i++ {
1785 
                if f := Atan2(vfatan2SC[i][0], vfatan2SC[i][1]); !alike(atan2SC[i], f) {
1786 
                        t.Errorf("Atan2(%g, %g) = %g, want %g", vfatan2SC[i][0], vfatan2SC[i][1], f, atan2SC[i])
1787 
                }
1788 
        }
1789 
}
1790 
 
1791 
func TestCbrt(t *testing.T) {
1792 
        for i := 0; i < len(vf); i++ {
1793 
                if f := Cbrt(vf[i]); !veryclose(cbrt[i], f) {
1794 
                        t.Errorf("Cbrt(%g) = %g, want %g", vf[i], f, cbrt[i])
1795 
                }
1796 
        }
1797 
        for i := 0; i < len(vfcbrtSC); i++ {
1798 
                if f := Cbrt(vfcbrtSC[i]); !alike(cbrtSC[i], f) {
1799 
                        t.Errorf("Cbrt(%g) = %g, want %g", vfcbrtSC[i], f, cbrtSC[i])
1800 
                }
1801 
        }
1802 
}
1803 
 
1804 
func TestCeil(t *testing.T) {
1805 
        for i := 0; i < len(vf); i++ {
1806 
                if f := Ceil(vf[i]); ceil[i] != f {
1807 
                        t.Errorf("Ceil(%g) = %g, want %g", vf[i], f, ceil[i])
1808 
                }
1809 
        }
1810 
        for i := 0; i < len(vfceilSC); i++ {
1811 
                if f := Ceil(vfceilSC[i]); !alike(ceilSC[i], f) {
1812 
                        t.Errorf("Ceil(%g) = %g, want %g", vfceilSC[i], f, ceilSC[i])
1813 
                }
1814 
        }
1815 
}
1816 
 
1817 
func TestCopysign(t *testing.T) {
1818 
        for i := 0; i < len(vf); i++ {
1819 
                if f := Copysign(vf[i], -1); copysign[i] != f {
1820 
                        t.Errorf("Copysign(%g, -1) = %g, want %g", vf[i], f, copysign[i])
1821 
                }
1822 
        }
1823 
        for i := 0; i < len(vf); i++ {
1824 
                if f := Copysign(vf[i], 1); -copysign[i] != f {
1825 
                        t.Errorf("Copysign(%g, 1) = %g, want %g", vf[i], f, -copysign[i])
1826 
                }
1827 
        }
1828 
        for i := 0; i < len(vfcopysignSC); i++ {
1829 
                if f := Copysign(vfcopysignSC[i], -1); !alike(copysignSC[i], f) {
1830 
                        t.Errorf("Copysign(%g, -1) = %g, want %g", vfcopysignSC[i], f, copysignSC[i])
1831 
                }
1832 
        }
1833 
}
1834 
 
1835 
func TestCos(t *testing.T) {
1836 
        for i := 0; i < len(vf); i++ {
1837 
                if f := Cos(vf[i]); !veryclose(cos[i], f) {
1838 
                        t.Errorf("Cos(%g) = %g, want %g", vf[i], f, cos[i])
1839 
                }
1840 
        }
1841 
        for i := 0; i < len(vfcosSC); i++ {
1842 
                if f := Cos(vfcosSC[i]); !alike(cosSC[i], f) {
1843 
                        t.Errorf("Cos(%g) = %g, want %g", vfcosSC[i], f, cosSC[i])
1844 
                }
1845 
        }
1846 
}
1847 
 
1848 
func TestCosh(t *testing.T) {
1849 
        for i := 0; i < len(vf); i++ {
1850 
                if f := Cosh(vf[i]); !close(cosh[i], f) {
1851 
                        t.Errorf("Cosh(%g) = %g, want %g", vf[i], f, cosh[i])
1852 
                }
1853 
        }
1854 
        for i := 0; i < len(vfcoshSC); i++ {
1855 
                if f := Cosh(vfcoshSC[i]); !alike(coshSC[i], f) {
1856 
                        t.Errorf("Cosh(%g) = %g, want %g", vfcoshSC[i], f, coshSC[i])
1857 
                }
1858 
        }
1859 
}
1860 
 
1861 
func TestErf(t *testing.T) {
1862 
        for i := 0; i < len(vf); i++ {
1863 
                a := vf[i] / 10
1864 
                if f := Erf(a); !veryclose(erf[i], f) {
1865 
                        t.Errorf("Erf(%g) = %g, want %g", a, f, erf[i])
1866 
                }
1867 
        }
1868 
        for i := 0; i < len(vferfSC); i++ {
1869 
                if f := Erf(vferfSC[i]); !alike(erfSC[i], f) {
1870 
                        t.Errorf("Erf(%g) = %g, want %g", vferfSC[i], f, erfSC[i])
1871 
                }
1872 
        }
1873 
}
1874 
 
1875 
func TestErfc(t *testing.T) {
1876 
        for i := 0; i < len(vf); i++ {
1877 
                a := vf[i] / 10
1878 
                if f := Erfc(a); !veryclose(erfc[i], f) {
1879 
                        t.Errorf("Erfc(%g) = %g, want %g", a, f, erfc[i])
1880 
                }
1881 
        }
1882 
        for i := 0; i < len(vferfcSC); i++ {
1883 
                if f := Erfc(vferfcSC[i]); !alike(erfcSC[i], f) {
1884 
                        t.Errorf("Erfc(%g) = %g, want %g", vferfcSC[i], f, erfcSC[i])
1885 
                }
1886 
        }
1887 
}
1888 
 
1889 
func TestExp(t *testing.T) {
1890 
        testExp(t, Exp, "Exp")
1891 
        testExp(t, ExpGo, "ExpGo")
1892 
}
1893 
 
1894 
func testExp(t *testing.T, Exp func(float64) float64, name string) {
1895 
        for i := 0; i < len(vf); i++ {
1896 
                if f := Exp(vf[i]); !close(exp[i], f) {
1897 
                        t.Errorf("%s(%g) = %g, want %g", name, vf[i], f, exp[i])
1898 
                }
1899 
        }
1900 
        for i := 0; i < len(vfexpSC); i++ {
1901 
                if f := Exp(vfexpSC[i]); !alike(expSC[i], f) {
1902 
                        t.Errorf("%s(%g) = %g, want %g", name, vfexpSC[i], f, expSC[i])
1903 
                }
1904 
        }
1905 
}
1906 
 
1907 
func TestExpm1(t *testing.T) {
1908 
        for i := 0; i < len(vf); i++ {
1909 
                a := vf[i] / 100
1910 
                if f := Expm1(a); !veryclose(expm1[i], f) {
1911 
                        t.Errorf("Expm1(%g) = %g, want %g", a, f, expm1[i])
1912 
                }
1913 
        }
1914 
        for i := 0; i < len(vfexpm1SC); i++ {
1915 
                if f := Expm1(vfexpm1SC[i]); !alike(expm1SC[i], f) {
1916 
                        t.Errorf("Expm1(%g) = %g, want %g", vfexpm1SC[i], f, expm1SC[i])
1917 
                }
1918 
        }
1919 
}
1920 
 
1921 
func TestExp2(t *testing.T) {
1922 
        testExp2(t, Exp2, "Exp2")
1923 
        testExp2(t, Exp2Go, "Exp2Go")
1924 
}
1925 
 
1926 
func testExp2(t *testing.T, Exp2 func(float64) float64, name string) {
1927 
        for i := 0; i < len(vf); i++ {
1928 
                if f := Exp2(vf[i]); !close(exp2[i], f) {
1929 
                        t.Errorf("%s(%g) = %g, want %g", name, vf[i], f, exp2[i])
1930 
                }
1931 
        }
1932 
        for i := 0; i < len(vfexpSC); i++ {
1933 
                if f := Exp2(vfexpSC[i]); !alike(expSC[i], f) {
1934 
                        t.Errorf("%s(%g) = %g, want %g", name, vfexpSC[i], f, expSC[i])
1935 
                }
1936 
        }
1937 
        for n := -1074; n < 1024; n++ {
1938 
                f := Exp2(float64(n))
1939 
                vf := Ldexp(1, n)
1940 
                if f != vf {
1941 
                        t.Errorf("%s(%d) = %g, want %g", name, n, f, vf)
1942 
                }
1943 
        }
1944 
}
1945 
 
1946 
func TestAbs(t *testing.T) {
1947 
        for i := 0; i < len(vf); i++ {
1948 
                if f := Abs(vf[i]); fabs[i] != f {
1949 
                        t.Errorf("Abs(%g) = %g, want %g", vf[i], f, fabs[i])
1950 
                }
1951 
        }
1952 
        for i := 0; i < len(vffabsSC); i++ {
1953 
                if f := Abs(vffabsSC[i]); !alike(fabsSC[i], f) {
1954 
                        t.Errorf("Abs(%g) = %g, want %g", vffabsSC[i], f, fabsSC[i])
1955 
                }
1956 
        }
1957 
}
1958 
 
1959 
func TestDim(t *testing.T) {
1960 
        for i := 0; i < len(vf); i++ {
1961 
                if f := Dim(vf[i], 0); fdim[i] != f {
1962 
                        t.Errorf("Dim(%g, %g) = %g, want %g", vf[i], 0.0, f, fdim[i])
1963 
                }
1964 
        }
1965 
        for i := 0; i < len(vffdimSC); i++ {
1966 
                if f := Dim(vffdimSC[i][0], vffdimSC[i][1]); !alike(fdimSC[i], f) {
1967 
                        t.Errorf("Dim(%g, %g) = %g, want %g", vffdimSC[i][0], vffdimSC[i][1], f, fdimSC[i])
1968 
                }
1969 
        }
1970 
}
1971 
 
1972 
func TestFloor(t *testing.T) {
1973 
        for i := 0; i < len(vf); i++ {
1974 
                if f := Floor(vf[i]); floor[i] != f {
1975 
                        t.Errorf("Floor(%g) = %g, want %g", vf[i], f, floor[i])
1976 
                }
1977 
        }
1978 
        for i := 0; i < len(vfceilSC); i++ {
1979 
                if f := Floor(vfceilSC[i]); !alike(ceilSC[i], f) {
1980 
                        t.Errorf("Floor(%g) = %g, want %g", vfceilSC[i], f, ceilSC[i])
1981 
                }
1982 
        }
1983 
}
1984 
 
1985 
func TestMax(t *testing.T) {
1986 
        for i := 0; i < len(vf); i++ {
1987 
                if f := Max(vf[i], ceil[i]); ceil[i] != f {
1988 
                        t.Errorf("Max(%g, %g) = %g, want %g", vf[i], ceil[i], f, ceil[i])
1989 
                }
1990 
        }
1991 
        for i := 0; i < len(vffdimSC); i++ {
1992 
                if f := Max(vffdimSC[i][0], vffdimSC[i][1]); !alike(fmaxSC[i], f) {
1993 
                        t.Errorf("Max(%g, %g) = %g, want %g", vffdimSC[i][0], vffdimSC[i][1], f, fmaxSC[i])
1994 
                }
1995 
        }
1996 
}
1997 
 
1998 
func TestMin(t *testing.T) {
1999 
        for i := 0; i < len(vf); i++ {
2000 
                if f := Min(vf[i], floor[i]); floor[i] != f {
2001 
                        t.Errorf("Min(%g, %g) = %g, want %g", vf[i], floor[i], f, floor[i])
2002 
                }
2003 
        }
2004 
        for i := 0; i < len(vffdimSC); i++ {
2005 
                if f := Min(vffdimSC[i][0], vffdimSC[i][1]); !alike(fminSC[i], f) {
2006 
                        t.Errorf("Min(%g, %g) = %g, want %g", vffdimSC[i][0], vffdimSC[i][1], f, fminSC[i])
2007 
                }
2008 
        }
2009 
}
2010 
 
2011 
func TestMod(t *testing.T) {
2012 
        for i := 0; i < len(vf); i++ {
2013 
                if f := Mod(10, vf[i]); fmod[i] != f {
2014 
                        t.Errorf("Mod(10, %g) = %g, want %g", vf[i], f, fmod[i])
2015 
                }
2016 
        }
2017 
        for i := 0; i < len(vffmodSC); i++ {
2018 
                if f := Mod(vffmodSC[i][0], vffmodSC[i][1]); !alike(fmodSC[i], f) {
2019 
                        t.Errorf("Mod(%g, %g) = %g, want %g", vffmodSC[i][0], vffmodSC[i][1], f, fmodSC[i])
2020 
                }
2021 
        }
2022 
}
2023 
 
2024 
func TestFrexp(t *testing.T) {
2025 
        for i := 0; i < len(vf); i++ {
2026 
                if f, j := Frexp(vf[i]); !veryclose(frexp[i].f, f) || frexp[i].i != j {
2027 
                        t.Errorf("Frexp(%g) = %g, %d, want %g, %d", vf[i], f, j, frexp[i].f, frexp[i].i)
2028 
                }
2029 
        }
2030 
        for i := 0; i < len(vffrexpSC); i++ {
2031 
                if f, j := Frexp(vffrexpSC[i]); !alike(frexpSC[i].f, f) || frexpSC[i].i != j {
2032 
                        t.Errorf("Frexp(%g) = %g, %d, want %g, %d", vffrexpSC[i], f, j, frexpSC[i].f, frexpSC[i].i)
2033 
                }
2034 
        }
2035 
        for i := 0; i < len(vffrexpBC); i++ {
2036 
                if f, j := Frexp(vffrexpBC[i]); !alike(frexpBC[i].f, f) || frexpBC[i].i != j {
2037 
                        t.Errorf("Frexp(%g) = %g, %d, want %g, %d", vffrexpBC[i], f, j, frexpBC[i].f, frexpBC[i].i)
2038 
                }
2039 
        }
2040 
}
2041 
 
2042 
func TestGamma(t *testing.T) {
2043 
        for i := 0; i < len(vf); i++ {
2044 
                if f := Gamma(vf[i]); !close(gamma[i], f) {
2045 
                        t.Errorf("Gamma(%g) = %g, want %g", vf[i], f, gamma[i])
2046 
                }
2047 
        }
2048 
        for i := 0; i < len(vfgammaSC); i++ {
2049 
                if f := Gamma(vfgammaSC[i]); !alike(gammaSC[i], f) {
2050 
                        t.Errorf("Gamma(%g) = %g, want %g", vfgammaSC[i], f, gammaSC[i])
2051 
                }
2052 
        }
2053 
}
2054 
 
2055 
func TestHypot(t *testing.T) {
2056 
        for i := 0; i < len(vf); i++ {
2057 
                a := Abs(1e200 * tanh[i] * Sqrt(2))
2058 
                if f := Hypot(1e200*tanh[i], 1e200*tanh[i]); !veryclose(a, f) {
2059 
                        t.Errorf("Hypot(%g, %g) = %g, want %g", 1e200*tanh[i], 1e200*tanh[i], f, a)
2060 
                }
2061 
        }
2062 
        for i := 0; i < len(vfhypotSC); i++ {
2063 
                if f := Hypot(vfhypotSC[i][0], vfhypotSC[i][1]); !alike(hypotSC[i], f) {
2064 
                        t.Errorf("Hypot(%g, %g) = %g, want %g", vfhypotSC[i][0], vfhypotSC[i][1], f, hypotSC[i])
2065 
                }
2066 
        }
2067 
}
2068 
 
2069 
func TestHypotGo(t *testing.T) {
2070 
        for i := 0; i < len(vf); i++ {
2071 
                a := Abs(1e200 * tanh[i] * Sqrt(2))
2072 
                if f := HypotGo(1e200*tanh[i], 1e200*tanh[i]); !veryclose(a, f) {
2073 
                        t.Errorf("HypotGo(%g, %g) = %g, want %g", 1e200*tanh[i], 1e200*tanh[i], f, a)
2074 
                }
2075 
        }
2076 
        for i := 0; i < len(vfhypotSC); i++ {
2077 
                if f := HypotGo(vfhypotSC[i][0], vfhypotSC[i][1]); !alike(hypotSC[i], f) {
2078 
                        t.Errorf("HypotGo(%g, %g) = %g, want %g", vfhypotSC[i][0], vfhypotSC[i][1], f, hypotSC[i])
2079 
                }
2080 
        }
2081 
}
2082 
 
2083 
func TestIlogb(t *testing.T) {
2084 
        for i := 0; i < len(vf); i++ {
2085 
                a := frexp[i].i - 1 // adjust because fr in the interval [½, 1)
2086 
                if e := Ilogb(vf[i]); a != e {
2087 
                        t.Errorf("Ilogb(%g) = %d, want %d", vf[i], e, a)
2088 
                }
2089 
        }
2090 
        for i := 0; i < len(vflogbSC); i++ {
2091 
                if e := Ilogb(vflogbSC[i]); ilogbSC[i] != e {
2092 
                        t.Errorf("Ilogb(%g) = %d, want %d", vflogbSC[i], e, ilogbSC[i])
2093 
                }
2094 
        }
2095 
        for i := 0; i < len(vffrexpBC); i++ {
2096 
                if e := Ilogb(vffrexpBC[i]); int(logbBC[i]) != e {
2097 
                        t.Errorf("Ilogb(%g) = %d, want %d", vffrexpBC[i], e, int(logbBC[i]))
2098 
                }
2099 
        }
2100 
}
2101 
 
2102 
func TestJ0(t *testing.T) {
2103 
        for i := 0; i < len(vf); i++ {
2104 
                if f := J0(vf[i]); !soclose(j0[i], f, 4e-14) {
2105 
                        t.Errorf("J0(%g) = %g, want %g", vf[i], f, j0[i])
2106 
                }
2107 
        }
2108 
        for i := 0; i < len(vfj0SC); i++ {
2109 
                if f := J0(vfj0SC[i]); !alike(j0SC[i], f) {
2110 
                        t.Errorf("J0(%g) = %g, want %g", vfj0SC[i], f, j0SC[i])
2111 
                }
2112 
        }
2113 
}
2114 
 
2115 
func TestJ1(t *testing.T) {
2116 
        for i := 0; i < len(vf); i++ {
2117 
                if f := J1(vf[i]); !close(j1[i], f) {
2118 
                        t.Errorf("J1(%g) = %g, want %g", vf[i], f, j1[i])
2119 
                }
2120 
        }
2121 
        for i := 0; i < len(vfj0SC); i++ {
2122 
                if f := J1(vfj0SC[i]); !alike(j1SC[i], f) {
2123 
                        t.Errorf("J1(%g) = %g, want %g", vfj0SC[i], f, j1SC[i])
2124 
                }
2125 
        }
2126 
}
2127 
 
2128 
func TestJn(t *testing.T) {
2129 
        for i := 0; i < len(vf); i++ {
2130 
                if f := Jn(2, vf[i]); !close(j2[i], f) {
2131 
                        t.Errorf("Jn(2, %g) = %g, want %g", vf[i], f, j2[i])
2132 
                }
2133 
                if f := Jn(-3, vf[i]); !close(jM3[i], f) {
2134 
                        t.Errorf("Jn(-3, %g) = %g, want %g", vf[i], f, jM3[i])
2135 
                }
2136 
        }
2137 
        for i := 0; i < len(vfj0SC); i++ {
2138 
                if f := Jn(2, vfj0SC[i]); !alike(j2SC[i], f) {
2139 
                        t.Errorf("Jn(2, %g) = %g, want %g", vfj0SC[i], f, j2SC[i])
2140 
                }
2141 
                if f := Jn(-3, vfj0SC[i]); !alike(jM3SC[i], f) {
2142 
                        t.Errorf("Jn(-3, %g) = %g, want %g", vfj0SC[i], f, jM3SC[i])
2143 
                }
2144 
        }
2145 
}
2146 
 
2147 
func TestLdexp(t *testing.T) {
2148 
        for i := 0; i < len(vf); i++ {
2149 
                if f := Ldexp(frexp[i].f, frexp[i].i); !veryclose(vf[i], f) {
2150 
                        t.Errorf("Ldexp(%g, %d) = %g, want %g", frexp[i].f, frexp[i].i, f, vf[i])
2151 
                }
2152 
        }
2153 
        for i := 0; i < len(vffrexpSC); i++ {
2154 
                if f := Ldexp(frexpSC[i].f, frexpSC[i].i); !alike(vffrexpSC[i], f) {
2155 
                        t.Errorf("Ldexp(%g, %d) = %g, want %g", frexpSC[i].f, frexpSC[i].i, f, vffrexpSC[i])
2156 
                }
2157 
        }
2158 
        for i := 0; i < len(vfldexpSC); i++ {
2159 
                if f := Ldexp(vfldexpSC[i].f, vfldexpSC[i].i); !alike(ldexpSC[i], f) {
2160 
                        t.Errorf("Ldexp(%g, %d) = %g, want %g", vfldexpSC[i].f, vfldexpSC[i].i, f, ldexpSC[i])
2161 
                }
2162 
        }
2163 
        for i := 0; i < len(vffrexpBC); i++ {
2164 
                if f := Ldexp(frexpBC[i].f, frexpBC[i].i); !alike(vffrexpBC[i], f) {
2165 
                        t.Errorf("Ldexp(%g, %d) = %g, want %g", frexpBC[i].f, frexpBC[i].i, f, vffrexpBC[i])
2166 
                }
2167 
        }
2168 
        for i := 0; i < len(vfldexpBC); i++ {
2169 
                if f := Ldexp(vfldexpBC[i].f, vfldexpBC[i].i); !alike(ldexpBC[i], f) {
2170 
                        t.Errorf("Ldexp(%g, %d) = %g, want %g", vfldexpBC[i].f, vfldexpBC[i].i, f, ldexpBC[i])
2171 
                }
2172 
        }
2173 
}
2174 
 
2175 
func TestLgamma(t *testing.T) {
2176 
        for i := 0; i < len(vf); i++ {
2177 
                if f, s := Lgamma(vf[i]); !close(lgamma[i].f, f) || lgamma[i].i != s {
2178 
                        t.Errorf("Lgamma(%g) = %g, %d, want %g, %d", vf[i], f, s, lgamma[i].f, lgamma[i].i)
2179 
                }
2180 
        }
2181 
        for i := 0; i < len(vflgammaSC); i++ {
2182 
                if f, s := Lgamma(vflgammaSC[i]); !alike(lgammaSC[i].f, f) || lgammaSC[i].i != s {
2183 
                        t.Errorf("Lgamma(%g) = %g, %d, want %g, %d", vflgammaSC[i], f, s, lgammaSC[i].f, lgammaSC[i].i)
2184 
                }
2185 
        }
2186 
}
2187 
 
2188 
func TestLog(t *testing.T) {
2189 
        for i := 0; i < len(vf); i++ {
2190 
                a := Abs(vf[i])
2191 
                if f := Log(a); log[i] != f {
2192 
                        t.Errorf("Log(%g) = %g, want %g", a, f, log[i])
2193 
                }
2194 
        }
2195 
        if f := Log(10); f != Ln10 {
2196 
                t.Errorf("Log(%g) = %g, want %g", 10.0, f, Ln10)
2197 
        }
2198 
        for i := 0; i < len(vflogSC); i++ {
2199 
                if f := Log(vflogSC[i]); !alike(logSC[i], f) {
2200 
                        t.Errorf("Log(%g) = %g, want %g", vflogSC[i], f, logSC[i])
2201 
                }
2202 
        }
2203 
}
2204 
 
2205 
func TestLogb(t *testing.T) {
2206 
        for i := 0; i < len(vf); i++ {
2207 
                if f := Logb(vf[i]); logb[i] != f {
2208 
                        t.Errorf("Logb(%g) = %g, want %g", vf[i], f, logb[i])
2209 
                }
2210 
        }
2211 
        for i := 0; i < len(vflogbSC); i++ {
2212 
                if f := Logb(vflogbSC[i]); !alike(logbSC[i], f) {
2213 
                        t.Errorf("Logb(%g) = %g, want %g", vflogbSC[i], f, logbSC[i])
2214 
                }
2215 
        }
2216 
        for i := 0; i < len(vffrexpBC); i++ {
2217 
                if f := Logb(vffrexpBC[i]); !alike(logbBC[i], f) {
2218 
                        t.Errorf("Logb(%g) = %g, want %g", vffrexpBC[i], f, logbBC[i])
2219 
                }
2220 
        }
2221 
}
2222 
 
2223 
func TestLog10(t *testing.T) {
2224 
        for i := 0; i < len(vf); i++ {
2225 
                a := Abs(vf[i])
2226 
                if f := Log10(a); !veryclose(log10[i], f) {
2227 
                        t.Errorf("Log10(%g) = %g, want %g", a, f, log10[i])
2228 
                }
2229 
        }
2230 
        if f := Log10(E); f != Log10E {
2231 
                t.Errorf("Log10(%g) = %g, want %g", E, f, Log10E)
2232 
        }
2233 
        for i := 0; i < len(vflogSC); i++ {
2234 
                if f := Log10(vflogSC[i]); !alike(logSC[i], f) {
2235 
                        t.Errorf("Log10(%g) = %g, want %g", vflogSC[i], f, logSC[i])
2236 
                }
2237 
        }
2238 
}
2239 
 
2240 
func TestLog1p(t *testing.T) {
2241 
        for i := 0; i < len(vf); i++ {
2242 
                a := vf[i] / 100
2243 
                if f := Log1p(a); !veryclose(log1p[i], f) {
2244 
                        t.Errorf("Log1p(%g) = %g, want %g", a, f, log1p[i])
2245 
                }
2246 
        }
2247 
        a := 9.0
2248 
        if f := Log1p(a); f != Ln10 {
2249 
                t.Errorf("Log1p(%g) = %g, want %g", a, f, Ln10)
2250 
        }
2251 
        for i := 0; i < len(vflogSC); i++ {
2252 
                if f := Log1p(vflog1pSC[i]); !alike(log1pSC[i], f) {
2253 
                        t.Errorf("Log1p(%g) = %g, want %g", vflog1pSC[i], f, log1pSC[i])
2254 
                }
2255 
        }
2256 
}
2257 
 
2258 
func TestLog2(t *testing.T) {
2259 
        for i := 0; i < len(vf); i++ {
2260 
                a := Abs(vf[i])
2261 
                if f := Log2(a); !veryclose(log2[i], f) {
2262 
                        t.Errorf("Log2(%g) = %g, want %g", a, f, log2[i])
2263 
                }
2264 
        }
2265 
        if f := Log2(E); f != Log2E {
2266 
                t.Errorf("Log2(%g) = %g, want %g", E, f, Log2E)
2267 
        }
2268 
        for i := 0; i < len(vflogSC); i++ {
2269 
                if f := Log2(vflogSC[i]); !alike(logSC[i], f) {
2270 
                        t.Errorf("Log2(%g) = %g, want %g", vflogSC[i], f, logSC[i])
2271 
                }
2272 
        }
2273 
}
2274 
 
2275 
func TestModf(t *testing.T) {
2276 
        for i := 0; i < len(vf); i++ {
2277 
                if f, g := Modf(vf[i]); !veryclose(modf[i][0], f) || !veryclose(modf[i][1], g) {
2278 
                        t.Errorf("Modf(%g) = %g, %g, want %g, %g", vf[i], f, g, modf[i][0], modf[i][1])
2279 
                }
2280 
        }
2281 
        for i := 0; i < len(vfmodfSC); i++ {
2282 
                if f, g := Modf(vfmodfSC[i]); !alike(modfSC[i][0], f) || !alike(modfSC[i][1], g) {
2283 
                        t.Errorf("Modf(%g) = %g, %g, want %g, %g", vfmodfSC[i], f, g, modfSC[i][0], modfSC[i][1])
2284 
                }
2285 
        }
2286 
}
2287 
 
2288 
func TestNextafter(t *testing.T) {
2289 
        for i := 0; i < len(vf); i++ {
2290 
                if f := Nextafter(vf[i], 10); nextafter[i] != f {
2291 
                        t.Errorf("Nextafter(%g, %g) = %g want %g", vf[i], 10.0, f, nextafter[i])
2292 
                }
2293 
        }
2294 
        for i := 0; i < len(vfnextafterSC); i++ {
2295 
                if f := Nextafter(vfnextafterSC[i][0], vfnextafterSC[i][1]); !alike(nextafterSC[i], f) {
2296 
                        t.Errorf("Nextafter(%g, %g) = %g want %g", vfnextafterSC[i][0], vfnextafterSC[i][1], f, nextafterSC[i])
2297 
                }
2298 
        }
2299 
}
2300 
 
2301 
func TestPow(t *testing.T) {
2302 
        for i := 0; i < len(vf); i++ {
2303 
                if f := Pow(10, vf[i]); !close(pow[i], f) {
2304 
                        t.Errorf("Pow(10, %g) = %g, want %g", vf[i], f, pow[i])
2305 
                }
2306 
        }
2307 
        for i := 0; i < len(vfpowSC); i++ {
2308 
                if f := Pow(vfpowSC[i][0], vfpowSC[i][1]); !alike(powSC[i], f) {
2309 
                        t.Errorf("Pow(%g, %g) = %g, want %g", vfpowSC[i][0], vfpowSC[i][1], f, powSC[i])
2310 
                }
2311 
        }
2312 
}
2313 
 
2314 
func TestPow10(t *testing.T) {
2315 
        for i := 0; i < len(vfpow10SC); i++ {
2316 
                if f := Pow10(vfpow10SC[i]); !alike(pow10SC[i], f) {
2317 
                        t.Errorf("Pow10(%d) = %g, want %g", vfpow10SC[i], f, pow10SC[i])
2318 
                }
2319 
        }
2320 
}
2321 
 
2322 
func TestRemainder(t *testing.T) {
2323 
        for i := 0; i < len(vf); i++ {
2324 
                if f := Remainder(10, vf[i]); remainder[i] != f {
2325 
                        t.Errorf("Remainder(10, %g) = %g, want %g", vf[i], f, remainder[i])
2326 
                }
2327 
        }
2328 
        for i := 0; i < len(vffmodSC); i++ {
2329 
                if f := Remainder(vffmodSC[i][0], vffmodSC[i][1]); !alike(fmodSC[i], f) {
2330 
                        t.Errorf("Remainder(%g, %g) = %g, want %g", vffmodSC[i][0], vffmodSC[i][1], f, fmodSC[i])
2331 
                }
2332 
        }
2333 
}
2334 
 
2335 
func TestSignbit(t *testing.T) {
2336 
        for i := 0; i < len(vf); i++ {
2337 
                if f := Signbit(vf[i]); signbit[i] != f {
2338 
                        t.Errorf("Signbit(%g) = %t, want %t", vf[i], f, signbit[i])
2339 
                }
2340 
        }
2341 
        for i := 0; i < len(vfsignbitSC); i++ {
2342 
                if f := Signbit(vfsignbitSC[i]); signbitSC[i] != f {
2343 
                        t.Errorf("Signbit(%g) = %t, want %t", vfsignbitSC[i], f, signbitSC[i])
2344 
                }
2345 
        }
2346 
}
2347 
func TestSin(t *testing.T) {
2348 
        for i := 0; i < len(vf); i++ {
2349 
                if f := Sin(vf[i]); !veryclose(sin[i], f) {
2350 
                        t.Errorf("Sin(%g) = %g, want %g", vf[i], f, sin[i])
2351 
                }
2352 
        }
2353 
        for i := 0; i < len(vfsinSC); i++ {
2354 
                if f := Sin(vfsinSC[i]); !alike(sinSC[i], f) {
2355 
                        t.Errorf("Sin(%g) = %g, want %g", vfsinSC[i], f, sinSC[i])
2356 
                }
2357 
        }
2358 
}
2359 
 
2360 
func TestSincos(t *testing.T) {
2361 
        for i := 0; i < len(vf); i++ {
2362 
                if s, c := Sincos(vf[i]); !veryclose(sin[i], s) || !veryclose(cos[i], c) {
2363 
                        t.Errorf("Sincos(%g) = %g, %g want %g, %g", vf[i], s, c, sin[i], cos[i])
2364 
                }
2365 
        }
2366 
}
2367 
 
2368 
func TestSinh(t *testing.T) {
2369 
        for i := 0; i < len(vf); i++ {
2370 
                if f := Sinh(vf[i]); !close(sinh[i], f) {
2371 
                        t.Errorf("Sinh(%g) = %g, want %g", vf[i], f, sinh[i])
2372 
                }
2373 
        }
2374 
        for i := 0; i < len(vfsinhSC); i++ {
2375 
                if f := Sinh(vfsinhSC[i]); !alike(sinhSC[i], f) {
2376 
                        t.Errorf("Sinh(%g) = %g, want %g", vfsinhSC[i], f, sinhSC[i])
2377 
                }
2378 
        }
2379 
}
2380 
 
2381 
func TestSqrt(t *testing.T) {
2382 
        for i := 0; i < len(vf); i++ {
2383 
                a := Abs(vf[i])
2384 
                if f := SqrtGo(a); sqrt[i] != f {
2385 
                        t.Errorf("SqrtGo(%g) = %g, want %g", a, f, sqrt[i])
2386 
                }
2387 
                a = Abs(vf[i])
2388 
                if f := Sqrt(a); sqrt[i] != f {
2389 
                        t.Errorf("Sqrt(%g) = %g, want %g", a, f, sqrt[i])
2390 
                }
2391 
        }
2392 
        for i := 0; i < len(vfsqrtSC); i++ {
2393 
                if f := SqrtGo(vfsqrtSC[i]); !alike(sqrtSC[i], f) {
2394 
                        t.Errorf("SqrtGo(%g) = %g, want %g", vfsqrtSC[i], f, sqrtSC[i])
2395 
                }
2396 
                if f := Sqrt(vfsqrtSC[i]); !alike(sqrtSC[i], f) {
2397 
                        t.Errorf("Sqrt(%g) = %g, want %g", vfsqrtSC[i], f, sqrtSC[i])
2398 
                }
2399 
        }
2400 
}
2401 
 
2402 
func TestTan(t *testing.T) {
2403 
        for i := 0; i < len(vf); i++ {
2404 
                if f := Tan(vf[i]); !veryclose(tan[i], f) {
2405 
                        t.Errorf("Tan(%g) = %g, want %g", vf[i], f, tan[i])
2406 
                }
2407 
        }
2408 
        // same special cases as Sin
2409 
        for i := 0; i < len(vfsinSC); i++ {
2410 
                if f := Tan(vfsinSC[i]); !alike(sinSC[i], f) {
2411 
                        t.Errorf("Tan(%g) = %g, want %g", vfsinSC[i], f, sinSC[i])
2412 
                }
2413 
        }
2414 
}
2415 
 
2416 
func TestTanh(t *testing.T) {
2417 
        for i := 0; i < len(vf); i++ {
2418 
                if f := Tanh(vf[i]); !veryclose(tanh[i], f) {
2419 
                        t.Errorf("Tanh(%g) = %g, want %g", vf[i], f, tanh[i])
2420 
                }
2421 
        }
2422 
        for i := 0; i < len(vftanhSC); i++ {
2423 
                if f := Tanh(vftanhSC[i]); !alike(tanhSC[i], f) {
2424 
                        t.Errorf("Tanh(%g) = %g, want %g", vftanhSC[i], f, tanhSC[i])
2425 
                }
2426 
        }
2427 
}
2428 
 
2429 
func TestTrunc(t *testing.T) {
2430 
        for i := 0; i < len(vf); i++ {
2431 
                if f := Trunc(vf[i]); trunc[i] != f {
2432 
                        t.Errorf("Trunc(%g) = %g, want %g", vf[i], f, trunc[i])
2433 
                }
2434 
        }
2435 
        for i := 0; i < len(vfceilSC); i++ {
2436 
                if f := Trunc(vfceilSC[i]); !alike(ceilSC[i], f) {
2437 
                        t.Errorf("Trunc(%g) = %g, want %g", vfceilSC[i], f, ceilSC[i])
2438 
                }
2439 
        }
2440 
}
2441 
 
2442 
func TestY0(t *testing.T) {
2443 
        for i := 0; i < len(vf); i++ {
2444 
                a := Abs(vf[i])
2445 
                if f := Y0(a); !close(y0[i], f) {
2446 
                        t.Errorf("Y0(%g) = %g, want %g", a, f, y0[i])
2447 
                }
2448 
        }
2449 
        for i := 0; i < len(vfy0SC); i++ {
2450 
                if f := Y0(vfy0SC[i]); !alike(y0SC[i], f) {
2451 
                        t.Errorf("Y0(%g) = %g, want %g", vfy0SC[i], f, y0SC[i])
2452 
                }
2453 
        }
2454 
}
2455 
 
2456 
func TestY1(t *testing.T) {
2457 
        for i := 0; i < len(vf); i++ {
2458 
                a := Abs(vf[i])
2459 
                if f := Y1(a); !soclose(y1[i], f, 2e-14) {
2460 
                        t.Errorf("Y1(%g) = %g, want %g", a, f, y1[i])
2461 
                }
2462 
        }
2463 
        for i := 0; i < len(vfy0SC); i++ {
2464 
                if f := Y1(vfy0SC[i]); !alike(y1SC[i], f) {
2465 
                        t.Errorf("Y1(%g) = %g, want %g", vfy0SC[i], f, y1SC[i])
2466 
                }
2467 
        }
2468 
}
2469 
 
2470 
func TestYn(t *testing.T) {
2471 
        for i := 0; i < len(vf); i++ {
2472 
                a := Abs(vf[i])
2473 
                if f := Yn(2, a); !close(y2[i], f) {
2474 
                        t.Errorf("Yn(2, %g) = %g, want %g", a, f, y2[i])
2475 
                }
2476 
                if f := Yn(-3, a); !close(yM3[i], f) {
2477 
                        t.Errorf("Yn(-3, %g) = %g, want %g", a, f, yM3[i])
2478 
                }
2479 
        }
2480 
        for i := 0; i < len(vfy0SC); i++ {
2481 
                if f := Yn(2, vfy0SC[i]); !alike(y2SC[i], f) {
2482 
                        t.Errorf("Yn(2, %g) = %g, want %g", vfy0SC[i], f, y2SC[i])
2483 
                }
2484 
                if f := Yn(-3, vfy0SC[i]); !alike(yM3SC[i], f) {
2485 
                        t.Errorf("Yn(-3, %g) = %g, want %g", vfy0SC[i], f, yM3SC[i])
2486 
                }
2487 
        }
2488 
}
2489 
 
2490 
// Check that math functions of high angle values
2491 
// return accurate results. [Since (vf[i] + large) - large != vf[i],
2492 
// testing for Trig(vf[i] + large) == Trig(vf[i]), where large is
2493 
// a multiple of 2*Pi, is misleading.]
2494 
func TestLargeCos(t *testing.T) {
2495 
        large := float64(100000 * Pi)
2496 
        for i := 0; i < len(vf); i++ {
2497 
                f1 := cosLarge[i]
2498 
                f2 := Cos(vf[i] + large)
2499 
                if !close(f1, f2) {
2500 
                        t.Errorf("Cos(%g) = %g, want %g", vf[i]+large, f2, f1)
2501 
                }
2502 
        }
2503 
}
2504 
 
2505 
func TestLargeSin(t *testing.T) {
2506 
        large := float64(100000 * Pi)
2507 
        for i := 0; i < len(vf); i++ {
2508 
                f1 := sinLarge[i]
2509 
                f2 := Sin(vf[i] + large)
2510 
                if !close(f1, f2) {
2511 
                        t.Errorf("Sin(%g) = %g, want %g", vf[i]+large, f2, f1)
2512 
                }
2513 
        }
2514 
}
2515 
 
2516 
func TestLargeSincos(t *testing.T) {
2517 
        large := float64(100000 * Pi)
2518 
        for i := 0; i < len(vf); i++ {
2519 
                f1, g1 := sinLarge[i], cosLarge[i]
2520 
                f2, g2 := Sincos(vf[i] + large)
2521 
                if !close(f1, f2) || !close(g1, g2) {
2522 
                        t.Errorf("Sincos(%g) = %g, %g, want %g, %g", vf[i]+large, f2, g2, f1, g1)
2523 
                }
2524 
        }
2525 
}
2526 
 
2527 
func TestLargeTan(t *testing.T) {
2528 
        large := float64(100000 * Pi)
2529 
        for i := 0; i < len(vf); i++ {
2530 
                f1 := tanLarge[i]
2531 
                f2 := Tan(vf[i] + large)
2532 
                if !close(f1, f2) {
2533 
                        t.Errorf("Tan(%g) = %g, want %g", vf[i]+large, f2, f1)
2534 
                }
2535 
        }
2536 
}
2537 
 
2538 
// Check that math constants are accepted by compiler
2539 
// and have right value (assumes strconv.ParseFloat works).
2540 
// http://code.google.com/p/go/issues/detail?id=201
2541 
 
2542 
type floatTest struct {
2543 
        val  interface{}
2544 
        name string
2545 
        str  string
2546 
}
2547 
 
2548 
var floatTests = []floatTest{
2549 
        {float64(MaxFloat64), "MaxFloat64", "1.7976931348623157e+308"},
2550 
        {float64(SmallestNonzeroFloat64), "SmallestNonzeroFloat64", "5e-324"},
2551 
        {float32(MaxFloat32), "MaxFloat32", "3.4028235e+38"},
2552 
        {float32(SmallestNonzeroFloat32), "SmallestNonzeroFloat32", "1e-45"},
2553 
}
2554 
 
2555 
func TestFloatMinMax(t *testing.T) {
2556 
        for _, tt := range floatTests {
2557 
                s := fmt.Sprint(tt.val)
2558 
                if s != tt.str {
2559 
                        t.Errorf("Sprint(%v) = %s, want %s", tt.name, s, tt.str)
2560 
                }
2561 
        }
2562 
}
2563 
 
2564 
// Benchmarks
2565 
 
2566 
func BenchmarkAcos(b *testing.B) {
2567 
        for i := 0; i < b.N; i++ {
2568 
                Acos(.5)
2569 
        }
2570 
}
2571 
 
2572 
func BenchmarkAcosh(b *testing.B) {
2573 
        for i := 0; i < b.N; i++ {
2574 
                Acosh(1.5)
2575 
        }
2576 
}
2577 
 
2578 
func BenchmarkAsin(b *testing.B) {
2579 
        for i := 0; i < b.N; i++ {
2580 
                Asin(.5)
2581 
        }
2582 
}
2583 
 
2584 
func BenchmarkAsinh(b *testing.B) {
2585 
        for i := 0; i < b.N; i++ {
2586 
                Asinh(.5)
2587 
        }
2588 
}
2589 
 
2590 
func BenchmarkAtan(b *testing.B) {
2591 
        for i := 0; i < b.N; i++ {
2592 
                Atan(.5)
2593 
        }
2594 
}
2595 
 
2596 
func BenchmarkAtanh(b *testing.B) {
2597 
        for i := 0; i < b.N; i++ {
2598 
                Atanh(.5)
2599 
        }
2600 
}
2601 
 
2602 
func BenchmarkAtan2(b *testing.B) {
2603 
        for i := 0; i < b.N; i++ {
2604 
                Atan2(.5, 1)
2605 
        }
2606 
}
2607 
 
2608 
func BenchmarkCbrt(b *testing.B) {
2609 
        for i := 0; i < b.N; i++ {
2610 
                Cbrt(10)
2611 
        }
2612 
}
2613 
 
2614 
func BenchmarkCeil(b *testing.B) {
2615 
        for i := 0; i < b.N; i++ {
2616 
                Ceil(.5)
2617 
        }
2618 
}
2619 
 
2620 
func BenchmarkCopysign(b *testing.B) {
2621 
        for i := 0; i < b.N; i++ {
2622 
                Copysign(.5, -1)
2623 
        }
2624 
}
2625 
 
2626 
func BenchmarkCos(b *testing.B) {
2627 
        for i := 0; i < b.N; i++ {
2628 
                Cos(.5)
2629 
        }
2630 
}
2631 
 
2632 
func BenchmarkCosh(b *testing.B) {
2633 
        for i := 0; i < b.N; i++ {
2634 
                Cosh(2.5)
2635 
        }
2636 
}
2637 
 
2638 
func BenchmarkErf(b *testing.B) {
2639 
        for i := 0; i < b.N; i++ {
2640 
                Erf(.5)
2641 
        }
2642 
}
2643 
 
2644 
func BenchmarkErfc(b *testing.B) {
2645 
        for i := 0; i < b.N; i++ {
2646 
                Erfc(.5)
2647 
        }
2648 
}
2649 
 
2650 
func BenchmarkExp(b *testing.B) {
2651 
        for i := 0; i < b.N; i++ {
2652 
                Exp(.5)
2653 
        }
2654 
}
2655 
 
2656 
func BenchmarkExpGo(b *testing.B) {
2657 
        for i := 0; i < b.N; i++ {
2658 
                ExpGo(.5)
2659 
        }
2660 
}
2661 
 
2662 
func BenchmarkExpm1(b *testing.B) {
2663 
        for i := 0; i < b.N; i++ {
2664 
                Expm1(.5)
2665 
        }
2666 
}
2667 
 
2668 
func BenchmarkExp2(b *testing.B) {
2669 
        for i := 0; i < b.N; i++ {
2670 
                Exp2(.5)
2671 
        }
2672 
}
2673 
 
2674 
func BenchmarkExp2Go(b *testing.B) {
2675 
        for i := 0; i < b.N; i++ {
2676 
                Exp2Go(.5)
2677 
        }
2678 
}
2679 
 
2680 
func BenchmarkAbs(b *testing.B) {
2681 
        for i := 0; i < b.N; i++ {
2682 
                Abs(.5)
2683 
        }
2684 
}
2685 
 
2686 
func BenchmarkDim(b *testing.B) {
2687 
        for i := 0; i < b.N; i++ {
2688 
                Dim(10, 3)
2689 
        }
2690 
}
2691 
 
2692 
func BenchmarkFloor(b *testing.B) {
2693 
        for i := 0; i < b.N; i++ {
2694 
                Floor(.5)
2695 
        }
2696 
}
2697 
 
2698 
func BenchmarkMax(b *testing.B) {
2699 
        for i := 0; i < b.N; i++ {
2700 
                Max(10, 3)
2701 
        }
2702 
}
2703 
 
2704 
func BenchmarkMin(b *testing.B) {
2705 
        for i := 0; i < b.N; i++ {
2706 
                Min(10, 3)
2707 
        }
2708 
}
2709 
 
2710 
func BenchmarkMod(b *testing.B) {
2711 
        for i := 0; i < b.N; i++ {
2712 
                Mod(10, 3)
2713 
        }
2714 
}
2715 
 
2716 
func BenchmarkFrexp(b *testing.B) {
2717 
        for i := 0; i < b.N; i++ {
2718 
                Frexp(8)
2719 
        }
2720 
}
2721 
 
2722 
func BenchmarkGamma(b *testing.B) {
2723 
        for i := 0; i < b.N; i++ {
2724 
                Gamma(2.5)
2725 
        }
2726 
}
2727 
 
2728 
func BenchmarkHypot(b *testing.B) {
2729 
        for i := 0; i < b.N; i++ {
2730 
                Hypot(3, 4)
2731 
        }
2732 
}
2733 
 
2734 
func BenchmarkHypotGo(b *testing.B) {
2735 
        for i := 0; i < b.N; i++ {
2736 
                HypotGo(3, 4)
2737 
        }
2738 
}
2739 
 
2740 
func BenchmarkIlogb(b *testing.B) {
2741 
        for i := 0; i < b.N; i++ {
2742 
                Ilogb(.5)
2743 
        }
2744 
}
2745 
 
2746 
func BenchmarkJ0(b *testing.B) {
2747 
        for i := 0; i < b.N; i++ {
2748 
                J0(2.5)
2749 
        }
2750 
}
2751 
 
2752 
func BenchmarkJ1(b *testing.B) {
2753 
        for i := 0; i < b.N; i++ {
2754 
                J1(2.5)
2755 
        }
2756 
}
2757 
 
2758 
func BenchmarkJn(b *testing.B) {
2759 
        for i := 0; i < b.N; i++ {
2760 
                Jn(2, 2.5)
2761 
        }
2762 
}
2763 
 
2764 
func BenchmarkLdexp(b *testing.B) {
2765 
        for i := 0; i < b.N; i++ {
2766 
                Ldexp(.5, 2)
2767 
        }
2768 
}
2769 
 
2770 
func BenchmarkLgamma(b *testing.B) {
2771 
        for i := 0; i < b.N; i++ {
2772 
                Lgamma(2.5)
2773 
        }
2774 
}
2775 
 
2776 
func BenchmarkLog(b *testing.B) {
2777 
        for i := 0; i < b.N; i++ {
2778 
                Log(.5)
2779 
        }
2780 
}
2781 
 
2782 
func BenchmarkLogb(b *testing.B) {
2783 
        for i := 0; i < b.N; i++ {
2784 
                Logb(.5)
2785 
        }
2786 
}
2787 
 
2788 
func BenchmarkLog1p(b *testing.B) {
2789 
        for i := 0; i < b.N; i++ {
2790 
                Log1p(.5)
2791 
        }
2792 
}
2793 
 
2794 
func BenchmarkLog10(b *testing.B) {
2795 
        for i := 0; i < b.N; i++ {
2796 
                Log10(.5)
2797 
        }
2798 
}
2799 
 
2800 
func BenchmarkLog2(b *testing.B) {
2801 
        for i := 0; i < b.N; i++ {
2802 
                Log2(.5)
2803 
        }
2804 
}
2805 
 
2806 
func BenchmarkModf(b *testing.B) {
2807 
        for i := 0; i < b.N; i++ {
2808 
                Modf(1.5)
2809 
        }
2810 
}
2811 
 
2812 
func BenchmarkNextafter(b *testing.B) {
2813 
        for i := 0; i < b.N; i++ {
2814 
                Nextafter(.5, 1)
2815 
        }
2816 
}
2817 
 
2818 
func BenchmarkPowInt(b *testing.B) {
2819 
        for i := 0; i < b.N; i++ {
2820 
                Pow(2, 2)
2821 
        }
2822 
}
2823 
 
2824 
func BenchmarkPowFrac(b *testing.B) {
2825 
        for i := 0; i < b.N; i++ {
2826 
                Pow(2.5, 1.5)
2827 
        }
2828 
}
2829 
 
2830 
func BenchmarkPow10Pos(b *testing.B) {
2831 
        for i := 0; i < b.N; i++ {
2832 
                Pow10(300)
2833 
        }
2834 
}
2835 
 
2836 
func BenchmarkPow10Neg(b *testing.B) {
2837 
        for i := 0; i < b.N; i++ {
2838 
                Pow10(-300)
2839 
        }
2840 
}
2841 
 
2842 
func BenchmarkRemainder(b *testing.B) {
2843 
        for i := 0; i < b.N; i++ {
2844 
                Remainder(10, 3)
2845 
        }
2846 
}
2847 
 
2848 
func BenchmarkSignbit(b *testing.B) {
2849 
        for i := 0; i < b.N; i++ {
2850 
                Signbit(2.5)
2851 
        }
2852 
}
2853 
 
2854 
func BenchmarkSin(b *testing.B) {
2855 
        for i := 0; i < b.N; i++ {
2856 
                Sin(.5)
2857 
        }
2858 
}
2859 
 
2860 
func BenchmarkSincos(b *testing.B) {
2861 
        for i := 0; i < b.N; i++ {
2862 
                Sincos(.5)
2863 
        }
2864 
}
2865 
 
2866 
func BenchmarkSinh(b *testing.B) {
2867 
        for i := 0; i < b.N; i++ {
2868 
                Sinh(2.5)
2869 
        }
2870 
}
2871 
 
2872 
func BenchmarkSqrt(b *testing.B) {
2873 
        for i := 0; i < b.N; i++ {
2874 
                Sqrt(10)
2875 
        }
2876 
}
2877 
 
2878 
func BenchmarkSqrtGo(b *testing.B) {
2879 
        for i := 0; i < b.N; i++ {
2880 
                SqrtGo(10)
2881 
        }
2882 
}
2883 
 
2884 
func BenchmarkTan(b *testing.B) {
2885 
        for i := 0; i < b.N; i++ {
2886 
                Tan(.5)
2887 
        }
2888 
}
2889 
 
2890 
func BenchmarkTanh(b *testing.B) {
2891 
        for i := 0; i < b.N; i++ {
2892 
                Tanh(2.5)
2893 
        }
2894 
}
2895 
func BenchmarkTrunc(b *testing.B) {
2896 
        for i := 0; i < b.N; i++ {
2897 
                Trunc(.5)
2898 
        }
2899 
}
2900 
 
2901 
func BenchmarkY0(b *testing.B) {
2902 
        for i := 0; i < b.N; i++ {
2903 
                Y0(2.5)
2904 
        }
2905 
}
2906 
 
2907 
func BenchmarkY1(b *testing.B) {
2908 
        for i := 0; i < b.N; i++ {
2909 
                Y1(2.5)
2910 
        }
2911 
}
2912 
 
2913 
func BenchmarkYn(b *testing.B) {
2914 
        for i := 0; i < b.N; i++ {
2915 
                Yn(2, 2.5)
2916 
        }
2917 
}

powered by: WebSVN 2.1.0

© copyright 1999-2024 OpenCores.org, equivalent to Oliscience, all rights reserved. OpenCores®, registered trademark.