![[BACK]](/cvsweb/icons/back.gif){kind=icon}
Return to paranoia.s CVS log ![[TXT]](/cvsweb/icons/text.gif){kind=icon}
![[DIR]](/cvsweb/icons/dir.gif){kind=icon}
Up to [Research Unix] / researchv10no / cmd / lcc / gen2 / mc-sun / tst|
Return to paranoia.s CVS log | Up to [Research Unix] / researchv10no / cmd / lcc / gen2 / mc-sun / tst |
1.1 root 1: .data
2: .globl _Zero
3: .align 2; _Zero:.long 0x0,0x0
4: .globl _Half
5: .align 2; _Half:.long 0x3fe00000,0x0
6: .globl _One
7: .align 2; _One:.long 0x3ff00000,0x0
8: .globl _Two
9: .align 2; _Two:.long 0x40000000,0x0
10: .globl _Three
11: .align 2; _Three:.long 0x40080000,0x0
12: .globl _Four
13: .align 2; _Four:.long 0x40100000,0x0
14: .globl _Five
15: .align 2; _Five:.long 0x40140000,0x0
16: .globl _Eight
17: .align 2; _Eight:.long 0x40200000,0x0
18: .globl _Nine
19: .align 2; _Nine:.long 0x40220000,0x0
20: .globl _TwentySeven
21: .align 2; _TwentySeven:.long 0x403b0000,0x0
22: .globl _ThirtyTwo
23: .align 2; _ThirtyTwo:.long 0x40400000,0x0
24: .globl _TwoForty
25: .align 2; _TwoForty:.long 0x406e0000,0x0
26: .globl _MinusOne
27: .align 2; _MinusOne:.long 0xbff00000,0x0
28: .globl _OneAndHalf
29: .align 2; _OneAndHalf:.long 0x3ff80000,0x0
30: .globl _NoTrials
31: .align 2; _NoTrials:.long 20
32: .text
33: .globl _sigfpe
34: .align 2
35: _sigfpe:link a6,#-8
36: addql #1,_fpecount:l
37: movl #L2,sp@
38: jbsr _printf:l
39: movl #__iob+20,sp@
40: jbsr _fflush:l
41: tstl _sigsave:l; jeq L4
42: movl #8,sp@
43: movl _sigsave:l,sp@(4)
44: jbsr _signal:l
45: clrl _sigsave:l
46: movl #_ovfl_buf,sp@
47: movl #1,sp@(4)
48: jbsr _longjmp:l
49: L4:
50: jbsr _abort:l
51: L1:
52: unlk a6
53: rts
54: .data
55: .align 2; L762:.long L763
56: .long L764
57: .long L765
58: .long L766
59: .text
60: .globl _main
61: .align 2
62: _main:link a6,#-144
63: fmovemx #0xff,a6@(-100)
64: moveml #0x498,a6@(-116)
65: movl #L7,sp@
66: movl #L8,sp@(4)
67: movl #L9,sp@(8)
68: lea a6@(-4),a2
69: movl a2,sp@(12)
70: jbsr _ieee_flags:l
71: fmoved L10:l,fp3
72: fmoved fp3,_Zero:l
73: fmoved L11:l,fp3
74: fmoved fp3,_One:l
75: fmoved _One:l,fp3
76: fmovex fp3,fp2; faddx fp3,fp2
77: fmoved fp2,_Two:l
78: fmoved _Two:l,fp2
79: fmovex fp2,fp7
80: fmovex fp2,fp1; faddx fp3,fp1
81: fmoved fp1,_Three:l
82: fmoved _Three:l,fp1
83: fmovex fp1,fp0; faddx fp3,fp0
84: fmoved fp0,_Four:l
85: fmoved _Four:l,fp0
86: fmovex fp0,fp2; faddx fp3,fp2
87: fmoved fp2,_Five:l
88: fmovex fp0,fp2; faddx fp0,fp2
89: fmoved fp2,_Eight:l
90: fmovex fp1,fp2; fmulx fp1,fp2
91: fmoved fp2,_Nine:l
92: fmoved _Nine:l,fp2
93: fmulx fp1,fp2
94: fmoved fp2,_TwentySeven:l
95: fmovex fp0,fp2; fmuld _Eight:l,fp2
96: fmoved fp2,_ThirtyTwo:l
97: fmovex fp0,fp2; fmuld _Five:l,fp2
98: fmulx fp1,fp2
99: fmulx fp0,fp2
100: fmoved fp2,_TwoForty:l
101: fnegx fp3,fp2
102: fmoved fp2,_MinusOne:l
103: fmovex fp3,fp2; fdivx fp7,fp2
104: fmoved fp2,_Half:l
105: faddd _Half:l,fp3
106: fmoved fp3,_OneAndHalf:l
107: clrl _ErrCnt:l
108: clrl _ErrCnt+4:l
109: clrl _ErrCnt+8:l
110: clrl _ErrCnt+12:l
111: movl #1,_PageNo:l
112: clrl _Milestone:l
113: movl #8,sp@
114: movl #_sigfpe,sp@(4)
115: jbsr _signal:l
116: jbsr _Instructions:l
117: jbsr _Pause:l
118: jbsr _Heading:l
119: jbsr _Pause:l
120: jbsr _Characteristics:l
121: jbsr _Pause:l
122: jbsr _History:l
123: jbsr _Pause:l
124: movl #7,_Milestone:l
125: movl #L15,sp@
126: jbsr _printf:l
127: movl #0,sp@
128: fmoved _Zero:l,fp3
129: fmovex fp3,fp2; faddx fp3,fp2
130: fcmpx fp3,fp2; fjne L18
131: fmoved _One:l,fp2
132: fmovex fp2,fp1; fsubx fp2,fp1
133: fcmpx fp3,fp1; fjne L18
134: fcmpx fp3,fp2; fjle L18
135: fmovex fp2,fp3; faddx fp2,fp3
136: fcmpd _Two:l,fp3; fjne L18
137: movl #1,d7
138: jra L19
139: L18:
140: clrl d7
141: L19:
142: movl d7,sp@(4)
143: movl #L17,sp@(8)
144: jbsr _TstCond:l
145: fnegd _Zero:l,fp3
146: fmoved fp3,_Z:l
147: fmoved _Z:l,fp3
148: fcmpd L10:l,fp3; fjeq L20
149: addql #1,_ErrCnt:l
150: movl #L22,sp@
151: jbsr _printf:l
152: fmoved L23:l,fp3
153: fmoved fp3,_U1:l
154: fmoved L11:l,fp3
155: fmoved fp3,_Radix:l
156: jbsr _TstPtUf:l
157: L20:
158: movl #0,sp@
159: fmoved _Three:l,fp3
160: fmoved _Two:l,fp2
161: fmovex fp2,fp7
162: fmoved _One:l,fp1
163: fmovex fp1,fp6
164: fmovex fp2,fp0; faddx fp1,fp0
165: fcmpx fp0,fp3; fjne L26
166: fmoved _Four:l,fp0
167: fmovex fp3,fp2; faddx fp1,fp2
168: fcmpx fp2,fp0; fjne L26
169: fnegx fp7,fp2
170: fmovex fp7,fp1; fmulx fp2,fp1
171: fmovex fp0,fp2; faddx fp1,fp2
172: fmoved _Zero:l,fp1
173: fcmpx fp1,fp2; fjne L26
174: fmovex fp0,fp2; fsubx fp3,fp2
175: fmovex fp2,fp3; fsubx fp6,fp3
176: fcmpx fp1,fp3; fjne L26
177: movl #1,d7
178: jra L27
179: L26:
180: clrl d7
181: L27:
182: movl d7,sp@(4)
183: movl #L25,sp@(8)
184: jbsr _TstCond:l
185: fmoved _MinusOne:l,fp3
186: fmoved L10:l,fp2
187: fmoved _One:l,fp1
188: fsubx fp1,fp2
189: fcmpx fp2,fp3; fjne L30
190: fmovex fp3,fp2; faddx fp1,fp2
191: fmoved _Zero:l,fp0
192: fcmpx fp0,fp2; fjne L30
193: fmovex fp1,fp2; faddx fp3,fp2
194: fcmpx fp0,fp2; fjne L30
195: fmoved fp1,sp@(0)
196: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
197: faddx fp0,fp3
198: fmoved _Zero:l,fp2
199: fcmpx fp2,fp3; fjne L30
200: fmoved _MinusOne:l,fp3
201: fmovex fp3,fp1; fmulx fp3,fp1
202: faddx fp1,fp3
203: fcmpx fp2,fp3; fjne L30
204: movl #1,d7
205: jra L31
206: L30:
207: clrl d7
208: L31:
209: movl #0,sp@
210: movl d7,sp@(4)
211: movl #L29,sp@(8)
212: jbsr _TstCond:l
213: movl #0,sp@
214: fmoved _Half:l,fp3
215: fmovex fp3,fp2; faddd _MinusOne:l,fp2
216: faddx fp3,fp2
217: fcmpd _Zero:l,fp2; fjne L34
218: movl #1,d7
219: jra L35
220: L34:
221: clrl d7
222: L35:
223: movl d7,sp@(4)
224: movl #L33,sp@(8)
225: jbsr _TstCond:l
226: movl #10,_Milestone:l
227: movl #0,sp@
228: fmoved _Nine:l,fp3
229: fmoved _Three:l,fp2
230: fmovex fp2,fp1; fmulx fp2,fp1
231: fcmpx fp1,fp3; fjne L38
232: fmoved _TwentySeven:l,fp1
233: fmulx fp2,fp3
234: fcmpx fp3,fp1; fjne L38
235: fmoved _Eight:l,fp3
236: fmoved _Four:l,fp2
237: fmovex fp2,fp0; faddx fp2,fp0
238: fcmpx fp0,fp3; fjne L38
239: fmoved _ThirtyTwo:l,fp0
240: fmulx fp2,fp3
241: fcmpx fp3,fp0; fjne L38
242: fmovex fp0,fp3; fsubx fp1,fp3
243: fsubx fp2,fp3
244: fsubd _One:l,fp3
245: fcmpd _Zero:l,fp3; fjne L38
246: movl #1,d7
247: jra L39
248: L38:
249: clrl d7
250: L39:
251: movl d7,sp@(4)
252: movl #L37,sp@(8)
253: jbsr _TstCond:l
254: movl #0,sp@
255: fmoved _Five:l,fp3
256: fmovex fp3,fp7
257: fmoved _Four:l,fp2
258: fmovex fp2,fp1; faddd _One:l,fp1
259: fcmpx fp1,fp3; fjne L42
260: fmoved _TwoForty:l,fp1
261: fmovex fp1,fp5
262: fmovex fp2,fp0; fmulx fp3,fp0
263: fmoved _Three:l,fp3
264: fmovex fp3,fp6
265: fmulx fp3,fp0
266: fmulx fp2,fp0
267: fcmpx fp0,fp1; fjne L42
268: fmovex fp1,fp0; fdivx fp3,fp0
269: fmovex fp2,fp3; fmulx fp2,fp3
270: fmulx fp7,fp3
271: fsubx fp3,fp0
272: fmoved _Zero:l,fp3
273: fcmpx fp3,fp0; fjne L42
274: fmovex fp1,fp0; fdivx fp2,fp0
275: fmovex fp7,fp1; fmulx fp6,fp1
276: fmulx fp2,fp1
277: fsubx fp1,fp0
278: fcmpx fp3,fp0; fjne L42
279: fmovex fp5,fp1; fdivx fp7,fp1
280: fmovex fp2,fp0; fmulx fp6,fp0
281: fmulx fp2,fp0
282: fmovex fp1,fp2; fsubx fp0,fp2
283: fcmpx fp3,fp2; fjne L42
284: movl #1,d7
285: jra L43
286: L42:
287: clrl d7
288: L43:
289: movl d7,sp@(4)
290: movl #L41,sp@(8)
291: jbsr _TstCond:l
292: tstl _ErrCnt:l; jne L44
293: movl #L46,sp@
294: jbsr _printf:l
295: movl #L47,sp@
296: jbsr _printf:l
297: L44:
298: movl #L48,sp@
299: jbsr _printf:l
300: fmoved _One:l,fp3
301: fmoved fp3,_W:l
302: L49:
303: fmoved _W:l,fp3
304: fmovex fp3,fp2; faddx fp3,fp2
305: fmoved fp2,_W:l
306: fmoved _W:l,fp3
307: fmoved _One:l,fp2
308: fmovex fp3,fp1; faddx fp2,fp1
309: fmoved fp1,_Y:l
310: fmoved _Y:l,fp1
311: fsubx fp3,fp1
312: fmoved fp1,_Z:l
313: fmoved _Z:l,fp3
314: fsubx fp2,fp3
315: fmoved fp3,_Y:l
316: L50:
317: fmoved _Y:l,fp3
318: fmoved fp3,sp@(0)
319: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
320: fmoved _MinusOne:l,fp3
321: faddx fp0,fp3
322: fcmpd _Zero:l,fp3; fjlt L49
323: fmoved _Zero:l,fp3
324: fmoved fp3,_Precision:l
325: fmoved _One:l,fp3
326: fmoved fp3,_Y:l
327: L52:
328: fmoved _W:l,fp3
329: fmoved _Y:l,fp2
330: fmovex fp3,fp1; faddx fp2,fp1
331: fmoved fp1,_Radix:l
332: fmovex fp2,fp1; faddx fp2,fp1
333: fmoved fp1,_Y:l
334: fmoved _Radix:l,fp2
335: fsubx fp3,fp2
336: fmoved fp2,_Radix:l
337: L53:
338: fmoved _Radix:l,fp3
339: fcmpd _Zero:l,fp3; fjeq L52
340: fmoved _Radix:l,fp3
341: fcmpd _Two:l,fp3; fjge L55
342: fmoved _One:l,fp3
343: fmoved fp3,_Radix:l
344: L55:
345: movl #L57,sp@
346: fmoved _Radix:l,fp3
347: fmoved fp3,sp@(4)
348: jbsr _printf:l
349: fmoved _Radix:l,fp3
350: fcmpd L11:l,fp3; fjeq L58
351: fmoved _One:l,fp3
352: fmoved fp3,_W:l
353: L60:
354: fmoved _Precision:l,fp3
355: fmoved _One:l,fp2
356: faddx fp2,fp3
357: fmoved fp3,_Precision:l
358: fmoved _W:l,fp3
359: fmuld _Radix:l,fp3
360: fmoved fp3,_W:l
361: fmoved _W:l,fp3
362: faddx fp2,fp3
363: fmoved fp3,_Y:l
364: L61:
365: fmoved _Y:l,fp3
366: fsubd _W:l,fp3
367: fcmpd _One:l,fp3; fjeq L60
368: L58:
369: fmoved _One:l,fp3
370: fdivd _W:l,fp3
371: fmoved fp3,_U1:l
372: fmoved _Radix:l,fp3
373: fmoved _U1:l,fp2
374: fmulx fp2,fp3
375: fmoved fp3,_U2:l
376: movl #L63,sp@
377: fmoved fp2,sp@(4)
378: jbsr _printf:l
379: movl #L64,sp@
380: jbsr _printf:l
381: fmoved _Radix:l,fp3
382: fmoved fp3,_E0:l
383: fmoved _U1:l,fp3
384: fmoved fp3,_E1:l
385: fmoved _U2:l,fp3
386: fmoved fp3,_E9:l
387: fmoved _Precision:l,fp3
388: fmoved fp3,_E3:l
389: fmoved _Four:l,fp3
390: fdivd _Three:l,fp3
391: fmoved fp3,_X:l
392: fmoved _X:l,fp3
393: fsubd _One:l,fp3
394: fmoved fp3,_Third:l
395: fmoved _Half:l,fp3
396: fmoved _Third:l,fp2
397: fsubx fp2,fp3
398: fmoved fp3,_F6:l
399: fmoved _F6:l,fp3
400: fmovex fp3,fp1; faddx fp3,fp1
401: fmoved fp1,_X:l
402: fmoved _X:l,fp3
403: fsubx fp2,fp3
404: fmoved fp3,sp@(0)
405: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
406: fmoved fp0,_X:l
407: fmoved _X:l,fp3
408: fcmpd _U2:l,fp3; fjge L65
409: fmoved _U2:l,fp3
410: fmoved fp3,_X:l
411: L65:
412: L67:
413: fmoved _X:l,fp3
414: fmoved fp3,_U2:l
415: fmoved _Half:l,fp3
416: fmoved _U2:l,fp2
417: fmulx fp2,fp3
418: fmoved _ThirtyTwo:l,fp1
419: fmulx fp2,fp1
420: fmulx fp2,fp1
421: faddx fp1,fp3
422: fmoved fp3,_Y:l
423: fmoved _One:l,fp3
424: fmovex fp3,fp2; faddd _Y:l,fp2
425: fmoved fp2,_Y:l
426: fmoved _Y:l,fp2
427: fsubx fp3,fp2
428: fmoved fp2,_X:l
429: L68:
430: fmoved _U2:l,fp3
431: fmoved _X:l,fp2
432: fcmpx fp2,fp3; fjle L70
433: fcmpd _Zero:l,fp2; fjgt L67
434: L70:
435: fmoved _Two:l,fp3
436: fdivd _Three:l,fp3
437: fmoved fp3,_X:l
438: fmoved _X:l,fp3
439: fmoved _Half:l,fp2
440: fsubx fp2,fp3
441: fmoved fp3,_F6:l
442: fmoved _F6:l,fp3
443: fmovex fp3,fp1; faddx fp3,fp1
444: fmoved fp1,_Third:l
445: fmoved _Third:l,fp1
446: fsubx fp2,fp1
447: fmoved fp1,_X:l
448: fmoved _X:l,fp2
449: faddx fp3,fp2
450: fmoved fp2,sp@(0)
451: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
452: fmoved fp0,_X:l
453: fmoved _X:l,fp3
454: fcmpd _U1:l,fp3; fjge L71
455: fmoved _U1:l,fp3
456: fmoved fp3,_X:l
457: L71:
458: L73:
459: fmoved _X:l,fp3
460: fmoved fp3,_U1:l
461: fmoved _Half:l,fp3
462: fmoved _U1:l,fp2
463: fmovex fp3,fp1; fmulx fp2,fp1
464: fmoved _ThirtyTwo:l,fp0
465: fmulx fp2,fp0
466: fmulx fp2,fp0
467: fmovex fp1,fp2; faddx fp0,fp2
468: fmoved fp2,_Y:l
469: fmovex fp3,fp2; fsubd _Y:l,fp2
470: fmoved fp2,_Y:l
471: fmovex fp3,fp2; faddd _Y:l,fp2
472: fmoved fp2,_X:l
473: fmovex fp3,fp2; fsubd _X:l,fp2
474: fmoved fp2,_Y:l
475: faddd _Y:l,fp3
476: fmoved fp3,_X:l
477: L74:
478: fmoved _U1:l,fp3
479: fmoved _X:l,fp2
480: fcmpx fp2,fp3; fjle L76
481: fcmpd _Zero:l,fp2; fjgt L73
482: L76:
483: fmoved _U1:l,fp3
484: fcmpd _E1:l,fp3; fjne L77
485: movl #L79,sp@
486: jbsr _printf:l
487: jra L78
488: L77:
489: movl #L80,sp@
490: fmoved _U1:l,fp3
491: fmoved fp3,sp@(4)
492: jbsr _printf:l
493: L78:
494: fmoved _One:l,fp3
495: fmoved _U1:l,fp2
496: fdivx fp2,fp3
497: fmoved fp3,_W:l
498: fmoved _Half:l,fp3
499: fmovex fp3,fp1; fsubx fp2,fp1
500: faddx fp3,fp1
501: fmoved fp1,_F9:l
502: fmoved _U2:l,fp3
503: fdivx fp2,fp3
504: faddd L81:l,fp3
505: fmoved fp3,sp@(0)
506: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
507: fmoved fp0,_Radix:l
508: fmoved _Radix:l,fp3
509: fcmpd _E0:l,fp3; fjne L82
510: movl #L84,sp@
511: jbsr _printf:l
512: jra L83
513: L82:
514: movl #L85,sp@
515: fmoved _Radix:l,fp3
516: fmoved fp3,sp@(4)
517: jbsr _printf:l
518: L83:
519: movl #2,sp@
520: fmoved _Radix:l,fp3
521: fmoved _Eight:l,fp2
522: fmovex fp2,fp1; faddx fp2,fp1
523: fcmpx fp1,fp3; fjgt L88
524: movl #1,d7
525: jra L89
526: L88:
527: clrl d7
528: L89:
529: movl d7,sp@(4)
530: movl #L87,sp@(8)
531: jbsr _TstCond:l
532: movl #3,sp@
533: fmoved _Radix:l,fp3
534: fcmpd _Two:l,fp3; fjeq L95
535: fcmpd L96:l,fp3; fjeq L95
536: fcmpd _One:l,fp3; fjne L92
537: L95:
538: movl #1,d7
539: jra L93
540: L92:
541: clrl d7
542: L93:
543: movl d7,sp@(4)
544: movl #L91,sp@(8)
545: jbsr _TstCond:l
546: movl #20,_Milestone:l
547: movl #0,sp@
548: fmoved _F9:l,fp3
549: fmoved _Half:l,fp2
550: fsubx fp2,fp3
551: fcmpx fp2,fp3; fjge L99
552: movl #1,d7
553: jra L100
554: L99:
555: clrl d7
556: L100:
557: movl d7,sp@(4)
558: movl #L98,sp@(8)
559: jbsr _TstCond:l
560: fmoved _F9:l,fp3
561: fmoved fp3,_X:l
562: movl #1,_I:l
563: fmoved _X:l,fp3
564: fmoved _Half:l,fp2
565: fsubx fp2,fp3
566: fmoved fp3,_Y:l
567: fmoved _Y:l,fp3
568: fsubx fp2,fp3
569: fmoved fp3,_Z:l
570: movl #0,sp@
571: fmoved _X:l,fp3
572: fcmpd _One:l,fp3; fjne L105
573: fmoved _Z:l,fp3
574: fcmpd _Zero:l,fp3; fjne L103
575: L105:
576: movl #1,d7
577: jra L104
578: L103:
579: clrl d7
580: L104:
581: movl d7,sp@(4)
582: movl #L102,sp@(8)
583: jbsr _TstCond:l
584: fmoved _One:l,fp3
585: fmoved _U2:l,fp2
586: fmovex fp3,fp1; faddx fp2,fp1
587: fmoved fp1,_X:l
588: clrl _I:l
589: movl #25,_Milestone:l
590: fmoved _Radix:l,fp1
591: fmovex fp1,fp0; fsubx fp3,fp0
592: fmoved fp0,_BMinusU2:l
593: fmoved _BMinusU2:l,fp0
594: fsubx fp2,fp0
595: fmovex fp0,fp2; faddx fp3,fp2
596: fmoved fp2,_BMinusU2:l
597: fcmpx fp3,fp1; fjeq L106
598: fmoved _U1:l,fp3
599: fmoved fp3,sp@(0)
600: jbsr _log:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
601: fmovex fp0,fp7
602: fmoved _Radix:l,fp3
603: fmoved fp3,sp@(0)
604: jbsr _log:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
605: fnegd _TwoForty:l,fp3
606: fmulx fp7,fp3
607: fdivx fp0,fp3
608: fmoved fp3,_X:l
609: fmoved _Half:l,fp3
610: faddd _X:l,fp3
611: fmoved fp3,sp@(0)
612: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
613: fmoved fp0,_Y:l
614: fmoved _X:l,fp3
615: fsubd _Y:l,fp3
616: fmoved fp3,sp@(0)
617: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
618: fmovex fp0,fp3; fmuld _Four:l,fp3
619: fcmpd _One:l,fp3; fjge L108
620: fmoved _Y:l,fp3
621: fmoved fp3,_X:l
622: L108:
623: fmoved _X:l,fp3
624: fdivd _TwoForty:l,fp3
625: fmoved fp3,_Precision:l
626: fmoved _Half:l,fp3
627: faddd _Precision:l,fp3
628: fmoved fp3,sp@(0)
629: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
630: fmoved fp0,_Y:l
631: fmoved _Precision:l,fp3
632: fsubd _Y:l,fp3
633: fmoved fp3,sp@(0)
634: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
635: fmovex fp0,fp3; fmuld _TwoForty:l,fp3
636: fcmpd _Half:l,fp3; fjge L110
637: fmoved _Y:l,fp3
638: fmoved fp3,_Precision:l
639: L110:
640: L106:
641: fmoved _Precision:l,fp3
642: fmoved fp3,sp@(0)
643: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
644: fcmpx fp0,fp3; fjne L114
645: fmoved _Radix:l,fp3
646: fcmpd _One:l,fp3; fjne L112
647: L114:
648: movl #L115,sp@
649: jbsr _printf:l
650: movl #L116,sp@
651: jbsr _printf:l
652: L112:
653: fmoved _Radix:l,fp3
654: fcmpd _One:l,fp3; fjne L117
655: movl #L119,sp@
656: jbsr _printf:l
657: jra L118
658: L117:
659: movl #L120,sp@
660: fmoved _Precision:l,fp3
661: fmoved fp3,sp@(4)
662: jbsr _printf:l
663: L118:
664: movl #1,sp@
665: fmoved _U2:l,fp3
666: fmoved _Nine:l,fp2
667: fmulx fp2,fp3
668: fmulx fp2,fp3
669: fmuld _TwoForty:l,fp3
670: fcmpd _One:l,fp3; fjge L123
671: movl #1,d7
672: jra L124
673: L123:
674: clrl d7
675: L124:
676: movl d7,sp@(4)
677: movl #L122,sp@(8)
678: jbsr _TstCond:l
679: movl #30,_Milestone:l
680: fmoved _Four:l,fp3
681: fmoved _Three:l,fp2
682: fmovex fp3,fp1; fdivx fp2,fp1
683: fmoved _One:l,fp0
684: fsubx fp0,fp1
685: fdivx fp3,fp0
686: fmovex fp1,fp3; fsubx fp0,fp3
687: fmulx fp2,fp3
688: fsubx fp0,fp3
689: fmoved fp3,sp@(0)
690: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
691: fmoved fp0,_X:l
692: L125:
693: fmoved _X:l,fp3
694: fmoved fp3,_Z2:l
695: fmoved _One:l,fp3
696: fmoved _Half:l,fp2
697: fmoved _Z2:l,fp1
698: fmulx fp1,fp2
699: fmoved _ThirtyTwo:l,fp0
700: fmulx fp1,fp0
701: fmulx fp1,fp0
702: faddx fp0,fp2
703: fmovex fp3,fp1; faddx fp2,fp1
704: fmovex fp1,fp2; fsubx fp3,fp2
705: fmoved fp2,_X:l
706: L126:
707: fmoved _Z2:l,fp3
708: fmoved _X:l,fp2
709: fcmpx fp2,fp3; fjle L128
710: fcmpd _Zero:l,fp2; fjgt L125
711: L128:
712: fmoved _Three:l,fp3
713: fmoved _Four:l,fp2
714: fmovex fp3,fp1; fdivx fp2,fp1
715: fmoved _Two:l,fp0
716: fdivx fp3,fp0
717: fsubx fp0,fp1
718: fmulx fp3,fp1
719: fmoved _One:l,fp3
720: fdivx fp2,fp3
721: fmovex fp1,fp2; fsubx fp3,fp2
722: fmoved fp2,sp@(0)
723: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
724: fmoved fp0,_Z:l
725: fmoved fp0,_Y:l
726: fmoved fp0,_X:l
727: L129:
728: fmoved _Z:l,fp3
729: fmoved fp3,_Z1:l
730: fmoved _One:l,fp3
731: fdivd _Two:l,fp3
732: fmoved _Half:l,fp2
733: fmoved _Z1:l,fp1
734: fmulx fp1,fp2
735: fmoved _ThirtyTwo:l,fp0
736: fmulx fp1,fp0
737: fmulx fp1,fp0
738: faddx fp0,fp2
739: fmovex fp3,fp1; fsubx fp2,fp1
740: fmovex fp1,fp2; faddx fp3,fp2
741: fmovex fp3,fp1; fsubx fp2,fp1
742: fmovex fp1,fp2; faddx fp3,fp2
743: fmoved fp2,_Z:l
744: L130:
745: fmoved _Z1:l,fp3
746: fmoved _Z:l,fp2
747: fcmpx fp2,fp3; fjle L132
748: fcmpd _Zero:l,fp2; fjgt L129
749: L132:
750: L133:
751: L136:
752: fmoved _Y:l,fp3
753: fmoved fp3,_Y1:l
754: fmoved _Half:l,fp3
755: fmoved _Y1:l,fp2
756: fmovex fp3,fp1; fmulx fp2,fp1
757: fmoved _ThirtyTwo:l,fp0
758: fmulx fp2,fp0
759: fmulx fp2,fp0
760: fmovex fp1,fp2; faddx fp0,fp2
761: fmovex fp3,fp1; fsubx fp2,fp1
762: fmovex fp1,fp2; faddx fp3,fp2
763: fmovex fp3,fp1; fsubx fp2,fp1
764: fmovex fp1,fp2; faddx fp3,fp2
765: fmoved fp2,_Y:l
766: L137:
767: fmoved _Y1:l,fp3
768: fmoved _Y:l,fp2
769: fcmpx fp2,fp3; fjle L139
770: fcmpd _Zero:l,fp2; fjgt L136
771: L139:
772: fmoved _X:l,fp3
773: fmoved fp3,_X1:l
774: fmoved _Half:l,fp3
775: fmoved _X1:l,fp2
776: fmulx fp2,fp3
777: fmoved _ThirtyTwo:l,fp1
778: fmulx fp2,fp1
779: fmulx fp2,fp1
780: faddx fp1,fp3
781: fmoved _F9:l,fp2
782: fsubx fp2,fp3
783: faddx fp2,fp3
784: fmoved fp3,_X:l
785: L134:
786: fmoved _X1:l,fp3
787: fmoved _X:l,fp2
788: fcmpx fp2,fp3; fjle L140
789: fcmpd _Zero:l,fp2; fjgt L133
790: L140:
791: fmoved _X1:l,fp3
792: fcmpd _Y1:l,fp3; fjne L143
793: fcmpd _Z1:l,fp3; fjeq L141
794: L143:
795: movl #1,sp@
796: movl #L144,sp@(4)
797: jbsr _BadCond:l
798: movl #L145,sp@
799: fmoved _X1:l,fp3
800: fmoved fp3,sp@(4)
801: fmoved _Y1:l,fp3
802: fmoved fp3,sp@(12)
803: fmoved _Z1:l,fp3
804: fmoved fp3,sp@(20)
805: jbsr _printf:l
806: movl #L146,sp@
807: jbsr _printf:l
808: movl #L147,sp@
809: jbsr _printf:l
810: movl #L148,sp@
811: jbsr _notify:l
812: fmoved _X1:l,fp3
813: fmoved _U1:l,fp2
814: fcmpx fp2,fp3; fjeq L152
815: fmoved _Y1:l,fp3
816: fcmpx fp2,fp3; fjeq L152
817: fmoved _Z1:l,fp3
818: fcmpx fp2,fp3; fjne L142
819: L152:
820: movl #L153,sp@
821: jbsr _printf:l
822: jra L142
823: L141:
824: fmoved _Z1:l,fp3
825: fcmpd _U1:l,fp3; fjne L156
826: fmoved _Z2:l,fp3
827: fcmpd _U2:l,fp3; fjeq L154
828: L156:
829: fmoved _Z1:l,fp3
830: fcmpd _U1:l,fp3; fjge L159
831: fmoved _Z2:l,fp3
832: fcmpd _U2:l,fp3; fjlt L157
833: L159:
834: movl #0,sp@
835: movl #L160,sp@(4)
836: jbsr _BadCond:l
837: movl #L161,sp@
838: jbsr _notify:l
839: movl #L162,sp@
840: fmoved _U1:l,fp3
841: fmoved fp3,sp@(4)
842: fmoved _Z1:l,fp2
843: fsubx fp3,fp2
844: fmoved fp2,sp@(12)
845: jbsr _printf:l
846: movl #L163,sp@
847: fmoved _U2:l,fp3
848: fmoved fp3,sp@(4)
849: fmoved _Z2:l,fp2
850: fsubx fp3,fp2
851: fmoved fp2,sp@(12)
852: jbsr _printf:l
853: jra L158
854: L157:
855: fmoved _Z1:l,fp3
856: fmoved _Zero:l,fp2
857: fcmpx fp2,fp3; fjle L166
858: fmoved _Z2:l,fp3
859: fcmpx fp2,fp3; fjgt L164
860: L166:
861: movl #L167,sp@
862: fmoved _Radix:l,fp3
863: fmoved fp3,sp@(4)
864: jbsr _printf:l
865: movl #L168,sp@
866: jbsr _printf:l
867: movl #L169,sp@
868: fmoved _Z1:l,fp3
869: fmoved fp3,sp@(4)
870: fmoved _Z2:l,fp3
871: fmoved fp3,sp@(12)
872: jbsr _printf:l
873: movl #L170,sp@
874: jbsr _notify:l
875: L164:
876: fmoved _Z1:l,fp3
877: fcmpd _Z2:l,fp3; fjne L173
878: fcmpd _Zero:l,fp3; fjle L171
879: L173:
880: fmoved _Z1:l,fp3
881: fdivd _U1:l,fp3
882: fmoved fp3,_X:l
883: fmoved _Z2:l,fp3
884: fdivd _U2:l,fp3
885: fmoved fp3,_Y:l
886: fmoved _Y:l,fp3
887: fcmpd _X:l,fp3; fjle L174
888: fmoved _Y:l,fp3
889: fmoved fp3,_X:l
890: L174:
891: fmoved _X:l,fp3
892: fmoved fp3,sp@(0)
893: jbsr _log:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
894: fnegx fp0,fp3
895: fmoved fp3,_Q:l
896: movl #L176,sp@
897: jbsr _printf:l
898: fmoved _Radix:l,fp3
899: fmoved fp3,sp@(0)
900: jbsr _log:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
901: movl #L177,sp@
902: fmoved _Q:l,fp3
903: fdivx fp0,fp3
904: fmoved fp3,sp@(4)
905: jbsr _printf:l
906: fmoved L96:l,fp3
907: fmoved fp3,sp@(0)
908: jbsr _log:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
909: movl #L178,sp@
910: fmoved _Q:l,fp3
911: fdivx fp0,fp3
912: fmoved fp3,sp@(4)
913: jbsr _printf:l
914: L171:
915: movl #L153,sp@
916: jbsr _printf:l
917: L158:
918: L154:
919: L142:
920: jbsr _Pause:l
921: movl #35,_Milestone:l
922: fmoved _Radix:l,fp3
923: fcmpd _Two:l,fp3; fjlt L179
924: fmoved _W:l,fp3
925: fmoved _Radix:l,fp2
926: fmovex fp2,fp1; fmulx fp2,fp1
927: fdivx fp1,fp3
928: fmoved fp3,_X:l
929: fmoved _X:l,fp3
930: fmovex fp3,fp2; faddd _One:l,fp2
931: fmoved fp2,_Y:l
932: fmoved _Y:l,fp2
933: fsubx fp3,fp2
934: fmoved fp2,_Z:l
935: fmoved _Z:l,fp3
936: fmovex fp3,fp2; faddd _U2:l,fp2
937: fmoved fp2,_T:l
938: fmoved _T:l,fp2
939: fsubx fp3,fp2
940: fmoved fp2,_X:l
941: movl #0,sp@
942: fmoved _X:l,fp3
943: fcmpd _U2:l,fp3; fjne L183
944: movl #1,d7
945: jra L184
946: L183:
947: clrl d7
948: L184:
949: movl d7,sp@(4)
950: movl #L182,sp@(8)
951: jbsr _TstCond:l
952: fmoved _X:l,fp3
953: fcmpd _U2:l,fp3; fjne L185
954: movl #L187,sp@
955: jbsr _printf:l
956: L185:
957: L179:
958: movl #L188,sp@
959: jbsr _printf:l
960: fmoved _F9:l,fp3
961: fmoved _One:l,fp2
962: fmovex fp2,fp7
963: fmovex fp3,fp1; fmulx fp2,fp1
964: fmoved fp1,_Y:l
965: fmovex fp2,fp1; fmulx fp3,fp1
966: fmoved fp1,_Z:l
967: fmoved _Half:l,fp1
968: fsubx fp1,fp3
969: fmoved fp3,_X:l
970: fmoved _Y:l,fp3
971: fsubx fp1,fp3
972: fmoved _X:l,fp0
973: fsubx fp0,fp3
974: fmoved fp3,_Y:l
975: fmoved _Z:l,fp3
976: fsubx fp1,fp3
977: fsubx fp0,fp3
978: fmoved fp3,_Z:l
979: fmoved _U2:l,fp3
980: fmovex fp2,fp1; faddx fp3,fp1
981: fmoved fp1,_X:l
982: fmoved _X:l,fp1
983: fmoved _Radix:l,fp0
984: fmovex fp1,fp2; fmulx fp0,fp2
985: fmoved fp2,_T:l
986: fmovex fp0,fp2; fmulx fp1,fp2
987: fmoved fp2,_R:l
988: fmoved _T:l,fp2
989: fsubx fp0,fp2
990: fmoved fp2,_X:l
991: fmoved _X:l,fp2
992: fmovex fp0,fp1; fmulx fp3,fp1
993: fmovex fp2,fp3; fsubx fp1,fp3
994: fmoved fp3,_X:l
995: fmoved _R:l,fp3
996: fsubx fp0,fp3
997: fmoved fp3,_T:l
998: fmoved _T:l,fp3
999: fsubx fp1,fp3
1000: fmoved fp3,_T:l
1001: fmoved _X:l,fp3
1002: fmovex fp0,fp2; fsubx fp7,fp2
1003: fmulx fp2,fp3
1004: fmoved fp3,_X:l
1005: fmoved _T:l,fp3
1006: fmulx fp2,fp3
1007: fmoved fp3,_T:l
1008: fmoved _X:l,fp3
1009: fmoved _Zero:l,fp2
1010: fcmpx fp2,fp3; fjne L189
1011: fmoved _Y:l,fp3
1012: fcmpx fp2,fp3; fjne L189
1013: fmoved _Z:l,fp3
1014: fcmpx fp2,fp3; fjne L189
1015: fmoved _T:l,fp3
1016: fcmpx fp2,fp3; fjne L189
1017: movl #1,_GMult:l
1018: jra L190
1019: L189:
1020: clrl _GMult:l
1021: movl #1,sp@
1022: clrl sp@(4)
1023: movl #L191,sp@(8)
1024: jbsr _TstCond:l
1025: L190:
1026: fmoved _Radix:l,fp3
1027: fmuld _U2:l,fp3
1028: fmoved fp3,_Z:l
1029: fmoved _One:l,fp3
1030: fmoved _Z:l,fp2
1031: faddx fp2,fp3
1032: fmoved fp3,_X:l
1033: fmoved _X:l,fp3
1034: fmovex fp3,fp1; faddx fp2,fp1
1035: fmovex fp3,fp2; fmulx fp3,fp2
1036: fmovex fp1,fp3; fsubx fp2,fp3
1037: fmoved fp3,sp@(0)
1038: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
1039: fmoved _U2:l,fp3
1040: fmovex fp0,fp2; fsubx fp3,fp2
1041: fmoved fp2,_Y:l
1042: fmoved _One:l,fp2
1043: fsubx fp3,fp2
1044: fmoved fp2,_X:l
1045: fmoved _X:l,fp2
1046: fmovex fp2,fp1; fsubx fp3,fp1
1047: fmovex fp2,fp3; fmulx fp2,fp3
1048: fmovex fp1,fp2; fsubx fp3,fp2
1049: fmoved fp2,sp@(0)
1050: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
1051: fmovex fp0,fp3; fsubd _U1:l,fp3
1052: fmoved fp3,_Z:l
1053: movl #0,sp@
1054: fmoved _Y:l,fp3
1055: fmoved _Zero:l,fp2
1056: fcmpx fp2,fp3; fjgt L194
1057: fmoved _Z:l,fp3
1058: fcmpx fp2,fp3; fjgt L194
1059: movl #1,d7
1060: jra L195
1061: L194:
1062: clrl d7
1063: L195:
1064: movl d7,sp@(4)
1065: movl #L193,sp@(8)
1066: jbsr _TstCond:l
1067: fmoved _One:l,fp3
1068: fmoved _U2:l,fp2
1069: fmovex fp3,fp1; fsubx fp2,fp1
1070: fmoved fp1,_Y:l
1071: fmovex fp3,fp1; faddx fp2,fp1
1072: fmoved fp1,_X:l
1073: fmovex fp3,fp2; fdivd _Y:l,fp2
1074: fmoved fp2,_Z:l
1075: fmoved _Z:l,fp2
1076: fsubd _X:l,fp2
1077: fmoved fp2,_Y:l
1078: fmoved _Three:l,fp2
1079: fdivx fp2,fp3
1080: fmoved fp3,_X:l
1081: fmoved _Nine:l,fp3
1082: fdivx fp3,fp2
1083: fmoved fp2,_Z:l
1084: fmoved _X:l,fp2
1085: fmoved _Z:l,fp1
1086: fsubx fp1,fp2
1087: fmoved fp2,_X:l
1088: fdivd _TwentySeven:l,fp3
1089: fmoved fp3,_T:l
1090: fmovex fp1,fp3; fsubd _T:l,fp3
1091: fmoved fp3,_Z:l
1092: movl #2,sp@
1093: fmoved _X:l,fp3
1094: fmoved _Zero:l,fp2
1095: fcmpx fp2,fp3; fjne L198
1096: fmoved _Y:l,fp3
1097: fcmpx fp2,fp3; fjne L198
1098: fmoved _Z:l,fp3
1099: fcmpx fp2,fp3; fjne L198
1100: movl #1,d7
1101: jra L199
1102: L198:
1103: clrl d7
1104: L199:
1105: movl d7,sp@(4)
1106: movl #L197,sp@(8)
1107: jbsr _TstCond:l
1108: fmoved _F9:l,fp3
1109: fmoved _One:l,fp2
1110: fmovex fp3,fp1; fdivx fp2,fp1
1111: fmoved fp1,_Y:l
1112: fmoved _Half:l,fp1
1113: fsubx fp1,fp3
1114: fmoved fp3,_X:l
1115: fmoved _Y:l,fp3
1116: fsubx fp1,fp3
1117: fsubd _X:l,fp3
1118: fmoved fp3,_Y:l
1119: fmovex fp2,fp3; faddd _U2:l,fp3
1120: fmoved fp3,_X:l
1121: fmoved _X:l,fp3
1122: fmovex fp3,fp1; fdivx fp2,fp1
1123: fmoved fp1,_T:l
1124: fmoved _T:l,fp2
1125: fsubx fp3,fp2
1126: fmoved fp2,_X:l
1127: fmoved _X:l,fp3
1128: fmoved _Zero:l,fp2
1129: fcmpx fp2,fp3; fjne L200
1130: fmoved _Y:l,fp3
1131: fcmpx fp2,fp3; fjne L200
1132: fmoved _Z:l,fp3
1133: fcmpx fp2,fp3; fjne L200
1134: movl #1,_GDiv:l
1135: jra L201
1136: L200:
1137: clrl _GDiv:l
1138: movl #1,sp@
1139: clrl sp@(4)
1140: movl #L202,sp@(8)
1141: jbsr _TstCond:l
1142: L201:
1143: fmoved _One:l,fp3
1144: fmovex fp3,fp2; faddd _U2:l,fp2
1145: fdivx fp2,fp3
1146: fmoved fp3,_X:l
1147: fmoved _X:l,fp3
1148: fmoved _Half:l,fp2
1149: fsubx fp2,fp3
1150: fsubx fp2,fp3
1151: fmoved fp3,_Y:l
1152: movl #1,sp@
1153: fmoved _Y:l,fp3
1154: fcmpd _Zero:l,fp3; fjge L205
1155: movl #1,d7
1156: jra L206
1157: L205:
1158: clrl d7
1159: L206:
1160: movl d7,sp@(4)
1161: movl #L204,sp@(8)
1162: jbsr _TstCond:l
1163: fmoved _One:l,fp3
1164: fmoved _U2:l,fp2
1165: fmovex fp3,fp1; fsubx fp2,fp1
1166: fmoved fp1,_X:l
1167: fmoved _Radix:l,fp1
1168: fmovex fp1,fp0; fmulx fp2,fp0
1169: faddx fp0,fp3
1170: fmoved fp3,_Y:l
1171: fmoved _X:l,fp3
1172: fmovex fp3,fp2; fmulx fp1,fp2
1173: fmoved fp2,_Z:l
1174: fmoved _Y:l,fp2
1175: fmovex fp2,fp0; fmulx fp1,fp0
1176: fmoved fp0,_T:l
1177: fmoved _Z:l,fp0
1178: fdivx fp1,fp0
1179: fmoved fp0,_R:l
1180: fmoved _T:l,fp0
1181: fdivx fp1,fp0
1182: fmoved fp0,_StickyBit:l
1183: fmoved _R:l,fp1
1184: fsubx fp3,fp1
1185: fmoved fp1,_X:l
1186: fmoved _StickyBit:l,fp3
1187: fsubx fp2,fp3
1188: fmoved fp3,_Y:l
1189: movl #0,sp@
1190: fmoved _X:l,fp3
1191: fmoved _Zero:l,fp2
1192: fcmpx fp2,fp3; fjne L209
1193: fmoved _Y:l,fp3
1194: fcmpx fp2,fp3; fjne L209
1195: movl #1,d7
1196: jra L210
1197: L209:
1198: clrl d7
1199: L210:
1200: movl d7,sp@(4)
1201: movl #L208,sp@(8)
1202: jbsr _TstCond:l
1203: fmoved _One:l,fp3
1204: fmoved _U1:l,fp2
1205: fmovex fp3,fp1; fsubx fp2,fp1
1206: fmoved fp1,_Y:l
1207: fmovex fp3,fp1; fsubd _F9:l,fp1
1208: fmoved fp1,_X:l
1209: fsubd _Y:l,fp3
1210: fmoved fp3,_Y:l
1211: fmoved _Radix:l,fp3
1212: fmoved _U2:l,fp1
1213: fmovex fp3,fp0; fsubx fp1,fp0
1214: fmoved fp0,_T:l
1215: fmovex fp3,fp0; fsubd _BMinusU2:l,fp0
1216: fmoved fp0,_Z:l
1217: fsubd _T:l,fp3
1218: fmoved fp3,_T:l
1219: fmoved _X:l,fp3
1220: fcmpx fp2,fp3; fjne L211
1221: fmoved _Y:l,fp3
1222: fcmpx fp2,fp3; fjne L211
1223: fmoved _Z:l,fp3
1224: fcmpx fp1,fp3; fjne L211
1225: fmoved _T:l,fp3
1226: fcmpx fp1,fp3; fjne L211
1227: movl #1,_GAddSub:l
1228: jra L212
1229: L211:
1230: clrl _GAddSub:l
1231: movl #1,sp@
1232: clrl sp@(4)
1233: movl #L213,sp@(8)
1234: jbsr _TstCond:l
1235: L212:
1236: fmoved _F9:l,fp3
1237: fmoved _One:l,fp2
1238: fcmpx fp2,fp3; fjeq L214
1239: fsubx fp2,fp3
1240: fcmpd _Zero:l,fp3; fjlt L214
1241: movl #1,sp@
1242: movl #L216,sp@(4)
1243: jbsr _BadCond:l
1244: movl #L217,sp@
1245: jbsr _printf:l
1246: movl #L218,sp@
1247: jbsr _printf:l
1248: movl #L219,sp@
1249: jbsr _printf:l
1250: L214:
1251: cmpl #1,_GMult:l; jne L220
1252: cmpl #1,_GDiv:l; jne L220
1253: cmpl #1,_GAddSub:l; jne L220
1254: movl #L222,sp@
1255: jbsr _printf:l
1256: L220:
1257: movl #40,_Milestone:l
1258: jbsr _Pause:l
1259: movl #L223,sp@
1260: jbsr _printf:l
1261: clrl _RMult:l
1262: clrl _RDiv:l
1263: clrl _RAddSub:l
1264: fmoved _Radix:l,fp3
1265: fmoved _Two:l,fp2
1266: fdivx fp2,fp3
1267: fmoved fp3,_RadixD2:l
1268: fmoved fp2,_A1:l
1269: clrl _Done:l
1270: L224:
1271: fmoved _Radix:l,fp3
1272: fmoved fp3,_AInvrse:l
1273: L227:
1274: fmoved _AInvrse:l,fp3
1275: fmoved fp3,_X:l
1276: fdivd _A1:l,fp3
1277: fmoved fp3,_AInvrse:l
1278: L228:
1279: fmoved _AInvrse:l,fp3
1280: fmoved fp3,sp@(0)
1281: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
1282: fcmpd _AInvrse:l,fp0; fjeq L227
1283: fmoved _X:l,fp3
1284: fcmpd _One:l,fp3; fjeq L233
1285: fmoved _A1:l,fp3
1286: fcmpd _Three:l,fp3; fjle L231
1287: L233:
1288: movl #1,d7
1289: jra L232
1290: L231:
1291: clrl d7
1292: L232:
1293: movl d7,_Done:l
1294: tstl _Done:l; jne L234
1295: fmoved _Nine:l,fp3
1296: faddd _One:l,fp3
1297: fmoved fp3,_A1:l
1298: L234:
1299: L225:
1300: tstl _Done:l; jeq L224
1301: fmoved _X:l,fp3
1302: fcmpd _One:l,fp3; fjne L236
1303: fmoved _Radix:l,fp3
1304: fmoved fp3,_A1:l
1305: L236:
1306: fmoved _One:l,fp3
1307: fmoved _A1:l,fp2
1308: fdivx fp2,fp3
1309: fmoved fp3,_AInvrse:l
1310: fmoved fp2,_X:l
1311: fmoved _AInvrse:l,fp3
1312: fmoved fp3,_Y:l
1313: clrl _Done:l
1314: L238:
1315: fmoved _X:l,fp3
1316: fmuld _Y:l,fp3
1317: fsubd _Half:l,fp3
1318: fmoved fp3,_Z:l
1319: movl #0,sp@
1320: fmoved _Z:l,fp3
1321: fcmpd _Half:l,fp3; fjne L243
1322: movl #1,d7
1323: jra L244
1324: L243:
1325: clrl d7
1326: L244:
1327: movl d7,sp@(4)
1328: movl #L242,sp@(8)
1329: jbsr _TstCond:l
1330: fmoved _X:l,fp3
1331: fcmpd _Radix:l,fp3; fjne L246
1332: movl #1,d7
1333: jra L247
1334: L246:
1335: clrl d7
1336: L247:
1337: movl d7,_Done:l
1338: fmoved _Radix:l,fp3
1339: fmoved fp3,_X:l
1340: fmoved _One:l,fp3
1341: fdivd _X:l,fp3
1342: fmoved fp3,_Y:l
1343: L239:
1344: tstl _Done:l; jeq L238
1345: fmoved _One:l,fp3
1346: fmoved _U2:l,fp2
1347: fmovex fp2,fp6
1348: fmovex fp3,fp1; faddx fp2,fp1
1349: fmoved fp1,_Y2:l
1350: fsubx fp2,fp3
1351: fmoved fp3,_Y1:l
1352: fmoved _OneAndHalf:l,fp3
1353: fmovex fp3,fp7
1354: fmovex fp3,fp1; fsubx fp2,fp1
1355: fmoved fp1,_X:l
1356: fmovex fp3,fp1; faddx fp2,fp1
1357: fmoved fp1,_Y:l
1358: fmoved _X:l,fp1
1359: fmovex fp1,fp0; fsubx fp2,fp0
1360: fmoved _Y2:l,fp3
1361: fmovex fp3,fp5
1362: fmulx fp3,fp0
1363: fmoved fp0,_Z:l
1364: fmoved _Y:l,fp0
1365: fmoved _Y1:l,fp2
1366: fmovex fp0,fp3; fmulx fp2,fp3
1367: fmoved fp3,_T:l
1368: fmoved _Z:l,fp3
1369: fsubx fp1,fp3
1370: fmoved fp3,_Z:l
1371: fmoved _T:l,fp3
1372: fsubx fp1,fp3
1373: fmoved fp3,_T:l
1374: fmovex fp1,fp3; fmulx fp5,fp3
1375: fmoved fp3,_X:l
1376: fmovex fp0,fp3; faddx fp6,fp3
1377: fmulx fp2,fp3
1378: fmoved fp3,_Y:l
1379: fmoved _X:l,fp3
1380: fsubx fp7,fp3
1381: fmoved fp3,_X:l
1382: fmoved _Y:l,fp3
1383: fsubx fp7,fp3
1384: fmoved fp3,_Y:l
1385: fmoved _X:l,fp3
1386: fmoved _Zero:l,fp2
1387: fcmpx fp2,fp3; fjne L248
1388: fmoved _Y:l,fp3
1389: fcmpx fp2,fp3; fjne L248
1390: fmoved _Z:l,fp3
1391: fcmpx fp2,fp3; fjne L248
1392: fmoved _T:l,fp3
1393: fcmpx fp2,fp3; fjgt L248
1394: fmoved _OneAndHalf:l,fp3
1395: fmovex fp3,fp7
1396: fmoved _U2:l,fp2
1397: fmovex fp3,fp1; faddx fp2,fp1
1398: fmoved _Y2:l,fp0
1399: fmovex fp0,fp6
1400: fmovex fp1,fp3; fmulx fp0,fp3
1401: fmoved fp3,_X:l
1402: fmovex fp7,fp3; fsubx fp2,fp3
1403: fmovex fp3,fp0; fsubx fp2,fp0
1404: fmoved fp0,_Y:l
1405: faddx fp2,fp1
1406: fmoved fp1,_Z:l
1407: fmoved _Y1:l,fp1
1408: fmovex fp1,fp5
1409: fmulx fp1,fp3
1410: fmoved fp3,_T:l
1411: fmoved _X:l,fp3
1412: fmoved _Z:l,fp0
1413: fmovex fp0,fp1; faddx fp2,fp1
1414: fmovex fp1,fp4
1415: fsubx fp1,fp3
1416: fmoved fp3,_X:l
1417: fmoved _Y:l,fp3
1418: fmovex fp3,fp1; fmulx fp5,fp1
1419: fmoved fp1,_StickyBit:l
1420: fmovex fp0,fp1; fmulx fp6,fp1
1421: fmoved fp1,_S:l
1422: fmoved _T:l,fp1
1423: fsubx fp3,fp1
1424: fmoved fp1,_T:l
1425: fmovex fp2,fp1; fsubx fp3,fp1
1426: fmovex fp1,fp3; faddd _StickyBit:l,fp3
1427: fmoved fp3,_Y:l
1428: fmoved _S:l,fp3
1429: fmovex fp4,fp1; faddx fp2,fp1
1430: fsubx fp1,fp3
1431: fmoved fp3,_Z:l
1432: fmovex fp6,fp3; faddx fp2,fp3
1433: fmulx fp5,fp3
1434: fmoved fp3,_StickyBit:l
1435: fmovex fp6,fp3; fmulx fp5,fp3
1436: fmoved fp3,_Y1:l
1437: fmoved _StickyBit:l,fp3
1438: fsubx fp6,fp3
1439: fmoved fp3,_StickyBit:l
1440: fmoved _Y1:l,fp3
1441: fmoved _Half:l,fp2
1442: fsubx fp2,fp3
1443: fmoved fp3,_Y1:l
1444: fmoved _X:l,fp3
1445: fmoved _Zero:l,fp1
1446: fcmpx fp1,fp3; fjne L250
1447: fmoved _Y:l,fp3
1448: fcmpx fp1,fp3; fjne L250
1449: fmoved _Z:l,fp3
1450: fcmpx fp1,fp3; fjne L250
1451: fmoved _T:l,fp3
1452: fcmpx fp1,fp3; fjne L250
1453: fmoved _StickyBit:l,fp3
1454: fcmpx fp1,fp3; fjne L250
1455: fmoved _Y1:l,fp3
1456: fcmpx fp2,fp3; fjne L250
1457: movl #1,_RMult:l
1458: movl #L252,sp@
1459: jbsr _printf:l
1460: jra L251
1461: L250:
1462: fmoved _X:l,fp3
1463: fmoved _U2:l,fp2
1464: faddx fp2,fp3
1465: fmoved _Zero:l,fp1
1466: fcmpx fp1,fp3; fjne L253
1467: fmoved _Y:l,fp3
1468: fcmpx fp1,fp3; fjge L253
1469: fmoved _Z:l,fp3
1470: faddx fp2,fp3
1471: fcmpx fp1,fp3; fjne L253
1472: fmoved _T:l,fp3
1473: fcmpx fp1,fp3; fjge L253
1474: fmoved _StickyBit:l,fp3
1475: faddx fp2,fp3
1476: fcmpx fp1,fp3; fjne L253
1477: fmoved _Y1:l,fp3
1478: fcmpd _Half:l,fp3; fjge L253
1479: movl #2,_RMult:l
1480: movl #L255,sp@
1481: jbsr _printf:l
1482: jra L254
1483: L253:
1484: movl #L256,sp@
1485: jbsr _printf:l
1486: L254:
1487: L251:
1488: cmpl #1,_RMult:l; jne L249
1489: tstl _GMult:l; jne L249
1490: movl #L259,sp@
1491: jbsr _notify:l
1492: jra L249
1493: L248:
1494: movl #L256,sp@
1495: jbsr _printf:l
1496: L249:
1497: movl #45,_Milestone:l
1498: fmoved _One:l,fp3
1499: fmoved _U2:l,fp2
1500: fmovex fp3,fp1; faddx fp2,fp1
1501: fmoved fp1,_Y2:l
1502: fsubx fp2,fp3
1503: fmoved fp3,_Y1:l
1504: fmoved _OneAndHalf:l,fp3
1505: fmovex fp3,fp6
1506: fmovex fp3,fp1; faddx fp2,fp1
1507: fmovex fp1,fp7
1508: fmovex fp1,fp0; faddx fp2,fp0
1509: fmoved fp0,_Z:l
1510: fmoved _Z:l,fp0
1511: fmoved _Y2:l,fp1
1512: fmovex fp1,fp5
1513: fmovex fp0,fp3; fdivx fp1,fp3
1514: fmoved fp3,_X:l
1515: fmovex fp6,fp3; fsubx fp2,fp3
1516: fsubx fp2,fp3
1517: fmoved fp3,_T:l
1518: fmoved _T:l,fp3
1519: fmovex fp3,fp4
1520: fmovex fp3,fp1; fsubx fp2,fp1
1521: fmoved _Y1:l,fp3
1522: fdivx fp3,fp1
1523: fmoved fp1,_Y:l
1524: fmovex fp0,fp1; faddx fp2,fp1
1525: fdivx fp5,fp1
1526: fmoved fp1,_Z:l
1527: fmoved _X:l,fp1
1528: fsubx fp6,fp1
1529: fmoved fp1,_X:l
1530: fmoved _Y:l,fp1
1531: fsubx fp4,fp1
1532: fmoved fp1,_Y:l
1533: fmovex fp4,fp1; fdivx fp3,fp1
1534: fmoved fp1,_T:l
1535: fmoved _Z:l,fp3
1536: fsubx fp7,fp3
1537: fmoved fp3,_Z:l
1538: fmovex fp2,fp3; fsubx fp6,fp3
1539: faddd _T:l,fp3
1540: fmoved fp3,_T:l
1541: fmoved _X:l,fp3
1542: fmoved _Zero:l,fp2
1543: fcmpx fp2,fp3; fjgt L260
1544: fmoved _Y:l,fp3
1545: fcmpx fp2,fp3; fjgt L260
1546: fmoved _Z:l,fp3
1547: fcmpx fp2,fp3; fjgt L260
1548: fmoved _T:l,fp3
1549: fcmpx fp2,fp3; fjgt L260
1550: fmoved _OneAndHalf:l,fp3
1551: fmoved _Y2:l,fp2
1552: fmovex fp2,fp7
1553: fmovex fp3,fp1; fdivx fp2,fp1
1554: fmoved fp1,_X:l
1555: fmoved _U2:l,fp1
1556: fmovex fp1,fp6
1557: fmovex fp3,fp0; fsubx fp1,fp0
1558: fmoved fp0,_Y:l
1559: fmovex fp3,fp0; faddx fp1,fp0
1560: fmoved fp0,_Z:l
1561: fmoved _X:l,fp0
1562: fmoved _Y:l,fp2
1563: fsubx fp2,fp0
1564: fmoved fp0,_X:l
1565: fmoved _Y1:l,fp0
1566: fmovex fp3,fp1; fdivx fp0,fp1
1567: fmoved fp1,_T:l
1568: fdivx fp0,fp2
1569: fmoved fp2,_Y:l
1570: fmoved _T:l,fp2
1571: fmoved _Z:l,fp1
1572: fmovex fp1,fp0; faddx fp6,fp0
1573: fsubx fp0,fp2
1574: fmoved fp2,_T:l
1575: fmoved _Y:l,fp2
1576: fsubx fp1,fp2
1577: fmoved fp2,_Y:l
1578: fmovex fp1,fp2; fdivx fp7,fp2
1579: fmoved fp2,_Z:l
1580: fmovex fp7,fp2; faddx fp6,fp2
1581: fdivx fp7,fp2
1582: fmoved fp2,_Y1:l
1583: fmoved _Z:l,fp2
1584: fsubx fp3,fp2
1585: fmoved fp2,_Z:l
1586: fmoved _Y1:l,fp3
1587: fsubx fp7,fp3
1588: fmoved fp3,_Y2:l
1589: fmoved _F9:l,fp3
1590: fmovex fp3,fp2; fsubd _U1:l,fp2
1591: fdivx fp3,fp2
1592: fmoved fp2,_Y1:l
1593: fmoved _X:l,fp2
1594: fmoved _Zero:l,fp1
1595: fcmpx fp1,fp2; fjne L262
1596: fmoved _Y:l,fp2
1597: fcmpx fp1,fp2; fjne L262
1598: fmoved _Z:l,fp2
1599: fcmpx fp1,fp2; fjne L262
1600: fmoved _T:l,fp2
1601: fcmpx fp1,fp2; fjne L262
1602: fmoved _Y2:l,fp2
1603: fcmpx fp1,fp2; fjne L262
1604: fcmpx fp1,fp2; fjne L262
1605: fmoved _Y1:l,fp2
1606: fmoved _Half:l,fp1
1607: fsubx fp1,fp2
1608: fsubx fp1,fp3
1609: fcmpx fp3,fp2; fjne L262
1610: movl #1,_RDiv:l
1611: movl #L264,sp@
1612: jbsr _printf:l
1613: tstl _GDiv:l; jne L263
1614: movl #L267,sp@
1615: jbsr _notify:l
1616: jra L263
1617: L262:
1618: fmoved _X:l,fp3
1619: fmoved _Zero:l,fp2
1620: fcmpx fp2,fp3; fjge L268
1621: fmoved _Y:l,fp3
1622: fcmpx fp2,fp3; fjge L268
1623: fmoved _Z:l,fp3
1624: fcmpx fp2,fp3; fjge L268
1625: fmoved _T:l,fp3
1626: fcmpx fp2,fp3; fjge L268
1627: fmoved _Y2:l,fp3
1628: fcmpx fp2,fp3; fjge L268
1629: fmoved _Y1:l,fp3
1630: fmoved _Half:l,fp2
1631: fsubx fp2,fp3
1632: fmoved _F9:l,fp1
1633: fsubx fp2,fp1
1634: fcmpx fp1,fp3; fjge L268
1635: movl #2,_RDiv:l
1636: movl #L270,sp@
1637: jbsr _printf:l
1638: L268:
1639: L263:
1640: L260:
1641: tstl _RDiv:l; jne L271
1642: movl #L273,sp@
1643: jbsr _printf:l
1644: L271:
1645: fmoved _One:l,fp3
1646: fdivd _Radix:l,fp3
1647: fmoved fp3,_BInvrse:l
1648: movl #0,sp@
1649: fmoved _BInvrse:l,fp3
1650: fmuld _Radix:l,fp3
1651: fmoved _Half:l,fp2
1652: fsubx fp2,fp3
1653: fcmpx fp2,fp3; fjne L276
1654: movl #1,d7
1655: jra L277
1656: L276:
1657: clrl d7
1658: L277:
1659: movl d7,sp@(4)
1660: movl #L275,sp@(8)
1661: jbsr _TstCond:l
1662: movl #50,_Milestone:l
1663: movl #0,sp@
1664: fmoved _F9:l,fp3
1665: faddd _U1:l,fp3
1666: fmoved _Half:l,fp2
1667: fsubx fp2,fp3
1668: fcmpx fp2,fp3; fjne L280
1669: fmoved _BMinusU2:l,fp3
1670: faddd _U2:l,fp3
1671: fmoved _One:l,fp2
1672: fsubx fp2,fp3
1673: fmoved _Radix:l,fp1
1674: fsubx fp2,fp1
1675: fcmpx fp1,fp3; fjne L280
1676: movl #1,d7
1677: jra L281
1678: L280:
1679: clrl d7
1680: L281:
1681: movl d7,sp@(4)
1682: movl #L279,sp@(8)
1683: jbsr _TstCond:l
1684: fmoved _One:l,fp3
1685: fmoved _U1:l,fp2
1686: fmovex fp2,fp1; fmulx fp2,fp1
1687: fmovex fp3,fp2; fsubx fp1,fp2
1688: fmoved fp2,_X:l
1689: fmoved _U2:l,fp2
1690: fmovex fp3,fp1; fsubx fp2,fp1
1691: fmulx fp1,fp2
1692: fmovex fp3,fp1; faddx fp2,fp1
1693: fmoved fp1,_Y:l
1694: fmoved _F9:l,fp2
1695: fmoved _Half:l,fp1
1696: fsubx fp1,fp2
1697: fmoved fp2,_Z:l
1698: fmoved _X:l,fp2
1699: fsubx fp1,fp2
1700: fsubd _Z:l,fp2
1701: fmoved fp2,_X:l
1702: fmoved _Y:l,fp2
1703: fsubx fp3,fp2
1704: fmoved fp2,_Y:l
1705: fmoved _X:l,fp3
1706: fmoved _Zero:l,fp2
1707: fcmpx fp2,fp3; fjne L282
1708: fmoved _Y:l,fp3
1709: fcmpx fp2,fp3; fjne L282
1710: movl #2,_RAddSub:l
1711: movl #L284,sp@
1712: jbsr _printf:l
1713: L282:
1714: cmpl #1,_GAddSub:l; jne L285
1715: fmoved _Half:l,fp3
1716: fmoved _U2:l,fp2
1717: fmovex fp3,fp1; faddx fp2,fp1
1718: fmulx fp2,fp1
1719: fmoved fp1,_X:l
1720: fsubx fp2,fp3
1721: fmulx fp2,fp3
1722: fmoved fp3,_Y:l
1723: fmoved _One:l,fp3
1724: fmovex fp3,fp1; faddd _X:l,fp1
1725: fmoved fp1,_X:l
1726: fmovex fp3,fp1; faddd _Y:l,fp1
1727: fmoved fp1,_Y:l
1728: fmovex fp3,fp1; faddx fp2,fp1
1729: fmovex fp1,fp2; fsubd _X:l,fp2
1730: fmoved fp2,_X:l
1731: fsubd _Y:l,fp3
1732: fmoved fp3,_Y:l
1733: fmoved _X:l,fp3
1734: fmoved _Zero:l,fp2
1735: fcmpx fp2,fp3; fjne L287
1736: fmoved _Y:l,fp3
1737: fcmpx fp2,fp3; fjne L287
1738: fmoved _Half:l,fp3
1739: fmoved _U2:l,fp2
1740: fmovex fp3,fp1; faddx fp2,fp1
1741: fmoved _U1:l,fp0
1742: fmulx fp0,fp1
1743: fmoved fp1,_X:l
1744: fsubx fp2,fp3
1745: fmulx fp0,fp3
1746: fmoved fp3,_Y:l
1747: fmoved _One:l,fp3
1748: fmovex fp3,fp2; fsubd _X:l,fp2
1749: fmoved fp2,_X:l
1750: fmovex fp3,fp2; fsubd _Y:l,fp2
1751: fmoved fp2,_Y:l
1752: fmoved _F9:l,fp2
1753: fsubd _X:l,fp2
1754: fmoved fp2,_X:l
1755: fsubd _Y:l,fp3
1756: fmoved fp3,_Y:l
1757: fmoved _X:l,fp3
1758: fmoved _Zero:l,fp2
1759: fcmpx fp2,fp3; fjne L289
1760: fmoved _Y:l,fp3
1761: fcmpx fp2,fp3; fjne L289
1762: movl #1,_RAddSub:l
1763: movl #L291,sp@
1764: jbsr _printf:l
1765: tstl _GAddSub:l; jne L286
1766: movl #L294,sp@
1767: jbsr _notify:l
1768: jra L286
1769: L289:
1770: movl #L295,sp@
1771: jbsr _printf:l
1772: jra L286
1773: L287:
1774: movl #L295,sp@
1775: jbsr _printf:l
1776: jra L286
1777: L285:
1778: movl #L295,sp@
1779: jbsr _printf:l
1780: L286:
1781: fmoved _One:l,fp3
1782: fmoved fp3,_S:l
1783: fmoved _Half:l,fp2
1784: fmovex fp3,fp1; faddx fp2,fp1
1785: fmovex fp2,fp0; fmulx fp1,fp0
1786: fmovex fp3,fp1; faddx fp0,fp1
1787: fmoved fp1,_X:l
1788: faddd _U2:l,fp3
1789: fmulx fp2,fp3
1790: fmoved fp3,_Y:l
1791: fmoved _X:l,fp3
1792: fmoved _Y:l,fp2
1793: fmovex fp3,fp1; fsubx fp2,fp1
1794: fmoved fp1,_Z:l
1795: fsubx fp3,fp2
1796: fmoved fp2,_T:l
1797: fmoved _Z:l,fp3
1798: faddd _T:l,fp3
1799: fmoved fp3,_StickyBit:l
1800: fmoved _StickyBit:l,fp3
1801: fcmpd _Zero:l,fp3; fjeq L296
1802: fmoved _Zero:l,fp3
1803: fmoved fp3,_S:l
1804: movl #3,sp@
1805: movl #L298,sp@(4)
1806: jbsr _BadCond:l
1807: L296:
1808: fmoved _Zero:l,fp3
1809: fmoved fp3,_StickyBit:l
1810: cmpl #1,_GMult:l; jne L299
1811: cmpl #1,_GDiv:l; jne L299
1812: cmpl #1,_GAddSub:l; jne L299
1813: cmpl #1,_RMult:l; jne L299
1814: cmpl #1,_RDiv:l; jne L299
1815: cmpl #1,_RAddSub:l; jne L299
1816: fmoved _RadixD2:l,fp3
1817: fmoved fp3,sp@(0)
1818: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
1819: fcmpd _RadixD2:l,fp0; fjne L299
1820: movl #L301,sp@
1821: jbsr _printf:l
1822: fmoved _Half:l,fp3
1823: fmovex fp3,fp2; faddd _U1:l,fp2
1824: fmoved _U2:l,fp1
1825: fmulx fp1,fp2
1826: fmoved fp2,_X:l
1827: fmulx fp1,fp3
1828: fmoved fp3,_Y:l
1829: fmoved _One:l,fp3
1830: fmovex fp3,fp2; faddd _Y:l,fp2
1831: fmoved fp2,_Z:l
1832: fmovex fp3,fp2; faddd _X:l,fp2
1833: fmoved fp2,_T:l
1834: fmoved _Z:l,fp2
1835: fsubx fp3,fp2
1836: fcmpd _Zero:l,fp2; fjgt L302
1837: fmoved _T:l,fp2
1838: fsubx fp3,fp2
1839: fcmpx fp1,fp2; fjlt L302
1840: fmoved _T:l,fp3
1841: fmovex fp3,fp2; faddd _Y:l,fp2
1842: fmoved fp2,_Z:l
1843: fmoved _Z:l,fp2
1844: fmovex fp2,fp1; fsubd _X:l,fp1
1845: fmoved fp1,_Y:l
1846: fsubx fp3,fp2
1847: fcmpd _U2:l,fp2; fjlt L304
1848: fmoved _Y:l,fp2
1849: fsubx fp3,fp2
1850: fcmpd _Zero:l,fp2; fjne L304
1851: fmoved _Half:l,fp3
1852: fmoved _U1:l,fp2
1853: fmovex fp3,fp1; faddx fp2,fp1
1854: fmulx fp2,fp1
1855: fmoved fp1,_X:l
1856: fmulx fp2,fp3
1857: fmoved fp3,_Y:l
1858: fmoved _One:l,fp3
1859: fmovex fp3,fp2; fsubd _Y:l,fp2
1860: fmoved fp2,_Z:l
1861: fmovex fp3,fp2; fsubd _X:l,fp2
1862: fmoved fp2,_T:l
1863: fmoved _Z:l,fp2
1864: fsubx fp3,fp2
1865: fmoved _Zero:l,fp3
1866: fcmpx fp3,fp2; fjne L306
1867: fmoved _T:l,fp2
1868: fsubd _F9:l,fp2
1869: fcmpx fp3,fp2; fjne L306
1870: fmoved _Half:l,fp3
1871: fmoved _U1:l,fp2
1872: fsubx fp2,fp3
1873: fmulx fp2,fp3
1874: fmoved fp3,_Z:l
1875: fmoved _F9:l,fp3
1876: fmovex fp3,fp1; fsubd _Z:l,fp1
1877: fmoved fp1,_T:l
1878: fmovex fp3,fp1; fsubd _Y:l,fp1
1879: fmoved fp1,_Q:l
1880: fmoved _T:l,fp1
1881: fsubx fp3,fp1
1882: fmoved _Zero:l,fp0
1883: fcmpx fp0,fp1; fjne L308
1884: fsubx fp2,fp3
1885: fsubd _Q:l,fp3
1886: fcmpx fp0,fp3; fjne L308
1887: fmoved _One:l,fp3
1888: fmoved _U2:l,fp2
1889: fmovex fp2,fp7
1890: fmovex fp3,fp1; faddx fp2,fp1
1891: fmoved _OneAndHalf:l,fp0
1892: fmulx fp0,fp1
1893: fmoved fp1,_Z:l
1894: fmovex fp0,fp1; faddx fp2,fp1
1895: fsubd _Z:l,fp1
1896: faddx fp2,fp1
1897: fmoved fp1,_T:l
1898: fmoved _Half:l,fp1
1899: fmoved _Radix:l,fp0
1900: fdivx fp0,fp1
1901: fmovex fp3,fp2; faddx fp1,fp2
1902: fmoved fp2,_X:l
1903: fmovex fp0,fp2; fmulx fp7,fp2
1904: faddx fp2,fp3
1905: fmoved fp3,_Y:l
1906: fmoved _X:l,fp3
1907: fmovex fp3,fp1; fmuld _Y:l,fp1
1908: fmoved fp1,_Z:l
1909: fmoved _T:l,fp1
1910: fmoved _Zero:l,fp0
1911: fcmpx fp0,fp1; fjne L310
1912: faddx fp2,fp3
1913: fsubd _Z:l,fp3
1914: fcmpx fp0,fp3; fjne L310
1915: fmoved _Radix:l,fp3
1916: fcmpd _Two:l,fp3; fjeq L312
1917: fmoved _Two:l,fp3
1918: fmovex fp3,fp2; faddd _U2:l,fp2
1919: fmoved fp2,_X:l
1920: fmoved _X:l,fp2
1921: fdivx fp3,fp2
1922: fmoved fp2,_Y:l
1923: fmoved _Y:l,fp3
1924: fsubd _One:l,fp3
1925: fcmpd _Zero:l,fp3; fjne L313
1926: fmoved _S:l,fp3
1927: fmoved fp3,_StickyBit:l
1928: jra L313
1929: L312:
1930: fmoved _S:l,fp3
1931: fmoved fp3,_StickyBit:l
1932: L313:
1933: L310:
1934: L308:
1935: L306:
1936: L304:
1937: L302:
1938: L299:
1939: fmoved _StickyBit:l,fp3
1940: fcmpd _One:l,fp3; fjne L316
1941: movl #L318,sp@
1942: jbsr _printf:l
1943: jra L317
1944: L316:
1945: movl #L319,sp@
1946: jbsr _printf:l
1947: L317:
1948: movl #3,sp@
1949: tstl _GMult:l; jeq L322
1950: tstl _GDiv:l; jeq L322
1951: tstl _GAddSub:l; jeq L322
1952: tstl _RMult:l; jeq L322
1953: tstl _RDiv:l; jeq L322
1954: tstl _RAddSub:l; jeq L322
1955: movl #1,d7
1956: jra L323
1957: L322:
1958: clrl d7
1959: L323:
1960: movl d7,sp@(4)
1961: movl #L321,sp@(8)
1962: jbsr _TstCond:l
1963: movl #60,_Milestone:l
1964: movl #L47,sp@
1965: jbsr _printf:l
1966: movl #L324,sp@
1967: jbsr _printf:l
1968: movl #L325,sp@
1969: movl _NoTrials:l,sp@(4)
1970: jbsr _printf:l
1971: fmoved L326:l,fp3
1972: fmoved fp3,sp@(0)
1973: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
1974: fmoved fp0,_Random9:l
1975: fmoved _Third:l,fp3
1976: fmoved fp3,_Random1:l
1977: movl #1,_I:l
1978: L327:
1979: jbsr _Random:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
1980: fmoved fp0,_X:l
1981: jbsr _Random:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
1982: fmoved fp0,_Y:l
1983: fmoved _Y:l,fp3
1984: fmoved _X:l,fp2
1985: fmovex fp3,fp1; fmulx fp2,fp1
1986: fmoved fp1,_Z9:l
1987: fmulx fp3,fp2
1988: fmoved fp2,_Z:l
1989: fmoved _Z:l,fp3
1990: fsubd _Z9:l,fp3
1991: fmoved fp3,_Z9:l
1992: addql #1,_I:l
1993: L328:
1994: movl _I:l,d4
1995: cmpl _NoTrials:l,d4; jgt L330
1996: fmoved _Z9:l,fp3
1997: fcmpd _Zero:l,fp3; fjeq L327
1998: L330:
1999: movl _I:l,d4
2000: cmpl _NoTrials:l,d4; jne L331
2001: fmoved _One:l,fp3
2002: fmoved _Half:l,fp2
2003: fdivd _Three:l,fp2
2004: fmovex fp3,fp1; faddx fp2,fp1
2005: fmoved fp1,_Random1:l
2006: fmoved _U2:l,fp2
2007: faddd _U1:l,fp2
2008: faddx fp3,fp2
2009: fmovex fp2,fp7
2010: fmoved fp2,_Random2:l
2011: fmoved _Random1:l,fp3
2012: fmoved _Random2:l,fp0
2013: fmovex fp3,fp2; fmulx fp0,fp2
2014: fmoved fp2,_Z:l
2015: fmovex fp0,fp2; fmulx fp3,fp2
2016: fmoved fp2,_Y:l
2017: fmovex fp1,fp3; fmulx fp7,fp3
2018: fmovex fp3,fp2; fsubx fp3,fp2
2019: fmoved fp2,_Z9:l
2020: L331:
2021: movl _I:l,d4
2022: cmpl _NoTrials:l,d4; jeq L333
2023: fmoved _Z9:l,fp3
2024: fcmpd _Zero:l,fp3; fjeq L333
2025: movl #2,sp@
2026: movl #L335,sp@(4)
2027: jbsr _BadCond:l
2028: jra L334
2029: L333:
2030: movl #L336,sp@
2031: movl _NoTrials:l,sp@(4)
2032: jbsr _printf:l
2033: L334:
2034: movl #70,_Milestone:l
2035: movl #L337,sp@
2036: jbsr _printf:l
2037: fmoved _Zero:l,fp3
2038: fmoved fp3,sp@(0)
2039: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2040: fcmpx fp0,fp3; fjne L340
2041: fnegd _Zero:l,fp3
2042: fmoved fp3,sp@(0)
2043: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2044: fcmpx fp0,fp3; fjne L340
2045: fmoved _One:l,fp3
2046: fmoved fp3,sp@(0)
2047: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2048: fcmpx fp0,fp3; fjne L340
2049: movl #1,d7
2050: jra L341
2051: L340:
2052: clrl d7
2053: L341:
2054: movl #0,sp@
2055: movl d7,sp@(4)
2056: movl #L339,sp@(8)
2057: jbsr _TstCond:l
2058: fmoved _Zero:l,fp3
2059: fmoved fp3,_MinSqEr:l
2060: fmoved fp3,_MaxSqEr:l
2061: fmoved fp3,_J:l
2062: fmoved _Radix:l,fp3
2063: fmoved fp3,_X:l
2064: fmoved _U2:l,fp3
2065: fmoved fp3,_OneUlp:l
2066: movl #1,sp@
2067: jbsr _SqXMinX:l
2068: fmoved _BInvrse:l,fp3
2069: fmoved fp3,_X:l
2070: fmuld _U1:l,fp3
2071: fmoved fp3,_OneUlp:l
2072: movl #1,sp@
2073: jbsr _SqXMinX:l
2074: fmoved _U1:l,fp3
2075: fmoved fp3,_X:l
2076: fmovex fp3,fp2; fmulx fp3,fp2
2077: fmoved fp2,_OneUlp:l
2078: movl #1,sp@
2079: jbsr _SqXMinX:l
2080: fmoved _J:l,fp3
2081: fcmpd _Zero:l,fp3; fjeq L342
2082: jbsr _Pause:l
2083: L342:
2084: movl #L344,sp@
2085: movl _NoTrials:l,sp@(4)
2086: jbsr _printf:l
2087: fmoved _Zero:l,fp3
2088: fmoved fp3,_J:l
2089: fmoved _Two:l,fp3
2090: fmoved fp3,_X:l
2091: fmoved _Radix:l,fp3
2092: fmoved fp3,_Y:l
2093: fcmpd _One:l,fp3; fjeq L345
2094: L347:
2095: fmoved _Y:l,fp3
2096: fmoved fp3,_X:l
2097: fmoved _Radix:l,fp2
2098: fmulx fp3,fp2
2099: fmoved fp2,_Y:l
2100: L348:
2101: fmoved _Y:l,fp3
2102: fsubd _X:l,fp3
2103: movl _NoTrials:l,d4
2104: fcmpl d4,fp3; fjlt L347
2105: L345:
2106: fmoved _X:l,fp3
2107: fmuld _U2:l,fp3
2108: fmoved fp3,_OneUlp:l
2109: movl #1,_I:l
2110: jra L351
2111: L350:
2112: fmoved _X:l,fp3
2113: faddd _One:l,fp3
2114: fmoved fp3,_X:l
2115: movl #2,sp@
2116: jbsr _SqXMinX:l
2117: fmoved _J:l,fp3
2118: fcmpd _Zero:l,fp3; fjle L353
2119: jra L352
2120: L353:
2121: addql #1,_I:l
2122: L351:
2123: movl _I:l,d4
2124: cmpl _NoTrials:l,d4; jle L350
2125: L352:
2126: movl #L355,sp@
2127: jbsr _printf:l
2128: movl #-1,_I:l
2129: fmoved _BMinusU2:l,fp3
2130: fmoved fp3,_X:l
2131: fmoved _Radix:l,fp3
2132: fmoved fp3,_Y:l
2133: fmovex fp3,fp2; fmuld _U2:l,fp2
2134: faddx fp2,fp3
2135: fmoved fp3,_Z:l
2136: clrl _NotMonot:l
2137: clrl _Monot:l
2138: jra L357
2139: L356:
2140: addql #1,_I:l
2141: fmoved _X:l,fp3
2142: fmoved fp3,sp@(0)
2143: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2144: fmoved fp0,_X:l
2145: fmoved _Y:l,fp3
2146: fmoved fp3,sp@(0)
2147: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2148: fmoved fp0,_Q:l
2149: fmoved _Z:l,fp3
2150: fmoved fp3,sp@(0)
2151: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2152: fmoved fp0,_Z:l
2153: fmoved _X:l,fp3
2154: fmoved _Q:l,fp2
2155: fcmpx fp2,fp3; fjgt L361
2156: fcmpd _Z:l,fp2; fjle L359
2157: L361:
2158: movl #1,_NotMonot:l
2159: jra L360
2160: L359:
2161: fmoved _Q:l,fp3
2162: faddd _Half:l,fp3
2163: fmoved fp3,sp@(0)
2164: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2165: fmoved fp0,_Q:l
2166: tstl _I:l; jgt L364
2167: fmoved _Radix:l,fp3
2168: fmoved _Q:l,fp2
2169: fmovex fp2,fp1; fmulx fp2,fp1
2170: fcmpx fp1,fp3; fjne L362
2171: L364:
2172: movl #1,_Monot:l
2173: jra L363
2174: L362:
2175: tstl _I:l; jle L365
2176: cmpl #1,_I:l; jle L367
2177: movl #1,_Monot:l
2178: jra L366
2179: L367:
2180: fmoved _Y:l,fp3
2181: fmuld _BInvrse:l,fp3
2182: fmoved fp3,_Y:l
2183: fmoved _Y:l,fp3
2184: fmoved _U1:l,fp2
2185: fmovex fp3,fp1; fsubx fp2,fp1
2186: fmoved fp1,_X:l
2187: faddx fp2,fp3
2188: fmoved fp3,_Z:l
2189: jra L366
2190: L365:
2191: fmoved _Q:l,fp3
2192: fmoved fp3,_Y:l
2193: fmoved _Y:l,fp3
2194: fmoved _U2:l,fp2
2195: fmovex fp3,fp1; fsubx fp2,fp1
2196: fmoved fp1,_X:l
2197: faddx fp2,fp3
2198: fmoved fp3,_Z:l
2199: L366:
2200: L363:
2201: L360:
2202: L357:
2203: tstl _NotMonot:l; jne L369
2204: tstl _Monot:l; jeq L356
2205: L369:
2206: tstl _Monot:l; jeq L370
2207: movl #L372,sp@
2208: jbsr _printf:l
2209: jra L371
2210: L370:
2211: movl #2,sp@
2212: movl #L160,sp@(4)
2213: jbsr _BadCond:l
2214: movl #L373,sp@
2215: fmoved _Y:l,fp3
2216: fmoved fp3,sp@(4)
2217: jbsr _printf:l
2218: L371:
2219: movl #80,_Milestone:l
2220: fmoved _MinSqEr:l,fp3
2221: fmoved _Half:l,fp2
2222: faddx fp2,fp3
2223: fmoved fp3,_MinSqEr:l
2224: fmoved _MaxSqEr:l,fp3
2225: fsubx fp2,fp3
2226: fmoved fp3,_MaxSqEr:l
2227: fmoved _One:l,fp3
2228: faddd _U2:l,fp3
2229: fmoved fp3,sp@(0)
2230: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2231: fmoved _One:l,fp3
2232: fmovex fp0,fp2; fsubx fp3,fp2
2233: fmoved _U2:l,fp1
2234: fdivx fp1,fp2
2235: fmoved fp2,_Y:l
2236: fmoved _Y:l,fp2
2237: fsubx fp3,fp2
2238: fmovex fp1,fp3; fdivd _Eight:l,fp3
2239: faddx fp3,fp2
2240: fmoved fp2,_SqEr:l
2241: fmoved _SqEr:l,fp3
2242: fcmpd _MaxSqEr:l,fp3; fjle L374
2243: fmoved _SqEr:l,fp3
2244: fmoved fp3,_MaxSqEr:l
2245: L374:
2246: fmoved _Y:l,fp3
2247: fmoved _U2:l,fp2
2248: fdivd _Eight:l,fp2
2249: faddx fp2,fp3
2250: fmoved fp3,_SqEr:l
2251: fmoved _SqEr:l,fp3
2252: fcmpd _MinSqEr:l,fp3; fjge L376
2253: fmoved _SqEr:l,fp3
2254: fmoved fp3,_MinSqEr:l
2255: L376:
2256: fmoved _F9:l,fp3
2257: fmoved fp3,sp@(0)
2258: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2259: fmoved _U2:l,fp3
2260: fmovex fp0,fp2; fsubx fp3,fp2
2261: fmoved _One:l,fp1
2262: fsubx fp3,fp1
2263: fmovex fp2,fp3; fsubx fp1,fp3
2264: fmoved _U1:l,fp2
2265: fdivx fp2,fp3
2266: fmoved fp3,_Y:l
2267: fmoved _Y:l,fp3
2268: fdivd _Eight:l,fp2
2269: faddx fp2,fp3
2270: fmoved fp3,_SqEr:l
2271: fmoved _SqEr:l,fp3
2272: fcmpd _MaxSqEr:l,fp3; fjle L378
2273: fmoved _SqEr:l,fp3
2274: fmoved fp3,_MaxSqEr:l
2275: L378:
2276: fmoved _Y:l,fp3
2277: faddd _One:l,fp3
2278: fmoved _U1:l,fp2
2279: fdivd _Eight:l,fp2
2280: faddx fp2,fp3
2281: fmoved fp3,_SqEr:l
2282: fmoved _SqEr:l,fp3
2283: fcmpd _MinSqEr:l,fp3; fjge L380
2284: fmoved _SqEr:l,fp3
2285: fmoved fp3,_MinSqEr:l
2286: L380:
2287: fmoved _U2:l,fp3
2288: fmoved fp3,_OneUlp:l
2289: fmoved _OneUlp:l,fp3
2290: fmoved fp3,_X:l
2291: movl #1,_Indx:l
2292: L382:
2293: fmoved _X:l,fp3
2294: fmovex fp3,fp2; faddd _U1:l,fp2
2295: faddx fp3,fp2
2296: fmovex fp2,fp3; faddd _F9:l,fp3
2297: fmoved fp3,sp@(0)
2298: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2299: fmoved fp0,_Y:l
2300: fmoved _Y:l,fp3
2301: fmoved _U2:l,fp2
2302: fsubx fp2,fp3
2303: fmoved _One:l,fp1
2304: fsubx fp2,fp1
2305: fmoved _X:l,fp2
2306: faddx fp2,fp1
2307: fsubx fp1,fp3
2308: fmoved _OneUlp:l,fp1
2309: fdivx fp1,fp3
2310: fmoved fp3,_Y:l
2311: fmoved _U1:l,fp3
2312: fsubx fp2,fp3
2313: faddd _F9:l,fp3
2314: fmoved _Half:l,fp0
2315: fmulx fp0,fp3
2316: fmulx fp2,fp3
2317: fmulx fp2,fp3
2318: fdivx fp1,fp3
2319: fmoved fp3,_Z:l
2320: fmoved _Y:l,fp3
2321: faddx fp0,fp3
2322: faddd _Z:l,fp3
2323: fmoved fp3,_SqEr:l
2324: fmoved _SqEr:l,fp3
2325: fcmpd _MinSqEr:l,fp3; fjge L386
2326: fmoved _SqEr:l,fp3
2327: fmoved fp3,_MinSqEr:l
2328: L386:
2329: fmoved _Y:l,fp3
2330: fsubd _Half:l,fp3
2331: faddd _Z:l,fp3
2332: fmoved fp3,_SqEr:l
2333: fmoved _SqEr:l,fp3
2334: fcmpd _MaxSqEr:l,fp3; fjle L388
2335: fmoved _SqEr:l,fp3
2336: fmoved fp3,_MaxSqEr:l
2337: L388:
2338: movl _Indx:l,d4
2339: cmpl #1,d4; jeq L392
2340: cmpl #3,d4; jne L390
2341: L392:
2342: fmoved _X:l,fp3
2343: fmoved fp3,sp@(0)
2344: jbsr _Sign:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2345: fmovex fp0,fp7
2346: fmoved _OneUlp:l,fp3
2347: fmoved fp3,sp@(0)
2348: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2349: fmoved _Eight:l,fp3
2350: fmoved _Nine:l,fp2
2351: fmulx fp0,fp2
2352: fdivx fp2,fp3
2353: fmoved fp3,sp@(0)
2354: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2355: fmoved _OneUlp:l,fp3
2356: fmulx fp7,fp3
2357: fmulx fp0,fp3
2358: fmoved fp3,_X:l
2359: jra L391
2360: L390:
2361: fmoved _U1:l,fp3
2362: fmoved fp3,_OneUlp:l
2363: fnegd _OneUlp:l,fp3
2364: fmoved fp3,_X:l
2365: L391:
2366: L383:
2367: addql #1,_Indx:l
2368: cmpl #3,_Indx:l; jle L382
2369: movl #85,_Milestone:l
2370: clrl _SqRWrng:l
2371: clrl _Anomaly:l
2372: clrl _RSqrt:l
2373: fmoved _Radix:l,fp3
2374: fcmpd _One:l,fp3; fjeq L393
2375: movl #L395,sp@
2376: jbsr _printf:l
2377: fmoved _Precision:l,fp3
2378: fmoved fp3,sp@(0)
2379: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2380: fmoved _Radix:l,fp2
2381: fmoved fp2,sp@(0)
2382: fmoved _One:l,fp2
2383: faddx fp3,fp2
2384: fmovex fp2,fp3; fsubx fp0,fp3
2385: fmoved fp3,sp@(8)
2386: jbsr _pow:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2387: fmoved _Half:l,fp3
2388: faddx fp0,fp3
2389: fmoved fp3,sp@(0)
2390: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2391: fmoved fp0,_D:l
2392: fmoved _D:l,fp3
2393: fmovex fp3,fp2; fdivd _Radix:l,fp2
2394: fmoved fp2,_X:l
2395: fdivd _A1:l,fp3
2396: fmoved fp3,_Y:l
2397: fmoved _X:l,fp3
2398: fmoved fp3,sp@(0)
2399: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2400: fcmpx fp0,fp3; fjne L398
2401: fmoved _Y:l,fp3
2402: fmoved fp3,sp@(0)
2403: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2404: fcmpx fp0,fp3; fjeq L396
2405: L398:
2406: movl #1,_Anomaly:l
2407: jra L397
2408: L396:
2409: fmoved _Zero:l,fp3
2410: fmoved fp3,_X:l
2411: fmoved _X:l,fp3
2412: fmoved fp3,_Z2:l
2413: fmoved _One:l,fp3
2414: fmoved fp3,_Y:l
2415: fmoved _Y:l,fp2
2416: fmoved fp2,_Y2:l
2417: fmoved _Radix:l,fp2
2418: fsubx fp3,fp2
2419: fmoved fp2,_Z1:l
2420: fmoved _Four:l,fp3
2421: fmuld _D:l,fp3
2422: fmoved fp3,_FourD:l
2423: L399:
2424: fmoved _Y2:l,fp3
2425: fcmpd _Z2:l,fp3; fjle L402
2426: fmoved _Radix:l,fp3
2427: fmoved fp3,_Q:l
2428: fmoved _Y:l,fp3
2429: fmoved fp3,_Y1:l
2430: L404:
2431: fmoved _Half:l,fp3
2432: fmoved _Q:l,fp2
2433: fmovex fp2,fp1; fdivd _Y1:l,fp1
2434: fsubx fp1,fp3
2435: fmoved fp3,sp@(0)
2436: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2437: fmovex fp0,fp3; fmuld _Y1:l,fp3
2438: faddx fp3,fp2
2439: fmoved fp2,sp@(0)
2440: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2441: fmoved fp0,_X1:l
2442: fmoved _Y1:l,fp3
2443: fmoved fp3,_Q:l
2444: fmoved _X1:l,fp3
2445: fmoved fp3,_Y1:l
2446: L405:
2447: fmoved _X1:l,fp3
2448: fcmpd _Zero:l,fp3; fjgt L404
2449: fmoved _Q:l,fp3
2450: fcmpd _One:l,fp3; fjgt L407
2451: fmoved _Y2:l,fp3
2452: fmoved fp3,_Z2:l
2453: fmoved _Y:l,fp3
2454: fmoved fp3,_Z:l
2455: L407:
2456: L402:
2457: fmoved _Y:l,fp3
2458: faddd _Two:l,fp3
2459: fmoved fp3,_Y:l
2460: fmoved _X:l,fp3
2461: faddd _Eight:l,fp3
2462: fmoved fp3,_X:l
2463: fmoved _Y2:l,fp3
2464: faddd _X:l,fp3
2465: fmoved fp3,_Y2:l
2466: fmoved _Y2:l,fp3
2467: fcmpd _FourD:l,fp3; fjlt L409
2468: fmoved _Y2:l,fp3
2469: fsubd _FourD:l,fp3
2470: fmoved fp3,_Y2:l
2471: L409:
2472: L400:
2473: fmoved _Y:l,fp3
2474: fcmpd _D:l,fp3; fjlt L399
2475: fmoved _FourD:l,fp3
2476: fmovex fp3,fp2; fsubd _Z2:l,fp2
2477: fmoved fp2,_X8:l
2478: fmoved _X8:l,fp2
2479: fmoved _Z:l,fp1
2480: fmovex fp1,fp0; fmulx fp1,fp0
2481: fmovex fp2,fp1; faddx fp0,fp1
2482: fdivx fp3,fp1
2483: fmoved fp1,_Q:l
2484: fmovex fp2,fp3; fdivd _Eight:l,fp3
2485: fmoved fp3,_X8:l
2486: fmoved _Q:l,fp3
2487: fmoved fp3,sp@(0)
2488: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2489: fcmpx fp0,fp3; fjeq L411
2490: movl #1,_Anomaly:l
2491: jra L412
2492: L411:
2493: clrl _Break:l
2494: L413:
2495: fmoved _Z1:l,fp3
2496: fmuld _Z:l,fp3
2497: fmoved fp3,_X:l
2498: fmoved _X:l,fp3
2499: fmovex fp3,fp2; fdivd _Radix:l,fp2
2500: fmoved fp2,sp@(0)
2501: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2502: fmovex fp0,fp2; fmuld _Radix:l,fp2
2503: fsubx fp2,fp3
2504: fmoved fp3,_X:l
2505: fmoved _X:l,fp3
2506: fcmpd _One:l,fp3; fjne L416
2507: movl #1,_Break:l
2508: jra L417
2509: L416:
2510: fmoved _Z1:l,fp3
2511: fsubd _One:l,fp3
2512: fmoved fp3,_Z1:l
2513: L417:
2514: L414:
2515: tstl _Break:l; jne L418
2516: fmoved _Z1:l,fp3
2517: fcmpd _Zero:l,fp3; fjgt L413
2518: L418:
2519: fmoved _Z1:l,fp3
2520: fcmpd _Zero:l,fp3; fjgt L419
2521: tstl _Break:l; jne L419
2522: movl #1,_Anomaly:l
2523: jra L420
2524: L419:
2525: fmoved _Z1:l,fp3
2526: fcmpd _RadixD2:l,fp3; fjle L421
2527: fmoved _Z1:l,fp3
2528: fsubd _Radix:l,fp3
2529: fmoved fp3,_Z1:l
2530: L421:
2531: L423:
2532: jbsr _NewD:l
2533: L424:
2534: fmoved _U2:l,fp3
2535: fmuld _D:l,fp3
2536: fcmpd _F9:l,fp3; fjlt L423
2537: fmoved _D:l,fp3
2538: fmovex fp3,fp2; fmuld _Radix:l,fp2
2539: fsubx fp3,fp2
2540: fmoved _W:l,fp1
2541: fsubx fp3,fp1
2542: fcmpx fp1,fp2; fjeq L426
2543: movl #1,_Anomaly:l
2544: jra L427
2545: L426:
2546: fmoved _D:l,fp3
2547: fmoved fp3,_Z2:l
2548: clrl _I:l
2549: fmoved _One:l,fp2
2550: fmoved _Z:l,fp1
2551: faddx fp1,fp2
2552: fmuld _Half:l,fp2
2553: fmovex fp3,fp0; faddx fp2,fp0
2554: fmoved fp0,_Y:l
2555: faddx fp1,fp3
2556: faddd _Q:l,fp3
2557: fmoved fp3,_X:l
2558: jbsr _SR3750:l
2559: fmoved _D:l,fp3
2560: fmoved _One:l,fp2
2561: fmoved _Z:l,fp1
2562: fsubx fp1,fp2
2563: fmuld _Half:l,fp2
2564: fmovex fp3,fp0; faddx fp2,fp0
2565: fmovex fp0,fp2; faddx fp3,fp2
2566: fmoved fp2,_Y:l
2567: fmovex fp3,fp2; fsubx fp1,fp2
2568: faddx fp3,fp2
2569: fmoved fp2,_X:l
2570: fmoved _X:l,fp3
2571: fmovex fp3,fp2; faddd _Q:l,fp2
2572: faddx fp3,fp2
2573: fmoved fp2,_X:l
2574: jbsr _SR3750:l
2575: jbsr _NewD:l
2576: fmoved _D:l,fp3
2577: fmoved _Z2:l,fp2
2578: fsubx fp2,fp3
2579: fmoved _W:l,fp1
2580: fsubx fp2,fp1
2581: fcmpx fp1,fp3; fjeq L428
2582: movl #1,_Anomaly:l
2583: jra L429
2584: L428:
2585: fmoved _D:l,fp3
2586: fmoved _Z2:l,fp2
2587: fsubx fp2,fp3
2588: fmoved _One:l,fp1
2589: fmoved _Z:l,fp0
2590: fmovex fp0,fp7
2591: fsubx fp0,fp1
2592: fmuld _Half:l,fp1
2593: fmovex fp2,fp0; faddx fp1,fp0
2594: fmovex fp3,fp1; faddx fp0,fp1
2595: fmoved fp1,_Y:l
2596: fsubx fp7,fp2
2597: faddd _Q:l,fp2
2598: faddx fp2,fp3
2599: fmoved fp3,_X:l
2600: jbsr _SR3750:l
2601: fmoved _One:l,fp3
2602: faddd _Z:l,fp3
2603: fmuld _Half:l,fp3
2604: fmoved fp3,_Y:l
2605: fmoved _Q:l,fp3
2606: fmoved fp3,_X:l
2607: jbsr _SR3750:l
2608: tstl _I:l; jne L430
2609: movl #1,_Anomaly:l
2610: L430:
2611: L429:
2612: L427:
2613: L420:
2614: L412:
2615: L397:
2616: tstl _I:l; jeq L434
2617: tstl _Anomaly:l; jeq L432
2618: L434:
2619: movl #0,sp@
2620: movl #L435,sp@(4)
2621: jbsr _BadCond:l
2622: movl #L436,sp@
2623: fmoved _W:l,fp3
2624: fmoved fp3,sp@(4)
2625: jbsr _printf:l
2626: movl #L437,sp@
2627: jbsr _printf:l
2628: movl #1,_SqRWrng:l
2629: L432:
2630: L393:
2631: tstl _Anomaly:l; jne L438
2632: fmoved _MinSqEr:l,fp3
2633: fmoved _Zero:l,fp2
2634: fcmpx fp2,fp3; fjlt L440
2635: fmoved _MaxSqEr:l,fp3
2636: fcmpx fp2,fp3; fjgt L440
2637: movl #1,_RSqrt:l
2638: movl #L442,sp@
2639: jbsr _printf:l
2640: jra L441
2641: L440:
2642: fmoved _MaxSqEr:l,fp3
2643: fmoved _U2:l,fp2
2644: faddx fp2,fp3
2645: fmoved _Half:l,fp1
2646: fsubx fp1,fp2
2647: fcmpx fp2,fp3; fjgt L446
2648: fmoved _MinSqEr:l,fp3
2649: fcmpx fp1,fp3; fjgt L446
2650: faddd _Radix:l,fp3
2651: fcmpx fp1,fp3; fjge L443
2652: L446:
2653: movl #1,_SqRWrng:l
2654: jra L444
2655: L443:
2656: movl #2,_RSqrt:l
2657: movl #L447,sp@
2658: jbsr _printf:l
2659: L444:
2660: L441:
2661: L438:
2662: tstl _SqRWrng:l; jeq L448
2663: movl #L450,sp@
2664: jbsr _printf:l
2665: movl #L451,sp@
2666: fmoved _MinSqEr:l,fp3
2667: fsubd _Half:l,fp3
2668: fmoved fp3,sp@(4)
2669: jbsr _printf:l
2670: movl #L452,sp@
2671: fmoved _Half:l,fp3
2672: faddd _MaxSqEr:l,fp3
2673: fmoved fp3,sp@(4)
2674: jbsr _printf:l
2675: movl #1,sp@
2676: fmoved _MaxSqEr:l,fp3
2677: fsubd _MinSqEr:l,fp3
2678: fmoved _Radix:l,fp2
2679: fmovex fp2,fp1; fmulx fp2,fp1
2680: fcmpx fp1,fp3; fjge L455
2681: movl #1,d7
2682: jra L456
2683: L455:
2684: clrl d7
2685: L456:
2686: movl d7,sp@(4)
2687: movl #L454,sp@(8)
2688: jbsr _TstCond:l
2689: L448:
2690: movl #90,_Milestone:l
2691: jbsr _Pause:l
2692: movl #L457,sp@
2693: jbsr _printf:l
2694: clrl _N:l
2695: clrl _I:l
2696: fnegd _Zero:l,fp3
2697: fmoved fp3,_Z:l
2698: movl #3,_M:l
2699: clrl _Break:l
2700: L458:
2701: fmoved _One:l,fp3
2702: fmoved fp3,_X:l
2703: jbsr _SR3980:l
2704: cmpl #10,_I:l; jgt L461
2705: movl #1023,_I:l
2706: jbsr _SR3980:l
2707: L461:
2708: fmoved _Z:l,fp3
2709: fcmpd _MinusOne:l,fp3; fjne L463
2710: movl #1,_Break:l
2711: jra L464
2712: L463:
2713: fmoved _MinusOne:l,fp3
2714: fmoved fp3,_Z:l
2715: jbsr _PrintIfNPositive:l
2716: clrl _N:l
2717: movl #-4,_I:l
2718: L464:
2719: L459:
2720: tstl _Break:l; jeq L458
2721: jbsr _PrintIfNPositive:l
2722: movl _N:l,_N1:l
2723: clrl _N:l
2724: fmoved _A1:l,fp3
2725: fmoved fp3,_Z:l
2726: fmoved _W:l,fp3
2727: fmoved fp3,sp@(0)
2728: jbsr _log:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2729: fmovex fp0,fp7
2730: fmoved _A1:l,fp3
2731: fmoved fp3,sp@(0)
2732: jbsr _log:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2733: fmoved _Two:l,fp3
2734: fmulx fp7,fp3
2735: fdivx fp0,fp3
2736: fmoved fp3,sp@(0)
2737: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2738: fintrzx fp0,fp0; fmovel fp0,d4
2739: movl d4,_M:l
2740: clrl _Break:l
2741: L465:
2742: fmoved _Z:l,fp3
2743: fmoved fp3,_X:l
2744: movl #1,_I:l
2745: jbsr _SR3980:l
2746: fmoved _Z:l,fp3
2747: fcmpd _AInvrse:l,fp3; fjne L468
2748: movl #1,_Break:l
2749: jra L469
2750: L468:
2751: fmoved _AInvrse:l,fp3
2752: fmoved fp3,_Z:l
2753: L469:
2754: L466:
2755: tstl _Break:l; jeq L465
2756: movl #100,_Milestone:l
2757: movl _NoTrials:l,_M:l
2758: fmoved _Three:l,fp3
2759: fmoved fp3,_Z:l
2760: L470:
2761: fmoved _Z:l,fp3
2762: fmoved fp3,_X:l
2763: movl #1,_I:l
2764: jbsr _SR3980:l
2765: L473:
2766: fmoved _Z:l,fp3
2767: faddd _Two:l,fp3
2768: fmoved fp3,_Z:l
2769: L474:
2770: fmoved _Z:l,fp3
2771: fmoved _Three:l,fp2
2772: fdivx fp2,fp3
2773: fmoved fp3,sp@(0)
2774: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2775: fmovex fp2,fp3; fmulx fp0,fp3
2776: fcmpd _Z:l,fp3; fjeq L473
2777: L471:
2778: fmoved _Z:l,fp3
2779: fmoved _Eight:l,fp2
2780: fmuld _Three:l,fp2
2781: fcmpx fp2,fp3; fjlt L470
2782: tstl _N:l; jle L476
2783: movl #L478,sp@
2784: jbsr _printf:l
2785: movl #L479,sp@
2786: jbsr _printf:l
2787: L476:
2788: jbsr _PrintIfNPositive:l
2789: movl _N:l,d4
2790: addl _N1:l,d4
2791: movl d4,_N:l
2792: tstl _N:l; jne L480
2793: movl #L482,sp@
2794: jbsr _printf:l
2795: L480:
2796: tstl _N:l; jle L483
2797: jbsr _Pause:l
2798: jra L484
2799: L483:
2800: movl #L47,sp@
2801: jbsr _printf:l
2802: L484:
2803: movl #110,_Milestone:l
2804: movl #L485,sp@
2805: jbsr _printf:l
2806: fmoved _U1:l,fp3
2807: fmoved fp3,_D:l
2808: fmoved _Precision:l,fp3
2809: fmoved fp3,sp@(0)
2810: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2811: fcmpx fp0,fp3; fjeq L486
2812: fmoved _BInvrse:l,fp3
2813: fmoved fp3,_D:l
2814: fmoved _Precision:l,fp3
2815: fmoved fp3,_X:l
2816: L488:
2817: fmoved _D:l,fp3
2818: fmuld _BInvrse:l,fp3
2819: fmoved fp3,_D:l
2820: fmoved _X:l,fp3
2821: fsubd _One:l,fp3
2822: fmoved fp3,_X:l
2823: L489:
2824: fmoved _X:l,fp3
2825: fcmpd _Zero:l,fp3; fjgt L488
2826: L486:
2827: fmoved _One:l,fp3
2828: fmoved fp3,_Y:l
2829: fmoved _D:l,fp3
2830: fmoved fp3,_Z:l
2831: L491:
2832: fmoved _Y:l,fp3
2833: fmoved fp3,_C:l
2834: fmoved _Z:l,fp3
2835: fmoved fp3,_Y:l
2836: fmoved _Y:l,fp3
2837: fmovex fp3,fp2; fmulx fp3,fp2
2838: fmoved fp2,_Z:l
2839: L492:
2840: fmoved _Y:l,fp3
2841: fmoved _Z:l,fp2
2842: fcmpx fp2,fp3; fjle L494
2843: fmovex fp2,fp3; faddx fp2,fp3
2844: fcmpx fp2,fp3; fjgt L491
2845: L494:
2846: fmoved _C:l,fp3
2847: fmoved fp3,_Y:l
2848: fmoved _Y:l,fp3
2849: fmuld _D:l,fp3
2850: fmoved fp3,_Z:l
2851: L495:
2852: fmoved _Y:l,fp3
2853: fmoved fp3,_C:l
2854: fmoved _Z:l,fp3
2855: fmoved fp3,_Y:l
2856: fmoved _Y:l,fp3
2857: fmuld _D:l,fp3
2858: fmoved fp3,_Z:l
2859: L496:
2860: fmoved _Y:l,fp3
2861: fmoved _Z:l,fp2
2862: fcmpx fp2,fp3; fjle L498
2863: fmovex fp2,fp3; faddx fp2,fp3
2864: fcmpx fp2,fp3; fjgt L495
2865: L498:
2866: fmoved _Radix:l,fp3
2867: fcmpd _Two:l,fp3; fjge L499
2868: fmoved _Two:l,fp3
2869: fmoved fp3,_HInvrse:l
2870: jra L500
2871: L499:
2872: fmoved _Radix:l,fp3
2873: fmoved fp3,_HInvrse:l
2874: L500:
2875: fmoved _One:l,fp3
2876: fmovex fp3,fp2; fdivd _HInvrse:l,fp2
2877: fmoved fp2,_H:l
2878: fmoved _C:l,fp2
2879: fdivx fp2,fp3
2880: fmoved fp3,_CInvrse:l
2881: fmoved fp2,_E0:l
2882: fmoved _E0:l,fp3
2883: fmuld _H:l,fp3
2884: fmoved fp3,_Z:l
2885: L501:
2886: fmoved _E0:l,fp3
2887: fmoved fp3,_Y:l
2888: fmoved _Z:l,fp3
2889: fmoved fp3,_E0:l
2890: fmoved _E0:l,fp3
2891: fmuld _H:l,fp3
2892: fmoved fp3,_Z:l
2893: L502:
2894: fmoved _E0:l,fp3
2895: fmoved _Z:l,fp2
2896: fcmpx fp2,fp3; fjle L504
2897: fmovex fp2,fp3; faddx fp2,fp3
2898: fcmpx fp2,fp3; fjgt L501
2899: L504:
2900: fmoved _E0:l,fp3
2901: fmoved fp3,_UfThold:l
2902: fmoved _Zero:l,fp3
2903: fmoved fp3,_E1:l
2904: fmoved fp3,_Q:l
2905: fmoved _U2:l,fp3
2906: fmoved fp3,_E9:l
2907: fmoved _One:l,fp3
2908: faddd _E9:l,fp3
2909: fmoved fp3,_S:l
2910: fmoved _C:l,fp3
2911: fmovex fp3,fp2; fmuld _S:l,fp2
2912: fmoved fp2,_D:l
2913: fmoved _D:l,fp2
2914: fcmpx fp3,fp2; fjgt L505
2915: fmoved _Radix:l,fp3
2916: fmuld _U2:l,fp3
2917: fmoved fp3,_E9:l
2918: fmoved _One:l,fp3
2919: faddd _E9:l,fp3
2920: fmoved fp3,_S:l
2921: fmoved _C:l,fp3
2922: fmovex fp3,fp2; fmuld _S:l,fp2
2923: fmoved fp2,_D:l
2924: fmoved _D:l,fp2
2925: fcmpx fp3,fp2; fjgt L506
2926: movl #0,sp@
2927: movl #L509,sp@(4)
2928: jbsr _BadCond:l
2929: fmoved _E0:l,fp3
2930: fmoved fp3,_Underflow:l
2931: fmoved _Zero:l,fp3
2932: fmoved fp3,_Y1:l
2933: fmoved _Z:l,fp3
2934: fmoved fp3,_PseudoZero:l
2935: jbsr _Pause:l
2936: jra L506
2937: L505:
2938: fmoved _D:l,fp3
2939: fmoved fp3,_Underflow:l
2940: fmoved _Underflow:l,fp3
2941: fmuld _H:l,fp3
2942: fmoved fp3,_PseudoZero:l
2943: fmoved _Zero:l,fp3
2944: fmoved fp3,_UfThold:l
2945: L510:
2946: fmoved _Underflow:l,fp3
2947: fmoved fp3,_Y1:l
2948: fmoved _PseudoZero:l,fp3
2949: fmoved fp3,_Underflow:l
2950: fmoved _E1:l,fp3
2951: fmovex fp3,fp2; faddx fp3,fp2
2952: fcmpx fp3,fp2; fjgt L513
2953: fmoved _Underflow:l,fp3
2954: fmuld _HInvrse:l,fp3
2955: fmoved fp3,_Y2:l
2956: fmoved _Y1:l,fp3
2957: fsubd _Y2:l,fp3
2958: fmoved fp3,sp@(0)
2959: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
2960: fmoved fp0,_E1:l
2961: fmoved _Y1:l,fp3
2962: fmoved fp3,_Q:l
2963: fmoved _UfThold:l,fp2
2964: fcmpd _Zero:l,fp2; fjne L515
2965: fcmpd _Y2:l,fp3; fjeq L515
2966: fmoved _Y1:l,fp3
2967: fmoved fp3,_UfThold:l
2968: L515:
2969: L513:
2970: fmoved _PseudoZero:l,fp3
2971: fmuld _H:l,fp3
2972: fmoved fp3,_PseudoZero:l
2973: L511:
2974: fmoved _Underflow:l,fp3
2975: fmoved _PseudoZero:l,fp2
2976: fcmpx fp2,fp3; fjle L517
2977: fmovex fp2,fp3; faddx fp2,fp3
2978: fcmpx fp2,fp3; fjgt L510
2979: L517:
2980: L506:
2981: fmoved _PseudoZero:l,fp3
2982: fcmpd _Zero:l,fp3; fjeq L518
2983: movl #L47,sp@
2984: jbsr _printf:l
2985: fmoved _PseudoZero:l,fp3
2986: fmoved fp3,_Z:l
2987: fcmpd _Zero:l,fp3; fjgt L520
2988: movl #0,sp@
2989: movl #L522,sp@(4)
2990: jbsr _BadCond:l
2991: movl #L523,sp@
2992: jbsr _printf:l
2993: movl #L524,sp@
2994: fmoved _PseudoZero:l,fp3
2995: fmoved fp3,sp@(4)
2996: jbsr _printf:l
2997: fnegd _PseudoZero:l,fp3
2998: fmoved fp3,_X:l
2999: fmoved _X:l,fp3
3000: fcmpd _Zero:l,fp3; fjgt L521
3001: movl #L527,sp@
3002: jbsr _printf:l
3003: movl #L528,sp@
3004: fmoved _X:l,fp3
3005: fmoved fp3,sp@(4)
3006: jbsr _printf:l
3007: jra L521
3008: L520:
3009: movl #3,sp@
3010: movl #L529,sp@(4)
3011: jbsr _BadCond:l
3012: movl #L530,sp@
3013: fmoved _PseudoZero:l,fp3
3014: fmoved fp3,sp@(4)
3015: jbsr _printf:l
3016: L521:
3017: jbsr _TstPtUf:l
3018: L518:
3019: movl #120,_Milestone:l
3020: fmoved _CInvrse:l,fp3
3021: fmovex fp3,fp2; fmuld _Y:l,fp2
3022: fmuld _Y1:l,fp3
3023: fcmpx fp3,fp2; fjle L531
3024: fmoved _H:l,fp3
3025: fmuld _S:l,fp3
3026: fmoved fp3,_S:l
3027: fmoved _Underflow:l,fp3
3028: fmoved fp3,_E0:l
3029: L531:
3030: fmoved _E1:l,fp3
3031: fcmpd _Zero:l,fp3; fjeq L533
3032: fcmpd _E0:l,fp3; fjeq L533
3033: movl #2,sp@
3034: movl #L160,sp@(4)
3035: jbsr _BadCond:l
3036: fmoved _E1:l,fp3
3037: fcmpd _E0:l,fp3; fjge L535
3038: movl #L537,sp@
3039: jbsr _printf:l
3040: movl #L538,sp@
3041: jbsr _printf:l
3042: fmoved _PseudoZero:l,fp3
3043: fcmpd _Zero:l,fp3; fjne L536
3044: fmoved _E1:l,fp3
3045: fmoved fp3,_E0:l
3046: jra L536
3047: L535:
3048: movl #L541,sp@
3049: jbsr _printf:l
3050: movl #L542,sp@
3051: jbsr _printf:l
3052: L536:
3053: L533:
3054: movl #L543,sp@
3055: fmoved _E0:l,fp3
3056: fmoved fp3,sp@(4)
3057: jbsr _printf:l
3058: fmoved _E0:l,fp3
3059: fmoved fp3,_Z:l
3060: jbsr _TstPtUf:l
3061: fmoved _E0:l,fp3
3062: fmoved fp3,_Underflow:l
3063: cmpl #1,_N:l; jne L544
3064: fmoved _Y:l,fp3
3065: fmoved fp3,_Underflow:l
3066: L544:
3067: movl #4,_I:l
3068: fmoved _E1:l,fp3
3069: fcmpd _Zero:l,fp3; fjne L546
3070: movl #3,_I:l
3071: L546:
3072: fmoved _UfThold:l,fp3
3073: fcmpd _Zero:l,fp3; fjne L548
3074: subql #2,_I:l
3075: L548:
3076: movl #1,_UfNGrad:l
3077: movl _I:l,d4
3078: cmpl #1,d4; jlt L550
3079: cmpl #4,d4; jgt L550
3080: movl @(L574-4:l,d4:l:4),a2
3081: jra a2@
3082: L574:
3083: .align 2; .long L552
3084: .align 2; .long L559
3085: .align 2; .long L564
3086: .align 2; .long L565
3087: L552:
3088: fmoved _Underflow:l,fp3
3089: fmoved fp3,_UfThold:l
3090: fmoved _CInvrse:l,fp3
3091: fmovex fp3,fp2; fmuld _Q:l,fp2
3092: fmuld _Y:l,fp3
3093: fmuld _S:l,fp3
3094: fcmpx fp3,fp2; fjeq L553
3095: fmoved _Y:l,fp3
3096: fmoved fp3,_UfThold:l
3097: movl #0,sp@
3098: movl #L555,sp@(4)
3099: jbsr _BadCond:l
3100: movl #L556,sp@
3101: fmoved _UfThold:l,fp3
3102: fmoved fp3,sp@(4)
3103: jbsr _printf:l
3104: movl #L557,sp@
3105: fmoved _C:l,fp3
3106: fmoved fp3,sp@(4)
3107: jbsr _printf:l
3108: movl #L558,sp@
3109: jbsr _printf:l
3110: L553:
3111: jbsr _Pause:l
3112: jra L551
3113: L559:
3114: movl #0,sp@
3115: movl #L560,sp@(4)
3116: jbsr _BadCond:l
3117: movl #L561,sp@
3118: jbsr _printf:l
3119: movl #L562,sp@
3120: fmoved _Q:l,fp3
3121: fmoved fp3,sp@(4)
3122: fmoved _Y2:l,fp3
3123: fmoved fp3,sp@(12)
3124: jbsr _printf:l
3125: fmoved _Q:l,fp3
3126: fsubd _Y2:l,fp3
3127: fmoved fp3,sp@(0)
3128: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3129: movl #L563,sp@
3130: fmoved fp0,sp@(4)
3131: jbsr _printf:l
3132: fmoved _Q:l,fp3
3133: fmoved fp3,_UfThold:l
3134: jra L551
3135: L564:
3136: fmoved _X:l,fp3
3137: fmoved fp3,_X:l
3138: jra L551
3139: L565:
3140: fmoved _Q:l,fp3
3141: fmoved _UfThold:l,fp2
3142: fcmpx fp2,fp3; fjne L566
3143: fmoved _E1:l,fp3
3144: fcmpd _E0:l,fp3; fjne L566
3145: fdivd _E9:l,fp3
3146: fsubx fp3,fp2
3147: fmoved fp2,sp@(0)
3148: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3149: fcmpd _E1:l,fp0; fjgt L566
3150: clrl _UfNGrad:l
3151: movl #L568,sp@
3152: jbsr _printf:l
3153: movl #L569,sp@
3154: jbsr _printf:l
3155: fmoved _E0:l,fp3
3156: fmoved _CInvrse:l,fp2
3157: fmulx fp2,fp3
3158: fmoved fp3,_Y:l
3159: fmoved _Y:l,fp3
3160: fmoved _OneAndHalf:l,fp1
3161: fmoved _U2:l,fp0
3162: faddx fp0,fp1
3163: fmulx fp1,fp3
3164: fmoved fp3,_Y:l
3165: fmoved _One:l,fp3
3166: faddx fp0,fp3
3167: fmulx fp3,fp2
3168: fmoved fp2,_X:l
3169: fmoved _Y:l,fp3
3170: fdivd _X:l,fp3
3171: fmoved fp3,_Y:l
3172: fmoved _Y:l,fp3
3173: fcmpd _E0:l,fp3; fjne L571
3174: movl #1,d7
3175: jra L572
3176: L571:
3177: clrl d7
3178: L572:
3179: movl d7,_IEEE:l
3180: L566:
3181: L550:
3182: L551:
3183: tstl _UfNGrad:l; jeq L576
3184: movl #L47,sp@
3185: jbsr _printf:l
3186: movl #_sigfpe,_sigsave:l
3187: movl #_ovfl_buf,sp@
3188: jbsr _setjmp:l
3189: tstl d0; jeq L578
3190: movl #L580,sp@
3191: jbsr _printf:l
3192: fmoved _H:l,fp3
3193: fmovex fp3,fp2; faddx fp3,fp2
3194: fmoved fp2,_R:l
3195: jra L579
3196: L578:
3197: fmoved _Underflow:l,fp3
3198: fdivd _UfThold:l,fp3
3199: fmoved fp3,sp@(0)
3200: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3201: fmoved fp0,_R:l
3202: L579:
3203: clrl _sigsave:l
3204: fmoved _R:l,fp3
3205: fcmpd _H:l,fp3; fjgt L581
3206: fmoved _R:l,fp3
3207: fmovex fp3,fp2; fmuld _UfThold:l,fp2
3208: fmoved fp2,_Z:l
3209: fmoved _Z:l,fp2
3210: fmovex fp2,fp7
3211: fmoved _One:l,fp1
3212: fmoved _H:l,fp0
3213: fmulx fp0,fp3
3214: fmovex fp1,fp2; faddx fp0,fp2
3215: fmulx fp2,fp3
3216: fmovex fp1,fp2; faddx fp3,fp2
3217: fmovex fp7,fp3; fmulx fp2,fp3
3218: fmoved fp3,_X:l
3219: jra L582
3220: L581:
3221: fmoved _UfThold:l,fp3
3222: fmoved fp3,_Z:l
3223: fmoved _Z:l,fp3
3224: fmovex fp3,fp7
3225: fmoved _One:l,fp2
3226: fmoved _H:l,fp1
3227: fmovex fp1,fp0; fmulx fp1,fp0
3228: fmovex fp2,fp3; faddx fp1,fp3
3229: fmovex fp0,fp1; fmulx fp3,fp1
3230: fmovex fp2,fp3; faddx fp1,fp3
3231: fmovex fp7,fp2; fmulx fp3,fp2
3232: fmoved fp2,_X:l
3233: L582:
3234: fmoved _X:l,fp3
3235: fmoved _Z:l,fp2
3236: fcmpx fp2,fp3; fjeq L583
3237: fsubx fp2,fp3
3238: fcmpd _Zero:l,fp3; fjne L583
3239: movl #3,sp@
3240: movl #L160,sp@(4)
3241: jbsr _BadCond:l
3242: movl #L585,sp@
3243: fmoved _X:l,fp3
3244: fmoved fp3,sp@(4)
3245: fmoved _Z:l,fp3
3246: fmoved fp3,sp@(12)
3247: jbsr _printf:l
3248: fmoved _X:l,fp3
3249: fsubd _Z:l,fp3
3250: fmoved fp3,_Z9:l
3251: movl #L586,sp@
3252: fmoved _Z9:l,fp3
3253: fmoved fp3,sp@(4)
3254: jbsr _printf:l
3255: movl #L587,sp@
3256: jbsr _printf:l
3257: movl #L588,sp@
3258: jbsr _printf:l
3259: movl #L589,sp@
3260: jbsr _printf:l
3261: movl #L590,sp@
3262: jbsr _printf:l
3263: movl #L591,sp@
3264: jbsr _printf:l
3265: movl #L592,sp@
3266: jbsr _printf:l
3267: movl #_sigfpe,_sigsave:l
3268: movl #_ovfl_buf,sp@
3269: jbsr _setjmp:l
3270: tstl d0; jeq L593
3271: movl #L595,sp@
3272: jbsr _printf:l
3273: jra L594
3274: L593:
3275: movl #L596,sp@
3276: fmoved _X:l,fp3
3277: fdivd _Z:l,fp3
3278: fmoved _Half:l,fp2
3279: fsubx fp2,fp3
3280: fsubx fp2,fp3
3281: fmoved fp3,sp@(4)
3282: jbsr _printf:l
3283: L594:
3284: clrl _sigsave:l
3285: L583:
3286: L576:
3287: movl #L597,sp@
3288: fmoved _UfThold:l,fp3
3289: fmoved fp3,sp@(4)
3290: movl #L598,sp@(12)
3291: jbsr _printf:l
3292: movl #L599,sp@
3293: jbsr _printf:l
3294: movl #L600,sp@
3295: jbsr _printf:l
3296: fmoved _U1:l,fp3
3297: fmovex fp3,fp2; fmulx fp3,fp2
3298: fmoved fp2,_Y2:l
3299: fmoved _Y2:l,fp2
3300: fmovex fp2,fp1; fmulx fp2,fp1
3301: fmoved fp1,_Y:l
3302: fmoved _Y:l,fp2
3303: fmulx fp3,fp2
3304: fmoved fp2,_Y2:l
3305: fmoved _Y2:l,fp3
3306: fcmpd _UfThold:l,fp3; fjgt L601
3307: fmoved _Y:l,fp3
3308: fcmpd _E0:l,fp3; fjle L603
3309: movl #2,sp@
3310: movl #L160,sp@(4)
3311: jbsr _BadCond:l
3312: movl #5,_I:l
3313: jra L604
3314: L603:
3315: movl #1,sp@
3316: movl #L160,sp@(4)
3317: jbsr _BadCond:l
3318: movl #4,_I:l
3319: L604:
3320: movl #L605,sp@
3321: movl _I:l,sp@(4)
3322: jbsr _printf:l
3323: L601:
3324: movl #130,_Milestone:l
3325: fmoved _UfThold:l,fp3
3326: fmoved fp3,sp@(0)
3327: jbsr _log:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3328: fmovex fp0,fp7
3329: fmoved _HInvrse:l,fp3
3330: fmoved fp3,sp@(0)
3331: jbsr _log:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3332: fmoved _Half:l,fp3
3333: fmoved _TwoForty:l,fp2
3334: fmulx fp7,fp2
3335: fdivx fp0,fp2
3336: fsubx fp2,fp3
3337: fmoved fp3,sp@(0)
3338: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3339: fnegx fp0,fp3
3340: fdivd _TwoForty:l,fp3
3341: fmoved fp3,_Y:l
3342: fmoved _Y:l,fp3
3343: fmovex fp3,fp2; faddx fp3,fp2
3344: fmoved fp2,_Y2:l
3345: movl #L606,sp@
3346: jbsr _printf:l
3347: movl #L607,sp@
3348: fmoved _HInvrse:l,fp3
3349: fmoved fp3,sp@(4)
3350: fmoved _Y:l,fp3
3351: fmoved fp3,sp@(12)
3352: jbsr _printf:l
3353: movl #L608,sp@
3354: fmoved _HInvrse:l,fp3
3355: fmoved fp3,sp@(4)
3356: fmoved _Y:l,fp3
3357: fmoved fp3,sp@(12)
3358: jbsr _printf:l
3359: fmoved _HInvrse:l,fp3
3360: fmoved fp3,sp@(0)
3361: fmoved _Y2:l,fp3
3362: fmoved fp3,sp@(8)
3363: jbsr _pow:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3364: fmoved fp0,_V9:l
3365: movl #L609,sp@
3366: fmoved _V9:l,fp3
3367: fmoved fp3,sp@(4)
3368: jbsr _printf:l
3369: fmoved _V9:l,fp3
3370: fcmpd _Zero:l,fp3; fjlt L612
3371: fmoved _Radix:l,fp2
3372: fmovex fp2,fp1; faddx fp2,fp1
3373: fmovex fp1,fp2; faddd _E9:l,fp2
3374: fmuld _UfThold:l,fp2
3375: fcmpx fp2,fp3; fjle L610
3376: L612:
3377: movl #1,sp@
3378: movl #L613,sp@(4)
3379: jbsr _BadCond:l
3380: movl #L614,sp@
3381: fmoved _UfThold:l,fp3
3382: fmoved fp3,sp@(4)
3383: jbsr _printf:l
3384: jra L611
3385: L610:
3386: fmoved _V9:l,fp3
3387: fmoved _UfThold:l,fp2
3388: fmoved _One:l,fp1
3389: faddd _E9:l,fp1
3390: fmulx fp1,fp2
3391: fcmpx fp2,fp3; fjgt L615
3392: movl #L617,sp@
3393: jbsr _printf:l
3394: jra L616
3395: L615:
3396: movl #2,sp@
3397: movl #L613,sp@(4)
3398: jbsr _BadCond:l
3399: movl #L614,sp@
3400: fmoved _UfThold:l,fp3
3401: fmoved fp3,sp@(4)
3402: jbsr _printf:l
3403: L616:
3404: L611:
3405: movl #140,_Milestone:l
3406: movl #L47,sp@
3407: jbsr _printf:l
3408: fmoved _Zero:l,fp3
3409: fmoved fp3,_X:l
3410: movl #2,_I:l
3411: fmoved _Two:l,fp2
3412: fmuld _Three:l,fp2
3413: fmoved fp2,_Y:l
3414: fmoved fp3,_Q:l
3415: clrl _N:l
3416: L618:
3417: fmoved _X:l,fp3
3418: fmoved fp3,_Z:l
3419: addql #1,_I:l
3420: fmoved _Y:l,fp3
3421: movl _I:l,d4
3422: movl d4,d3; addl d4,d3
3423: fdivl d3,fp3
3424: fmoved fp3,_Y:l
3425: fmoved _Y:l,fp3
3426: faddd _Q:l,fp3
3427: fmoved fp3,_R:l
3428: fmoved _Z:l,fp3
3429: fmoved _R:l,fp2
3430: fmovex fp3,fp1; faddx fp2,fp1
3431: fmoved fp1,_X:l
3432: fsubd _X:l,fp3
3433: faddx fp2,fp3
3434: fmoved fp3,_Q:l
3435: L619:
3436: fmoved _X:l,fp3
3437: fcmpd _Z:l,fp3; fjgt L618
3438: fmoved _OneAndHalf:l,fp3
3439: fmoved _One:l,fp2
3440: fdivd _Eight:l,fp2
3441: fmovex fp3,fp1; faddx fp2,fp1
3442: fmoved _X:l,fp2
3443: fmuld _ThirtyTwo:l,fp3
3444: fdivx fp3,fp2
3445: fmovex fp1,fp3; faddx fp2,fp3
3446: fmoved fp3,_Z:l
3447: fmoved _Z:l,fp3
3448: fmovex fp3,fp2; fmulx fp3,fp2
3449: fmoved fp2,_X:l
3450: fmoved _X:l,fp3
3451: fmovex fp3,fp2; fmulx fp3,fp2
3452: fmoved fp2,_Exp2:l
3453: fmoved _F9:l,fp3
3454: fmoved fp3,_X:l
3455: fmoved _X:l,fp3
3456: fsubd _U1:l,fp3
3457: fmoved fp3,_Y:l
3458: movl #L621,sp@
3459: fmoved _Exp2:l,fp3
3460: fmoved fp3,sp@(4)
3461: jbsr _printf:l
3462: movl #1,_I:l
3463: L622:
3464: fmoved _X:l,fp3
3465: fmovex fp3,fp7
3466: fmoved _BInvrse:l,fp2
3467: fmovex fp3,fp1; fsubx fp2,fp1
3468: fmoved fp1,_Z:l
3469: fmoved _One:l,fp1
3470: fmovex fp3,fp0; faddx fp1,fp0
3471: fmoved _Z:l,fp3
3472: fsubx fp2,fp1
3473: fsubx fp1,fp3
3474: fmovex fp0,fp2; fdivx fp3,fp2
3475: fmoved fp2,_Z:l
3476: fmoved fp7,sp@(0)
3477: fmoved _Z:l,fp3
3478: fmoved fp3,sp@(8)
3479: jbsr _pow:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3480: fmovex fp0,fp3; fsubd _Exp2:l,fp3
3481: fmoved fp3,_Q:l
3482: fmoved _Q:l,fp3
3483: fmoved fp3,sp@(0)
3484: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3485: fmoved _TwoForty:l,fp3
3486: fmuld _U2:l,fp3
3487: fcmpx fp3,fp0; fjle L626
3488: movl #1,_N:l
3489: fmoved _X:l,fp3
3490: fmoved _BInvrse:l,fp2
3491: fsubx fp2,fp3
3492: fmoved _One:l,fp1
3493: fsubx fp2,fp1
3494: fsubx fp1,fp3
3495: fmoved fp3,_V9:l
3496: movl #2,sp@
3497: movl #L628,sp@(4)
3498: jbsr _BadCond:l
3499: fmoved _X:l,fp3
3500: fmoved fp3,sp@(0)
3501: fmoved _Z:l,fp3
3502: fmoved fp3,sp@(8)
3503: jbsr _pow:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3504: movl #L629,sp@
3505: fmoved fp0,sp@(4)
3506: jbsr _printf:l
3507: movl #L630,sp@
3508: fmoved _V9:l,fp3
3509: fmoved fp3,sp@(4)
3510: fmoved _Z:l,fp3
3511: fmoved fp3,sp@(12)
3512: jbsr _printf:l
3513: movl #L631,sp@
3514: fmoved _Q:l,fp3
3515: fmoved fp3,sp@(4)
3516: jbsr _printf:l
3517: movl #L632,sp@
3518: jbsr _printf:l
3519: movl #L633,sp@
3520: jbsr _printf:l
3521: jra L624
3522: L626:
3523: fmoved _Y:l,fp3
3524: fmovex fp3,fp2; fsubd _X:l,fp2
3525: fmuld _Two:l,fp2
3526: faddx fp3,fp2
3527: fmoved fp2,_Z:l
3528: fmoved fp3,_X:l
3529: fmoved _Z:l,fp3
3530: fmoved fp3,_Y:l
3531: fmoved _One:l,fp3
3532: fmoved _X:l,fp2
3533: fsubd _F9:l,fp2
3534: fmovex fp2,fp1; fmulx fp2,fp1
3535: fmovex fp3,fp2; faddx fp1,fp2
3536: fmoved fp2,_Z:l
3537: fmoved _Z:l,fp2
3538: fcmpx fp3,fp2; fjle L634
3539: movl _I:l,d4
3540: cmpl _NoTrials:l,d4; jge L634
3541: addql #1,_I:l
3542: jra L622
3543: L634:
3544: fmoved _X:l,fp3
3545: fcmpd _One:l,fp3; fjle L636
3546: tstl _N:l; jne L624
3547: movl #L640,sp@
3548: jbsr _printf:l
3549: jra L624
3550: L636:
3551: fmoved _One:l,fp3
3552: fmoved _U2:l,fp2
3553: faddx fp2,fp3
3554: fmoved fp3,_X:l
3555: fmovex fp2,fp3; faddx fp2,fp3
3556: fmoved fp3,_Y:l
3557: fmoved _Y:l,fp3
3558: faddd _X:l,fp3
3559: fmoved fp3,_Y:l
3560: movl #1,_I:l
3561: jra L622
3562: L624:
3563: movl #150,_Milestone:l
3564: movl #L641,sp@
3565: jbsr _printf:l
3566: clrl _N:l
3567: fmoved _A1:l,fp3
3568: fmoved fp3,_Z:l
3569: fmoved _C:l,fp3
3570: fmoved fp3,sp@(0)
3571: jbsr _log:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3572: fmovex fp0,fp7
3573: fmoved _A1:l,fp3
3574: fmoved fp3,sp@(0)
3575: jbsr _log:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3576: fmoved _Half:l,fp3
3577: fmovex fp7,fp2; fdivx fp0,fp2
3578: fsubx fp2,fp3
3579: fmoved fp3,sp@(0)
3580: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3581: fmoved fp0,_Q:l
3582: clrl _Break:l
3583: L642:
3584: fmoved _CInvrse:l,fp3
3585: fmoved fp3,_X:l
3586: fmoved _Z:l,fp3
3587: fmoved fp3,sp@(0)
3588: fmoved _Q:l,fp3
3589: fmoved fp3,sp@(8)
3590: jbsr _pow:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3591: fmoved fp0,_Y:l
3592: jbsr _IsYeqX:l
3593: fnegd _Q:l,fp3
3594: fmoved fp3,_Q:l
3595: fmoved _C:l,fp3
3596: fmoved fp3,_X:l
3597: fmoved _Z:l,fp3
3598: fmoved fp3,sp@(0)
3599: fmoved _Q:l,fp3
3600: fmoved fp3,sp@(8)
3601: jbsr _pow:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3602: fmoved fp0,_Y:l
3603: jbsr _IsYeqX:l
3604: fmoved _Z:l,fp3
3605: fcmpd _One:l,fp3; fjge L645
3606: movl #1,_Break:l
3607: jra L646
3608: L645:
3609: fmoved _AInvrse:l,fp3
3610: fmoved fp3,_Z:l
3611: L646:
3612: L643:
3613: tstl _Break:l; jeq L642
3614: jbsr _PrintIfNPositive:l
3615: tstl _N:l; jne L647
3616: movl #L649,sp@
3617: jbsr _printf:l
3618: L647:
3619: movl #L47,sp@
3620: jbsr _printf:l
3621: movl #160,_Milestone:l
3622: jbsr _Pause:l
3623: movl #L650,sp@
3624: jbsr _printf:l
3625: movl #L651,sp@
3626: jbsr _printf:l
3627: fnegd _CInvrse:l,fp3
3628: fmoved fp3,_Y:l
3629: fmoved _HInvrse:l,fp3
3630: fmuld _Y:l,fp3
3631: fmoved fp3,_V9:l
3632: movl #_sigfpe,_sigsave:l
3633: movl #_ovfl_buf,sp@
3634: jbsr _setjmp:l
3635: tstl d0; jeq L652
3636: clrl _I:l
3637: fmoved _Y:l,fp3
3638: fmoved fp3,_V9:l
3639: jra L654
3640: L652:
3641: L655:
3642: fmoved _Y:l,fp3
3643: fmoved fp3,_V:l
3644: fmoved _V9:l,fp3
3645: fmoved fp3,_Y:l
3646: fmoved _HInvrse:l,fp3
3647: fmuld _Y:l,fp3
3648: fmoved fp3,_V9:l
3649: L656:
3650: fmoved _V9:l,fp3
3651: fcmpd _Y:l,fp3; fjlt L655
3652: movl #1,_I:l
3653: L654:
3654: clrl _sigsave:l
3655: fmoved _V9:l,fp3
3656: fmoved fp3,_Z:l
3657: movl #L658,sp@
3658: jbsr _printf:l
3659: movl #L659,sp@
3660: fmoved _Y:l,fp3
3661: fmoved fp3,sp@(4)
3662: jbsr _printf:l
3663: fmoved _Y:l,fp3
3664: fnegx fp3,fp2
3665: fmoved fp2,_V9:l
3666: fmoved _V9:l,fp2
3667: fmoved fp2,_V0:l
3668: fmoved _V:l,fp2
3669: fmovex fp2,fp1; fsubx fp3,fp1
3670: fmovex fp2,fp3; faddd _V0:l,fp3
3671: fcmpx fp3,fp1; fjne L660
3672: movl #L662,sp@
3673: jbsr _printf:l
3674: jra L661
3675: L660:
3676: movl #L663,sp@
3677: jbsr _printf:l
3678: movl #3,sp@
3679: movl #L664,sp@(4)
3680: jbsr _BadCond:l
3681: L661:
3682: fmoved _Z:l,fp3
3683: fcmpd _Y:l,fp3; fjeq L665
3684: movl #1,sp@
3685: movl #L160,sp@(4)
3686: jbsr _BadCond:l
3687: movl #L667,sp@
3688: fmoved _Y:l,fp3
3689: fmoved fp3,sp@(4)
3690: fmoved _Z:l,fp3
3691: fmoved fp3,sp@(12)
3692: jbsr _printf:l
3693: L665:
3694: tstl _I:l; jeq L668
3695: fmoved _V:l,fp3
3696: fmoved _HInvrse:l,fp2
3697: fmoved _U2:l,fp1
3698: fmovex fp1,fp7
3699: fmovex fp2,fp0; fmulx fp1,fp0
3700: fsubx fp2,fp0
3701: fmovex fp3,fp1; fmulx fp0,fp1
3702: fmoved fp1,_Y:l
3703: fmoved _Y:l,fp1
3704: fmoved _One:l,fp0
3705: fsubx fp2,fp0
3706: fmovex fp0,fp2; fmulx fp7,fp2
3707: fmulx fp3,fp2
3708: fmovex fp1,fp3; faddx fp2,fp3
3709: fmoved fp3,_Z:l
3710: fmoved _Z:l,fp3
3711: fcmpd _V0:l,fp3; fjge L670
3712: fmoved _Z:l,fp3
3713: fmoved fp3,_Y:l
3714: L670:
3715: fmoved _Y:l,fp3
3716: fcmpd _V0:l,fp3; fjge L672
3717: fmoved _Y:l,fp3
3718: fmoved fp3,_V:l
3719: L672:
3720: fmoved _V0:l,fp3
3721: fmovex fp3,fp2; fsubd _V:l,fp2
3722: fcmpx fp3,fp2; fjge L669
3723: fmoved _V0:l,fp3
3724: fmoved fp3,_V:l
3725: jra L669
3726: L668:
3727: fmoved _Y:l,fp3
3728: fmoved _HInvrse:l,fp2
3729: fmoved _U2:l,fp1
3730: fmovex fp1,fp7
3731: fmovex fp2,fp0; fmulx fp1,fp0
3732: fsubx fp2,fp0
3733: fmovex fp3,fp1; fmulx fp0,fp1
3734: fmoved fp1,_V:l
3735: fmoved _V:l,fp1
3736: fmoved _One:l,fp0
3737: fsubx fp2,fp0
3738: fmovex fp0,fp2; fmulx fp7,fp2
3739: fmulx fp3,fp2
3740: fmovex fp1,fp3; faddx fp2,fp3
3741: fmoved fp3,_V:l
3742: L669:
3743: movl #L676,sp@
3744: fmoved _V:l,fp3
3745: fmoved fp3,sp@(4)
3746: jbsr _printf:l
3747: tstl _I:l; jeq L677
3748: movl #L679,sp@
3749: fmoved _V0:l,fp3
3750: fmoved fp3,sp@(4)
3751: jbsr _printf:l
3752: jra L678
3753: L677:
3754: movl #L680,sp@
3755: jbsr _printf:l
3756: L678:
3757: fmoved _V:l,fp3
3758: fmuld _One:l,fp3
3759: fmoved fp3,_V9:l
3760: movl #L681,sp@
3761: fmoved _V9:l,fp3
3762: fmoved fp3,sp@(4)
3763: jbsr _printf:l
3764: fmoved _V:l,fp3
3765: fdivd _One:l,fp3
3766: fmoved fp3,_V9:l
3767: movl #L682,sp@
3768: fmoved _V9:l,fp3
3769: fmoved fp3,sp@(4)
3770: jbsr _printf:l
3771: movl #L683,sp@
3772: jbsr _printf:l
3773: movl #L684,sp@
3774: jbsr _printf:l
3775: movl #170,_Milestone:l
3776: fmoved _V:l,fp3
3777: fnegx fp3,fp2
3778: fcmpx fp3,fp2; fjge L689
3779: fmoved _V0:l,fp2
3780: fnegx fp2,fp1
3781: fcmpx fp2,fp1; fjge L689
3782: fmoved _UfThold:l,fp2
3783: fnegx fp2,fp1
3784: fcmpx fp3,fp1; fjge L689
3785: fcmpx fp3,fp2; fjlt L685
3786: L689:
3787: movl #0,sp@
3788: movl #L690,sp@(4)
3789: jbsr _BadCond:l
3790: movl #L691,sp@
3791: fmoved _V:l,fp3
3792: fmoved fp3,sp@(4)
3793: fmoved _V0:l,fp3
3794: fmoved fp3,sp@(12)
3795: fmoved _UfThold:l,fp3
3796: fmoved fp3,sp@(20)
3797: jbsr _printf:l
3798: L685:
3799: movl #175,_Milestone:l
3800: movl #L47,sp@
3801: jbsr _printf:l
3802: movl #1,_Indx:l
3803: L692:
3804: movl _Indx:l,d4
3805: cmpl #1,d4; jeq L698
3806: cmpl #2,d4; jeq L699
3807: cmpl #3,d4; jeq L700
3808: jra L696
3809: L698:
3810: fmoved _UfThold:l,fp3
3811: fmoved fp3,_Z:l
3812: jra L697
3813: L699:
3814: fmoved _E0:l,fp3
3815: fmoved fp3,_Z:l
3816: jra L697
3817: L700:
3818: fmoved _PseudoZero:l,fp3
3819: fmoved fp3,_Z:l
3820: L696:
3821: L697:
3822: fmoved _Z:l,fp3
3823: fcmpd _Zero:l,fp3; fjeq L701
3824: fmoved _Z:l,fp3
3825: fmoved fp3,sp@(0)
3826: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3827: fmoved fp0,_V9:l
3828: fmoved _V9:l,fp3
3829: fmovex fp3,fp2; fmulx fp3,fp2
3830: fmoved fp2,_Y:l
3831: fmoved _Y:l,fp3
3832: fmoved _One:l,fp2
3833: fmovex fp2,fp7
3834: fmoved _Radix:l,fp1
3835: fmuld _E9:l,fp1
3836: fmovex fp2,fp0; fsubx fp1,fp0
3837: fmovex fp3,fp2; fdivx fp0,fp2
3838: fmoved _Z:l,fp0
3839: fcmpx fp0,fp2; fjlt L705
3840: fmovex fp7,fp2; faddx fp1,fp2
3841: fmulx fp0,fp2
3842: fcmpx fp2,fp3; fjle L703
3843: L705:
3844: fmoved _V9:l,fp3
3845: fcmpd _U1:l,fp3; fjle L706
3846: movl #1,sp@
3847: movl #L160,sp@(4)
3848: jbsr _BadCond:l
3849: jra L707
3850: L706:
3851: movl #2,sp@
3852: movl #L160,sp@(4)
3853: jbsr _BadCond:l
3854: L707:
3855: movl #L708,sp@
3856: fmoved _Z:l,fp3
3857: fmoved fp3,sp@(4)
3858: jbsr _printf:l
3859: movl #L709,sp@
3860: fmoved _Y:l,fp3
3861: fmoved fp3,sp@(4)
3862: jbsr _printf:l
3863: L703:
3864: L701:
3865: L693:
3866: addql #1,_Indx:l
3867: cmpl #3,_Indx:l; jle L692
3868: movl #180,_Milestone:l
3869: movl #1,_Indx:l
3870: L710:
3871: cmpl #1,_Indx:l; jne L714
3872: fmoved _V:l,fp3
3873: fmoved fp3,_Z:l
3874: jra L715
3875: L714:
3876: fmoved _V0:l,fp3
3877: fmoved fp3,_Z:l
3878: L715:
3879: fmoved _Z:l,fp3
3880: fmoved fp3,sp@(0)
3881: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
3882: fmoved fp0,_V9:l
3883: fmoved _One:l,fp3
3884: fmoved _Radix:l,fp2
3885: fmoved _E9:l,fp1
3886: fmovex fp1,fp7
3887: fmovex fp2,fp0; fmulx fp1,fp0
3888: fmovex fp3,fp1; fsubx fp0,fp1
3889: fmoved _V9:l,fp0
3890: fmulx fp0,fp1
3891: fmoved fp1,_X:l
3892: fmovex fp0,fp1; fmuld _X:l,fp1
3893: fmoved fp1,_V9:l
3894: fmoved _V9:l,fp1
3895: fmoved _Two:l,fp0
3896: fmulx fp2,fp0
3897: fmovex fp0,fp2; fmulx fp7,fp2
3898: fsubx fp2,fp3
3899: fmoved _Z:l,fp2
3900: fmulx fp2,fp3
3901: fcmpx fp3,fp1; fjlt L718
3902: fcmpx fp2,fp1; fjle L716
3903: L718:
3904: fmoved _V9:l,fp3
3905: fmoved fp3,_Y:l
3906: fmoved _X:l,fp3
3907: fcmpd _W:l,fp3; fjge L719
3908: movl #1,sp@
3909: movl #L160,sp@(4)
3910: jbsr _BadCond:l
3911: jra L720
3912: L719:
3913: movl #2,sp@
3914: movl #L160,sp@(4)
3915: jbsr _BadCond:l
3916: L720:
3917: movl #L721,sp@
3918: fmoved _Z:l,fp3
3919: fmoved fp3,sp@(4)
3920: jbsr _printf:l
3921: movl #L722,sp@
3922: fmoved _Y:l,fp3
3923: fmoved fp3,sp@(4)
3924: jbsr _printf:l
3925: L716:
3926: L711:
3927: addql #1,_Indx:l
3928: cmpl #2,_Indx:l; jle L710
3929: movl #190,_Milestone:l
3930: jbsr _Pause:l
3931: fmoved _UfThold:l,fp3
3932: fmuld _V:l,fp3
3933: fmoved fp3,_X:l
3934: fmoved _Radix:l,fp3
3935: fmovex fp3,fp2; fmulx fp3,fp2
3936: fmoved fp2,_Y:l
3937: fmoved _X:l,fp3
3938: fmoved _Y:l,fp2
3939: fmovex fp3,fp1; fmulx fp2,fp1
3940: fcmpd _One:l,fp1; fjlt L725
3941: fcmpx fp2,fp3; fjle L723
3942: L725:
3943: fmoved _X:l,fp3
3944: fmoved _Y:l,fp2
3945: fmovex fp3,fp1; fmulx fp2,fp1
3946: fmoved _U1:l,fp0
3947: fcmpx fp0,fp1; fjlt L728
3948: fdivx fp0,fp2
3949: fcmpx fp2,fp3; fjle L726
3950: L728:
3951: movl #2,sp@
3952: movl #L729,sp@(4)
3953: jbsr _BadCond:l
3954: jra L727
3955: L726:
3956: movl #3,sp@
3957: movl #L160,sp@(4)
3958: jbsr _BadCond:l
3959: L727:
3960: movl #L730,sp@
3961: fmoved _X:l,fp3
3962: fmoved fp3,sp@(4)
3963: movl #L731,sp@(12)
3964: jbsr _printf:l
3965: L723:
3966: movl #200,_Milestone:l
3967: movl #1,_Indx:l
3968: L732:
3969: fmoved _F9:l,fp3
3970: fmoved fp3,_X:l
3971: movl _Indx:l,d4
3972: cmpl #2,d4; jlt L736
3973: cmpl #5,d4; jgt L736
3974: movl @(L743-8:l,d4:l:4),a2
3975: jra a2@
3976: L743:
3977: .align 2; .long L738
3978: .align 2; .long L739
3979: .align 2; .long L740
3980: .align 2; .long L741
3981: L738:
3982: fmoved _One:l,fp3
3983: faddd _U2:l,fp3
3984: fmoved fp3,_X:l
3985: jra L737
3986: L739:
3987: fmoved _V:l,fp3
3988: fmoved fp3,_X:l
3989: jra L737
3990: L740:
3991: fmoved _UfThold:l,fp3
3992: fmoved fp3,_X:l
3993: jra L737
3994: L741:
3995: fmoved _Radix:l,fp3
3996: fmoved fp3,_X:l
3997: L736:
3998: L737:
3999: fmoved _X:l,fp3
4000: fmoved fp3,_Y:l
4001: movl #_sigfpe,_sigsave:l
4002: movl #_ovfl_buf,sp@
4003: jbsr _setjmp:l
4004: tstl d0; jeq L745
4005: movl #L747,sp@
4006: fmoved _X:l,fp3
4007: fmoved fp3,sp@(4)
4008: jbsr _printf:l
4009: jra L746
4010: L745:
4011: fmoved _Y:l,fp3
4012: fdivd _X:l,fp3
4013: fmoved _Half:l,fp2
4014: fsubx fp2,fp3
4015: fsubx fp2,fp3
4016: fmoved fp3,_V9:l
4017: fmoved _V9:l,fp3
4018: fcmpd _Zero:l,fp3; fjne L748
4019: jra L733
4020: L748:
4021: fmoved _V9:l,fp3
4022: fnegd _U1:l,fp2
4023: fcmpx fp2,fp3; fjne L750
4024: cmpl #5,_Indx:l; jge L750
4025: movl #3,sp@
4026: movl #L160,sp@(4)
4027: jbsr _BadCond:l
4028: jra L751
4029: L750:
4030: movl #1,sp@
4031: movl #L160,sp@(4)
4032: jbsr _BadCond:l
4033: L751:
4034: movl #L752,sp@
4035: fmoved _X:l,fp3
4036: fmoved fp3,sp@(4)
4037: jbsr _printf:l
4038: movl #L753,sp@
4039: fmoved _V9:l,fp3
4040: fmoved fp3,sp@(4)
4041: jbsr _printf:l
4042: L746:
4043: clrl _sigsave:l
4044: L733:
4045: addql #1,_Indx:l
4046: cmpl #5,_Indx:l; jle L732
4047: movl #210,_Milestone:l
4048: fmoved _Zero:l,fp3
4049: fmoved fp3,_MyZero:l
4050: movl #L47,sp@
4051: jbsr _printf:l
4052: movl #L754,sp@
4053: jbsr _printf:l
4054: movl #_sigfpe,_sigsave:l
4055: movl #L755,sp@
4056: jbsr _printf:l
4057: movl #_ovfl_buf,sp@
4058: jbsr _setjmp:l
4059: tstl d0; jne L756
4060: movl #L758,sp@
4061: fmoved _One:l,fp3
4062: fdivd _MyZero:l,fp3
4063: fmoved fp3,sp@(4)
4064: jbsr _printf:l
4065: L756:
4066: clrl _sigsave:l
4067: movl #_sigfpe,_sigsave:l
4068: movl #L759,sp@
4069: jbsr _printf:l
4070: movl #_ovfl_buf,sp@
4071: jbsr _setjmp:l
4072: tstl d0; jne L760
4073: movl #L758,sp@
4074: fmoved _Zero:l,fp3
4075: fdivd _MyZero:l,fp3
4076: fmoved fp3,sp@(4)
4077: jbsr _printf:l
4078: L760:
4079: clrl _sigsave:l
4080: movl #220,_Milestone:l
4081: jbsr _Pause:l
4082: movl #L47,sp@
4083: jbsr _printf:l
4084: clrl d7
4085: L767:
4086: tstl @(_ErrCnt:l,d7:l:4); jeq L771
4087: movl #L773,sp@
4088: movl @(L762:l,d7:l:4),sp@(4)
4089: movl @(_ErrCnt:l,d7:l:4),sp@(8)
4090: jbsr _printf:l
4091: L771:
4092: L768:
4093: addql #1,d7
4094: cmpl #4,d7; jlt L767
4095: movl #L47,sp@
4096: jbsr _printf:l
4097: movl _ErrCnt:l,d4
4098: addl _ErrCnt+4:l,d4
4099: addl _ErrCnt+8:l,d4
4100: addl _ErrCnt+12:l,d4
4101: tstl d4; jle L774
4102: movl _ErrCnt:l,d4
4103: addl _ErrCnt+4:l,d4
4104: addl _ErrCnt+8:l,d4
4105: tstl d4; jne L779
4106: tstl _ErrCnt+12:l; jle L779
4107: movl #L784,sp@
4108: jbsr _printf:l
4109: movl #L785,sp@
4110: jbsr _printf:l
4111: L779:
4112: movl _ErrCnt:l,d4
4113: addl _ErrCnt+4:l,d4
4114: tstl d4; jne L786
4115: tstl _ErrCnt+8:l; jle L786
4116: movl #L790,sp@
4117: jbsr _printf:l
4118: movl #L791,sp@
4119: jbsr _printf:l
4120: L786:
4121: movl _ErrCnt:l,d4
4122: addl _ErrCnt+4:l,d4
4123: tstl d4; jle L792
4124: movl #L795,sp@
4125: jbsr _printf:l
4126: movl #L796,sp@
4127: jbsr _printf:l
4128: L792:
4129: tstl _ErrCnt:l; jle L775
4130: movl #L799,sp@
4131: jbsr _printf:l
4132: movl #L800,sp@
4133: jbsr _printf:l
4134: jra L775
4135: L774:
4136: movl #L801,sp@
4137: jbsr _printf:l
4138: cmpl #1,_RMult:l; jne L806
4139: cmpl #1,_RDiv:l; jne L806
4140: cmpl #1,_RAddSub:l; jne L806
4141: cmpl #1,_RSqrt:l; jeq L802
4142: L806:
4143: movl #L807,sp@
4144: jbsr _printf:l
4145: jra L803
4146: L802:
4147: fmoved _StickyBit:l,fp3
4148: fmoved _One:l,fp2
4149: fcmpx fp2,fp3; fjlt L808
4150: fmoved _Radix:l,fp3
4151: fmovex fp3,fp1; fsubd _Two:l,fp1
4152: fsubd _Nine:l,fp3
4153: fsubx fp2,fp3
4154: fmovex fp1,fp2; fmulx fp3,fp2
4155: fcmpd _Zero:l,fp2; fjne L808
4156: movl #L810,sp@
4157: jbsr _printf:l
4158: movl #L811,sp@
4159: jbsr _printf:l
4160: fmoved _Radix:l,fp3
4161: fmoved _Two:l,fp2
4162: fcmpx fp2,fp3; fjne L812
4163: fmoved _Precision:l,fp3
4164: fmoved _Four:l,fp1
4165: fmuld _Three:l,fp1
4166: fmulx fp2,fp1
4167: fmovex fp3,fp2; fsubx fp1,fp2
4168: fmoved _TwentySeven:l,fp1
4169: fsubx fp1,fp3
4170: fsubx fp1,fp3
4171: faddd _One:l,fp3
4172: fmulx fp3,fp2
4173: fcmpd _Zero:l,fp2; fjne L812
4174: movl #L814,sp@
4175: jbsr _printf:l
4176: jra L813
4177: L812:
4178: movl #L815,sp@
4179: jbsr _printf:l
4180: L813:
4181: tstl _IEEE:l; jeq L816
4182: movl #L818,sp@
4183: jbsr _printf:l
4184: jra L817
4185: L816:
4186: movl #L819,sp@
4187: jbsr _printf:l
4188: movl #L820,sp@
4189: jbsr _printf:l
4190: L817:
4191: L808:
4192: movl #L821,sp@
4193: jbsr _printf:l
4194: L803:
4195: L775:
4196: tstl _fpecount:l; jeq L822
4197: movl #L824,sp@
4198: movl _fpecount:l,sp@(4)
4199: jbsr _printf:l
4200: L822:
4201: movl #L825,sp@
4202: jbsr _printf:l
4203: L6:
4204: moveml a6@(-116),#0x498
4205: fmovemx a6@(-100),#0xff
4206: unlk a6
4207: rts
4208: .globl _Sign
4209: .align 2
4210: _Sign:link a6,#-36
4211: fmovemx #0x91,a6@(-36)
4212: fmoved a6@(8),fp3
4213: fcmpd L10:l,fp3; fjlt L863
4214: fmoved L11:l,fp7
4215: jra L864
4216: L863:
4217: fmoved L865:l,fp7
4218: L864:
4219: fmovex fp7,fp0
4220: L861:
4221: fmoved fp0,sp@-; movl sp@+,d0; movl sp@+,d1
4222: fmovemx a6@(-36),#0x91
4223: unlk a6
4224: rts
4225: .globl _Pause
4226: .align 2
4227: _Pause:link a6,#-8
4228: movl #L867,sp@
4229: movl _Milestone:l,sp@(4)
4230: jbsr _printf:l
4231: movl #L868,sp@
4232: movl _PageNo:l,sp@(4)
4233: jbsr _printf:l
4234: addql #1,_Milestone:l
4235: addql #1,_PageNo:l
4236: L866:
4237: unlk a6
4238: rts
4239: .globl _TstCond
4240: .align 2
4241: _TstCond:link a6,#-8
4242: tstl a6@(12); jne L870
4243: movl a6@(8),sp@
4244: movl a6@(16),sp@(4)
4245: jbsr _BadCond:l
4246: movl #L818,sp@
4247: jbsr _printf:l
4248: L870:
4249: L869:
4250: unlk a6
4251: rts
4252: .data
4253: .align 2; L873:.long L874
4254: .long L875
4255: .long L876
4256: .long L877
4257: .text
4258: .globl _BadCond
4259: .align 2
4260: _BadCond:link a6,#-20
4261: moveml #0x480,a6@(-8)
4262: movl a6@(8),d7
4263: lea @(_ErrCnt:l,d7:l:4),a2
4264: addql #1,a2@
4265: movl #L878,sp@
4266: movl @(L873:l,d7:l:4),sp@(4)
4267: movl a6@(12),sp@(8)
4268: jbsr _printf:l
4269: L872:
4270: moveml a6@(-8),#0x480
4271: unlk a6
4272: rts
4273: .globl _Random
4274: .align 2
4275: _Random:link a6,#-68
4276: fmovemx #0xb3,a6@(-60)
4277: fmoved _Random1:l,fp3
4278: faddd _Random9:l,fp3
4279: fmovex fp3,fp7
4280: fmovex fp7,fp3; fmulx fp7,fp3
4281: fmovex fp3,fp6
4282: fmovex fp6,fp3; fmulx fp6,fp3
4283: fmovex fp3,fp6
4284: fmovex fp7,fp3; fmulx fp6,fp3
4285: fmovex fp3,fp7
4286: fmoved fp7,sp@(0)
4287: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
4288: fmovex fp7,fp3; fsubx fp0,fp3
4289: fmovex fp3,fp6
4290: fmoved L880:l,fp3
4291: fmulx fp7,fp3
4292: fmovex fp6,fp2; faddx fp3,fp2
4293: fmoved fp2,_Random1:l
4294: fmoved _Random1:l,fp0
4295: L879:
4296: fmoved fp0,sp@-; movl sp@+,d0; movl sp@+,d1
4297: fmovemx a6@(-60),#0xb3
4298: unlk a6
4299: rts
4300: .globl _SqXMinX
4301: .align 2
4302: _SqXMinX:link a6,#-84
4303: fmovemx #0xb1,a6@(-56)
4304: fmoved _X:l,fp3
4305: fmovex fp3,fp2; fmuld _BInvrse:l,fp2
4306: fmovex fp2,fp7
4307: fmovex fp3,fp2; fsubx fp7,fp2
4308: fmoved fp2,a6@(-8)
4309: fmovex fp3,fp2; fmulx fp3,fp2
4310: fmoved fp2,sp@(0)
4311: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
4312: fmovex fp0,fp3; fsubx fp7,fp3
4313: fsubd a6@(-8),fp3
4314: fdivd _OneUlp:l,fp3
4315: fmoved fp3,_SqEr:l
4316: fmoved _SqEr:l,fp3
4317: fcmpd _Zero:l,fp3; fjeq L882
4318: fmoved _SqEr:l,fp3
4319: fcmpd _MinSqEr:l,fp3; fjge L884
4320: fmoved _SqEr:l,fp3
4321: fmoved fp3,_MinSqEr:l
4322: L884:
4323: fmoved _SqEr:l,fp3
4324: fcmpd _MaxSqEr:l,fp3; fjle L886
4325: fmoved _SqEr:l,fp3
4326: fmoved fp3,_MaxSqEr:l
4327: L886:
4328: fmoved _J:l,fp3
4329: faddd L11:l,fp3
4330: fmoved fp3,_J:l
4331: movl a6@(8),sp@
4332: movl #L47,sp@(4)
4333: jbsr _BadCond:l
4334: movl #L888,sp@
4335: fmoved _X:l,fp3
4336: fmovex fp3,fp2; fmulx fp3,fp2
4337: fmoved fp2,sp@(4)
4338: fmoved fp3,sp@(12)
4339: fmoved _OneUlp:l,fp3
4340: fmuld _SqEr:l,fp3
4341: fmoved fp3,sp@(20)
4342: jbsr _printf:l
4343: movl #L889,sp@
4344: jbsr _printf:l
4345: L882:
4346: L881:
4347: fmovemx a6@(-56),#0xb1
4348: unlk a6
4349: rts
4350: .globl _NewD
4351: .align 2
4352: _NewD:link a6,#-92
4353: fmovemx #0xf7,a6@(-84)
4354: fmoved _Z1:l,fp3
4355: fmuld _Q:l,fp3
4356: fmoved fp3,_X:l
4357: fmoved _Half:l,fp3
4358: fmoved _X:l,fp2
4359: fdivd _Radix:l,fp2
4360: fsubx fp2,fp3
4361: fmoved fp3,sp@(0)
4362: jbsr _floor:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
4363: fmoved _Radix:l,fp3
4364: fmovex fp3,fp7
4365: fmovex fp0,fp2; fmulx fp3,fp2
4366: faddd _X:l,fp2
4367: fmoved fp2,_X:l
4368: fmoved _Q:l,fp2
4369: fmoved _X:l,fp1
4370: fmovex fp1,fp5
4371: fmoved _Z:l,fp0
4372: fmovex fp0,fp6
4373: fmovex fp1,fp3; fmulx fp0,fp3
4374: fsubx fp3,fp2
4375: fmovex fp2,fp3; fdivx fp7,fp3
4376: fmovex fp1,fp2; fmulx fp1,fp2
4377: fmoved _D:l,fp0
4378: fmovex fp0,fp1; fdivx fp7,fp1
4379: fmulx fp1,fp2
4380: faddx fp2,fp3
4381: fmoved fp3,_Q:l
4382: fmoved _Two:l,fp3
4383: fmulx fp5,fp3
4384: fmulx fp0,fp3
4385: fmovex fp6,fp2; fsubx fp3,fp2
4386: fmoved fp2,_Z:l
4387: fmoved _Z:l,fp3
4388: fcmpd _Zero:l,fp3; fjgt L891
4389: fnegd _Z:l,fp3
4390: fmoved fp3,_Z:l
4391: fnegd _Z1:l,fp3
4392: fmoved fp3,_Z1:l
4393: L891:
4394: fmoved _Radix:l,fp3
4395: fmuld _D:l,fp3
4396: fmoved fp3,_D:l
4397: L890:
4398: fmovemx a6@(-84),#0xf7
4399: unlk a6
4400: rts
4401: .globl _SR3750
4402: .align 2
4403: _SR3750:link a6,#-68
4404: fmovemx #0xf1,a6@(-60)
4405: fmoved _X:l,fp3
4406: fmovex fp3,fp7
4407: fmoved _Radix:l,fp2
4408: fmovex fp3,fp1; fsubx fp2,fp1
4409: fmoved _Z2:l,fp0
4410: fmovex fp0,fp3; fsubx fp2,fp3
4411: fcmpx fp3,fp1; fjlt L897
4412: fmovex fp7,fp3; fsubx fp0,fp3
4413: fmoved _W:l,fp2
4414: fsubx fp0,fp2
4415: fcmpx fp2,fp3; fjgt L897
4416: addql #1,_I:l
4417: fmoved _X:l,fp3
4418: fmuld _D:l,fp3
4419: fmoved fp3,sp@(0)
4420: jbsr _sqrt:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
4421: fmoved fp0,_X2:l
4422: fmoved _X2:l,fp3
4423: fmoved _Z2:l,fp2
4424: fsubx fp2,fp3
4425: fmoved _Y:l,fp1
4426: fmovex fp1,fp0; fsubx fp2,fp0
4427: fsubx fp0,fp3
4428: fmoved fp3,_Y2:l
4429: fmoved _X8:l,fp3
4430: fmoved _Half:l,fp2
4431: fsubx fp2,fp1
4432: fdivx fp1,fp3
4433: fmoved fp3,_X2:l
4434: fmoved _X2:l,fp3
4435: fmovex fp2,fp1; fmulx fp3,fp1
4436: fmulx fp3,fp1
4437: fsubx fp1,fp3
4438: fmoved fp3,_X2:l
4439: fmoved _Y2:l,fp3
4440: faddx fp2,fp3
4441: fsubd _X2:l,fp2
4442: faddx fp2,fp3
4443: fmoved fp3,_SqEr:l
4444: fmoved _SqEr:l,fp3
4445: fcmpd _MinSqEr:l,fp3; fjge L899
4446: fmoved _SqEr:l,fp3
4447: fmoved fp3,_MinSqEr:l
4448: L899:
4449: fmoved _Y2:l,fp3
4450: fsubd _X2:l,fp3
4451: fmoved fp3,_SqEr:l
4452: fmoved _SqEr:l,fp3
4453: fcmpd _MaxSqEr:l,fp3; fjle L901
4454: fmoved _SqEr:l,fp3
4455: fmoved fp3,_MaxSqEr:l
4456: L901:
4457: L897:
4458: L896:
4459: fmovemx a6@(-60),#0xf1
4460: unlk a6
4461: rts
4462: .globl _IsYeqX
4463: .align 2
4464: _IsYeqX:link a6,#-44
4465: fmovemx #0x30,a6@(-24)
4466: fmoved _Y:l,fp3
4467: fcmpd _X:l,fp3; fjeq L905
4468: tstl _N:l; jgt L907
4469: fmoved _Z:l,fp3
4470: fmoved _Zero:l,fp2
4471: fcmpx fp2,fp3; fjne L909
4472: fmoved _Q:l,fp3
4473: fcmpx fp2,fp3; fjgt L909
4474: movl #L911,sp@
4475: jbsr _printf:l
4476: jra L910
4477: L909:
4478: movl #2,sp@
4479: movl #L912,sp@(4)
4480: jbsr _BadCond:l
4481: L910:
4482: movl #L913,sp@
4483: fmoved _Z:l,fp3
4484: fmoved fp3,sp@(4)
4485: fmoved _Q:l,fp3
4486: fmoved fp3,sp@(12)
4487: jbsr _printf:l
4488: movl #L914,sp@
4489: fmoved _Y:l,fp3
4490: fmoved fp3,sp@(4)
4491: jbsr _printf:l
4492: movl #L915,sp@
4493: fmoved _X:l,fp3
4494: fmoved fp3,sp@(4)
4495: jbsr _printf:l
4496: movl #L916,sp@
4497: fmoved _Y:l,fp3
4498: fsubd _X:l,fp3
4499: fmoved fp3,sp@(4)
4500: jbsr _printf:l
4501: L907:
4502: addql #1,_N:l
4503: L905:
4504: L904:
4505: fmovemx a6@(-24),#0x30
4506: unlk a6
4507: rts
4508: .globl _SR3980
4509: .align 2
4510: _SR3980:link a6,#-44
4511: fmovemx #0x90,a6@(-24)
4512: moveml #0x10,a6@(-28)
4513: L918:
4514: movl _I:l,d4
4515: fmovel d4,fp3
4516: fmoved fp3,_Q:l
4517: fmoved _Z:l,fp3
4518: fmoved fp3,sp@(0)
4519: fmoved _Q:l,fp3
4520: fmoved fp3,sp@(8)
4521: jbsr _pow:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
4522: fmoved fp0,_Y:l
4523: jbsr _IsYeqX:l
4524: movl _I:l,d4
4525: addql #1,d4
4526: movl d4,_I:l
4527: cmpl _M:l,d4; jle L921
4528: jra L920
4529: L921:
4530: fmoved _Z:l,fp3
4531: fmuld _X:l,fp3
4532: fmoved fp3,_X:l
4533: L919:
4534: fmoved _X:l,fp3
4535: fcmpd _W:l,fp3; fjlt L918
4536: L920:
4537: L917:
4538: moveml a6@(-28),#0x10
4539: fmovemx a6@(-24),#0x90
4540: unlk a6
4541: rts
4542: .globl _PrintIfNPositive
4543: .align 2
4544: _PrintIfNPositive:link a6,#-8
4545: tstl _N:l; jle L924
4546: movl #L926,sp@
4547: movl _N:l,sp@(4)
4548: jbsr _printf:l
4549: L924:
4550: L923:
4551: unlk a6
4552: rts
4553: .globl _TstPtUf
4554: .align 2
4555: _TstPtUf:link a6,#-64
4556: fmovemx #0xf0,a6@(-48)
4557: moveml #0x10,a6@(-52)
4558: clrl _N:l
4559: fmoved _Z:l,fp3
4560: fcmpd _Zero:l,fp3; fjeq L928
4561: movl #L930,sp@
4562: jbsr _printf:l
4563: movl #L931,sp@
4564: jbsr _printf:l
4565: movl #_sigfpe,_sigsave:l
4566: movl #_ovfl_buf,sp@
4567: jbsr _setjmp:l
4568: tstl d0; jeq L932
4569: jra L934
4570: L932:
4571: fmoved _Z:l,fp3
4572: fmovex fp3,fp2; faddx fp3,fp2
4573: fdivx fp3,fp2
4574: fmoved fp2,_Q9:l
4575: movl #L935,sp@
4576: fmoved _Q9:l,fp3
4577: fmoved fp3,sp@(4)
4578: jbsr _printf:l
4579: fmoved _Q9:l,fp3
4580: fsubd _Two:l,fp3
4581: fmoved fp3,sp@(0)
4582: jbsr _fabs:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
4583: fmoved _Radix:l,fp3
4584: fmuld _U2:l,fp3
4585: fcmpx fp3,fp0; fjge L936
4586: movl #L938,sp@
4587: jbsr _printf:l
4588: movl #L939,sp@
4589: jbsr _printf:l
4590: jra L937
4591: L936:
4592: fmoved _Q9:l,fp3
4593: fcmpd _One:l,fp3; fjlt L942
4594: fcmpd _Two:l,fp3; fjle L940
4595: L942:
4596: L934:
4597: movl #1,_N:l
4598: movl _ErrCnt+4:l,d4
4599: addql #1,d4
4600: movl d4,_ErrCnt+4:l
4601: movl #L945,sp@
4602: jbsr _printf:l
4603: jra L941
4604: L940:
4605: movl #1,_N:l
4606: movl _ErrCnt+8:l,d4
4607: addql #1,d4
4608: movl d4,_ErrCnt+8:l
4609: movl #L948,sp@
4610: jbsr _printf:l
4611: L941:
4612: L937:
4613: clrl _sigsave:l
4614: fmoved _Z:l,fp3
4615: fmoved _One:l,fp2
4616: fmovex fp3,fp1; fmulx fp2,fp1
4617: fmoved fp1,_V9:l
4618: fmoved _V9:l,fp1
4619: fmoved fp1,_Random1:l
4620: fmovex fp2,fp1; fmulx fp3,fp1
4621: fmoved fp1,_V9:l
4622: fmoved _V9:l,fp1
4623: fmoved fp1,_Random2:l
4624: fmovex fp3,fp1; fdivx fp2,fp1
4625: fmoved fp1,_V9:l
4626: fcmpd _Random1:l,fp3; fjne L949
4627: fcmpd _Random2:l,fp3; fjne L949
4628: fcmpd _V9:l,fp3; fjne L949
4629: tstl _N:l; jle L950
4630: jbsr _Pause:l
4631: jra L950
4632: L949:
4633: movl #1,_N:l
4634: movl #2,sp@
4635: movl #L953,sp@(4)
4636: jbsr _BadCond:l
4637: movl #L954,sp@
4638: fmoved _Z:l,fp3
4639: fmoved fp3,sp@(4)
4640: jbsr _printf:l
4641: fmoved _Z:l,fp3
4642: fcmpd _Random1:l,fp3; fjeq L955
4643: movl #L957,sp@
4644: fmoved _Random1:l,fp3
4645: fmoved fp3,sp@(4)
4646: jbsr _printf:l
4647: L955:
4648: fmoved _Z:l,fp3
4649: fmoved _Random2:l,fp2
4650: fcmpx fp2,fp3; fjeq L958
4651: fcmpd _Random1:l,fp2; fjeq L958
4652: movl #L960,sp@
4653: fmoved _Random2:l,fp3
4654: fmoved fp3,sp@(4)
4655: jbsr _printf:l
4656: L958:
4657: fmoved _Z:l,fp3
4658: fcmpd _V9:l,fp3; fjeq L961
4659: movl #L963,sp@
4660: fmoved _V9:l,fp3
4661: fmoved fp3,sp@(4)
4662: jbsr _printf:l
4663: L961:
4664: fmoved _Random2:l,fp3
4665: fcmpd _Random1:l,fp3; fjeq L964
4666: movl _ErrCnt+8:l,d4
4667: addql #1,d4
4668: movl d4,_ErrCnt+8:l
4669: movl #2,sp@
4670: movl #L968,sp@(4)
4671: jbsr _BadCond:l
4672: movl #L969,sp@
4673: fmoved _Random2:l,fp3
4674: fmoved fp3,sp@(4)
4675: jbsr _printf:l
4676: movl #L970,sp@
4677: fmoved _Random1:l,fp3
4678: fmoved fp3,sp@(4)
4679: jbsr _printf:l
4680: L964:
4681: jbsr _Pause:l
4682: L950:
4683: L928:
4684: L927:
4685: moveml a6@(-52),#0x10
4686: fmovemx a6@(-48),#0xf0
4687: unlk a6
4688: rts
4689: .globl _notify
4690: .align 2
4691: _notify:link a6,#-8
4692: movl #L972,sp@
4693: movl a6@(8),sp@(4)
4694: jbsr _printf:l
4695: movl #L973,sp@
4696: jbsr _printf:l
4697: L971:
4698: unlk a6
4699: rts
4700: .globl _msglist
4701: .align 2
4702: _msglist:link a6,#-16
4703: moveml #0x2400,a6@(-8)
4704: movl a6@(8),a5
4705: jra L976
4706: L975:
4707: movl #L978,sp@
4708: movl a5,a2
4709: lea a2@(4),a5
4710: movl a2@,sp@(4)
4711: jbsr _printf:l
4712: L976:
4713: tstl a5@; jne L975
4714: L974:
4715: moveml a6@(-8),#0x2400
4716: unlk a6
4717: rts
4718: .data
4719: .align 2; L980:.long L981
4720: .long L982
4721: .long L983
4722: .long L984
4723: .long L985
4724: .long L986
4725: .long L987
4726: .long L988
4727: .long L989
4728: .long 0x0
4729: .text
4730: .globl _Instructions
4731: .align 2
4732: _Instructions:link a6,#-4
4733: movl #L980,sp@
4734: jbsr _msglist:l
4735: L979:
4736: unlk a6
4737: rts
4738: .data
4739: .align 2; L991:.long L992
4740: .long L993
4741: .long L994
4742: .long L995
4743: .long L996
4744: .long L997
4745: .long L998
4746: .long L999
4747: .long L1000
4748: .long L1001
4749: .long L1002
4750: .long L1003
4751: .long L1004
4752: .long L1005
4753: .long 0x0
4754: .text
4755: .globl _Heading
4756: .align 2
4757: _Heading:link a6,#-4
4758: movl #L991,sp@
4759: jbsr _msglist:l
4760: L990:
4761: unlk a6
4762: rts
4763: .data
4764: .align 2; L1007:.long L1008
4765: .long L1009
4766: .long L1010
4767: .long L1011
4768: .long L1012
4769: .long L1013
4770: .long L1014
4771: .long L1015
4772: .long L1016
4773: .long L1017
4774: .long L1018
4775: .long L1019
4776: .long L1020
4777: .long L1021
4778: .long L1022
4779: .long L1023
4780: .long L1024
4781: .long L1025
4782: .long L1026
4783: .long 0x0
4784: .text
4785: .globl _Characteristics
4786: .align 2
4787: _Characteristics:link a6,#-4
4788: movl #L1007,sp@
4789: jbsr _msglist:l
4790: L1006:
4791: unlk a6
4792: rts
4793: .data
4794: .align 2; L1028:.long L1029
4795: .long L1030
4796: .long L1031
4797: .long L1032
4798: .long L1033
4799: .long L1034
4800: .long L1035
4801: .long L1036
4802: .long L1037
4803: .long L1038
4804: .long L1039
4805: .long L1040
4806: .long L1041
4807: .long L1042
4808: .long L1043
4809: .long L1044
4810: .long L1045
4811: .long 0x0
4812: .text
4813: .globl _History
4814: .align 2
4815: _History:link a6,#-4
4816: movl #L1028,sp@
4817: jbsr _msglist:l
4818: L1027:
4819: unlk a6
4820: rts
4821: .globl _pow
4822: .align 2
4823: _pow:link a6,#-104
4824: fmovemx #0x97,a6@(-76)
4825: moveml #0x4d0,a6@(-92)
4826: fmoved a6@(8),fp7
4827: fmoved a6@(16),fp6
4828: clrl d6
4829: clrl a6@(-8)
4830: fcmpd L10:l,fp6; fjne L1047
4831: fmoved L11:l,fp0
4832: jra L1046
4833: L1047:
4834: fcmpd L1052:l,fp6; fjlt L1051
4835: fcmpd L1053:l,fp6; fjle L1049
4836: L1051:
4837: fcmpd L865:l,fp7; fjeq L1049
4838: fmoved fp7,sp@(0)
4839: jbsr _log:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
4840: fmovex fp6,fp3; fmulx fp0,fp3
4841: fmoved fp3,sp@(0)
4842: jbsr _exp:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
4843: jra L1046
4844: L1049:
4845: fcmpd L10:l,fp6; fjge L1054
4846: fnegx fp6,fp3
4847: fmovex fp3,fp6
4848: movl #1,a6@(-8)
4849: L1054:
4850: fmoved fp6,sp@(0)
4851: lea a6@(-16),a2
4852: movl a2,sp@(8)
4853: jbsr _modf:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
4854: fmovex fp0,fp6
4855: fcmpd L10:l,fp6; fjeq L1056
4856: fmoved fp7,sp@(0)
4857: jbsr _log:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
4858: fmovex fp6,fp3; fmulx fp0,fp3
4859: fmoved fp3,sp@(0)
4860: jbsr _exp:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
4861: fmovex fp0,fp5
4862: jra L1057
4863: L1056:
4864: fmoved L11:l,fp5
4865: L1057:
4866: fmoved fp7,sp@(0)
4867: lea a6@(-4),a2
4868: movl a2,sp@(8)
4869: jbsr _frexp:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
4870: fmovex fp0,fp7
4871: fmoved a6@(-16),fp3
4872: fintrzx fp3,fp0; fmovel fp0,d4
4873: movl d4,d7
4874: tstl d4; jeq L1058
4875: L1060:
4876: movl d7,d4; andl #1,d4
4877: tstl d4; jeq L1064
4878: fmovex fp5,fp3; fmulx fp7,fp3
4879: fmovex fp3,fp5
4880: movl d6,d4; addl a6@(-4),d4
4881: movl d4,d6
4882: L1064:
4883: movl d7,d4; asrl #1,d4
4884: movl d4,d7
4885: tstl d4; jne L1066
4886: jra L1062
4887: L1066:
4888: fmovex fp7,fp3; fmulx fp7,fp3
4889: fmovex fp3,fp7
4890: movl a6@(-4),d4
4891: asll #1,d4
4892: movl d4,a6@(-4)
4893: fcmpd L1070:l,fp7; fjge L1060
4894: fmoved L1071:l,fp3
4895: fmulx fp7,fp3
4896: fmovex fp3,fp7
4897: subql #1,a6@(-4)
4898: jra L1060
4899: L1062:
4900: L1058:
4901: tstl a6@(-8); jeq L1072
4902: fmoved L11:l,fp3
4903: fdivx fp5,fp3
4904: fmovex fp3,fp5
4905: movl d6,d4; negl d4
4906: movl d4,d6
4907: L1072:
4908: fmoved fp5,sp@(0)
4909: movl d6,sp@(8)
4910: jbsr _ldexp:l; movl d1,sp@-; movl d0,sp@-; fmoved sp@+,fp0
4911: L1046:
4912: fmoved fp0,sp@-; movl sp@+,d0; movl sp@+,d1
4913: moveml a6@(-92),#0x4d0
4914: fmovemx a6@(-76),#0x97
4915: unlk a6
4916: rts
4917: .bss
4918: .globl _UfNGrad
4919: .comm _UfNGrad,4
4920: .globl _SqRWrng
4921: .comm _SqRWrng,4
4922: .globl _IEEE
4923: .comm _IEEE,4
4924: .globl _Anomaly
4925: .comm _Anomaly,4
4926: .globl _Monot
4927: .comm _Monot,4
4928: .globl _NotMonot
4929: .comm _NotMonot,4
4930: .globl _Done
4931: .comm _Done,4
4932: .globl _Break
4933: .comm _Break,4
4934: .globl _RSqrt
4935: .comm _RSqrt,4
4936: .globl _RAddSub
4937: .comm _RAddSub,4
4938: .globl _RDiv
4939: .comm _RDiv,4
4940: .globl _RMult
4941: .comm _RMult,4
4942: .globl _GAddSub
4943: .comm _GAddSub,4
4944: .globl _GDiv
4945: .comm _GDiv,4
4946: .globl _GMult
4947: .comm _GMult,4
4948: .globl _N1
4949: .comm _N1,4
4950: .globl _N
4951: .comm _N,4
4952: .globl _M
4953: .comm _M,4
4954: .globl _PageNo
4955: .comm _PageNo,4
4956: .globl _Milestone
4957: .comm _Milestone,4
4958: .globl _fpecount
4959: .comm _fpecount,4
4960: .globl _ErrCnt
4961: .comm _ErrCnt,16
4962: .globl _Z9
4963: .comm _Z9,8
4964: .globl _Z2
4965: .comm _Z2,8
4966: .globl _Z1
4967: .comm _Z1,8
4968: .globl _PseudoZero
4969: .comm _PseudoZero,8
4970: .globl _Z
4971: .comm _Z,8
4972: .globl _Random2
4973: .comm _Random2,8
4974: .globl _Y2
4975: .comm _Y2,8
4976: .globl _Y1
4977: .comm _Y1,8
4978: .globl _Y
4979: .comm _Y,8
4980: .globl _Random1
4981: .comm _Random1,8
4982: .globl _X8
4983: .comm _X8,8
4984: .globl _X2
4985: .comm _X2,8
4986: .globl _X1
4987: .comm _X1,8
4988: .globl _X
4989: .comm _X,8
4990: .globl _W
4991: .comm _W,8
4992: .globl _V9
4993: .comm _V9,8
4994: .globl _V0
4995: .comm _V0,8
4996: .globl _V
4997: .comm _V,8
4998: .globl _U2
4999: .comm _U2,8
5000: .globl _U1
5001: .comm _U1,8
5002: .globl _UfThold
5003: .comm _UfThold,8
5004: .globl _OneUlp
5005: .comm _OneUlp,8
5006: .globl _S
5007: .comm _S,8
5008: .globl _Underflow
5009: .comm _Underflow,8
5010: .globl _T
5011: .comm _T,8
5012: .globl _Random9
5013: .comm _Random9,8
5014: .globl _R
5015: .comm _R,8
5016: .globl _Q9
5017: .comm _Q9,8
5018: .globl _Q
5019: .comm _Q,8
5020: .globl _Precision
5021: .comm _Precision,8
5022: .globl _MyZero
5023: .comm _MyZero,8
5024: .globl _J
5025: .comm _J,8
5026: .globl _StickyBit
5027: .comm _StickyBit,8
5028: .globl _I
5029: .comm _I,4
5030: .globl _HInvrse
5031: .comm _HInvrse,8
5032: .globl _H
5033: .comm _H,8
5034: .globl _F9
5035: .comm _F9,8
5036: .globl _F6
5037: .comm _F6,8
5038: .globl _Third
5039: .comm _Third,8
5040: .globl _E9
5041: .comm _E9,8
5042: .globl _MaxSqEr
5043: .comm _MaxSqEr,8
5044: .globl _SqEr
5045: .comm _SqEr,8
5046: .globl _MinSqEr
5047: .comm _MinSqEr,8
5048: .globl _E3
5049: .comm _E3,8
5050: .globl _Exp2
5051: .comm _Exp2,8
5052: .globl _E1
5053: .comm _E1,8
5054: .globl _E0
5055: .comm _E0,8
5056: .globl _FourD
5057: .comm _FourD,8
5058: .globl _D
5059: .comm _D,8
5060: .globl _CInvrse
5061: .comm _CInvrse,8
5062: .globl _C
5063: .comm _C,8
5064: .globl _A1
5065: .comm _A1,8
5066: .globl _AInvrse
5067: .comm _AInvrse,8
5068: .globl _ch
5069: .comm _ch,8
5070: .globl _Indx
5071: .comm _Indx,4
5072: .globl _BMinusU2
5073: .comm _BMinusU2,8
5074: .globl _RadixD2
5075: .comm _RadixD2,8
5076: .globl _BInvrse
5077: .comm _BInvrse,8
5078: .globl _Radix
5079: .comm _Radix,8
5080: .globl _sigsave
5081: .comm _sigsave,4
5082: .globl _ovfl_buf
5083: .comm _ovfl_buf,232
5084: .text
5085: .align 2; L1071:.long 0x40000000,0x0
5086: .align 2; L1070:.long 0x3fe00000,0x0
5087: .align 2; L1053:.long 0x40913000,0x0
5088: .align 2; L1052:.long 0xc0913000,0x0
5089: .align 2; L1045:.ascii "see source comments for more history.\000"
5090: .align 2; L1044:.ascii "BASIC version of this program (C) 1983 by Prof. W. M. Kahan;\000"
5091: .align 2; L1043:.ascii "as used by certain early WANG machines.\012\000"
5092: .align 2; L1042:.ascii "floating-point numbers, but also allows logarithmic encoding\000"
5093: .align 2; L1041:.ascii "\012The program is based upon a conventional radix representatio"
5094: .ascii "n for\000"
5095: .align 2; L1040:.ascii "of pathologies, and to say how well the arithmetic is implemente"
5096: .ascii "d.\000"
5097: .align 2; L1039:.ascii "of the arithmetic, this program tries to cope with a wider varie"
5098: .ascii "ty\000"
5099: .align 2; L1038:.ascii "the Radix, Precision and range (over/underflow thresholds)\000"
5100: .align 2; L1037:.ascii "W. J. Cody and W. Waite. Although both programs try to discover\000"
5101: .align 2; L1036:.ascii "book `Software Manual for the Elementary Functions' (1980) by\000"
5102: .align 2; L1035:.ascii "program called `MACHAR', which can be found at the end of the\000"
5103: .align 2; L1034:.ascii "The diagnostic capabilities of this program go beyond an earlier"
5104: .ascii "\000"
5105: .align 2; L1033:.ascii "Failures may confound subsequent diagnoses.\012\000"
5106: .align 2; L1032:.ascii " FAILUREs, like 2+2 == 5 .\000"
5107: .align 2; L1031:.ascii " Serious DEFECTs, like lack of a guard digit, and\000"
5108: .align 2; L1030:.ascii " FLAWs, like lack of a sticky bit,\000"
5109: .align 2; L1029:.ascii "The program attempts to discriminate among\000"
5110: .align 2; L1026:.ascii " Decimal-Binary conversion is NOT YET tested for accuracy.\000"
5111: .align 2; L1025:.ascii " Extra-precise subexpressions are revealed but NOT YET teste"
5112: .ascii "d.\000"
5113: .align 2; L1024:.ascii " Sqrt is tested. Y^X is not tested.\000"
5114: .align 2; L1023:.ascii "\011and for contamination with pseudo-zeros.\000"
5115: .align 2; L1022:.ascii " Comparisions are checked for consistency with subtraction\000"
5116: .align 2; L1021:.ascii " V0 tells, roughly, whether Infinity is represented.\000"
5117: .align 2; L1020:.ascii " V = an overflow threshold, roughly.\000"
5118: .align 2; L1019:.ascii " E0 and PseudoZero tell whether underflow is abrupt, gradual"
5119: .ascii ", or fuzzy.\000"
5120: .align 2; L1018:.ascii " UnderflowThreshold = an underflow threshold.\000"
5121: .align 2; L1017:.ascii " Whether a Sticky Bit used correctly for rounding.\000"
5122: .align 2; L1016:.ascii "\011for Mult., Div., Add/Subt. and Sqrt.\000"
5123: .align 2; L1015:.ascii " Whether arithmetic is chopped, correctly rounded, or someth"
5124: .ascii "ing else\000"
5125: .align 2; L1014:.ascii " Adequacy of guard digits for Mult., Div. and Subt.\000"
5126: .align 2; L1013:.ascii " U1 = 1/Radix^Precision = One Ulp of numbers a little less t"
5127: .ascii "han 1.0 .\000"
5128: .align 2; L1012:.ascii "\011(OneUlpnit in the Last Place) of 1.000xxx .\000"
5129: .align 2; L1011:.ascii " U2 = Radix/Radix^Precision = One Ulp\000"
5130: .align 2; L1010:.ascii " Precision = number of significant digits carried.\000"
5131: .align 2; L1009:.ascii " Radix = 1, 2, 4, 8, 10, 16, 100, 256 ...\000"
5132: .align 2; L1008:.ascii "Running this program should reveal these characteristics:\000"
5133: .align 2; L1005:.ascii "\011Other relevant compiler options:\000"
5134: .align 2; L1004:.ascii "\011Optimization level:\012\000"
5135: .align 2; L1003:.ascii "\011Compiler:\012\000"
5136: .align 2; L1002:.ascii "\011Computer:\012\000"
5137: .align 2; L1001:.ascii "\011Version:\01110 February 1989;\000"
5138: .align 2; L1000:.ascii "\011Precision:\011double;\000"
5139: .align 2; L999:.ascii "In doing so, please include the following information:\000"
5140: .align 2; L998:.ascii "\011San Francisco, CA 94143-0704, USA\012\000"
5141: .align 2; L997:.ascii "\011University of California\000"
5142: .align 2; L996:.ascii "\011Computer Center U-76\000"
5143: .align 2; L995:.ascii "\011Richard Karpinski\000"
5144: .align 2; L994:.ascii "Please send suggestions and interesting results to\000"
5145: .align 2; L993:.ascii "cope with unanticipated and newly uncovered arithmetic pathologi"
5146: .ascii "es.\012\000"
5147: .align 2; L992:.ascii "Users are invited to help debug and augment this program so it w"
5148: .ascii "ill\000"
5149: .align 2; L989:.ascii "Answer questions with Y, y, N or n (unless otherwise indicated)."
5150: .ascii "\012\000"
5151: .align 2; L988:.ascii "amend it to make further progress.\012\000"
5152: .align 2; L987:.ascii "program anyway to see how many milestones it passes, and then\000"
5153: .align 2; L986:.ascii "warning. If persuasion avails naught, don't despair but run thi"
5154: .ascii "s\000"
5155: .align 2; L985:.ascii "to persevere with a surrogate value after, perhaps, displaying s"
5156: .ascii "ome\000"
5157: .align 2; L984:.ascii "error like Over/Underflow or Division by Zero occurs, but rather"
5158: .ascii "\000"
5159: .align 2; L983:.ascii "try to persuade the computer NOT to terminate execution when an\000"
5160: .align 2; L982:.ascii " `END OF TEST',\012\000"
5161: .align 2; L981:.ascii "Lest this program stop prematurely, i.e. before displaying\012\000"
5162: .align 2; L978:.byte 37,115,10,0
5163: .align 2; L973:.ascii " PLEASE NOTIFY KARPINKSI!\012\000"
5164: .align 2; L972:.ascii "%s test appears to be inconsistent...\012\000"
5165: .align 2; L970:.ascii "\011differs from Z * 1 = %.17e\012\000"
5166: .align 2; L969:.ascii "\011Comparison alleges that 1 * Z = %.17e\012\000"
5167: .align 2; L968:.ascii "Multiplication does not commute!\012\000"
5168: .align 2; L963:.ascii "Z / 1 = %.17e\012\000"
5169: .align 2; L960:.ascii "1 * Z == %g\012\000"
5170: .align 2; L957:.ascii "Z * 1 = %.17e \000"
5171: .align 2; L954:.ascii "%.17e\012\011compares different from \000"
5172: .align 2; L953:.ascii "What prints as Z = \000"
5173: .align 2; L948:.ascii "This is a DEFECT!\012\000"
5174: .align 2; L945:.ascii "This is a VERY SERIOUS DEFECT!\012\000"
5175: .align 2; L939:.ascii " has NOT just been signaled.\012\000"
5176: .align 2; L938:.ascii "This is O.K., provided Over/Underflow\000"
5177: .align 2; L935:.ascii "What the machine gets for (Z + Z) / Z is %.17e .\012\000"
5178: .align 2; L931:.ascii "(Z + Z) / Z should be safe.\012\000"
5179: .align 2; L930:.ascii "Since comparison denies Z = 0, evaluating \000"
5180: .align 2; L926:.ascii "Similar discrepancies have occurred %d times.\012\000"
5181: .align 2; L916:.ascii "\011\011they differ by %.17e .\012\000"
5182: .align 2; L915:.ascii "\011which compared unequal to correct %.17e ;\012\000"
5183: .align 2; L914:.ascii "\011yielded %.17e;\012\000"
5184: .align 2; L913:.ascii "\011(%.17e) ^ (%.17e)\012\000"
5185: .align 2; L912:.ascii "computing\012\000"
5186: .align 2; L911:.ascii "WARNING: computing\012\000"
5187: .align 2; L889:.ascii "\011instead of correct value 0 .\012\000"
5188: .align 2; L888:.ascii "sqrt( %.17e) - %.17e = %.17e\012\000"
5189: .align 2; L880:.long 0x3ed4f8b5,0x88e368f0
5190: .align 2; L878:.ascii "%s: %s\000"
5191: .align 2; L877:.ascii "FLAW\000"
5192: .align 2; L876:.ascii "DEFECT\000"
5193: .align 2; L875:.ascii "SERIOUS DEFECT\000"
5194: .align 2; L874:.ascii "FAILURE\000"
5195: .align 2; L868:.ascii " Page: %d\012\012\000"
5196: .align 2; L867:.ascii "\012Diagnosis resumes after milestone Number %d\000"
5197: .align 2; L865:.long 0xbff00000,0x0
5198: .align 2; L825:.ascii "END OF TEST.\012\000"
5199: .align 2; L824:.ascii "\012A total of %d floating point exceptions were registered.\012"
5200: .ascii "\000"
5201: .align 2; L821:.ascii "The arithmetic diagnosed appears to be Excellent!\012\000"
5202: .align 2; L820:.ascii " during Gradual Underflow.\012\000"
5203: .align 2; L819:.ascii ",\012except for possibly Double Rounding\000"
5204: .align 2; L818:.byte 46,10,0
5205: .align 2; L815:.ascii "854\000"
5206: .align 2; L814:.ascii "754\000"
5207: .align 2; L811:.ascii "the proposed IEEE standard P\000"
5208: .align 2; L810:.ascii "Rounding appears to conform to \000"
5209: .align 2; L807:.ascii "The arithmetic diagnosed seems Satisfactory.\012\000"
5210: .align 2; L801:.ascii "No failures, defects nor flaws have been discovered.\012\000"
5211: .align 2; L800:.ascii " program's subsequent diagnoses.\012\000"
5212: .align 2; L799:.ascii "Potentially fatal FAILURE may have spoiled this\000"
5213: .align 2; L796:.ascii "unacceptable Serious Defects.\012\000"
5214: .align 2; L795:.ascii "The arithmetic diagnosed has \000"
5215: .align 2; L791:.ascii "despite inconvenient Defects.\012\000"
5216: .align 2; L790:.ascii "The arithmetic diagnosed may be Acceptable\012\000"
5217: .align 2; L785:.ascii "Satisfactory though flawed.\012\000"
5218: .align 2; L784:.ascii "The arithmetic diagnosed seems \000"
5219: .align 2; L773:.ascii "The number of %-29s %d.\012\000"
5220: .align 2; L766:.ascii "FLAWs discovered =\000"
5221: .align 2; L765:.ascii "DEFECTs discovered =\000"
5222: .align 2; L764:.ascii "SERIOUS DEFECTs discovered =\000"
5223: .align 2; L763:.ascii "FAILUREs encountered =\000"
5224: .align 2; L759:.ascii "\012 Trying to compute 0 / 0 produces ...\000"
5225: .align 2; L758:.ascii " %.7e .\012\000"
5226: .align 2; L755:.ascii " Trying to compute 1 / 0 produces ...\000"
5227: .align 2; L754:.ascii "What message and/or values does Division by Zero produce?\012\000"
5228: .align 2; L753:.ascii " instead, X / X - 1/2 - 1/2 = %.17e .\012\000"
5229: .align 2; L752:.ascii " X / X differs from 1 when X = %.17e\012\000"
5230: .align 2; L747:.ascii " X / X traps when X = %g\012\000"
5231: .align 2; L731:.ascii "is too far from 1.\012\000"
5232: .align 2; L730:.ascii " unbalanced range; UfThold * V = %.17e\012\011%s\012\000"
5233: .align 2; L729:.ascii "Badly\000"
5234: .align 2; L722:.ascii " is too far from sqrt(Z) ^ 2 (%.17e) .\012\000"
5235: .align 2; L721:.ascii "Comparison alleges that Z = %17e\012\000"
5236: .align 2; L709:.ascii " is too far from sqrt(Z) ^ 2 = %.17e .\012\000"
5237: .align 2; L708:.ascii "Comparison alleges that what prints as Z = %.17e\012\000"
5238: .align 2; L691:.ascii "+-%g, +-%g\012and +-%g are confused by Overflow.\000"
5239: .align 2; L690:.ascii "Comparisons involving \000"
5240: .align 2; L684:.ascii "above is a DEFECT.\012\000"
5241: .align 2; L683:.ascii "Any overflow signal separating this * from the one\012\000"
5242: .align 2; L682:.ascii " nor for V / 1 = %.17e .\012\000"
5243: .align 2; L681:.ascii "No Overflow should be signaled for V * 1 = %.17e\012\000"
5244: .align 2; L680:.ascii "There is no saturation value because the system traps on overflo"
5245: .ascii "w.\012\000"
5246: .align 2; L679:.ascii "Overflow saturates at V0 = %.17e .\012\000"
5247: .align 2; L676:.ascii "Overflow threshold is V = %.17e .\012\000"
5248: .align 2; L667:.ascii "overflow past %.17e\012\011shrinks to %.17e .\012\000"
5249: .align 2; L664:.ascii "-(-Y) differs from Y.\012\000"
5250: .align 2; L663:.ascii "finds a \000"
5251: .align 2; L662:.ascii "Seems O.K.\012\000"
5252: .align 2; L659:.ascii "Trying it on Y = %.17e .\012\000"
5253: .align 2; L658:.ascii "Can `Z = -Y' overflow?\012\000"
5254: .align 2; L651:.ascii "This may generate an error.\012\000"
5255: .align 2; L650:.ascii "Searching for Overflow threshold:\012\000"
5256: .align 2; L649:.ascii " ... no discrepancies found.\012\000"
5257: .align 2; L641:.ascii "Testing powers Z^Q at four nearly extreme values.\012\000"
5258: .align 2; L640:.ascii "Accuracy seems adequate.\012\000"
5259: .align 2; L633:.ascii "\011calculations involving tiny interest rates.\012\000"
5260: .align 2; L632:.ascii "\011This much error may spoil financial\012\000"
5261: .align 2; L631:.ascii "\011differs from correct value by %.17e .\012\000"
5262: .align 2; L630:.ascii "\011(1 + (%.17e) ^ (%.17e);\012\000"
5263: .align 2; L629:.ascii " %.17e for\012\000"
5264: .align 2; L628:.ascii "Calculated\000"
5265: .align 2; L621:.ascii "Testing X^((X + 1) / (X - 1)) vs. exp(2) = %.17e as X -> 1.\012\000"
5266: .align 2; L617:.ascii "This computed value is O.K.\012\000"
5267: .align 2; L614:.ascii " threshold = %.17e .\012\000"
5268: .align 2; L613:.ascii "this is not between 0 and underflow\012\000"
5269: .align 2; L609:.ascii "actually calculating yields: %.17e .\012\000"
5270: .align 2; L608:.ascii "should afflict the expression\012\011(%.17e) ^ (%.17e);\012\000"
5271: .align 2; L607:.ascii "UfThold = (%.17e) ^ (%.17e)\012only underflow \000"
5272: .align 2; L606:.ascii "Since underflow occurs below the threshold\012\000"
5273: .align 2; L605:.ascii "Range is too narrow; U1^%d Underflows.\012\000"
5274: .align 2; L600:.ascii "merely roundoff.\012\000"
5275: .align 2; L599:.ascii "calculation may suffer larger Relative error than \000"
5276: .align 2; L598:.ascii " below which\000"
5277: .align 2; L597:.ascii "The Underflow threshold is %.17e, %s\012\000"
5278: .align 2; L596:.ascii "X / Z = 1 + %g .\012\000"
5279: .align 2; L595:.ascii "X / Z fails!\012\000"
5280: .align 2; L592:.ascii "encounter Division by Zero although actually\012\000"
5281: .align 2; L591:.ascii " ... (f(X) - f(Z)) / (X - Z) ...\012\000"
5282: .align 2; L590:.ascii " if (X == Z) ... else\000"
5283: .align 2; L589:.ascii "confusion when innocent statements like\012\000"
5284: .align 2; L588:.ascii "this is a SERIOUS DEFECT\012that causes \000"
5285: .align 2; L587:.ascii " Should this NOT signal Underflow, \000"
5286: .align 2; L586:.ascii "yet X - Z yields %.17e .\012\000"
5287: .align 2; L585:.ascii "X = %.17e\012\011is not equal to Z = %.17e .\012\000"
5288: .align 2; L580:.ascii "Underflow / UfThold failed!\012\000"
5289: .align 2; L569:.ascii "(roundoff in UfThold) < E0.\012\000"
5290: .align 2; L568:.ascii "Underflow is gradual; it incurs Absolute Error =\012\000"
5291: .align 2; L563:.ascii "|Q - Y| = %.17e .\012\000"
5292: .align 2; L562:.ascii "print out as Q = %.17e, Y = %.17e .\012\000"
5293: .align 2; L561:.ascii "Q == Y while denying that |Q - Y| == 0; these values\012\000"
5294: .align 2; L560:.ascii "Underflow confuses Comparison, which alleges that\012\000"
5295: .align 2; L558:.ascii " or else multiplication gets too many last digits wrong.\012\000"
5296: .align 2; L557:.ascii " coming down from %.17e\012\000"
5297: .align 2; L556:.ascii "approach a threshold = %.17e\012\000"
5298: .align 2; L555:.ascii "Either accuracy deteriorates as numbers\012\000"
5299: .align 2; L543:.ascii "Smallest strictly positive number found is E0 = %g .\012\000"
5300: .align 2; L542:.ascii " threshold than products.\012\000"
5301: .align 2; L541:.ascii "Difference underflows at a higher\000"
5302: .align 2; L538:.ascii " threshold than differences.\012\000"
5303: .align 2; L537:.ascii "Products underflow at a higher\000"
5304: .align 2; L530:.ascii "value PseudoZero that prints out as %g .\012\000"
5305: .align 2; L529:.ascii "Underflow can stick at an allegedly positive\012\000"
5306: .align 2; L528:.ascii "positive, isn't; it prints out as %g .\012\000"
5307: .align 2; L527:.ascii "But -PseudoZero, which should be\012\000"
5308: .align 2; L524:.ascii "PseudoZero that prints out as: %g .\012\000"
5309: .align 2; L523:.ascii "allegedly negative value\012\000"
5310: .align 2; L522:.ascii "Positive expressions can underflow to an\012\000"
5311: .align 2; L509:.ascii "multiplication gets too many last digits wrong.\012\000"
5312: .align 2; L485:.ascii "Seeking Underflow thresholds UfThold and E0.\012\000"
5313: .align 2; L482:.ascii "... no discrepancis found.\012\000"
5314: .align 2; L479:.ascii "\011involving interest rates.\012\000"
5315: .align 2; L478:.ascii "Errors like this may invalidate financial calculations\012\000"
5316: .align 2; L457:.ascii "Testing powers Z^i for small Integers Z and i.\012\000"
5317: .align 2; L454:.ascii "sqrt gets too many last digits wrong\000"
5318: .align 2; L452:.ascii "to %.7e ulps.\012\000"
5319: .align 2; L451:.ascii "Observed errors run from %.7e \000"
5320: .align 2; L450:.ascii "Square root is neither chopped nor correctly rounded.\012\000"
5321: .align 2; L447:.ascii "Square root appears to be chopped.\012\000"
5322: .align 2; L442:.ascii "Square root appears to be correctly rounded.\012\000"
5323: .align 2; L437:.ascii " fails test whether sqrt rounds or chops.\012\000"
5324: .align 2; L436:.ascii "Radix^Precision = %.7e\012\000"
5325: .align 2; L435:.ascii "Anomalous arithmetic with Integer < \000"
5326: .align 2; L395:.ascii "Testing whether sqrt is rounded or chopped.\012\000"
5327: .align 2; L373:.ascii "sqrt(X) is non-monotonic for X near %.7e .\012\000"
5328: .align 2; L372:.ascii "sqrt has passed a test for Monotonicity.\012\000"
5329: .align 2; L355:.ascii "Test for sqrt monotonicity.\012\000"
5330: .align 2; L344:.ascii "Testing if sqrt(X * X) == X for %d Integers X.\012\000"
5331: .align 2; L339:.ascii "Square root of 0.0, -0.0 or 1.0 wrong\000"
5332: .align 2; L337:.ascii "\012Running test of square root(x).\012\000"
5333: .align 2; L336:.ascii " No failures found in %d integer pairs.\012\000"
5334: .align 2; L335:.ascii "X * Y == Y * X trial fails.\012\000"
5335: .align 2; L326:.long 0x40080000,0x0
5336: .align 2; L325:.ascii "Testing on %d random pairs.\012\000"
5337: .align 2; L324:.ascii "Does Multiplication commute? \000"
5338: .align 2; L321:.ascii "lack(s) of guard digits or failure(s) to correctly round or chop"
5339: .ascii "\012(noted above) count as one flaw in the final tally below\000"
5340: .align 2; L319:.ascii "Sticky bit used incorrectly or not at all.\012\000"
5341: .align 2; L318:.ascii "Sticky bit apparently used correctly.\012\000"
5342: .align 2; L301:.ascii "Checking for sticky bit.\012\000"
5343: .align 2; L298:.ascii "(X - Y) + (Y - X) is non zero!\012\000"
5344: .align 2; L295:.ascii "Addition/Subtraction neither rounds nor chops.\012\000"
5345: .align 2; L294:.ascii "Add/Subtract\000"
5346: .align 2; L291:.ascii "Addition/Subtraction appears to round correctly.\012\000"
5347: .align 2; L284:.ascii "Add/Subtract appears to be chopped.\012\000"
5348: .align 2; L279:.ascii "Incomplete carry-propagation in Addition\000"
5349: .align 2; L275:.ascii "Radix * ( 1 / Radix ) differs from 1\000"
5350: .align 2; L273:.ascii "/ is neither chopped nor correctly rounded.\012\000"
5351: .align 2; L270:.ascii "Division appears to chop.\012\000"
5352: .align 2; L267:.ascii "Division\000"
5353: .align 2; L264:.ascii "Division appears to round correctly.\012\000"
5354: .align 2; L259:.ascii "Multiplication\000"
5355: .align 2; L256:.ascii "* is neither chopped nor correctly rounded.\012\000"
5356: .align 2; L255:.ascii "Multiplication appears to chop.\012\000"
5357: .align 2; L252:.ascii "Multiplication appears to round correctly.\012\000"
5358: .align 2; L242:.ascii "X * (1/X) differs from 1\000"
5359: .align 2; L223:.ascii "Checking rounding on multiply, divide and add/subtract.\012\000"
5360: .align 2; L222:.ascii " *, /, and - appear to have guard digits, as they should.\012"
5361: .ascii "\000"
5362: .align 2; L219:.ascii " ... if (X == 1.0) {.....} else {.../(X-1.0)...}\012\000"
5363: .align 2; L218:.ascii " such precautions against division by zero as\012\000"
5364: .align 2; L217:.ascii " subtraction yields (1-U1) - 1 = 0 , thereby vitiating\012\000"
5365: .align 2; L216:.ascii "comparison alleges (1-U1) < 1 although\012\000"
5366: .align 2; L213:.ascii "- lacks Guard Digit, so cancellation is obscured\000"
5367: .align 2; L208:.ascii "* and/or / gets too many last digits wrong\000"
5368: .align 2; L204:.ascii "Computed value of 1/1.000..1 >= 1\000"
5369: .align 2; L202:.ascii "Division lacks a Guard Digit, so X/1 != X\000"
5370: .align 2; L197:.ascii "Division lacks a Guard Digit, so error can exceed 1 ulp\012or 1"
5371: .ascii "/3 and 3/9 and 9/27 may disagree\000"
5372: .align 2; L193:.ascii "* gets too many final digits wrong.\012\000"
5373: .align 2; L191:.ascii "* lacks a Guard Digit, so 1*X != X\000"
5374: .align 2; L188:.ascii "\012Checking for guard digit in *, /, and -.\012\000"
5375: .align 2; L187:.ascii "Subtraction appears to be normalized, as it should be.\000"
5376: .align 2; L182:.ascii "Subtraction is not normalized X=Y,X+Z != Y+Z!\000"
5377: .align 2; L178:.ascii "roughly %g extra significant decimals.\012\000"
5378: .align 2; L177:.ascii "precisely with about %g extra B-digits, i.e.\012\000"
5379: .align 2; L176:.ascii "Some subexpressions appear to be calculated extra\012\000"
5380: .align 2; L170:.ascii "of an\012extra-precision\000"
5381: .align 2; L169:.ascii "Z1 = %.7e, or Z2 = %.7e \000"
5382: .align 2; L168:.ascii ", or exact rational arithmetic a result\012\000"
5383: .align 2; L167:.ascii "Because of unusual Radix = %f\000"
5384: .align 2; L163:.ascii "\011U2 = %.7e, Z2 - U2 = %.7e\012\000"
5385: .align 2; L162:.ascii "\011U1 = %.7e, Z1 - U1 = %.7e\012\000"
5386: .align 2; L161:.ascii "Precision\000"
5387: .align 2; L160:.byte 0
5388: .align 2; L153:.ascii "That feature is not tested further by this program.\012\000"
5389: .align 2; L148:.ascii "Possibly some part of this\000"
5390: .align 2; L147:.ascii "by extra-precise evaluation of arithmetic subexpressions.\012\000"
5391: .align 2; L146:.ascii "are symptoms of inconsistencies introduced\012\000"
5392: .align 2; L145:.ascii "respectively %.7e, %.7e, %.7e,\012\000"
5393: .align 2; L144:.ascii "Disagreements among the values X1, Y1, Z1,\012\000"
5394: .align 2; L122:.ascii "Precision worse than 5 decimal figures \000"
5395: .align 2; L120:.ascii "The number of significant digits of the Radix is %f .\012\000"
5396: .align 2; L119:.ascii "logarithmic encoding has precision characterized solely by U1.\012"
5397: .ascii "\000"
5398: .align 2; L116:.ascii "of significant digits but, by itself, this is a minor flaw.\012\000"
5399: .align 2; L115:.ascii "Precision cannot be characterized by an Integer number\012\000"
5400: .align 2; L102:.ascii "Comparison is fuzzy,X=1 but X-1/2-1/2 != 0\000"
5401: .align 2; L98:.ascii "(1-U1)-1/2 < 1/2 is FALSE, prog. fails?\000"
5402: .align 2; L96:.long 0x40240000,0x0
5403: .align 2; L91:.ascii "Radix is not as good as 2 or 10\000"
5404: .align 2; L87:.ascii "Radix is too big: roundoff problems\000"
5405: .align 2; L85:.ascii "MYSTERY: recalculated Radix = %.7e .\012\000"
5406: .align 2; L84:.ascii "Radix confirmed.\012\000"
5407: .align 2; L81:.long 0x3f847ae1,0x47ae147a
5408: .align 2; L80:.ascii "gets better closest relative separation U1 = %.7e .\012\000"
5409: .align 2; L79:.ascii "confirms closest relative separation U1 .\012\000"
5410: .align 2; L64:.ascii "Recalculating radix and precision\012 \000"
5411: .align 2; L63:.ascii "Closest relative separation found is U1 = %.7e .\012\012\000"
5412: .align 2; L57:.ascii "Radix = %f .\012\000"
5413: .align 2; L48:.ascii "Searching for Radix and Precision.\012\000"
5414: .align 2; L47:.byte 10,0
5415: .align 2; L46:.ascii "-1, 0, 1/2, 1, 2, 3, 4, 5, 9, 27, 32 & 240 are O.K.\012\000"
5416: .align 2; L41:.ascii "5 != 4+1, 240/3 != 80, 240/4 != 60, or 240/5 != 48\000"
5417: .align 2; L37:.ascii "9 != 3*3, 27 != 9*3, 32 != 8*4, or 32-27-4-1 != 0\000"
5418: .align 2; L33:.ascii "1/2 + (-1) + 1/2 != 0\000"
5419: .align 2; L29:.ascii "-1+1 != 0, (-1)+abs(1) != 0, or -1+(-1)*(-1) != 0\000"
5420: .align 2; L25:.ascii "3 != 2+1, 4 != 3+1, 4+2*(-2) != 0, or 4-3-1 != 0\000"
5421: .align 2; L23:.long 0x3f50624d,0xd2f1a9fb
5422: .align 2; L22:.ascii "Comparison alleges that -0.0 is Non-zero!\012\000"
5423: .align 2; L17:.ascii "0+0 != 0, 1-1 != 0, 1 <= 0, or 1+1 != 2\000"
5424: .align 2; L15:.ascii "Program is now RUNNING tests on small integers:\012\000"
5425: .align 2; L11:.long 0x3ff00000,0x0
5426: .align 2; L10:.long 0x0,0x0
5427: .align 2; L9:.ascii "double\000"
5428: .align 2; L8:.ascii "precision\000"
5429: .align 2; L7:.ascii "set\000"
5430: .align 2; L2:.ascii "\012* * * FLOATING-POINT ERROR * * *\012\000"
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.