liyujie
2025-08-28 d9927380ed7c8366f762049be9f3fee225860833
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
// Copyright 2017 The Go Authors. All rights reserved.
// Use of this source code is governed by a BSD-style
// license that can be found in the LICENSE file.
 
#include "textflag.h"
 
#define Neg2p11 0xC000E147AE147AE1
#define Pos15   0x402E
 
// Minimax polynomial coefficients and other constants
DATA ·erfcrodataL38<> + 0(SB)/8, $.234875460637085087E-01
DATA ·erfcrodataL38<> + 8(SB)/8, $.234469449299256284E-01
DATA ·erfcrodataL38<> + 16(SB)/8, $-.606918710392844955E-04
DATA ·erfcrodataL38<> + 24(SB)/8, $-.198827088077636213E-04
DATA ·erfcrodataL38<> + 32(SB)/8, $.257805645845475331E-06
DATA ·erfcrodataL38<> + 40(SB)/8, $-.184427218110620284E-09
DATA ·erfcrodataL38<> + 48(SB)/8, $.122408098288933181E-10
DATA ·erfcrodataL38<> + 56(SB)/8, $.484691106751495392E-07
DATA ·erfcrodataL38<> + 64(SB)/8, $-.150147637632890281E-08
DATA ·erfcrodataL38<> + 72(SB)/8, $23.999999999973521625
DATA ·erfcrodataL38<> + 80(SB)/8, $27.226017111108365754
DATA ·erfcrodataL38<> + 88(SB)/8, $-2.0
DATA ·erfcrodataL38<> + 96(SB)/8, $0.100108802034478228E+00
DATA ·erfcrodataL38<> + 104(SB)/8, $0.244588413746558125E+00
DATA ·erfcrodataL38<> + 112(SB)/8, $-.669188879646637174E-01
DATA ·erfcrodataL38<> + 120(SB)/8, $0.151311447000953551E-01
DATA ·erfcrodataL38<> + 128(SB)/8, $-.284720833493302061E-02
DATA ·erfcrodataL38<> + 136(SB)/8, $0.455491239358743212E-03
DATA ·erfcrodataL38<> + 144(SB)/8, $-.631850539280720949E-04
DATA ·erfcrodataL38<> + 152(SB)/8, $0.772532660726086679E-05
DATA ·erfcrodataL38<> + 160(SB)/8, $-.843706007150936940E-06
DATA ·erfcrodataL38<> + 168(SB)/8, $-.735330214904227472E-08
DATA ·erfcrodataL38<> + 176(SB)/8, $0.753002008837084967E-09
DATA ·erfcrodataL38<> + 184(SB)/8, $0.832482036660624637E-07
DATA ·erfcrodataL38<> + 192(SB)/8, $-0.75
DATA ·erfcrodataL38<> + 200(SB)/8, $.927765678007128609E-01
DATA ·erfcrodataL38<> + 208(SB)/8, $.903621209344751506E-01
DATA ·erfcrodataL38<> + 216(SB)/8, $-.344203375025257265E-02
DATA ·erfcrodataL38<> + 224(SB)/8, $-.869243428221791329E-03
DATA ·erfcrodataL38<> + 232(SB)/8, $.174699813107105603E-03
DATA ·erfcrodataL38<> + 240(SB)/8, $.649481036316130000E-05
DATA ·erfcrodataL38<> + 248(SB)/8, $-.895265844897118382E-05
DATA ·erfcrodataL38<> + 256(SB)/8, $.135970046909529513E-05
DATA ·erfcrodataL38<> + 264(SB)/8, $.277617717014748015E-06
DATA ·erfcrodataL38<> + 272(SB)/8, $.810628018408232910E-08
DATA ·erfcrodataL38<> + 280(SB)/8, $.210430084693497985E-07
DATA ·erfcrodataL38<> + 288(SB)/8, $-.342138077525615091E-08
DATA ·erfcrodataL38<> + 296(SB)/8, $-.165467946798610800E-06
DATA ·erfcrodataL38<> + 304(SB)/8, $5.999999999988412824
DATA ·erfcrodataL38<> + 312(SB)/8, $.468542210149072159E-01
DATA ·erfcrodataL38<> + 320(SB)/8, $.465343528567604256E-01
DATA ·erfcrodataL38<> + 328(SB)/8, $-.473338083650201733E-03
DATA ·erfcrodataL38<> + 336(SB)/8, $-.147220659069079156E-03
DATA ·erfcrodataL38<> + 344(SB)/8, $.755284723554388339E-05
DATA ·erfcrodataL38<> + 352(SB)/8, $.116158570631428789E-05
DATA ·erfcrodataL38<> + 360(SB)/8, $-.155445501551602389E-06
DATA ·erfcrodataL38<> + 368(SB)/8, $-.616940119847805046E-10
DATA ·erfcrodataL38<> + 376(SB)/8, $-.728705590727563158E-10
DATA ·erfcrodataL38<> + 384(SB)/8, $-.983452460354586779E-08
DATA ·erfcrodataL38<> + 392(SB)/8, $.365156164194346316E-08
DATA ·erfcrodataL38<> + 400(SB)/8, $11.999999999996530775
DATA ·erfcrodataL38<> + 408(SB)/8, $0.467773498104726584E-02
DATA ·erfcrodataL38<> + 416(SB)/8, $0.206669853540920535E-01
DATA ·erfcrodataL38<> + 424(SB)/8, $0.413339707081841473E-01
DATA ·erfcrodataL38<> + 432(SB)/8, $0.482229658262131320E-01
DATA ·erfcrodataL38<> + 440(SB)/8, $0.344449755901841897E-01
DATA ·erfcrodataL38<> + 448(SB)/8, $0.130890907240765465E-01
DATA ·erfcrodataL38<> + 456(SB)/8, $-.459266344100642687E-03
DATA ·erfcrodataL38<> + 464(SB)/8, $-.337888800856913728E-02
DATA ·erfcrodataL38<> + 472(SB)/8, $-.159103061687062373E-02
DATA ·erfcrodataL38<> + 480(SB)/8, $-.501128905515922644E-04
DATA ·erfcrodataL38<> + 488(SB)/8, $0.262775855852903132E-03
DATA ·erfcrodataL38<> + 496(SB)/8, $0.103860982197462436E-03
DATA ·erfcrodataL38<> + 504(SB)/8, $-.548835785414200775E-05
DATA ·erfcrodataL38<> + 512(SB)/8, $-.157075054646618214E-04
DATA ·erfcrodataL38<> + 520(SB)/8, $-.480056366276045110E-05
DATA ·erfcrodataL38<> + 528(SB)/8, $0.198263013759701555E-05
DATA ·erfcrodataL38<> + 536(SB)/8, $-.224394262958888780E-06
DATA ·erfcrodataL38<> + 544(SB)/8, $-.321853693146683428E-06
DATA ·erfcrodataL38<> + 552(SB)/8, $0.445073894984683537E-07
DATA ·erfcrodataL38<> + 560(SB)/8, $0.660425940000555729E-06
DATA ·erfcrodataL38<> + 568(SB)/8, $2.0
DATA ·erfcrodataL38<> + 576(SB)/8, $8.63616855509444462538e-78
DATA ·erfcrodataL38<> + 584(SB)/8, $1.00000000000000222044
DATA ·erfcrodataL38<> + 592(SB)/8, $0.500000000000004237e+00
DATA ·erfcrodataL38<> + 600(SB)/8, $0.416666664838056960e-01
DATA ·erfcrodataL38<> + 608(SB)/8, $0.166666666630345592e+00
DATA ·erfcrodataL38<> + 616(SB)/8, $0.138926439368309441e-02
DATA ·erfcrodataL38<> + 624(SB)/8, $0.833349307718286047e-02
DATA ·erfcrodataL38<> + 632(SB)/8, $-.693147180558298714e+00
DATA ·erfcrodataL38<> + 640(SB)/8, $-.164659495826017651e-11
DATA ·erfcrodataL38<> + 648(SB)/8, $.179001151181866548E+00
DATA ·erfcrodataL38<> + 656(SB)/8, $-.144269504088896339e+01
DATA ·erfcrodataL38<> + 664(SB)/8, $+281475245147134.9375
DATA ·erfcrodataL38<> + 672(SB)/8, $.163116780021877404E+00
DATA ·erfcrodataL38<> + 680(SB)/8, $-.201574395828120710E-01
DATA ·erfcrodataL38<> + 688(SB)/8, $-.185726336009394125E-02
DATA ·erfcrodataL38<> + 696(SB)/8, $.199349204957273749E-02
DATA ·erfcrodataL38<> + 704(SB)/8, $-.554902415532606242E-03
DATA ·erfcrodataL38<> + 712(SB)/8, $-.638914789660242846E-05
DATA ·erfcrodataL38<> + 720(SB)/8, $-.424441522653742898E-04
DATA ·erfcrodataL38<> + 728(SB)/8, $.827967511921486190E-04
DATA ·erfcrodataL38<> + 736(SB)/8, $.913965446284062654E-05
DATA ·erfcrodataL38<> + 744(SB)/8, $.277344791076320853E-05
DATA ·erfcrodataL38<> + 752(SB)/8, $-.467239678927239526E-06
DATA ·erfcrodataL38<> + 760(SB)/8, $.344814065920419986E-07
DATA ·erfcrodataL38<> + 768(SB)/8, $-.366013491552527132E-05
DATA ·erfcrodataL38<> + 776(SB)/8, $.181242810023783439E-05
DATA ·erfcrodataL38<> + 784(SB)/8, $2.999999999991234567
DATA ·erfcrodataL38<> + 792(SB)/8, $1.0
GLOBL ·erfcrodataL38<> + 0(SB), RODATA, $800
 
