Use better mul/div routines

desktop/desktop_lc.s

@@ -254,53 +254,80 @@ op:     lda     dummy1234
 .endproc
 
 ;;; ============================================================
-;;; Input: numbers in A,X, Y
+;;; Input: numbers in A,X, Y (all unsigned)
-;;; Output: number in A,X
+;;; Output: number in A,X (unsigned)
 
 .proc Multiply_16_8_16
-        stax    num16
+        stax    num1
-        sty     num8
+        sty     num2
-        copy16  #0, accum
-        ldx     num8
-        beq     done
 
-loop:   add16   num16, accum, accum
+        ;; Accumulate directly into A,X
-        dex
+        lda     #0
+        tax
+        beq     test
+
+add:    clc
+        adc     num1
+        tay
+
+        txa
+        adc     num1+1
+        tax
+        tya
+
+loop:   asl     num1
+        rol     num1+1
+test:   lsr     num2
+        bcs     add
         bne     loop
 
-done:   ldax    accum
         rts
 
-num16:  .word   0
+num1:   .word   0
-num8:   .byte   0
+num2:   .byte   0
-accum:  .word   0
-
 .endproc
 
 ;;; ============================================================
-;;; Input: numerator in A,X, divisor in Y
+;;; Input: dividend in A,X, divisor in Y (all unsigned)
-;;; Output: result in A,X
+;;; Output: quotient in A,X (unsigned)
 
 .proc Divide_16_8_16
-        stax    numerator
+        result := dividend
+
+        stax    dividend
         sty     divisor
-        copy16  #0, result
+        lda     #0
+        sta     divisor+1
+        sta     remainder
+        sta     remainder+1
+        ldx     #16
 
-loop:   sub16_8 numerator, divisor, numerator
+loop:   asl     dividend
-        bmi     done
+        rol     dividend+1
-        inc16   result
+        rol     remainder
-        jmp     loop
+        rol     remainder+1
+        lda     remainder
+        sec
+        sbc     divisor
+        tay
+        lda     remainder+1
+        sbc     divisor+1
+        bcc     skip
+        sta     remainder+1
+        sty     remainder
+        inc     result
 
-done:   ldax    result
+skip:   dex
+        bne     loop
+        ldax    dividend
         rts
 
-numerator:
+dividend:
         .word   0
 divisor:
-        .byte   0
+        .word   0
-
+remainder:
-result: .word   0
+        .word   0
-
 .endproc
 
 ;;; ============================================================