home *** CD-ROM | disk | FTP | other *** search
/ GEMini Atari / GEMini_Atari_CD-ROM_Walnut_Creek_December_1993.iso / files / language / forst / maths.s < prev    next >
Encoding:
Text File  |  1993-10-23  |  3.9 KB  |  229 lines
  1. ; MATHS.S: Additional integer maths and logic module
  2. ; Copyright <C> John Redmond 1989, 1990
  3. ; Public domain for non-commercial use.
  4. ;
  5.     section    text
  6.     even
  7. ;
  8. ;_XDMOD:    (u64,u64--u64,u64)
  9. _xdmod: movem.l d2-d6,-(a7)
  10.     moveq    #63,d6        ;loop counter
  11.     pop    d4
  12.     pop    d5        ;divisor
  13.     pop    d2
  14.     pop    d3        ;dividend
  15.     moveq    #0,d0
  16.     moveq    #0,d1
  17. .xdmlp: add.l    d3,d3
  18.     addx.l    d2,d2
  19.     addx.l    d1,d1    
  20.     addx.l    d0,d0        ;shift dividend left
  21.     sub.l    d5,d1
  22.     subx.l    d4,d0
  23.     bcs    .xdm5        ;if cannot divide
  24.     addq.w    #1,d3
  25.     dbra    d6,.xdmlp
  26.     bra    .xdm6
  27. .xdm5:    add.l    d5,d1        ;cancel subtraction
  28.     addx.l    d4,d0
  29.     dbra    d6,.xdmlp
  30. .xdm6:    push    d1
  31.     push    d0        ;remainder
  32.     push    d3
  33.     push    d2        ;quotient
  34.     movem.l (a7)+,d2-d6
  35.     rts
  36. ;
  37. ;_DMOD: (n32,n32--rem32,quot32)
  38. _dmod:    pop    d0
  39.     move.l    d0,-(a7)    ;save divisor sign
  40.     move.l    d0,-(a7)    ;will be quotient sign
  41.     move.l    d0,-(a7)    ;divisor
  42.     move.l    (a6),d0        ;dividend
  43.     eor.l    d0,4(a7)    ;net sign
  44.     bsr    _abs        ;abs dividend
  45.     push    (a7)+
  46.     bsr    _abs        ;abs divisor
  47.     bsr    _udmod
  48.     move.l    (a7)+,d0    ;quotient sign back
  49.     bpl    .dmz
  50.     bsr    _negate        ;quotient
  51. .dmx:    move.l    (a7)+,d0    ;divisor sign back
  52.     bmi    .dmy
  53.     pop    -(a7)        ;quotient
  54.     bsr    _negate        ;remainder
  55.     push    (a7)+        ;quotient
  56. .dmy:    rts
  57. .dmz:    move.l    (a7)+,d0    ;divisor sign back
  58.     bpl    .dmy
  59.     pop    -(a7)        ;quotient
  60.     bsr    _negate        ;remainder
  61.     push    (a7)+        ;quotient
  62.     rts
  63. ;
  64. ; _DIV: (n32,n32--n32)
  65. _div:    bsr    _dmod
  66.     move.l    (a6)+,(a6)
  67.     rts
  68. ;
  69. ; _USSMOD (u32,u32,u32--u32,u32)
  70. _ussmod: pop    -(a7)        ;divisor
  71.     bsr    _uxmult
  72.     push    (a7)+
  73.     push    #0
  74.     bsr    _xdmod        ;64-bit division
  75.     pop    d0
  76.     pop    d1
  77.     move.l    d1,(a6)        ;trim results to 32 bits
  78.     rts
  79. ;
  80. _ssmod: pop    d0        ;divisor
  81.     move.l    d0,-(a7)    ;divisor sign
  82.     move.l    d0,-(a7)    ;will be overall sign
  83.     move    4(a6),d1
  84.     eor.l    d1,(a7)
  85.     move.l    8(a6),d1
  86.     eor.l    d1,(a7)        ;result sign
  87.     move.l    d0,-(a7)    ;save divisor
  88.     pop    -(a7)
  89.     bsr    _abs
  90.     push    (a7)+
  91.     bsr    _abs
  92.     bsr    _uxmult
  93.     push    (a7)+        ;divisor
  94.     bsr    _abs
  95.     clr.l    -(a6)        ;64 bits
  96.     bsr    _xdmod        ;64-bit division
  97.     addq.l    #4,a6        ;trim quotient to 32 bits
  98.     move.l    (a7)+,d0    ;quotient sign back
  99.     bpl    .dmz
  100.     bsr    _negate        ;quotient
  101. .dmx:    move.l    (a6)+,(a6)    ;trim remainder
  102.     move.l    (a7)+,d0    ;divisor sign back
  103.     bmi    .dmy
  104.     pop    -(a7)        ;quotient
  105.     bsr    _negate        ;remainder
  106.     push    (a7)+        ;quotient
  107. .dmy:    rts
  108. .dmz:    move.l    (a6)+,(a6)    ;trim remainder
  109.     move.l    (a7)+,d0    ;divisor sign back
  110.     bpl    .dmy
  111.     pop    -(a7)        ;quotient
  112.     bsr    _negate        ;remainder
  113.     push    (a7)+        ;quotient
  114.     rts
  115. ;
  116. _ssl:    bsr    _ssmod
  117.     move.l    (a6)+,(a6)    ;drop remainder
  118.     rts
  119. ;
  120. ; _MOD: (n32,n32--n32)
  121. _mod:    bsr    _dmod
  122.     addq.l    #4,a6
  123.     rts
  124. ;
  125. _mult:    pop    d0
  126.     move.l    d0,-(a7)    ;for product sign
  127.     move.l    d0,-(a7)    ;multiplier
  128.     pop    d0
  129.     eor.l    d0,4(a7)    ;overall sign
  130.     push    d0
  131.     bsr    _abs
  132.     push    (a7)+        ;multiplier back
  133.     bsr    _abs
  134.     bsr    _uxmult
  135.     addq.l    #4,a6        ;drop top 32 bits
  136.     tst.l    (a7)+
  137.     bpl    .multx
  138.     bsr    _negate
  139. .multx:
  140.     rts
  141. ;
  142. ; _XMULT: (n32,n32--n64)
  143. _xmult: move.l    d5,-(a7)
  144.     move.l    (a6),d5
  145.     move.l    4(a6),d0    ;get net sign
  146.     eor.l    d0,d5
  147.     move.l    d5,-(a7)
  148.     bsr    _abs
  149.     move.l    (a6)+,-(a7)
  150.     bsr    _abs
  151.     move.l    (a7)+,-(a6)
  152.     bsr    _uxmult
  153.     tst.l    (a7)+        ;get sign back
  154.     bpl    .xmultx
  155.     bsr    _xnegate
  156. .xmultx: move.l (a7)+,d5
  157.     rts
  158. ;
  159. _abs:    pop    d0
  160.     bpl    .absx
  161.     neg.l    d0
  162. .absx    push    d0
  163.     rts
  164. ;
  165. _xabs:    tst.w    (a6)        ;64-bit
  166.     bpl    .xabsx
  167.     neg.l    4(a6)
  168.     negx.l    (a6)
  169. .xabsx    rts
  170. ;
  171. _min:    pop    d0
  172.     pop    d1
  173.     cmp.l    d1,d0
  174.     ble    .mx
  175.     move.l    d1,d0
  176. .mx:    push    d0
  177.     rts
  178. ;
  179. _max:    pop    d0
  180.     pop    d1
  181.     cmp.l    d1,d0
  182.     bge    .mx
  183.     move.l    d1,d0
  184. .mx:    push    d0
  185.     rts
  186. ;
  187.     section    data
  188.     even
  189. ;
  190. ;    dc.b    $81,'*'!$80
  191. ;    ptrs    _mult,14
  192. ;
  193.     dc.b    $82,'D*',$a0
  194.     ptrs    _xmult,16
  195. ;
  196.     dc.b    $84,'/MOD',$A0
  197.     ptrs    _dmod,18
  198. ;
  199.     dc.b    $82,'*/',$a0
  200.     ptrs    _ssl,16
  201. ;
  202.     dc.b    $86,'U*/MOD',$a0
  203.     ptrs    _ussmod,20
  204. ;
  205.     dc.b    $85,'*/MO','D'!$80
  206.     ptrs    _ssmod,18
  207. ;
  208. ;    dc.b    $81,'/'!$80
  209. ;    ptrs    _div,14
  210. ;
  211. ;    dc.b    $83,'MO','D'!$80
  212. ;    ptrs    _mod,16
  213. ;
  214.     dc.b    $86,'UD/MOD',$a0
  215.     ptrs    _xdmod,20
  216. ;
  217.     dc.b    $83,'MI','N'!$80
  218.     ptrs    _min,16
  219. ;
  220.     dc.b    $83,'MA','X'!$80
  221.     ptrs    _max,16
  222. ;
  223. ;    dc.b    $83,'AB','S'!$80
  224. ;    ptrs    _abs,16
  225. ;
  226.     dc.b    $84,'DABS',$a0
  227.     ptrs    _xabs,18
  228. ;
  229.