// Table of log correction terms
DATA ·erfctab2069<> + 0(SB)/8, $0.442737824274138381e-01
DATA ·erfctab2069<> + 8(SB)/8, $0.263602189790660309e-01
DATA ·erfctab2069<> + 16(SB)/8, $0.122565642281703586e-01
DATA ·erfctab2069<> + 24(SB)/8, $0.143757052860721398e-02
DATA ·erfctab2069<> + 32(SB)/8, $-.651375034121276075e-02
DATA ·erfctab2069<> + 40(SB)/8, $-.119317678849450159e-01
DATA ·erfctab2069<> + 48(SB)/8, $-.150868749549871069e-01
DATA ·erfctab2069<> + 56(SB)/8, $-.161992609578469234e-01
DATA ·erfctab2069<> + 64(SB)/8, $-.154492360403337917e-01
DATA ·erfctab2069<> + 72(SB)/8, $-.129850717389178721e-01
DATA ·erfctab2069<> + 80(SB)/8, $-.892902649276657891e-02
DATA ·erfctab2069<> + 88(SB)/8, $-.338202636596794887e-02
DATA ·erfctab2069<> + 96(SB)/8, $0.357266307045684762e-02
DATA ·erfctab2069<> + 104(SB)/8, $0.118665304327406698e-01
DATA ·erfctab2069<> + 112(SB)/8, $0.214434994118118914e-01
DATA ·erfctab2069<> + 120(SB)/8, $0.322580645161290314e-01
GLOBL ·erfctab2069<> + 0(SB), RODATA, $128
 
