aulrem.s,v

来自「用TCPIP进行网际互联源码」· S,V 代码 · 共 73 行

head	1.1;
access;
symbols;
locks
	dls:1.1; strict;
comment	@# @;


1.1
date	97.09.21.19.27.40;	author dls;	state Dist;
branches;
next	;


desc
@@


1.1
log
@pre-3e code
@
text
@|J. Test	1/81
|addressed unsigned long remainder: *dividend = *dividend % divisor

	.globl	aulrem
	.text

aulrem:	link	a6,#0
	moveml	#0x3000,sp@@-	|need d2,d3 registers
	movl	a6@@(8),a0	|a0 = dividend pointer
	movl	a0@@,d0		|d0 = dividend
	movl	d0,d2		|save dividend
	movl	a6@@(12),d1	|divisor

	cmpl	#0x10000,d1	|divisor < 2 ** 16?
	bge	1$		|no, divisor must be < 2 ** 16
	clrw	d0		|d0 =
	swap	d0		|   dividend high
	divu	d1,d0		|yes, divide
	movw	d2,d0		|d0 = remainder high + quotient low
	divu	d1,d0		|divide
	clrw	d0		|d0 = 
	swap	d0		|   remainder
	bra	4$		|return

1$:	movl	d1,d3		|save divisor
2$:	asrl	#1,d0		|shift dividend
	asrl	#1,d1		|shift divisor
	andl	#0x7FFFFFFF,d0	|assure positive
	andl	#0x7FFFFFFF,d1	|  sign bit
	cmpl	#0x10000,d1	|divisor < 2 ** 16?
	bge	2$		|no, continue shift
	divu	d1,d0		|yes, divide
	andl	#0xFFFF,d0	|erase remainder
	movl	d0,sp@@-		|call ulmul with quotient
	movl	d3,sp@@-		|  and saved divisor on stack
	jsr	ulmul		|  as arguments
	addql	#8,sp		|restore sp
	cmpl	d0,d2		|original dividend >= lmul result?
	jge	3$		|yes, quotient should be correct
	subl	d3,d0		|no, fixup 
3$:	subl	d2,d0		|calculate
	negl	d0		|  remainder

4$:	movl	d0,a0@@		|store result through pointer
	moveml	sp@@+,#0xC	|restore registers
	unlk	a6
	rts

@