home *** CD-ROM | disk | FTP | other *** search
/ Source Code 1992 March / Source_Code_CD-ROM_Walnut_Creek_March_1992.iso / usenet / altsrcs / 3 / 3137 / sflonum.s < prev    next >
Encoding:
Text File  |  1991-03-27  |  15.0 KB  |  816 lines
  1. #
  2. #    Floating point support code.  What a crock!
  3. #
  4. #    A float looks like:
  5. #
  6. #    |S|E.x.p ... |M.a.n.t.i.s.s.a ... |
  7. #    +-+----------+--------------------+
  8. #
  9. #    where s is the sign bit, Exp is 8 bits of exponent, interpreted
  10. #    as E + 126, and Mantissa is 23 bits of fraction, with a
  11. #    hidden bit.  The point is to the left of the hidden bit.
  12. #    Doubles have another word of mantissa following.
  13. #
  14. #    All these routines have calling sequences like c routines,
  15. #    ie args on stack in backwards order, return values in d0
  16. #
  17.  
  18. # union double_di { double d; int i[2]; };
  19. # union flt_or_int { int i; float f; };
  20. #
  21. ## #ifdef L_divdf3
  22. ## double
  23. ## _divdf3 (a, b)
  24. ##      double a, b;
  25. ## {
  26. ##   return a / b;
  27. ## }
  28. ## #endif
  29. #.text
  30. #    .even
  31. #.globl ___divdf3
  32. #__divdf3:
  33. #    movel    sp@(12),sp@-
  34. #    movel    sp@(8),sp@-
  35. #    jsr    __divsf3
  36. #    addql    #8,sp
  37. #    clrl    d1        | kludge!!!
  38. #    rts            | sigh
  39. ## 
  40. ## #ifdef L_muldf3
  41. ## double
  42. ## _muldf3 (a, b)
  43. ##      double a, b;
  44. ## {
  45. ##   return a * b;
  46. ## }
  47. ## #endif
  48. #
  49. #.text
  50. #    .even
  51. #.globl __muldf3
  52. #__muldf3:
  53. #    movel    sp@(12),sp@-
  54. #    movel    sp@(8),sp@-
  55. #    jsr    __mulsf3
  56. #    addql    #8,sp
  57. #    clrl    d1        | kludge!!!
  58. #    rts            | sigh
  59. ## 
  60. ## #ifdef L_negdf2
  61. ## double
  62. ## _negdf2 (a)
  63. ##      double a;
  64. ## {
  65. ##   return -a;
  66. ## }
  67. ## #endif
  68. #
  69. #.text
  70. #    .even
  71. #.globl __negdf2
  72. #__negdf2:
  73. #    movel    sp@(8),d1        | get a lo
  74. #    movel    sp@(4),d0        | get a hi
  75. #    beq    negdf2_z        | zero, leave it
  76. #    eorl    #0x80000000,d0        | twiddle sign
  77. #negdf2_z:
  78. #    rts
  79. ## 
  80. ## #ifdef L_adddf3
  81. ## double
  82. ## _adddf3 (a, b)
  83. ##      double a, b;
  84. ## {
  85. ##   return a + b;
  86. ## }
  87. ## #endif
  88. #
  89. #.text
  90. #    .even
  91. #.globl __adddf3
  92. #__adddf3:
  93. #    movel    sp@(12),sp@-
  94. #    movel    sp@(8),sp@-
  95. #    jsr    __addsf3
  96. #    addql    #8,sp
  97. #    clrl    d1        | kludge!!!
  98. #    rts            | sigh
  99. ## 
  100. ## #ifdef L_subdf3
  101. ## double
  102. ## _subdf3 (a, b)
  103. ##      double a, b;
  104. ## {
  105. ##   return a - b;
  106. ## }
  107. ## #endif
  108. #
  109. #.text
  110. #    .even
  111. #.globl __subdf3
  112. #__subdf3:
  113. #    movel    sp@(12),sp@-
  114. #    movel    sp@(8),sp@-
  115. #    jsr    __subsf3
  116. #    addql    #8,sp
  117. #    clrl    d1        | kludge!!!
  118. #    rts            | sigh
  119. ## 
  120. ## #ifdef L_cmpdf2
  121. ## int
  122. ## _cmpdf2 (a, b)
  123. ##      double a, b;
  124. ## {
  125. ##   if (a > b)
  126. ##     return 1;
  127. ##   else if (a < b)
  128. ##     return -1;
  129. ##   return 0;
  130. ## }
  131. ## #endif
  132. #
  133. #.text
  134. #    .even
  135. #.globl __cmpdf2
  136. #__cmpdf2:
  137. #    movel    sp@(4),d0    | get b hi
  138. #    movel    sp@(12),d1    | get a hi
  139. #|
  140. #| crockery.  If both neg and not equal, this algorithm lose.  find a better one!
  141. #|
  142. #    bpl    cmpdf2_p
  143. #    tstl    d0
  144. #    bpl    cmpdf2_p
  145. #    cmpl    d1,d0
  146. #    bgt    cmpdf2_m
  147. #    blt    cmpdf2_1
  148. #    beq    cmpdf2_z
  149. #cmpdf2_p:
  150. #    cmpl    d1,d0        | get a hi
  151. #    beq    cmpdf2_z    | if eq, return 0
  152. #    bgt    cmpdf2_1    | if greater, return 1
  153. #cmpdf2_m:
  154. #    movel    #-1,d0        | else return -1
  155. #    rts
  156. #cmpdf2_z:
  157. #    clrl    d0
  158. #    rts
  159. #cmpdf2_1:
  160. #    movel    #1,d0
  161. #    rts            | sigh
  162. ## 
  163. ## #ifdef L_fixunsdfsi
  164. ## _fixunsdfsi (a)
  165. ##      double a;
  166. ## {
  167. ##   return (unsigned int) a;
  168. ## }
  169. ## #endif
  170. #
  171. #.text
  172. #    .even
  173. #.globl __fixunsdfsi
  174. #__fixunsdfsi:
  175. #    clrl d0
  176. #    clrl d1
  177. #    rts            | sigh
  178. ## 
  179. ## #ifdef L_fixunsdfdi
  180. ## double
  181. ## _fixunsdfdi (a)
  182. ##      double a;
  183. ## {
  184. ##   union double_di u;
  185. ##   u.i[LOW] = (unsigned int) a;
  186. ##   u.i[HIGH] = 0;
  187. ##   return u.d;
  188. ## }
  189. ## #endif
  190. #
  191. #.text
  192. #    .even
  193. #.globl __fixunsdfdi
  194. #__fixunsdfdi:
  195. #    clrl d0
  196. #    clrl d1
  197. #    rts            | sigh
  198. ## 
  199. ## #ifdef L_fixdfsi
  200. ## _fixdfsi (a)
  201. ##      double a;
  202. ## {
  203. ##   return (int) a;
  204. ## }
  205. ## #endif
  206. #
  207. #.text
  208. #    .even
  209. #.globl __fixdfsi
  210. #__fixdfsi:
  211. #    link    a6,#0
  212. #    movel    d2,sp@-        | save reg
  213. #    clrl    d2        | sign flag
  214. #    movel    a6@(8),d0    | get the float
  215. #    beq    fixdfsi_ret
  216. #    bpl    fixdfsi_1
  217. #    addql    #1,d2
  218. #fixdfsi_1:
  219. #    movel    d0,d1        | snag the exp
  220. #    lsrl    #7,d1
  221. #    lsrl    #8,d1
  222. #    lsrl    #8,d1
  223. #    andl    #0xFF,d1
  224. #    subl    #126,d1
  225. #    andl    #0x7FFFFF,d0    | zap cruft
  226. #    orl    #0x800000,d0    | put back hidden bit
  227. #|
  228. #| at this point the mantissa looks like 2^24 * integer value.
  229. #| if Exp is 24, were done.  If its less, we shift right,
  230. #| else left
  231. #|
  232. #fixdfsi_2:
  233. #    cmpl    #24,d1
  234. #    beq    fixdfsi_4    | were done
  235. #    bmi    fixdfsi_3    | less, shift right
  236. #    lsll    #1,d0        | greater, shift it left one
  237. #    subql    #1,d1        | and dec exp
  238. #    bra    fixdfsi_2
  239. #fixdfsi_3:
  240. #    lsrl    #1,d0        | shift right one
  241. #    addql    #1,d1        | and inc exp
  242. #    bra    fixdfsi_2
  243. #fixdfsi_4:
  244. #    tstl    d2        | negative?
  245. #    beq    fixdfsi_ret
  246. #    negl    d0
  247. #fixdfsi_ret:
  248. #    movel    sp@+,d2        | get d2 back
  249. #    unlk    a6
  250. #    rts
  251. ## 
  252. ## #ifdef L_fixdfdi
  253. ## double
  254. ## _fixdfdi (a)
  255. ##      double a;
  256. ## {
  257. ##   union double_di u;
  258. ##   u.i[LOW] = (int) a;
  259. ##   u.i[HIGH] = (int) a < 0 ? -1 : 0;
  260. ##   return u.d;
  261. ## }
  262. ## #endif
  263. #
  264. #.text
  265. #    .even
  266. #.globl __fixdfdi
  267. #__fixdfdi:
  268. #    clrl d0
  269. #    clrl d1
  270. #    rts            | sigh
  271. ## 
  272. ## #ifdef L_floatsidf
  273. ## double
  274. ## _floatsidf (a)
  275. ##      int a;
  276. ## {
  277. ##   return (double) a;
  278. ## }
  279. ## #endif
  280. #
  281. #.text
  282. #    .even
  283. #.globl __floatsidf
  284. #__floatsidf:
  285. #    link    a6,#0        | set up frame
  286. #    movel    d2,sp@-        | save d2; sign flag
  287. #    clrl    d2        | not negative yet
  288. #    movel    #24,d1        | exponent so far
  289. #    movel    a6@(8),d0    | get the int
  290. #    beq    floatsidf_ret    | zero?
  291. #    bpl    floatsidf_1    | pos, its ok
  292. #    negl    d0        | negate it
  293. #    addql    #1,d2        | bump sign
  294. #floatsidf_1:
  295. #|
  296. #| normalize right if necessary
  297. #|
  298. #    cmpl    #0xFFFFFF,d0    | too big?
  299. #    ble    floatsidf_2    | nope, see if need to slide left
  300. #    addql    #1,d1        | bump exp
  301. #    lsrl    #1,d0        | and slide mantissa right one
  302. #    bra    floatsidf_1
  303. #floatsidf_2:
  304. #    btst    #23,d0        | got a bit up here yet?
  305. #    bne    floatsidf_3    | nope, go left
  306. #    subql    #1,d1        | dec exp
  307. #    lsll    #1,d0        | and slide mantissa left one
  308. #    bra    floatsidf_2
  309. #floatsidf_3:
  310. #|
  311. #| now put it all together
  312. #|
  313. #    andl    #0x7FFFFF,d0    | zap hidden bit
  314. #    addl    #126,d1        | offset exp
  315. #    andl    #0xFF,d1    | trim it
  316. #    lsll    #8,d1        | shift up
  317. #    lsll    #8,d1
  318. #    lsll    #7,d1
  319. #    orl    d1,d0        | stuff it in
  320. #    tstl    d2        | negative?
  321. #    beq    floatsidf_ret
  322. #    orl    #0x80000000,d0
  323. #floatsidf_ret:
  324. #    movel    sp@+,d2
  325. #    clrl    d1        | ???
  326. #    unlk    a6
  327. #    rts            | sigh
  328. ## 
  329. ## #ifdef L_floatdidf
  330. ## double
  331. ## _floatdidf (u)
  332. ##      union double_di u;
  333. ## {
  334. ##   register double hi
  335. ##     = ((double) u.i[HIGH]) * (double) 0x10000 * (double) 0x10000;
  336. ##   register double low = (unsigned int) u.i[LOW];
  337. ##   return hi + low;
  338. ## }
  339. ## #endif
  340. #
  341. #.text
  342. #    .even
  343. #.globl __floatdidf
  344. #__floatdidf:
  345. #    clrl d0
  346. #    clrl d1
  347. #    rts            | sigh
  348. ## 
  349. # #define INTIFY(FLOATVAL)  (intify.f = (FLOATVAL), intify.i)
  350. # 
  351. # #ifdef L_addsf3
  352. # int
  353. # _addsf3 (a, b)
  354. #      union flt_or_int a, b;
  355. # {
  356. #   union flt_or_int intify;
  357. #   return INTIFY (a.f + b.f);
  358. # }
  359. # #endif
  360.  
  361. .text
  362.     .even
  363. .globl ___addsf3
  364. .globl __addsf3
  365. ___addsf3:
  366. __addsf3:
  367.     link    a6,#0        | dont need any locals
  368.     moveml    #0x3F00,sp@-    | save all data registers
  369.     movel    a6@(8),d0    | get a
  370.     beq    addsf_ret_b    |  zero .. just return b
  371.     movel    #23,d6        | shift count
  372.     movel    d0,d2        | get the exponent
  373.     lsrl    d6,d2        | and shift right
  374.     andl    #0xFF,d2    | no sign bit
  375.     subl    #126,d2        | offset the exponent
  376.     movel    a6@(12),d1    | get b
  377.     beq    addsf_ret_a
  378.     movel    d1,d3        | get the exponent for b
  379.     lsrl    d6,d3        | and shift right, with implicit extend
  380.     andl    #0xFF,d3    | make sure we didnt get a sign bit
  381.     subl    #126,d3        | off set this one too
  382.  
  383.     andl    #0x7FFFFF,d0    | mask a for mantissa
  384.     orl    #0x800000,d0    | and put in hidden bit
  385.     tstl    a6@(8)        | test the original value
  386.     bpl    addsf_1        | pos, ok
  387.     negl    d0        | neg, negate the mantissa
  388. addsf_1:
  389.     andl    #0x7FFFFF,d1    | mask b for mantissa
  390.     orl    #0x800000,d1    | ditto
  391.     tstl    a6@(12)        | test ...
  392.     bpl    addsf_2
  393.     negl    d1        | negate this one
  394. addsf_2:
  395.     cmpl    d2,d3        | compare Ea to Eb
  396.     blt    addsf_3        | Ea > Eb
  397.  
  398.     movel    d3,d5        | get Eb
  399.     subl    d2,d5        | subtract Ea
  400.     asrl    d5,d0        |  yielding count to shift Ma right
  401.     movel    d3,d5        | use this as resultant exponent
  402.     bra    addsf_4        | and go rejoin common part
  403. addsf_3:
  404.     movel    d2,d5        | get Ea
  405.     subl    d3,d5        | subtract Eb
  406.     asrl    d5,d1        |  yielding count to shift Mb right
  407.     movel    d2,d5        | use this as resultant exponent
  408.  
  409. addsf_4:
  410.     clrl    d7        | zap sign flag
  411.     addl    d1,d0        | add Mb to Ma
  412.  
  413.  
  414.     beq    addsf_z        | zero? ok, go return zero
  415.     bpl    addsf_5        | positive? ok, go scale it
  416.     negl    d0        | negate Mr
  417.     movel    #1,d7        | remember sign
  418. addsf_5:
  419.     btst    #24,d0        | carry?
  420.     beq    addsf_6        | nope, its ok as is
  421.     asrl    #1,d0        | shift right one
  422.     addql    #1,d5        | inc exp
  423.  
  424. | zzz check for overflow in here someplace
  425.  
  426. addsf_6:
  427.     btst    #23,d0        | got a bit in the right place yet?
  428.     bne    addsf_7        | yes, were done
  429.     lsll    #1,d0        | shift left one
  430.     subql    #1,d5        | dec exponent
  431.     bra    addsf_6
  432. addsf_7:
  433.     andl    #0x7FFFFF,d0    | zap out hidden bit
  434.     addl    #126,d5        | add offset to exp
  435.     andl    #0xFF,d5    | zap to 8 bits
  436.     movel    #23,d6        | shift count
  437.     lsll    d6,d5        | shift the exp up
  438.     orl    d5,d0        | stick the exp in
  439.     tstl    d7        | negative?
  440.     beq    addsf_ret_a
  441.     orl    #0x80000000,d0    | yup, negate it
  442.     bra    addsf_ret_a
  443. addsf_z:
  444.     clrl    d0
  445.     bra    addsf_ret_a
  446. addsf_ret_b:
  447.     movel    a6@(12),d0
  448. addsf_ret_a:
  449.     moveml    sp@+,#0x00FC    | snarf back all regs
  450.     unlk    a6
  451.     rts            | sigh
  452. # 
  453. # #ifdef L_negsf2
  454. # int
  455. # _negsf2 (a)
  456. #      union flt_or_int a;
  457. # {
  458. #   union flt_or_int intify;
  459. #   return INTIFY (-a.f);
  460. # }
  461. # #endif
  462.  
  463. .text
  464.     .even
  465. .globl ___negsf2
  466. .globl __negsf2
  467. ___negsf2:
  468. __negsf2:
  469.     movel    sp@(4),d0
  470.     beq    negsf2_z
  471.     eorl    #0x80000000,d0
  472. negsf2_z:
  473.     rts            | sigh
  474. # 
  475. # #ifdef L_subsf3
  476. # int
  477. # _subsf3 (a, b)
  478. #      union flt_or_int a, b;
  479. # {
  480. #   union flt_or_int intify;
  481. #   return INTIFY (a.f - b.f);
  482. # }
  483. # #endif
  484.  
  485. .text
  486.     .even
  487. .globl ___subsf3
  488. .globl __subsf3
  489. ___subsf3:
  490. __subsf3:
  491.     tstl    sp@(8)        | kludge.  just negate b and add
  492.     beq    subsf_bz    | zero.  dont bother
  493.     eorl    #0x80000000,sp@(8)    | negate it
  494.     jmp    ___addsf3
  495. subsf_bz:
  496.     movel    sp@(4),d0
  497.     rts
  498.  
  499. # 
  500. # #ifdef L_cmpsf2
  501. # int
  502. # _cmpsf2 (a, b)
  503. #      union flt_or_int a, b;
  504. # {
  505. #   union flt_or_int intify;
  506. #   if (a.f > b.f)
  507. #     return 1;
  508. #   else if (a.f < b.f)
  509. #     return -1;
  510. #   return 0;
  511. # }
  512. # #endif
  513.  
  514. .text
  515.     .even
  516. .globl ___cmpsf2
  517. .globl __cmpsf2
  518. ___cmpsf2:
  519. __cmpsf2:
  520.     movel    sp@(4),d0
  521.     movel    sp@(12),d1    | get a hi
  522. |
  523. | crockery.  If both neg and not equal, this algorithm lose.  find a better one!
  524. |
  525.     bpl    cmpsf2_p
  526.     tstl    d0
  527.     bpl    cmpsf2_p
  528.     cmpl    d1,d0
  529.     bgt    cmpsf2_m
  530.     blt    cmpsf2_1
  531.     beq    cmpsf2_z
  532. cmpsf2_p:
  533.     cmpl    d1,d0
  534.     beq    cmpsf2_z
  535.     bgt    cmpsf2_1
  536. cmpsf2_m:
  537.     movel    #-1,d0
  538.     rts
  539. cmpsf2_z:
  540.     clrl    d0
  541.     rts
  542. cmpsf2_1:
  543.     movel    #1,d0
  544.     rts            | sigh
  545. # 
  546. # #ifdef L_mulsf3
  547. # int
  548. # _mulsf3 (a, b)
  549. #      union flt_or_int a, b;
  550. # {
  551. #   union flt_or_int intify;
  552. #   return INTIFY (a.f * b.f);
  553. # }
  554. # #endif
  555.  
  556. .text
  557.     .even
  558. .globl ___mulsf3
  559. .globl __mulsf3
  560. ___mulsf3:
  561. __mulsf3:
  562. |
  563. | multiply.  take the numbers apart.  shift each exponent down to
  564. | 16 bits.  unsigned multiply those.  shift that down to 24 bits.
  565. | exponent is Ea + Eb.
  566. |
  567.  
  568.     link    a6,#-8        | 64 bit accum for mult
  569.     moveml    #0x3F00,sp@-    | save all data registers
  570.     movel    a6@(8),d0    | get a
  571.     beq    mulsf3_z
  572.     movel    a6@(12),d1    | get b
  573.     beq    mulsf3_z
  574.     movel    #23,d6        | shift count
  575.     movel    d0,d2        | get the exponent
  576.     lsrl    d6,d2        | and shift right
  577.     andl    #0xFF,d2
  578.     subl    #126,d2        | offset the exponent
  579.     movel    d1,d3        | get the exponent for b
  580.     lsrl    d6,d3        | and shift right
  581.     andl    #0xFF,d2
  582.     subl    #126,d3        | off set this one too
  583.  
  584.     clrl    d7        | negative result flag
  585.     andl    #0x7FFFFF,d0    | mask a for mantissa
  586.     orl    #0x800000,d0    | and put in hidden bit
  587.     tstl    a6@(8)        | test the original value
  588.     bpl    mulsf3_1    | pos, ok
  589.     eorl    #1,d7        | remember negative
  590. mulsf3_1:
  591.     andl    #0x7FFFFF,d1    | mask b for mantissa
  592.     orl    #0x800000,d1    | ditto
  593.     tstl    a6@(12)        | test ...
  594.     bpl    mulsf3_2
  595.     eorl    #1,d7
  596. mulsf3_2:
  597. |    lsrl    #8,d1        | shift this one down
  598. |    lsrl    #8,d0        | this one too...
  599. |    mulu    d1,d0        | do the multiply
  600.  
  601. |    beq    mulsf3_ret    | zero? ok, just return
  602. |    lsrl    #8,d0        | shift right again
  603.  
  604. |
  605. | we have mantissas as follows:
  606. |
  607. |    |...ah...|...al...|    |...bh...|...bl...|
  608. |
  609. | product is composed as:
  610. |
  611. |            |....al * bl....|
  612. |        |....al * bh....|
  613. |        |....ah * bl....|
  614. |    |....ah * bh....|
  615. |
  616. | then take the 24 bit chunk thats 16 bits in.
  617.  
  618.     movel    d0,d4
  619.     andl    #0xFFFF,d4    | al
  620.     movel    d1,d5
  621.     andl    #0xFFFF,d5    | bl
  622.     mulu    d5,d4        | thats al * bl
  623.     movel    d4,a6@(-4)    | into the accum
  624.     clrl    a6@(-8)        | zap the top part
  625.  
  626.     movel    d0,d4
  627.     andl    #0xFFFF,d4    | al
  628.     movel    d1,d5
  629.     movel    #16,d6        | shift count
  630.     lsrl    d6,d5        | bh
  631.     mulu    d5,d4        | al * bh
  632.     addl    d4,a6@(-6)
  633.  
  634.     movel    d0,d4
  635.     lsrl    d6,d4        | ah
  636.     movel    d1,d5
  637.     andl    #0xFFFF,d5    | bl
  638.     mulu    d5,d4        | ah * bl
  639.     addl    d4,a6@(-6)
  640.  
  641.     movel    d0,d4
  642.     lsrl    d6,d4        | ah
  643.     movel    d1,d5
  644.     lsrl    d6,d5        | bh
  645.     mulu    d5,d4        | ah * bh
  646.     addl    d4,a6@(-8)
  647.  
  648.     movel    a6@(-6),d0    | get the relevant part
  649.     lsrl    #8,d0        | and shift it down
  650.  
  651. mulsf3_norm:
  652.     btst    #23,d0        | normalized?
  653.     bne    mulsf3_ok
  654.     lsll    #1,d0
  655.     subql    #1,d2
  656.     bra    mulsf3_norm
  657.  
  658. mulsf3_ok:
  659.     andl    #0x7FFFFF,d0    | zap hidden bit
  660.     addl    d3,d2        | add Eb to Ea
  661.     addl    #126,d2        | fix offset
  662.     andl    #0xFF,d2    | whack to 8 bits
  663.     movel    #23,d6        | shift count
  664.     lsll    d6,d2        | shift up to right place
  665.     orl    d2,d0        | shove it in
  666.     tstl    d7        | sign neg?
  667.     beq    mulsf3_ret
  668.     orl    #0x80000000,d0    | set sign bit
  669.     bra    mulsf3_ret
  670. mulsf3_z:
  671.     clrl    d0
  672. mulsf3_ret:
  673.     moveml    sp@+,#0x00FC    | snarf back all regs
  674.     unlk    a6
  675.     rts            | sigh
  676. # 
  677. # #ifdef L_divsf3
  678. # int
  679. # _divsf3 (a, b)
  680. #      union flt_or_int a, b;
  681. # {
  682. #   union flt_or_int intify;
  683. #   return INTIFY (a.f / b.f);
  684. # }
  685. # #endif
  686.  
  687. .text
  688.     .even
  689. .globl ___divsf3
  690. .globl __divsf3
  691. ___divsf3:
  692. __divsf3:
  693. |
  694. | divide.  sort of like mult, exc we do shifts and subtracts to
  695. | do the division of the mantissa.  resultant exponent is Ea - Eb.
  696. |
  697.  
  698.     link    a6,#0        | dont need any locals
  699.     moveml    #0x3F00,sp@-    | save all data registers
  700.     movel    a6@(8),d0    | get a
  701.     movel    a6@(12),d1    | get b
  702.     movel    #23,d6        | shift count
  703.     movel    d0,d2        | get the exponent
  704.     lsrl    d6,d2        | and shift right
  705.     andl    #0xFF,d2
  706.     subl    #127,d2        | offset the exponent
  707.     movel    d1,d3        | get the exponent for b
  708.     lsrl    d6,d3        | and shift right
  709.     andl    #0xFF,d3
  710.     subl    #127,d3        | off set this one too
  711.  
  712.     clrl    d7        | negative result flag
  713.     andl    #0x7FFFFF,d0    | mask a for mantissa
  714.     orl    #0x800000,d0    | and put in hidden bit
  715.     tstl    a6@(8)        | test the original value
  716.     bpl    divsf3_1    | pos, ok
  717.     eorl    #1,d7        | remember negative
  718. divsf3_1:
  719.     andl    #0x7FFFFF,d1    | mask b for mantissa
  720.     orl    #0x800000,d1    | ditto
  721.     tstl    a6@(12)        | test ...
  722.     bpl    divsf3_2
  723.     eorl    #1,d7
  724. divsf3_2:
  725. |
  726. | for now, kludge.  shift Ma left and Mb right, then do an unsigned divide
  727. | and shift the result left.  Blech
  728. |
  729.  
  730. |    lsrl    #8,d1        | shift this one down
  731. |    lsll    #7,d0        | this one up
  732. |    divu    d1,d0        | do the divide
  733. |    andl    #0xFFFF,d0    | and mask off cruft
  734.  
  735. |    beq    divsf3_ret    | zero? ok, just return
  736. |    lsll    #8,d0        | shift left again
  737.  
  738. | same sort of trick as long divide, exc its easier here, cause
  739. | the numbers (mantissas) are already bit-aligned.
  740.  
  741.     clrl    d4        | accumulator
  742.     movel    #0x800000,d5    | bit
  743.     lsll    #7,d0        | buy a little extra accuracy...
  744.     lsll    #7,d1
  745. divsf3_2a:
  746.     cmpl    d1,d0        | compare dividend to divisor
  747.     bmi    divsf3_2b    | nope, no bit here
  748.     orl    d5,d4        | put in the bit
  749.     subl    d1,d0        | and subtract
  750. divsf3_2b:
  751.     lsrl    #1,d1        | slide divisor down
  752.     lsrl    #1,d5        | slide bit down
  753.     bne    divsf3_2a    | and go round again
  754.     movel    d4,d0        | leave the result here
  755.  
  756. divsf3_3:
  757.     btst    #23,d0        | right place yet?
  758.     bne    divsf3_4
  759.     lsll    #1,d0
  760.     subql    #1,d2
  761.     bra    divsf3_3
  762. divsf3_4:
  763.     andl    #0x7FFFFF,d0    | zap hidden bit
  764.     subl    d3,d2        | sub Eb from Ea
  765.     addl    #127,d2        | fix offset
  766.     andl    #0xFF,d2    | whack to 8 bits
  767.     lsll    d6,d2        | shift up to right place
  768.     orl    d2,d0        | shove it in
  769.     tstl    d7        | sign neg?
  770.     beq    divsf3_ret
  771.     orl    #0x80000000,d0    | set sign bit
  772. divsf3_ret:
  773.     moveml    sp@+,#0x00FC    | snarf back all regs
  774.     unlk    a6
  775.  
  776.     rts            | sigh
  777. ## 
  778. ## #ifdef L_truncdfsf2
  779. ## int
  780. ## _truncdfsf2 (a)
  781. ##      double a;
  782. ## {
  783. ##   union flt_or_int intify;
  784. ##   return INTIFY (a);
  785. ## }
  786. ## #endif
  787. #
  788. #.text
  789. #    .even
  790. #.globl __truncdfsf2
  791. #__truncdfsf2:
  792. #    movel    sp@(4),d0
  793. #    rts
  794. #
  795. ## 
  796. ## #ifdef L_extendsfdf2
  797. ## double
  798. ## _extendsfdf2 (a)
  799. ##      union flt_or_int a;
  800. ## {
  801. ##   union flt_or_int intify;
  802. ##   return a.f;
  803. ## }
  804. ## #endif
  805. #
  806. #.text
  807. #    .even
  808. #.globl __extendsfdf2
  809. #__extendsfdf2:
  810. #    movel    sp@(4),d0
  811. #    clrl    d1
  812. #    rts            | sigh
  813. #    .even
  814. #.text
  815.  
  816.