// Erfc returns the complementary error function of the argument.
//
// Special cases are:
//      Erfc(+Inf) = 0
//      Erfc(-Inf) = 2
//      Erfc(NaN) = NaN
// The algorithm used is minimax polynomial approximation
// with coefficients determined with a Remez exchange algorithm.
// This assembly implementation handles inputs in the range [-2.11, +15].
// For all other inputs we call the generic Go implementation.
 
TEXT    ·erfcAsm(SB), NOSPLIT|NOFRAME, $0-16
   MOVD    x+0(FP), R1
   MOVD    $Neg2p11, R2
   CMPUBGT    R1, R2, usego
 
   FMOVD    x+0(FP), F0
   MOVD    $·erfcrodataL38<>+0(SB), R9
   FMOVD    F0, F2
   SRAD    $48, R1
   MOVH    R1, R2
   ANDW    $0x7FFF, R1
   MOVH    $Pos15, R3
   CMPW    R1, R3
   BGT    usego
   MOVH    $0x3FFF, R3
   MOVW    R1, R6
   MOVW    R3, R7
   CMPBGT    R6, R7, L2
   MOVH    $0x3FEF, R3
   MOVW    R3, R7
   CMPBGT    R6, R7, L3
   MOVH    $0x2FFF, R2
   MOVW    R2, R7
   CMPBGT    R6, R7, L4
   FMOVD    792(R9), F0
   WFSDB    V2, V0, V2
   FMOVD    F2, ret+8(FP)
   RET
 
L2:
   LTDBR    F0, F0
   MOVH    $0x0, R4
   BLTU    L3
   FMOVD    F0, F1
