mirror of
https://github.com/klzgrad/naiveproxy.git
synced 2024-11-24 22:36:09 +03:00
1261 lines
20 KiB
NASM
1261 lines
20 KiB
NASM
%ifidn __OUTPUT_FORMAT__,obj
|
|
section code use32 class=code align=64
|
|
%elifidn __OUTPUT_FORMAT__,win32
|
|
%ifdef __YASM_VERSION_ID__
|
|
%if __YASM_VERSION_ID__ < 01010000h
|
|
%error yasm version 1.1.0 or later needed.
|
|
%endif
|
|
; Yasm automatically includes .00 and complains about redefining it.
|
|
; https://www.tortall.net/projects/yasm/manual/html/objfmt-win32-safeseh.html
|
|
%else
|
|
$@feat.00 equ 1
|
|
%endif
|
|
section .text code align=64
|
|
%else
|
|
section .text code
|
|
%endif
|
|
global _bn_mul_comba8
|
|
align 16
|
|
_bn_mul_comba8:
|
|
L$_bn_mul_comba8_begin:
|
|
push esi
|
|
mov esi,DWORD [12+esp]
|
|
push edi
|
|
mov edi,DWORD [20+esp]
|
|
push ebp
|
|
push ebx
|
|
xor ebx,ebx
|
|
mov eax,DWORD [esi]
|
|
xor ecx,ecx
|
|
mov edx,DWORD [edi]
|
|
; ################## Calculate word 0
|
|
xor ebp,ebp
|
|
; mul a[0]*b[0]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ecx,edx
|
|
mov edx,DWORD [edi]
|
|
adc ebp,0
|
|
mov DWORD [eax],ebx
|
|
mov eax,DWORD [4+esi]
|
|
; saved r[0]
|
|
; ################## Calculate word 1
|
|
xor ebx,ebx
|
|
; mul a[1]*b[0]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [4+edi]
|
|
adc ebx,0
|
|
; mul a[0]*b[1]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebp,edx
|
|
mov edx,DWORD [edi]
|
|
adc ebx,0
|
|
mov DWORD [4+eax],ecx
|
|
mov eax,DWORD [8+esi]
|
|
; saved r[1]
|
|
; ################## Calculate word 2
|
|
xor ecx,ecx
|
|
; mul a[2]*b[0]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [4+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [4+edi]
|
|
adc ecx,0
|
|
; mul a[1]*b[1]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [8+edi]
|
|
adc ecx,0
|
|
; mul a[0]*b[2]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebx,edx
|
|
mov edx,DWORD [edi]
|
|
adc ecx,0
|
|
mov DWORD [8+eax],ebp
|
|
mov eax,DWORD [12+esi]
|
|
; saved r[2]
|
|
; ################## Calculate word 3
|
|
xor ebp,ebp
|
|
; mul a[3]*b[0]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [8+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [4+edi]
|
|
adc ebp,0
|
|
; mul a[2]*b[1]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [4+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [8+edi]
|
|
adc ebp,0
|
|
; mul a[1]*b[2]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [12+edi]
|
|
adc ebp,0
|
|
; mul a[0]*b[3]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ecx,edx
|
|
mov edx,DWORD [edi]
|
|
adc ebp,0
|
|
mov DWORD [12+eax],ebx
|
|
mov eax,DWORD [16+esi]
|
|
; saved r[3]
|
|
; ################## Calculate word 4
|
|
xor ebx,ebx
|
|
; mul a[4]*b[0]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [12+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [4+edi]
|
|
adc ebx,0
|
|
; mul a[3]*b[1]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [8+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [8+edi]
|
|
adc ebx,0
|
|
; mul a[2]*b[2]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [4+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [12+edi]
|
|
adc ebx,0
|
|
; mul a[1]*b[3]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [16+edi]
|
|
adc ebx,0
|
|
; mul a[0]*b[4]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebp,edx
|
|
mov edx,DWORD [edi]
|
|
adc ebx,0
|
|
mov DWORD [16+eax],ecx
|
|
mov eax,DWORD [20+esi]
|
|
; saved r[4]
|
|
; ################## Calculate word 5
|
|
xor ecx,ecx
|
|
; mul a[5]*b[0]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [16+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [4+edi]
|
|
adc ecx,0
|
|
; mul a[4]*b[1]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [12+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [8+edi]
|
|
adc ecx,0
|
|
; mul a[3]*b[2]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [8+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [12+edi]
|
|
adc ecx,0
|
|
; mul a[2]*b[3]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [4+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [16+edi]
|
|
adc ecx,0
|
|
; mul a[1]*b[4]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [20+edi]
|
|
adc ecx,0
|
|
; mul a[0]*b[5]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebx,edx
|
|
mov edx,DWORD [edi]
|
|
adc ecx,0
|
|
mov DWORD [20+eax],ebp
|
|
mov eax,DWORD [24+esi]
|
|
; saved r[5]
|
|
; ################## Calculate word 6
|
|
xor ebp,ebp
|
|
; mul a[6]*b[0]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [20+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [4+edi]
|
|
adc ebp,0
|
|
; mul a[5]*b[1]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [16+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [8+edi]
|
|
adc ebp,0
|
|
; mul a[4]*b[2]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [12+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [12+edi]
|
|
adc ebp,0
|
|
; mul a[3]*b[3]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [8+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [16+edi]
|
|
adc ebp,0
|
|
; mul a[2]*b[4]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [4+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [20+edi]
|
|
adc ebp,0
|
|
; mul a[1]*b[5]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [24+edi]
|
|
adc ebp,0
|
|
; mul a[0]*b[6]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ecx,edx
|
|
mov edx,DWORD [edi]
|
|
adc ebp,0
|
|
mov DWORD [24+eax],ebx
|
|
mov eax,DWORD [28+esi]
|
|
; saved r[6]
|
|
; ################## Calculate word 7
|
|
xor ebx,ebx
|
|
; mul a[7]*b[0]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [24+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [4+edi]
|
|
adc ebx,0
|
|
; mul a[6]*b[1]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [20+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [8+edi]
|
|
adc ebx,0
|
|
; mul a[5]*b[2]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [16+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [12+edi]
|
|
adc ebx,0
|
|
; mul a[4]*b[3]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [12+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [16+edi]
|
|
adc ebx,0
|
|
; mul a[3]*b[4]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [8+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [20+edi]
|
|
adc ebx,0
|
|
; mul a[2]*b[5]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [4+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [24+edi]
|
|
adc ebx,0
|
|
; mul a[1]*b[6]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [28+edi]
|
|
adc ebx,0
|
|
; mul a[0]*b[7]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebp,edx
|
|
mov edx,DWORD [4+edi]
|
|
adc ebx,0
|
|
mov DWORD [28+eax],ecx
|
|
mov eax,DWORD [28+esi]
|
|
; saved r[7]
|
|
; ################## Calculate word 8
|
|
xor ecx,ecx
|
|
; mul a[7]*b[1]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [24+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [8+edi]
|
|
adc ecx,0
|
|
; mul a[6]*b[2]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [20+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [12+edi]
|
|
adc ecx,0
|
|
; mul a[5]*b[3]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [16+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [16+edi]
|
|
adc ecx,0
|
|
; mul a[4]*b[4]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [12+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [20+edi]
|
|
adc ecx,0
|
|
; mul a[3]*b[5]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [8+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [24+edi]
|
|
adc ecx,0
|
|
; mul a[2]*b[6]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [4+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [28+edi]
|
|
adc ecx,0
|
|
; mul a[1]*b[7]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebx,edx
|
|
mov edx,DWORD [8+edi]
|
|
adc ecx,0
|
|
mov DWORD [32+eax],ebp
|
|
mov eax,DWORD [28+esi]
|
|
; saved r[8]
|
|
; ################## Calculate word 9
|
|
xor ebp,ebp
|
|
; mul a[7]*b[2]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [24+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [12+edi]
|
|
adc ebp,0
|
|
; mul a[6]*b[3]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [20+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [16+edi]
|
|
adc ebp,0
|
|
; mul a[5]*b[4]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [16+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [20+edi]
|
|
adc ebp,0
|
|
; mul a[4]*b[5]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [12+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [24+edi]
|
|
adc ebp,0
|
|
; mul a[3]*b[6]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [8+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [28+edi]
|
|
adc ebp,0
|
|
; mul a[2]*b[7]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ecx,edx
|
|
mov edx,DWORD [12+edi]
|
|
adc ebp,0
|
|
mov DWORD [36+eax],ebx
|
|
mov eax,DWORD [28+esi]
|
|
; saved r[9]
|
|
; ################## Calculate word 10
|
|
xor ebx,ebx
|
|
; mul a[7]*b[3]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [24+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [16+edi]
|
|
adc ebx,0
|
|
; mul a[6]*b[4]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [20+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [20+edi]
|
|
adc ebx,0
|
|
; mul a[5]*b[5]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [16+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [24+edi]
|
|
adc ebx,0
|
|
; mul a[4]*b[6]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [12+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [28+edi]
|
|
adc ebx,0
|
|
; mul a[3]*b[7]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebp,edx
|
|
mov edx,DWORD [16+edi]
|
|
adc ebx,0
|
|
mov DWORD [40+eax],ecx
|
|
mov eax,DWORD [28+esi]
|
|
; saved r[10]
|
|
; ################## Calculate word 11
|
|
xor ecx,ecx
|
|
; mul a[7]*b[4]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [24+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [20+edi]
|
|
adc ecx,0
|
|
; mul a[6]*b[5]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [20+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [24+edi]
|
|
adc ecx,0
|
|
; mul a[5]*b[6]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [16+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [28+edi]
|
|
adc ecx,0
|
|
; mul a[4]*b[7]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebx,edx
|
|
mov edx,DWORD [20+edi]
|
|
adc ecx,0
|
|
mov DWORD [44+eax],ebp
|
|
mov eax,DWORD [28+esi]
|
|
; saved r[11]
|
|
; ################## Calculate word 12
|
|
xor ebp,ebp
|
|
; mul a[7]*b[5]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [24+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [24+edi]
|
|
adc ebp,0
|
|
; mul a[6]*b[6]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [20+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [28+edi]
|
|
adc ebp,0
|
|
; mul a[5]*b[7]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ecx,edx
|
|
mov edx,DWORD [24+edi]
|
|
adc ebp,0
|
|
mov DWORD [48+eax],ebx
|
|
mov eax,DWORD [28+esi]
|
|
; saved r[12]
|
|
; ################## Calculate word 13
|
|
xor ebx,ebx
|
|
; mul a[7]*b[6]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [24+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [28+edi]
|
|
adc ebx,0
|
|
; mul a[6]*b[7]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebp,edx
|
|
mov edx,DWORD [28+edi]
|
|
adc ebx,0
|
|
mov DWORD [52+eax],ecx
|
|
mov eax,DWORD [28+esi]
|
|
; saved r[13]
|
|
; ################## Calculate word 14
|
|
xor ecx,ecx
|
|
; mul a[7]*b[7]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebx,edx
|
|
adc ecx,0
|
|
mov DWORD [56+eax],ebp
|
|
; saved r[14]
|
|
; save r[15]
|
|
mov DWORD [60+eax],ebx
|
|
pop ebx
|
|
pop ebp
|
|
pop edi
|
|
pop esi
|
|
ret
|
|
global _bn_mul_comba4
|
|
align 16
|
|
_bn_mul_comba4:
|
|
L$_bn_mul_comba4_begin:
|
|
push esi
|
|
mov esi,DWORD [12+esp]
|
|
push edi
|
|
mov edi,DWORD [20+esp]
|
|
push ebp
|
|
push ebx
|
|
xor ebx,ebx
|
|
mov eax,DWORD [esi]
|
|
xor ecx,ecx
|
|
mov edx,DWORD [edi]
|
|
; ################## Calculate word 0
|
|
xor ebp,ebp
|
|
; mul a[0]*b[0]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ecx,edx
|
|
mov edx,DWORD [edi]
|
|
adc ebp,0
|
|
mov DWORD [eax],ebx
|
|
mov eax,DWORD [4+esi]
|
|
; saved r[0]
|
|
; ################## Calculate word 1
|
|
xor ebx,ebx
|
|
; mul a[1]*b[0]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [4+edi]
|
|
adc ebx,0
|
|
; mul a[0]*b[1]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebp,edx
|
|
mov edx,DWORD [edi]
|
|
adc ebx,0
|
|
mov DWORD [4+eax],ecx
|
|
mov eax,DWORD [8+esi]
|
|
; saved r[1]
|
|
; ################## Calculate word 2
|
|
xor ecx,ecx
|
|
; mul a[2]*b[0]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [4+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [4+edi]
|
|
adc ecx,0
|
|
; mul a[1]*b[1]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [8+edi]
|
|
adc ecx,0
|
|
; mul a[0]*b[2]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebx,edx
|
|
mov edx,DWORD [edi]
|
|
adc ecx,0
|
|
mov DWORD [8+eax],ebp
|
|
mov eax,DWORD [12+esi]
|
|
; saved r[2]
|
|
; ################## Calculate word 3
|
|
xor ebp,ebp
|
|
; mul a[3]*b[0]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [8+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [4+edi]
|
|
adc ebp,0
|
|
; mul a[2]*b[1]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [4+esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [8+edi]
|
|
adc ebp,0
|
|
; mul a[1]*b[2]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [esi]
|
|
adc ecx,edx
|
|
mov edx,DWORD [12+edi]
|
|
adc ebp,0
|
|
; mul a[0]*b[3]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ecx,edx
|
|
mov edx,DWORD [4+edi]
|
|
adc ebp,0
|
|
mov DWORD [12+eax],ebx
|
|
mov eax,DWORD [12+esi]
|
|
; saved r[3]
|
|
; ################## Calculate word 4
|
|
xor ebx,ebx
|
|
; mul a[3]*b[1]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [8+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [8+edi]
|
|
adc ebx,0
|
|
; mul a[2]*b[2]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [4+esi]
|
|
adc ebp,edx
|
|
mov edx,DWORD [12+edi]
|
|
adc ebx,0
|
|
; mul a[1]*b[3]
|
|
mul edx
|
|
add ecx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebp,edx
|
|
mov edx,DWORD [8+edi]
|
|
adc ebx,0
|
|
mov DWORD [16+eax],ecx
|
|
mov eax,DWORD [12+esi]
|
|
; saved r[4]
|
|
; ################## Calculate word 5
|
|
xor ecx,ecx
|
|
; mul a[3]*b[2]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [8+esi]
|
|
adc ebx,edx
|
|
mov edx,DWORD [12+edi]
|
|
adc ecx,0
|
|
; mul a[2]*b[3]
|
|
mul edx
|
|
add ebp,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ebx,edx
|
|
mov edx,DWORD [12+edi]
|
|
adc ecx,0
|
|
mov DWORD [20+eax],ebp
|
|
mov eax,DWORD [12+esi]
|
|
; saved r[5]
|
|
; ################## Calculate word 6
|
|
xor ebp,ebp
|
|
; mul a[3]*b[3]
|
|
mul edx
|
|
add ebx,eax
|
|
mov eax,DWORD [20+esp]
|
|
adc ecx,edx
|
|
adc ebp,0
|
|
mov DWORD [24+eax],ebx
|
|
; saved r[6]
|
|
; save r[7]
|
|
mov DWORD [28+eax],ecx
|
|
pop ebx
|
|
pop ebp
|
|
pop edi
|
|
pop esi
|
|
ret
|
|
global _bn_sqr_comba8
|
|
align 16
|
|
_bn_sqr_comba8:
|
|
L$_bn_sqr_comba8_begin:
|
|
push esi
|
|
push edi
|
|
push ebp
|
|
push ebx
|
|
mov edi,DWORD [20+esp]
|
|
mov esi,DWORD [24+esp]
|
|
xor ebx,ebx
|
|
xor ecx,ecx
|
|
mov eax,DWORD [esi]
|
|
; ############### Calculate word 0
|
|
xor ebp,ebp
|
|
; sqr a[0]*a[0]
|
|
mul eax
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov edx,DWORD [esi]
|
|
adc ebp,0
|
|
mov DWORD [edi],ebx
|
|
mov eax,DWORD [4+esi]
|
|
; saved r[0]
|
|
; ############### Calculate word 1
|
|
xor ebx,ebx
|
|
; sqr a[1]*a[0]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebx,0
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov eax,DWORD [8+esi]
|
|
adc ebx,0
|
|
mov DWORD [4+edi],ecx
|
|
mov edx,DWORD [esi]
|
|
; saved r[1]
|
|
; ############### Calculate word 2
|
|
xor ecx,ecx
|
|
; sqr a[2]*a[0]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ecx,0
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov eax,DWORD [4+esi]
|
|
adc ecx,0
|
|
; sqr a[1]*a[1]
|
|
mul eax
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov edx,DWORD [esi]
|
|
adc ecx,0
|
|
mov DWORD [8+edi],ebp
|
|
mov eax,DWORD [12+esi]
|
|
; saved r[2]
|
|
; ############### Calculate word 3
|
|
xor ebp,ebp
|
|
; sqr a[3]*a[0]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebp,0
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov eax,DWORD [8+esi]
|
|
adc ebp,0
|
|
mov edx,DWORD [4+esi]
|
|
; sqr a[2]*a[1]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebp,0
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov eax,DWORD [16+esi]
|
|
adc ebp,0
|
|
mov DWORD [12+edi],ebx
|
|
mov edx,DWORD [esi]
|
|
; saved r[3]
|
|
; ############### Calculate word 4
|
|
xor ebx,ebx
|
|
; sqr a[4]*a[0]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebx,0
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov eax,DWORD [12+esi]
|
|
adc ebx,0
|
|
mov edx,DWORD [4+esi]
|
|
; sqr a[3]*a[1]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebx,0
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov eax,DWORD [8+esi]
|
|
adc ebx,0
|
|
; sqr a[2]*a[2]
|
|
mul eax
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov edx,DWORD [esi]
|
|
adc ebx,0
|
|
mov DWORD [16+edi],ecx
|
|
mov eax,DWORD [20+esi]
|
|
; saved r[4]
|
|
; ############### Calculate word 5
|
|
xor ecx,ecx
|
|
; sqr a[5]*a[0]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ecx,0
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov eax,DWORD [16+esi]
|
|
adc ecx,0
|
|
mov edx,DWORD [4+esi]
|
|
; sqr a[4]*a[1]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ecx,0
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov eax,DWORD [12+esi]
|
|
adc ecx,0
|
|
mov edx,DWORD [8+esi]
|
|
; sqr a[3]*a[2]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ecx,0
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov eax,DWORD [24+esi]
|
|
adc ecx,0
|
|
mov DWORD [20+edi],ebp
|
|
mov edx,DWORD [esi]
|
|
; saved r[5]
|
|
; ############### Calculate word 6
|
|
xor ebp,ebp
|
|
; sqr a[6]*a[0]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebp,0
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov eax,DWORD [20+esi]
|
|
adc ebp,0
|
|
mov edx,DWORD [4+esi]
|
|
; sqr a[5]*a[1]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebp,0
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov eax,DWORD [16+esi]
|
|
adc ebp,0
|
|
mov edx,DWORD [8+esi]
|
|
; sqr a[4]*a[2]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebp,0
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov eax,DWORD [12+esi]
|
|
adc ebp,0
|
|
; sqr a[3]*a[3]
|
|
mul eax
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov edx,DWORD [esi]
|
|
adc ebp,0
|
|
mov DWORD [24+edi],ebx
|
|
mov eax,DWORD [28+esi]
|
|
; saved r[6]
|
|
; ############### Calculate word 7
|
|
xor ebx,ebx
|
|
; sqr a[7]*a[0]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebx,0
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov eax,DWORD [24+esi]
|
|
adc ebx,0
|
|
mov edx,DWORD [4+esi]
|
|
; sqr a[6]*a[1]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebx,0
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov eax,DWORD [20+esi]
|
|
adc ebx,0
|
|
mov edx,DWORD [8+esi]
|
|
; sqr a[5]*a[2]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebx,0
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov eax,DWORD [16+esi]
|
|
adc ebx,0
|
|
mov edx,DWORD [12+esi]
|
|
; sqr a[4]*a[3]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebx,0
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov eax,DWORD [28+esi]
|
|
adc ebx,0
|
|
mov DWORD [28+edi],ecx
|
|
mov edx,DWORD [4+esi]
|
|
; saved r[7]
|
|
; ############### Calculate word 8
|
|
xor ecx,ecx
|
|
; sqr a[7]*a[1]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ecx,0
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov eax,DWORD [24+esi]
|
|
adc ecx,0
|
|
mov edx,DWORD [8+esi]
|
|
; sqr a[6]*a[2]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ecx,0
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov eax,DWORD [20+esi]
|
|
adc ecx,0
|
|
mov edx,DWORD [12+esi]
|
|
; sqr a[5]*a[3]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ecx,0
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov eax,DWORD [16+esi]
|
|
adc ecx,0
|
|
; sqr a[4]*a[4]
|
|
mul eax
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov edx,DWORD [8+esi]
|
|
adc ecx,0
|
|
mov DWORD [32+edi],ebp
|
|
mov eax,DWORD [28+esi]
|
|
; saved r[8]
|
|
; ############### Calculate word 9
|
|
xor ebp,ebp
|
|
; sqr a[7]*a[2]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebp,0
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov eax,DWORD [24+esi]
|
|
adc ebp,0
|
|
mov edx,DWORD [12+esi]
|
|
; sqr a[6]*a[3]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebp,0
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov eax,DWORD [20+esi]
|
|
adc ebp,0
|
|
mov edx,DWORD [16+esi]
|
|
; sqr a[5]*a[4]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebp,0
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov eax,DWORD [28+esi]
|
|
adc ebp,0
|
|
mov DWORD [36+edi],ebx
|
|
mov edx,DWORD [12+esi]
|
|
; saved r[9]
|
|
; ############### Calculate word 10
|
|
xor ebx,ebx
|
|
; sqr a[7]*a[3]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebx,0
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov eax,DWORD [24+esi]
|
|
adc ebx,0
|
|
mov edx,DWORD [16+esi]
|
|
; sqr a[6]*a[4]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebx,0
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov eax,DWORD [20+esi]
|
|
adc ebx,0
|
|
; sqr a[5]*a[5]
|
|
mul eax
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov edx,DWORD [16+esi]
|
|
adc ebx,0
|
|
mov DWORD [40+edi],ecx
|
|
mov eax,DWORD [28+esi]
|
|
; saved r[10]
|
|
; ############### Calculate word 11
|
|
xor ecx,ecx
|
|
; sqr a[7]*a[4]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ecx,0
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov eax,DWORD [24+esi]
|
|
adc ecx,0
|
|
mov edx,DWORD [20+esi]
|
|
; sqr a[6]*a[5]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ecx,0
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov eax,DWORD [28+esi]
|
|
adc ecx,0
|
|
mov DWORD [44+edi],ebp
|
|
mov edx,DWORD [20+esi]
|
|
; saved r[11]
|
|
; ############### Calculate word 12
|
|
xor ebp,ebp
|
|
; sqr a[7]*a[5]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebp,0
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov eax,DWORD [24+esi]
|
|
adc ebp,0
|
|
; sqr a[6]*a[6]
|
|
mul eax
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov edx,DWORD [24+esi]
|
|
adc ebp,0
|
|
mov DWORD [48+edi],ebx
|
|
mov eax,DWORD [28+esi]
|
|
; saved r[12]
|
|
; ############### Calculate word 13
|
|
xor ebx,ebx
|
|
; sqr a[7]*a[6]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebx,0
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov eax,DWORD [28+esi]
|
|
adc ebx,0
|
|
mov DWORD [52+edi],ecx
|
|
; saved r[13]
|
|
; ############### Calculate word 14
|
|
xor ecx,ecx
|
|
; sqr a[7]*a[7]
|
|
mul eax
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
adc ecx,0
|
|
mov DWORD [56+edi],ebp
|
|
; saved r[14]
|
|
mov DWORD [60+edi],ebx
|
|
pop ebx
|
|
pop ebp
|
|
pop edi
|
|
pop esi
|
|
ret
|
|
global _bn_sqr_comba4
|
|
align 16
|
|
_bn_sqr_comba4:
|
|
L$_bn_sqr_comba4_begin:
|
|
push esi
|
|
push edi
|
|
push ebp
|
|
push ebx
|
|
mov edi,DWORD [20+esp]
|
|
mov esi,DWORD [24+esp]
|
|
xor ebx,ebx
|
|
xor ecx,ecx
|
|
mov eax,DWORD [esi]
|
|
; ############### Calculate word 0
|
|
xor ebp,ebp
|
|
; sqr a[0]*a[0]
|
|
mul eax
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov edx,DWORD [esi]
|
|
adc ebp,0
|
|
mov DWORD [edi],ebx
|
|
mov eax,DWORD [4+esi]
|
|
; saved r[0]
|
|
; ############### Calculate word 1
|
|
xor ebx,ebx
|
|
; sqr a[1]*a[0]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebx,0
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov eax,DWORD [8+esi]
|
|
adc ebx,0
|
|
mov DWORD [4+edi],ecx
|
|
mov edx,DWORD [esi]
|
|
; saved r[1]
|
|
; ############### Calculate word 2
|
|
xor ecx,ecx
|
|
; sqr a[2]*a[0]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ecx,0
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov eax,DWORD [4+esi]
|
|
adc ecx,0
|
|
; sqr a[1]*a[1]
|
|
mul eax
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov edx,DWORD [esi]
|
|
adc ecx,0
|
|
mov DWORD [8+edi],ebp
|
|
mov eax,DWORD [12+esi]
|
|
; saved r[2]
|
|
; ############### Calculate word 3
|
|
xor ebp,ebp
|
|
; sqr a[3]*a[0]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebp,0
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov eax,DWORD [8+esi]
|
|
adc ebp,0
|
|
mov edx,DWORD [4+esi]
|
|
; sqr a[2]*a[1]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebp,0
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
mov eax,DWORD [12+esi]
|
|
adc ebp,0
|
|
mov DWORD [12+edi],ebx
|
|
mov edx,DWORD [4+esi]
|
|
; saved r[3]
|
|
; ############### Calculate word 4
|
|
xor ebx,ebx
|
|
; sqr a[3]*a[1]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ebx,0
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov eax,DWORD [8+esi]
|
|
adc ebx,0
|
|
; sqr a[2]*a[2]
|
|
mul eax
|
|
add ecx,eax
|
|
adc ebp,edx
|
|
mov edx,DWORD [8+esi]
|
|
adc ebx,0
|
|
mov DWORD [16+edi],ecx
|
|
mov eax,DWORD [12+esi]
|
|
; saved r[4]
|
|
; ############### Calculate word 5
|
|
xor ecx,ecx
|
|
; sqr a[3]*a[2]
|
|
mul edx
|
|
add eax,eax
|
|
adc edx,edx
|
|
adc ecx,0
|
|
add ebp,eax
|
|
adc ebx,edx
|
|
mov eax,DWORD [12+esi]
|
|
adc ecx,0
|
|
mov DWORD [20+edi],ebp
|
|
; saved r[5]
|
|
; ############### Calculate word 6
|
|
xor ebp,ebp
|
|
; sqr a[3]*a[3]
|
|
mul eax
|
|
add ebx,eax
|
|
adc ecx,edx
|
|
adc ebp,0
|
|
mov DWORD [24+edi],ebx
|
|
; saved r[6]
|
|
mov DWORD [28+edi],ecx
|
|
pop ebx
|
|
pop ebp
|
|
pop edi
|
|
pop esi
|
|
ret
|