| 1 |
15 |
hellwig |
\I\X47:Add nonzero numbers and \PB{\&{return}}\X
|
| 2 |
|
|
\U46.
|
| 3 |
|
|
\I\X49:Adjust for difference in exponents\X
|
| 4 |
|
|
\U47.
|
| 5 |
|
|
\I\X14:Check that \PB{$\|x<\|z$}; otherwise give trivial answer\X
|
| 6 |
|
|
\U13.
|
| 7 |
|
|
\I\X51:Compare two numbers with respect to epsilon and \PB{\&{return}}\X
|
| 8 |
|
|
\U50.
|
| 9 |
|
|
\I\X53:Compute the difference of fraction parts, \PB{\|o}\X
|
| 10 |
|
|
\U51.
|
| 11 |
|
|
\I\X65:Compute the significant digits in the large-exponent case\X
|
| 12 |
|
|
\U64.
|
| 13 |
|
|
\I\X64:Compute the significant digits \PB{\|s} and decimal exponent \PB{\|e}\X
|
| 14 |
|
|
\U54.
|
| 15 |
|
|
\I\X90:Convert to short float\X
|
| 16 |
|
|
\U89.
|
| 17 |
|
|
\I\X83:Determine the binary fraction and binary exponent\X
|
| 18 |
|
|
\U78.
|
| 19 |
|
|
\I\X16:Determine the number of significant places \PB{\|n} in the divisor \PB{%
|
| 20 |
|
|
\|v}\X
|
| 21 |
|
|
\U13.
|
| 22 |
|
|
\I\X20:Determine the quotient digit \PB{\|q[\|j]}\X
|
| 23 |
|
|
\U13.
|
| 24 |
|
|
\I\X45:Divide nonzero numbers and \PB{\&{return}}\X
|
| 25 |
|
|
\U44.
|
| 26 |
|
|
\I\X48:Exchange \PB{\|y} with \PB{\|z}\X
|
| 27 |
|
|
\Us47\ET51.
|
| 28 |
|
|
\I\X55:Extract the exponent \PB{\|e} and determine the fraction interval $[f%
|
| 29 |
|
|
\dts g]$ or $(f\dts g)$\X
|
| 30 |
|
|
\U54.
|
| 31 |
|
|
\I\X21:Find the trial quotient, $\hat q$\X
|
| 32 |
|
|
\U20.
|
| 33 |
|
|
\I\X4, 9, 30, 32, 69, 75:Global variables\X
|
| 34 |
|
|
\U1.
|
| 35 |
|
|
\I\X57:Handle the special case when the fraction part is zero\X
|
| 36 |
|
|
\U55.
|
| 37 |
|
|
\I\X23:If the result was negative, decrease $\hat q$ by 1\X
|
| 38 |
|
|
\U20.
|
| 39 |
|
|
\I\X87:Integerize and \PB{\&{return}}\X
|
| 40 |
|
|
\U86.
|
| 41 |
|
|
\I\X56, 66:Local variables for \PB{\\{print\_float}}\X
|
| 42 |
|
|
\U54.
|
| 43 |
|
|
\I\X70, 76, 81:Local variables for \PB{\\{scan\_const}}\X
|
| 44 |
|
|
\U68.
|
| 45 |
|
|
\I\X79:Move the digits from \PB{\\{buf}} to \PB{$\ff$}\X
|
| 46 |
|
|
\U78.
|
| 47 |
|
|
\I\X43:Multiply nonzero numbers and \PB{\&{return}}\X
|
| 48 |
|
|
\U41.
|
| 49 |
|
|
\I\X17:Normalize the divisor\X
|
| 50 |
|
|
\U13.
|
| 51 |
|
|
\I\X36, 59:Other type definitions\X
|
| 52 |
|
|
\U1.
|
| 53 |
|
|
\I\X84:Pack and round the answer\X
|
| 54 |
|
|
\U78.
|
| 55 |
|
|
\I\X19:Pack \PB{\|q} and \PB{\|u} to \PB{\\{acc}} and \PB{\\{aux}}\X
|
| 56 |
|
|
\U13.
|
| 57 |
|
|
\I\X11:Pack \PB{\|w} into the outputs \PB{\\{aux}} and \PB{\\{acc}}\X
|
| 58 |
|
|
\U8.
|
| 59 |
|
|
\I\X67:Print the significant digits with proper context\X
|
| 60 |
|
|
\U54.
|
| 61 |
|
|
\I\X80:Put the 9-digit number \PB{${*}\|p$}\thinspace\dots\thinspace\PB{${*}(%
|
| 62 |
|
|
\|p+\T{8})$} into \PB{$\ff.\\{dat}[\|k]$}\X
|
| 63 |
|
|
\U79.
|
| 64 |
|
|
\I\X95:Reduce \PB{$(\\{ye},\\{yf})$} by a multiple of \PB{\\{zf}}; \PB{\&{goto}
|
| 65 |
|
|
\\{zero\_out}} if the remainder is zero, \PB{\&{goto} \\{try\_complement}} if
|
| 66 |
|
|
appropriate\X
|
| 67 |
|
|
\U94.
|
| 68 |
|
|
\I\X94:Remainderize nonzero numbers and \PB{\&{return}}\X
|
| 69 |
|
|
\U93.
|
| 70 |
|
|
\I\X78:Return a floating point constant\X
|
| 71 |
|
|
\U73.
|
| 72 |
|
|
\I\X72:Return infinity\X
|
| 73 |
|
|
\U68.
|
| 74 |
|
|
\I\X71:Return the standard NaN\X
|
| 75 |
|
|
\U68.
|
| 76 |
|
|
\I\X33:Round and return the result\X
|
| 77 |
|
|
\U31.
|
| 78 |
|
|
\I\X35:Round and return the short result\X
|
| 79 |
|
|
\U34.
|
| 80 |
|
|
\I\X74:Scan a fraction part\X
|
| 81 |
|
|
\U73.
|
| 82 |
|
|
\I\X73:Scan a number and \PB{\&{return}}\X
|
| 83 |
|
|
\U68.
|
| 84 |
|
|
\I\X77:Scan an exponent\X
|
| 85 |
|
|
\U73.
|
| 86 |
|
|
\I\X63:Store $f$ and $g$ as multiprecise integers\X
|
| 87 |
|
|
\U54.
|
| 88 |
|
|
\I\X2:Stuff for \CEE/ preprocessor\X
|
| 89 |
|
|
\U1.
|
| 90 |
|
|
\I\X5, 6, 7, 8, 12, 13, 24, 25, 26, 27, 28, 29, 31, 34, 37, 38, 39, 40, 41, 44,
|
| 91 |
|
|
46, 50, 54, 60, 61, 62, 68, 82, 85, 86, 88, 89, 91, 93:Subroutines\X
|
| 92 |
|
|
\U1.
|
| 93 |
|
|
\I\X22:Subtract $b^j\hat q v$ from \PB{\|u}\X
|
| 94 |
|
|
\U20.
|
| 95 |
|
|
\I\X92:Take the square root and \PB{\&{return}}\X
|
| 96 |
|
|
\U91.
|
| 97 |
|
|
\I\X3:Tetrabyte and octabyte type definitions\X
|
| 98 |
|
|
\U1.
|
| 99 |
|
|
\I\X42:The usual NaN cases\X
|
| 100 |
|
|
\Us41, 44, 46\ETs93.
|
| 101 |
|
|
\I\X18:Unnormalize the remainder\X
|
| 102 |
|
|
\U13.
|
| 103 |
|
|
\I\X15:Unpack the dividend and divisor to \PB{\|u} and \PB{\|v}\X
|
| 104 |
|
|
\U13.
|
| 105 |
|
|
\I\X10:Unpack the multiplier and multiplicand to \PB{\|u} and \PB{\|v}\X
|
| 106 |
|
|
\U8.
|
| 107 |
|
|
\I\X52:Unsubnormalize \PB{\|y} and \PB{\|z}, if they are subnormal\X
|
| 108 |
|
|
\U51.
|
| 109 |
|
|
|