L9:
   MOVH    $0x400F, R3
   MOVW    R1, R6
   MOVW    R3, R7
   CMPBGT    R6, R7, L10
   FMOVD    784(R9), F3
   FSUB    F1, F3
   VLEG    $0, 776(R9), V20
   WFDDB    V1, V3, V6
   VLEG    $0, 768(R9), V18
   FMOVD    760(R9), F7
   FMOVD    752(R9), F5
   VLEG    $0, 744(R9), V16
   FMOVD    736(R9), F3
   FMOVD    728(R9), F2
   FMOVD    720(R9), F4
   WFMDB    V6, V6, V1
   FMUL    F0, F0
   MOVH    $0x0, R3
   WFMADB    V1, V7, V20, V7
   WFMADB    V1, V5, V18, V5
   WFMADB    V1, V7, V16, V7
   WFMADB    V1, V5, V3, V5
   WFMADB    V1, V7, V4, V7
   WFMADB    V1, V5, V2, V5
   FMOVD    712(R9), F2
   WFMADB    V1, V7, V2, V7
   FMOVD    704(R9), F2
   WFMADB    V1, V5, V2, V5
   FMOVD    696(R9), F2
   WFMADB    V1, V7, V2, V7
   FMOVD    688(R9), F2
   MOVH    $0x0, R1
   WFMADB    V1, V5, V2, V5
   FMOVD    680(R9), F2
   WFMADB    V1, V7, V2, V7
   FMOVD    672(R9), F2
   WFMADB    V1, V5, V2, V1
   FMOVD    664(R9), F3
   WFMADB    V6, V7, V1, V7
   FMOVD    656(R9), F5
   FMOVD    648(R9), F2
   WFMADB    V0, V5, V3, V5
   WFMADB    V6, V7, V2, V7
L11:
   LGDR    F5, R6
   WFSDB    V0, V0, V2
   WORD    $0xED509298    //sdb    %f5,.L55-.L38(%r9)
   BYTE    $0x00
   BYTE    $0x1B
   FMOVD    640(R9), F6
   FMOVD    632(R9), F4
   WFMSDB    V5, V6, V2, V6
   WFMSDB    V5, V4, V0, V4
   FMOVD    624(R9), F2
   FADD    F6, F4
   FMOVD    616(R9), F0
   FMOVD    608(R9), F6
   WFMADB    V4, V0, V2, V0
   FMOVD    600(R9), F3
   WFMDB    V4, V4, V2
   MOVH    R6,R6
   ADD    R6, R3
   WFMADB    V4, V3, V6, V3
   FMOVD    592(R9), F6
   WFMADB    V0, V2, V3, V0
   FMOVD    584(R9), F3
   WFMADB    V4, V6, V3, V6
   WORD    $0xECC339BC    //risbg    %r12,%r3,57,128+60,3
   BYTE    $0x03
   BYTE    $0x55
   WFMADB    V2, V0, V6, V0
   MOVD    $·erfctab2069<>+0(SB), R5
   WORD    $0x682C5000    //ld    %f2,0(%r12,%r5)
   FMADD    F2, F4, F4
   WORD    $0xEC43000F    //risbgn    %r4,%r3,64-64+0,64-64+0+16-1,64-0-16
   BYTE    $0x30
   BYTE    $0x59
   WFMADB    V4, V0, V2, V4
   LDGR    R4, F2
   FMADD    F4, F2, F2
   MOVW    R2, R6
   CMPBLE    R6, $0, L20
   MOVW    R1, R6
   CMPBEQ    R6, $0, L21
   WORD    $0xED709240    //mdb    %f7,.L66-.L38(%r9)
   BYTE    $0x00
   BYTE    $0x1C
L21:
   FMUL    F7, F2
L1:
   FMOVD    F2, ret+8(FP)
   RET
L3:
   LTDBR    F0, F0
   BLTU    L30
   FMOVD    568(R9), F2
   WFSDB    V0, V2, V0
