Optimize 48-bit multiplication

Author: calc84maniac

src/crt/i48mulu.src

@@ -14,12 +14,6 @@
 __i48mulu:
 	;backup af
 	push af
-
-	;backup interrupt
-	ld a, i ; P = IEF2
-	di
-	push af
-
 	push ix
 	ld ix, 0
 	add ix, sp
@@ -30,20 +24,13 @@ __i48mulu:
 	push iy
 	push bc
 
-	; load ude:uhl into ude':ubc, load uiy:ubc into uhl':ubc'
-	
-	ld bc, (ix-6) ; bc = hl, b[0], b[1], b[2]
-	exx
-	ld hl, (ix-3) ; hl' = iy, b[3], b[4], b[5]
-	ld de, (ix-9) ; de' = de, a[3], a[4], a[5]
-	ld bc, (ix-12) ; bc' = bc, a[0], a[1], a[2]
-	exx
-	
+	; bc = a[0], a[1]
+	ld a, l ; a = b[0]
+	ld iy, (ix-5) ; iy = b[1], b[2]
 
-	ld de, 0
-	push de ; sum overflow
-	push de ; upper bytes of sum at -15
-	push de ; lower bytes of sum at -18
+	or a, a
+	sbc hl, hl
+	push hl ; upper bytes of sum at -15
 	;Stack Use:
 	; ix-1 : deu b[5]
 	; ix-2 : d   b[4]
@@ -57,221 +44,172 @@ __i48mulu:
 	; ix-10: bcu a[2]
 	; ix-11: b   a[1]
 	; ix-12: c   a[0]
-	; ix-13: overflow
-	; ix-14: overflow
-	; ix-15: overflow
-	; ix-16:   sum[5]
-	; ix-17:   sum[4]
-	; ix-18:   sum[3]
-	; ix-19:   sum[2]
-	; ix-20:   sum[1]
-	; ix-21:   sum[0]
+	; ix-13:   sum[5]
+	; ix-14:   sum[4]
+	; ix-15:   sum[3]
+	; ix-16:   sum[2]
+	; ix-17:   sum[1]
+	; ix-18:   sum[0]
 
 	; sum[0-1] =============================================================================================================================
 
-	exx
-	ld iyl, c ; iyl = a[0]
-	ld iyh, b ; iyh = a[1]
-	exx
-
 	; a[0]*b[0]
-	ld d, iyl ; d = a[0]
-	ld e, c ; e = b[0]
+	ld d, c ; d = a[0]
+	ld e, a ; e = b[0]
 	mlt de
-	ld (ix-21), de ; multiplaction result in sum[0-1]
+	push de ; lower bytes of sum at -18
 
 	; sum[1-2] =============================================================================================================================
-	
-	; a[0]*b[1]
-	ld hl, (ix-20) ; hl will store current partial sum
-	ld d, iyl ; d = a[0]
-	ld e, b ; e = b[1]
+	ld l, d ; hl will store current partial sum
+
+	; a[1]*b[0]
+	ld d, b ; d = a[1]
+	ld e, a ; e = b[0]
 	mlt de
 	add hl, de
 
-	; a[1]*b[0]
-	ld d, iyh; d = a[1]
-	ld e, c ; e = b[0]
+	; a[0]*b[1]
+	ld d, c ; d = a[0]
+	ld e, iyl ; e = b[1]
 	mlt de
 	add hl, de
 
-	ld (ix-20), hl
+	ld (ix-17), hl
 
 	; sum[2-3] =============================================================================================================================
-	ld a, (ix-10) ; a = a[2]
-	ld hl, (ix-19) ; hl will store current partial sum
+	ld hl, (ix-16) ; hl will store current partial sum
 
-	; a[1]*b[1]
-	ld d, iyh ; d = a[1]
-	ld e, b ; e = b[1]
+	; a[0]*b[2]
+	ld d, c ; d = a[0]
+	ld e, iyh ; e = b[2]
 	mlt de
 	add hl, de
 
-	; a[2]*b[0]
-	ld d, a ; d = a[2]
-	ld e, c ; e = b[0]
+	; a[1]*b[1]
+	ld d, b ; d = a[1]
+	ld e, iyl ; e = b[1]
 	mlt de
 	add hl, de
 
-	; a[0]*b[2]
-	ld d, iyl ; d = a[0]
-	ld e, (ix-4) ; e = b[2]
+	; a[2]*b[0]
+	ld d, (ix-10) ; d = a[2]
+	ld e, a ; e = b[0]
 	mlt de
 	add hl, de
 
-	ld (ix-19), hl
+	ld (ix-16), hl
 
 	; sum[3-4] =============================================================================================================================
-	ld hl, (ix-18) ; hl will store current partial sum
+	ld hl, (ix-15) ; hl will store current partial sum
 
-	; a[2]*b[1]
-	ld d, a ; d = a[2]
-	ld e, b ; e = b[1]
+	; a[0]*b[3]
+	ld d, c ; d = a[0]
+	ld e, (ix-3) ; e = b[3]
 	mlt de
 	add hl, de
 
 	; a[1]*b[2]
-	ld d, iyh ; d = a[1]
-	ld e, (ix-4) ; e = b[2]
+	ld d, b ; d = a[1]
+	ld e, iyh ; e = b[2]
 	mlt de
 	add hl, de
 
