/
/ This Source Code Form is subject to the terms of the Mozilla Public
/ License, v. 2.0. If a copy of the MPL was not distributed with this
/ file, You can obtain one at http://mozilla.org/MPL/2.0/.

.text

 /  ebp - 36:	caller's esi
 /  ebp - 32:	caller's edi
 /  ebp - 28:	
 /  ebp - 24:	
 /  ebp - 20:	
 /  ebp - 16:	
 /  ebp - 12:	
 /  ebp - 8:	
 /  ebp - 4:	
 /  ebp + 0:	caller's ebp
 /  ebp + 4:	return address
 /  ebp + 8:	a	argument
 /  ebp + 12:	a_len	argument
 /  ebp + 16:	b	argument
 /  ebp + 20:	c	argument
 /  registers:
 / 	eax:
 /	ebx:	carry
 /	ecx:	a_len
 /	edx:
 /	esi:	a ptr
 /	edi:	c ptr
.globl	s_mpv_mul_d
.type	s_mpv_mul_d,@function
s_mpv_mul_d:
    push   %ebp
    mov    %esp,%ebp
    sub    $28,%esp
    push   %edi
    push   %esi
    push   %ebx
    movl   $0,%ebx		/ carry = 0
    mov    12(%ebp),%ecx	/ ecx = a_len
    mov    20(%ebp),%edi
    cmp    $0,%ecx
    je     L2			/ jmp if a_len == 0
    mov    8(%ebp),%esi		/ esi = a
    cld
L1:
    lodsl			/ eax = [ds:esi]; esi += 4
    mov    16(%ebp),%edx	/ edx = b
    mull   %edx			/ edx:eax = Phi:Plo = a_i * b

    add    %ebx,%eax		/ add carry (%ebx) to edx:eax
    adc    $0,%edx
    mov    %edx,%ebx		/ high half of product becomes next carry

    stosl			/ [es:edi] = ax; edi += 4;
    dec    %ecx			/ --a_len
    jnz    L1			/ jmp if a_len != 0
L2:
    mov    %ebx,0(%edi)		/ *c = carry
    pop    %ebx
    pop    %esi
    pop    %edi
    leave  
    ret    
    nop

 /  ebp - 36:	caller's esi
 /  ebp - 32:	caller's edi
 /  ebp - 28:	
 /  ebp - 24:	
 /  ebp - 20:	
 /  ebp - 16:	
 /  ebp - 12:	
 /  ebp - 8:	
 /  ebp - 4:	
 /  ebp + 0:	caller's ebp
 /  ebp + 4:	return address
 /  ebp + 8:	a	argument
 /  ebp + 12:	a_len	argument
 /  ebp + 16:	b	argument
 /  ebp + 20:	c	argument
 /  registers:
 / 	eax:
 /	ebx:	carry
 /	ecx:	a_len
 /	edx:
 /	esi:	a ptr
 /	edi:	c ptr
.globl	s_mpv_mul_d_add
.type	s_mpv_mul_d_add,@function
s_mpv_mul_d_add:
    push   %ebp
    mov    %esp,%ebp
    sub    $28,%esp
    push   %edi
    push   %esi
    push   %ebx
    movl   $0,%ebx		/ carry = 0
    mov    12(%ebp),%ecx	/ ecx = a_len
    mov    20(%ebp),%edi
    cmp    $0,%ecx
    je     L4			/ jmp if a_len == 0
    mov    8(%ebp),%esi		/ esi = a
    cld
L3:
    lodsl			/ eax = [ds:esi]; esi += 4
    mov    16(%ebp),%edx	/ edx = b
    mull   %edx			/ edx:eax = Phi:Plo = a_i * b

    add    %ebx,%eax		/ add carry (%ebx) to edx:eax
    adc    $0,%edx
    mov    0(%edi),%ebx		/ add in current word from *c
    add    %ebx,%eax		
    adc    $0,%edx
    mov    %edx,%ebx		/ high half of product becomes next carry

    stosl			/ [es:edi] = ax; edi += 4;
    dec    %ecx			/ --a_len
    jnz    L3			/ jmp if a_len != 0
L4:
    mov    %ebx,0(%edi)		/ *c = carry
    pop    %ebx
    pop    %esi
    pop    %edi
    leave  
    ret    
    nop

 /  ebp - 36:	caller's esi
 /  ebp - 32:	caller's edi
 /  ebp - 28:	
 /  ebp - 24:	
 /  ebp - 20:	
 /  ebp - 16:	
 /  ebp - 12:	
 /  ebp - 8:	
 /  ebp - 4:	
 /  ebp + 0:	caller's ebp
 /  ebp + 4:	return address
 /  ebp + 8:	a	argument
 /  ebp + 12:	a_len	argument
 /  ebp + 16:	b	argument
 /  ebp + 20:	c	argument
 /  registers:
 / 	eax:
 /	ebx:	carry
 /	ecx:	a_len
 /	edx:
 /	esi:	a ptr
 /	edi:	c ptr