L8:
   WFMDB    V0, V0, V4
   FMOVD    560(R9), F2
   FMOVD    552(R9), F6
   FMOVD    544(R9), F1
   WFMADB    V4, V6, V2, V6
   FMOVD    536(R9), F2
   WFMADB    V4, V1, V2, V1
   FMOVD    528(R9), F3
   FMOVD    520(R9), F2
   WFMADB    V4, V6, V3, V6
   WFMADB    V4, V1, V2, V1
   FMOVD    512(R9), F3
   FMOVD    504(R9), F2
   WFMADB    V4, V6, V3, V6
   WFMADB    V4, V1, V2, V1
   FMOVD    496(R9), F3
   FMOVD    488(R9), F2
   WFMADB    V4, V6, V3, V6
   WFMADB    V4, V1, V2, V1
   FMOVD    480(R9), F3
   FMOVD    472(R9), F2
   WFMADB    V4, V6, V3, V6
   WFMADB    V4, V1, V2, V1
   FMOVD    464(R9), F3
   FMOVD    456(R9), F2
   WFMADB    V4, V6, V3, V6
   WFMADB    V4, V1, V2, V1
   FMOVD    448(R9), F3
   FMOVD    440(R9), F2
   WFMADB    V4, V6, V3, V6
   WFMADB    V4, V1, V2, V1
   FMOVD    432(R9), F3
   FMOVD    424(R9), F2
   WFMADB    V4, V6, V3, V6
   WFMADB    V4, V1, V2, V1
   FMOVD    416(R9), F3
   FMOVD    408(R9), F2
   WFMADB    V4, V6, V3, V6
   FMADD    F1, F4, F2
   FMADD    F6, F0, F2
   MOVW    R2, R6
   CMPBGE    R6, $0, L1
   FMOVD    568(R9), F0
   WFSDB    V2, V0, V2
   BR    L1
L10:
   MOVH    $0x401F, R3
   MOVW    R1, R6
   MOVW    R3, R7
   CMPBLE    R6, R7, L36
   MOVH    $0x402F, R3
   MOVW    R3, R7
   CMPBGT    R6, R7, L13
   FMOVD    400(R9), F3
   FSUB    F1, F3
   VLEG    $0, 392(R9), V20
   WFDDB    V1, V3, V6
   VLEG    $0, 384(R9), V18
   FMOVD    376(R9), F2
   FMOVD    368(R9), F4
   VLEG    $0, 360(R9), V16
   FMOVD    352(R9), F7
   FMOVD    344(R9), F3
   FMUL    F0, F0
   WFMDB    V6, V6, V1
   FMOVD    656(R9), F5
   MOVH    $0x0, R3
   WFMADB    V1, V2, V20, V2
   WFMADB    V1, V4, V18, V4
   WFMADB    V1, V2, V16, V2
   WFMADB    V1, V4, V7, V4
   WFMADB    V1, V2, V3, V2
   FMOVD    336(R9), F3
   WFMADB    V1, V4, V3, V4
   FMOVD    328(R9), F3
   WFMADB    V1, V2, V3, V2
   FMOVD    320(R9), F3
   WFMADB    V1, V4, V3, V1
   FMOVD    312(R9), F7
   WFMADB    V6, V2, V1, V2
   MOVH    $0x0, R1
   FMOVD    664(R9), F3
   FMADD    F2, F6, F7
   WFMADB    V0, V5, V3, V5
   BR    L11
L35:
   WORD    $0xB3130010    //lcdbr    %f1,%f0
   BR    L9
L36:
   FMOVD    304(R9), F3
   FSUB    F1, F3
   VLEG    $0, 296(R9), V20
   WFDDB    V1, V3, V6
   FMOVD    288(R9), F5
   FMOVD    280(R9), F1
   FMOVD    272(R9), F2
   VLEG    $0, 264(R9), V18
   VLEG    $0, 256(R9), V16
   FMOVD    248(R9), F3
   FMOVD    240(R9), F4
   WFMDB    V6, V6, V7
   FMUL    F0, F0
   MOVH    $0x0, R3
   FMADD    F5, F7, F1
   WFMADB    V7, V2, V20, V2
   WFMADB    V7, V1, V18, V1
   WFMADB    V7, V2, V16, V2
   WFMADB    V7, V1, V3, V1
   WFMADB    V7, V2, V4, V2
   FMOVD    232(R9), F4
   WFMADB    V7, V1, V4, V1
   FMOVD    224(R9), F4
   WFMADB    V7, V2, V4, V2
   FMOVD    216(R9), F4
   WFMADB    V7, V1, V4, V1
   FMOVD    208(R9), F4
   MOVH    $0x0, R1
   WFMADB    V7, V2, V4, V7
   FMOVD    656(R9), F5
   WFMADB    V6, V1, V7, V1
   FMOVD    664(R9), F3
   FMOVD    200(R9), F7
   WFMADB    V0, V5, V3, V5
   FMADD    F1, F6, F7
   BR    L11