-	ld a, iyl ; a = a[0]
-
-	exx
-	ld iyl, e ; iyl = a[3]
-	ex de, hl
-	; hl' = a[3], a[4], a[5]
-	; de' = b[3], b[4], b[5]
-	ld iyh, e ; iyh = b[3]
-	exx
-
-	; a[3]*b[0]
-	ld d, iyl ; d = a[3]
-	ld e, c ; e = b[0]
+	; a[2]*b[1]
+	ld d, (ix-10) ; d = a[2]
+	ld e, iyl ; e = b[1]
 	mlt de
 	add hl, de
 
-	; a[0]*b[3]
-	ld d, a ; d = a[0]
-	ld e, iyh ; e = b[3]
+	; a[3]*b[0]
+	ld d, (ix-9) ; d = a[3]
+	ld e, a ; e = b[0]
 	mlt de
 	add hl, de
 
-	ld (ix-18), hl
+	ld (ix-15), hl
 
 	; sum [4-5] =============================================================================================================================
-	ld hl, (ix-17)
+	ld hl, (ix-14) ; hl will store current partial sum
 
-	; a[3]*b[1]
-	ld d, iyl ; d = a[3]
-	ld e, b ; e = b[1]
+	; a[0]*b[4]
+	ld d, c ; d = a[0]
+	ld e, (ix-2) ; e = b[4]
 	mlt de
 	add hl, de
 
-	; a[2]*b[2]
-	ld d, (ix-10) ; d = a[2]
-	ld e, (ix-4) ; e = b[2]
+	; a[1]*b[3]
+	ld d, b ; d = a[1]
+	ld e, (ix-3) ; e = b[3]
 	mlt de
 	add hl, de
 
-	; a[1]*b[3]
-	exx
-	ld iyl, b ; iyl = a[1]
-	ld iyh, e ; iyh = b[3]
-	exx
-	ld d, iyl ; d = a[1]
-	ld e, iyh ; e = b[3]
+	; a[2]*b[2]
+	ld d, (ix-10) ; d = a[2]
+	ld e, iyh ; e = b[2]
 	mlt de
 	add hl, de
 
-	; a[0]*b[4]
-	exx
-	ld iyh, d ; iyh = b[4]
-	exx
-	ld d, a ; d = a[0]
-	ld e, iyh ; e = b[4]
+	; a[3]*b[1]
+	ld d, (ix-9) ; d = a[3]
+	ld e, iyl ; e = b[1]
 	mlt de
 	add hl, de
 
 	; a[4]*b[0]
-	exx
-	ex de, hl
-	; hl' = b[3], b[4], b[5]
-	; de' = a[3], a[4], a[5]
-	ld iyl, d ; iyl = a[4]
-	ld iyh, e ; iyh = a[3], used later
-	exx
-	ld d, iyl ; d = a[4]
-	ld e, c ; e = b[0]
+	ld d, (ix-8) ; d = a[4]
+	ld e, a ; e = b[0]
 	mlt de
 	add hl, de
 
-	ld (ix-17), hl
+	ld (ix-14), l
 
 	; sum [5] (and overflow) ====================================================================================================================
-	ld hl, (ix-16)
-
-	; a[4]*b[1]
-	ld d, iyl ; d = a[4]
-	ld e, b ; e = b[1]
-	mlt de
-	add hl, de
-
-	; a[3]*b[2]
-	ld d, iyh ; d = a[3]
-	ld e, (ix-4) ; e = b[2]
-	mlt de
-	add hl, de
-
-	; a[5]*b[0]
-	ld d, (ix-7) ; d = a[5]
-	ld e, c ; e = b[0]
-	mlt de
-	add hl, de
+	ld l, h ; hl will store current partial sum
 
 	; a[0]*b[5]
-	ld d, a ; d = a[0]
+	ld d, c ; d = a[0]
 	ld e, (ix-1) ; e = b[5]
 	mlt de
 	add hl, de
 
+	; a[1]*b[4]
+	ld c, (ix-2) ; c = b[4]
+	mlt bc
+	add hl, bc
+
 	; a[2]*b[3]
-	exx
-	ex de, hl
-	; hl' = a[3], a[4], a[5]
-	; de' = b[3], b[4], b[5]
-	ld iyl, e ; iyl = b[3]
-	ld iyh, d ; iyh = b[4], used later
-	exx
-	ld d, (ix-10) ; d = a[2]
-	ld e, iyl ; e = b[3]
+	ld bc, (ix-10) ; bc = a[2], a[3]
+	ld d, c ; d = a[2]
+	ld e, (ix-3) ; e = b[3]
 	mlt de
 	add hl, de
 
-	; a[1]*b[4]
-	exx
-	ld iyl, b ; iyl = a[1]
-	exx
-	ld d, iyl ; d = a[1]
-	ld e, iyh ; e = b[4]
+	; a[3]*b[2]
+	ld c, iyh ; c = b[2]
+	mlt bc
+	add hl, bc
+
+	; a[4]*b[1]
+	ld bc, (ix-8) ; bc = a[4], a[5]
+	ld d, c ; d = a[4]
+	ld e, iyl ; e = b[1]
 	mlt de
 	add hl, de
 
-	ld (ix-16), hl
+	; a[5]*b[0]
+	ld c, a ; c = b[0]
+	mlt bc
+	add hl, bc
+
+	ld (ix-13), l
 
 	;clean up stack and restore registers
 	pop hl
 	pop de
 	pop bc
-	pop bc
 	pop iy
+
 	ld sp, ix
 	pop ix
-
-	pop af ; interrupt stuff back, P/V = IEF2
-	jp po, .skipEI
-	ei
-.skipEI:
 	pop af
 	ret