.globl	s_mpv_mul_d_add_prop
.type	s_mpv_mul_d_add_prop,@function
s_mpv_mul_d_add_prop:
    push   %ebp
    mov    %esp,%ebp
    sub    $28,%esp
    push   %edi
    push   %esi
    push   %ebx
    movl   $0,%ebx		/ carry = 0
    mov    12(%ebp),%ecx	/ ecx = a_len
    mov    20(%ebp),%edi
    cmp    $0,%ecx
    je     L6			/ jmp if a_len == 0
    cld
    mov    8(%ebp),%esi		/ esi = a
L5:
    lodsl			/ eax = [ds:esi]; esi += 4
    mov    16(%ebp),%edx	/ edx = b
    mull   %edx			/ edx:eax = Phi:Plo = a_i * b

    add    %ebx,%eax		/ add carry (%ebx) to edx:eax
    adc    $0,%edx
    mov    0(%edi),%ebx		/ add in current word from *c
    add    %ebx,%eax		
    adc    $0,%edx
    mov    %edx,%ebx		/ high half of product becomes next carry

    stosl			/ [es:edi] = ax; edi += 4;
    dec    %ecx			/ --a_len
    jnz    L5			/ jmp if a_len != 0
L6:
    cmp    $0,%ebx		/ is carry zero?
    jz     L8
    mov    0(%edi),%eax		/ add in current word from *c
    add	   %ebx,%eax
    stosl			/ [es:edi] = ax; edi += 4;
    jnc    L8
L7:
    mov    0(%edi),%eax		/ add in current word from *c
    adc	   $0,%eax
    stosl			/ [es:edi] = ax; edi += 4;
    jc     L7
L8:
    pop    %ebx
    pop    %esi
    pop    %edi
    leave  
    ret    
    nop

 /  ebp - 20:	caller's esi
 /  ebp - 16:	caller's edi
 /  ebp - 12:	
 /  ebp - 8:	carry
 /  ebp - 4:	a_len	local
 /  ebp + 0:	caller's ebp
 /  ebp + 4:	return address
 /  ebp + 8:	pa	argument
 /  ebp + 12:	a_len	argument
 /  ebp + 16:	ps	argument
 /  ebp + 20:	
 /  registers:
 / 	eax:
 /	ebx:	carry
 /	ecx:	a_len
 /	edx:
 /	esi:	a ptr
 /	edi:	c ptr

.globl	s_mpv_sqr_add_prop
.type	s_mpv_sqr_add_prop,@function
s_mpv_sqr_add_prop:
     push   %ebp
     mov    %esp,%ebp
     sub    $12,%esp
     push   %edi
     push   %esi
     push   %ebx
     movl   $0,%ebx		/ carry = 0
     mov    12(%ebp),%ecx	/ a_len
     mov    16(%ebp),%edi	/ edi = ps
     cmp    $0,%ecx
     je     L11			/ jump if a_len == 0
     cld
     mov    8(%ebp),%esi	/ esi = pa
L10:
     lodsl			/ %eax = [ds:si]; si += 4;
     mull   %eax

     add    %ebx,%eax		/ add "carry"
     adc    $0,%edx
     mov    0(%edi),%ebx
     add    %ebx,%eax		/ add low word from result
     mov    4(%edi),%ebx
     stosl			/ [es:di] = %eax; di += 4;
     adc    %ebx,%edx		/ add high word from result
     movl   $0,%ebx
     mov    %edx,%eax
     adc    $0,%ebx
     stosl			/ [es:di] = %eax; di += 4;
     dec    %ecx		/ --a_len
     jnz    L10			/ jmp if a_len != 0
L11:
    cmp    $0,%ebx		/ is carry zero?
    jz     L14
    mov    0(%edi),%eax		/ add in current word from *c
    add	   %ebx,%eax
    stosl			/ [es:edi] = ax; edi += 4;
    jnc    L14
L12:
    mov    0(%edi),%eax		/ add in current word from *c
    adc	   $0,%eax
    stosl			/ [es:edi] = ax; edi += 4;
    jc     L12
L14:
    pop    %ebx
    pop    %esi
    pop    %edi
    leave  
    ret    
    nop

 /
 / Divide 64-bit (Nhi,Nlo) by 32-bit divisor, which must be normalized
 / so its high bit is 1.   This code is from NSPR.
 /
 / mp_err s_mpv_div_2dx1d(mp_digit Nhi, mp_digit Nlo, mp_digit divisor,
 / 		          mp_digit *qp, mp_digit *rp)

 /  esp +  0:   Caller's ebx
 /  esp +  4:	return address
 /  esp +  8:	Nhi	argument
 /  esp + 12:	Nlo	argument
 /  esp + 16:	divisor	argument
 /  esp + 20:	qp	argument
 /  esp + 24:   rp	argument
 /  registers:
 / 	eax:
 /	ebx:	carry
 /	ecx:	a_len
 /	edx:
 /	esi:	a ptr
 /	edi:	c ptr
 / 

.globl	s_mpv_div_2dx1d
.type	s_mpv_div_2dx1d,@function
s_mpv_div_2dx1d:
       push   %ebx
       mov    8(%esp),%edx
       mov    12(%esp),%eax
       mov    16(%esp),%ebx
       div    %ebx
       mov    20(%esp),%ebx
       mov    %eax,0(%ebx)
       mov    24(%esp),%ebx
       mov    %edx,0(%ebx)
       xor    %eax,%eax		/ return zero
       pop    %ebx
       ret    
       nop