L4:
   FMOVD    192(R9), F1
   FMADD    F0, F0, F1
   FMOVD    184(R9), F3
   WFMDB    V1, V1, V0
   FMOVD    176(R9), F4
   FMOVD    168(R9), F6
   WFMADB    V0, V4, V3, V4
   FMOVD    160(R9), F3
   WFMADB    V0, V6, V3, V6
   FMOVD    152(R9), F5
   FMOVD    144(R9), F3
   WFMADB    V0, V4, V5, V4
   WFMADB    V0, V6, V3, V6
   FMOVD    136(R9), F5
   FMOVD    128(R9), F3
   WFMADB    V0, V4, V5, V4
   WFMADB    V0, V6, V3, V6
   FMOVD    120(R9), F5
   FMOVD    112(R9), F3
   WFMADB    V0, V4, V5, V4
   WFMADB    V0, V6, V3, V6
   FMOVD    104(R9), F5
   FMOVD    96(R9), F3
   WFMADB    V0, V4, V5, V4
   WFMADB    V0, V6, V3, V0
   FMOVD    F2, F6
   FMADD    F4, F1, F0
   WORD    $0xED609318    //sdb    %f6,.L39-.L38(%r9)
   BYTE    $0x00
   BYTE    $0x1B
   WFMSDB    V2, V0, V6, V2
   FMOVD    F2, ret+8(FP)
   RET
L30:
   WORD    $0xED009238    //adb    %f0,.L67-.L38(%r9)
   BYTE    $0x00
   BYTE    $0x1A
   BR    L8
L20:
   FMOVD    88(R9), F0
   WFMADB    V7, V2, V0, V2
   WORD    $0xB3130022    //lcdbr    %f2,%f2
   FMOVD    F2, ret+8(FP)
   RET
L13:
   MOVH    $0x403A, R3
   MOVW    R1, R6
   MOVW    R3, R7
   CMPBLE    R6, R7, L4
   WORD    $0xED109050    //cdb    %f1,.L128-.L38(%r9)
   BYTE    $0x00
   BYTE    $0x19
   BGE    L37
   BVS    L37
   FMOVD    72(R9), F6
   FSUB    F1, F6
   MOVH    $0x1000, R3
   FDIV    F1, F6
   MOVH    $0x1000, R1
L17:
   WFMDB    V6, V6, V1
   FMOVD    64(R9), F2
   FMOVD    56(R9), F4
   FMOVD    48(R9), F3
   WFMADB    V1, V3, V2, V3
   FMOVD    40(R9), F2
   WFMADB    V1, V2, V4, V2
   FMOVD    32(R9), F4
   WFMADB    V1, V3, V4, V3
   FMOVD    24(R9), F4
   WFMADB    V1, V2, V4, V2
   FMOVD    16(R9), F4
   WFMADB    V1, V3, V4, V3
   FMOVD    8(R9), F4
   WFMADB    V1, V2, V4, V1
   FMUL    F0, F0
   WFMADB    V3, V6, V1, V3
   FMOVD    656(R9), F5
   FMOVD    664(R9), F4
   FMOVD    0(R9), F7
   WFMADB    V0, V5, V4, V5
   FMADD    F6, F3, F7
   BR    L11
L14:
   FMOVD    72(R9), F6
   FSUB    F1, F6
   MOVH    $0x403A, R3
   FDIV    F1, F6
   MOVW    R1, R6
   MOVW    R3, R7
   CMPBEQ    R6, R7, L23
   MOVH    $0x0, R3
   MOVH    $0x0, R1
   BR    L17
L37:
   WFCEDBS    V0, V0, V0
   BVS    L1
   MOVW    R2, R6
   CMPBLE    R6, $0, L18
   MOVH    $0x7FEF, R2
   MOVW    R1, R6
   MOVW    R2, R7
   CMPBGT    R6, R7, L24
 
   WORD    $0xA5400010    //iihh    %r4,16
   LDGR    R4, F2
   FMUL    F2, F2
   BR    L1
L23:
   MOVH    $0x1000, R3
   MOVH    $0x1000, R1
   BR    L17
L24:
   FMOVD    $0, F2
   BR    L1
L18:
   MOVH    $0x7FEF, R2
   MOVW    R1, R6
   MOVW    R2, R7
   CMPBGT    R6, R7, L25
   WORD    $0xA5408010    //iihh    %r4,32784
   FMOVD    568(R9), F2
   LDGR    R4, F0
   FMADD    F2, F0, F2
   BR    L1
L25:
   FMOVD    568(R9), F2
   BR    L1
usego:
   BR    ·erfc(SB)