x86_mont.asm 6.8 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360
  1. %ifidn __OUTPUT_FORMAT__,obj
  2. section code use32 class=code align=256
  3. %elifidn __OUTPUT_FORMAT__,win32
  4. [email protected] equ 1
  5. section .text code align=256
  6. %else
  7. section .text code
  8. %endif
  9. global _bn_mul_mont
  10. align 16
  11. _bn_mul_mont:
  12. L$_bn_mul_mont_begin:
  13. push ebp
  14. push ebx
  15. push esi
  16. push edi
  17. xor eax,eax
  18. mov edi,DWORD [40+esp]
  19. cmp edi,4
  20. jl NEAR L$000just_leave
  21. lea esi,[20+esp]
  22. lea edx,[24+esp]
  23. add edi,2
  24. neg edi
  25. lea ebp,[edi*4+esp-32]
  26. neg edi
  27. mov eax,ebp
  28. sub eax,edx
  29. and eax,2047
  30. sub ebp,eax
  31. xor edx,ebp
  32. and edx,2048
  33. xor edx,2048
  34. sub ebp,edx
  35. and ebp,-64
  36. mov eax,esp
  37. sub eax,ebp
  38. and eax,-4096
  39. mov edx,esp
  40. lea esp,[eax*1+ebp]
  41. mov eax,DWORD [esp]
  42. cmp esp,ebp
  43. ja NEAR L$001page_walk
  44. jmp NEAR L$002page_walk_done
  45. align 16
  46. L$001page_walk:
  47. lea esp,[esp-4096]
  48. mov eax,DWORD [esp]
  49. cmp esp,ebp
  50. ja NEAR L$001page_walk
  51. L$002page_walk_done:
  52. mov eax,DWORD [esi]
  53. mov ebx,DWORD [4+esi]
  54. mov ecx,DWORD [8+esi]
  55. mov ebp,DWORD [12+esi]
  56. mov esi,DWORD [16+esi]
  57. mov esi,DWORD [esi]
  58. mov DWORD [4+esp],eax
  59. mov DWORD [8+esp],ebx
  60. mov DWORD [12+esp],ecx
  61. mov DWORD [16+esp],ebp
  62. mov DWORD [20+esp],esi
  63. lea ebx,[edi-3]
  64. mov DWORD [24+esp],edx
  65. mov esi,DWORD [8+esp]
  66. lea ebp,[1+ebx]
  67. mov edi,DWORD [12+esp]
  68. xor ecx,ecx
  69. mov edx,esi
  70. and ebp,1
  71. sub edx,edi
  72. lea eax,[4+ebx*4+edi]
  73. or ebp,edx
  74. mov edi,DWORD [edi]
  75. jz NEAR L$003bn_sqr_mont
  76. mov DWORD [28+esp],eax
  77. mov eax,DWORD [esi]
  78. xor edx,edx
  79. align 16
  80. L$004mull:
  81. mov ebp,edx
  82. mul edi
  83. add ebp,eax
  84. lea ecx,[1+ecx]
  85. adc edx,0
  86. mov eax,DWORD [ecx*4+esi]
  87. cmp ecx,ebx
  88. mov DWORD [28+ecx*4+esp],ebp
  89. jl NEAR L$004mull
  90. mov ebp,edx
  91. mul edi
  92. mov edi,DWORD [20+esp]
  93. add eax,ebp
  94. mov esi,DWORD [16+esp]
  95. adc edx,0
  96. imul edi,DWORD [32+esp]
  97. mov DWORD [32+ebx*4+esp],eax
  98. xor ecx,ecx
  99. mov DWORD [36+ebx*4+esp],edx
  100. mov DWORD [40+ebx*4+esp],ecx
  101. mov eax,DWORD [esi]
  102. mul edi
  103. add eax,DWORD [32+esp]
  104. mov eax,DWORD [4+esi]
  105. adc edx,0
  106. inc ecx
  107. jmp NEAR L$0052ndmadd
  108. align 16
  109. L$0061stmadd:
  110. mov ebp,edx
  111. mul edi
  112. add ebp,DWORD [32+ecx*4+esp]
  113. lea ecx,[1+ecx]
  114. adc edx,0
  115. add ebp,eax
  116. mov eax,DWORD [ecx*4+esi]
  117. adc edx,0
  118. cmp ecx,ebx
  119. mov DWORD [28+ecx*4+esp],ebp
  120. jl NEAR L$0061stmadd
  121. mov ebp,edx
  122. mul edi
  123. add eax,DWORD [32+ebx*4+esp]
  124. mov edi,DWORD [20+esp]
  125. adc edx,0
  126. mov esi,DWORD [16+esp]
  127. add ebp,eax
  128. adc edx,0
  129. imul edi,DWORD [32+esp]
  130. xor ecx,ecx
  131. add edx,DWORD [36+ebx*4+esp]
  132. mov DWORD [32+ebx*4+esp],ebp
  133. adc ecx,0
  134. mov eax,DWORD [esi]
  135. mov DWORD [36+ebx*4+esp],edx
  136. mov DWORD [40+ebx*4+esp],ecx
  137. mul edi
  138. add eax,DWORD [32+esp]
  139. mov eax,DWORD [4+esi]
  140. adc edx,0
  141. mov ecx,1
  142. align 16
  143. L$0052ndmadd:
  144. mov ebp,edx
  145. mul edi
  146. add ebp,DWORD [32+ecx*4+esp]
  147. lea ecx,[1+ecx]
  148. adc edx,0
  149. add ebp,eax
  150. mov eax,DWORD [ecx*4+esi]
  151. adc edx,0
  152. cmp ecx,ebx
  153. mov DWORD [24+ecx*4+esp],ebp
  154. jl NEAR L$0052ndmadd
  155. mov ebp,edx
  156. mul edi
  157. add ebp,DWORD [32+ebx*4+esp]
  158. adc edx,0
  159. add ebp,eax
  160. adc edx,0
  161. mov DWORD [28+ebx*4+esp],ebp
  162. xor eax,eax
  163. mov ecx,DWORD [12+esp]
  164. add edx,DWORD [36+ebx*4+esp]
  165. adc eax,DWORD [40+ebx*4+esp]
  166. lea ecx,[4+ecx]
  167. mov DWORD [32+ebx*4+esp],edx
  168. cmp ecx,DWORD [28+esp]
  169. mov DWORD [36+ebx*4+esp],eax
  170. je NEAR L$007common_tail
  171. mov edi,DWORD [ecx]
  172. mov esi,DWORD [8+esp]
  173. mov DWORD [12+esp],ecx
  174. xor ecx,ecx
  175. xor edx,edx
  176. mov eax,DWORD [esi]
  177. jmp NEAR L$0061stmadd
  178. align 16
  179. L$003bn_sqr_mont:
  180. mov DWORD [esp],ebx
  181. mov DWORD [12+esp],ecx
  182. mov eax,edi
  183. mul edi
  184. mov DWORD [32+esp],eax
  185. mov ebx,edx
  186. shr edx,1
  187. and ebx,1
  188. inc ecx
  189. align 16
  190. L$008sqr:
  191. mov eax,DWORD [ecx*4+esi]
  192. mov ebp,edx
  193. mul edi
  194. add eax,ebp
  195. lea ecx,[1+ecx]
  196. adc edx,0
  197. lea ebp,[eax*2+ebx]
  198. shr eax,31
  199. cmp ecx,DWORD [esp]
  200. mov ebx,eax
  201. mov DWORD [28+ecx*4+esp],ebp
  202. jl NEAR L$008sqr
  203. mov eax,DWORD [ecx*4+esi]
  204. mov ebp,edx
  205. mul edi
  206. add eax,ebp
  207. mov edi,DWORD [20+esp]
  208. adc edx,0
  209. mov esi,DWORD [16+esp]
  210. lea ebp,[eax*2+ebx]
  211. imul edi,DWORD [32+esp]
  212. shr eax,31
  213. mov DWORD [32+ecx*4+esp],ebp
  214. lea ebp,[edx*2+eax]
  215. mov eax,DWORD [esi]
  216. shr edx,31
  217. mov DWORD [36+ecx*4+esp],ebp
  218. mov DWORD [40+ecx*4+esp],edx
  219. mul edi
  220. add eax,DWORD [32+esp]
  221. mov ebx,ecx
  222. adc edx,0
  223. mov eax,DWORD [4+esi]
  224. mov ecx,1
  225. align 16
  226. L$0093rdmadd:
  227. mov ebp,edx
  228. mul edi
  229. add ebp,DWORD [32+ecx*4+esp]
  230. adc edx,0
  231. add ebp,eax
  232. mov eax,DWORD [4+ecx*4+esi]
  233. adc edx,0
  234. mov DWORD [28+ecx*4+esp],ebp
  235. mov ebp,edx
  236. mul edi
  237. add ebp,DWORD [36+ecx*4+esp]
  238. lea ecx,[2+ecx]
  239. adc edx,0
  240. add ebp,eax
  241. mov eax,DWORD [ecx*4+esi]
  242. adc edx,0
  243. cmp ecx,ebx
  244. mov DWORD [24+ecx*4+esp],ebp
  245. jl NEAR L$0093rdmadd
  246. mov ebp,edx
  247. mul edi
  248. add ebp,DWORD [32+ebx*4+esp]
  249. adc edx,0
  250. add ebp,eax
  251. adc edx,0
  252. mov DWORD [28+ebx*4+esp],ebp
  253. mov ecx,DWORD [12+esp]
  254. xor eax,eax
  255. mov esi,DWORD [8+esp]
  256. add edx,DWORD [36+ebx*4+esp]
  257. adc eax,DWORD [40+ebx*4+esp]
  258. mov DWORD [32+ebx*4+esp],edx
  259. cmp ecx,ebx
  260. mov DWORD [36+ebx*4+esp],eax
  261. je NEAR L$007common_tail
  262. mov edi,DWORD [4+ecx*4+esi]
  263. lea ecx,[1+ecx]
  264. mov eax,edi
  265. mov DWORD [12+esp],ecx
  266. mul edi
  267. add eax,DWORD [32+ecx*4+esp]
  268. adc edx,0
  269. mov DWORD [32+ecx*4+esp],eax
  270. xor ebp,ebp
  271. cmp ecx,ebx
  272. lea ecx,[1+ecx]
  273. je NEAR L$010sqrlast
  274. mov ebx,edx
  275. shr edx,1
  276. and ebx,1
  277. align 16
  278. L$011sqradd:
  279. mov eax,DWORD [ecx*4+esi]
  280. mov ebp,edx
  281. mul edi
  282. add eax,ebp
  283. lea ebp,[eax*1+eax]
  284. adc edx,0
  285. shr eax,31
  286. add ebp,DWORD [32+ecx*4+esp]
  287. lea ecx,[1+ecx]
  288. adc eax,0
  289. add ebp,ebx
  290. adc eax,0
  291. cmp ecx,DWORD [esp]
  292. mov DWORD [28+ecx*4+esp],ebp
  293. mov ebx,eax
  294. jle NEAR L$011sqradd
  295. mov ebp,edx
  296. add edx,edx
  297. shr ebp,31
  298. add edx,ebx
  299. adc ebp,0
  300. L$010sqrlast:
  301. mov edi,DWORD [20+esp]
  302. mov esi,DWORD [16+esp]
  303. imul edi,DWORD [32+esp]
  304. add edx,DWORD [32+ecx*4+esp]
  305. mov eax,DWORD [esi]
  306. adc ebp,0
  307. mov DWORD [32+ecx*4+esp],edx
  308. mov DWORD [36+ecx*4+esp],ebp
  309. mul edi
  310. add eax,DWORD [32+esp]
  311. lea ebx,[ecx-1]
  312. adc edx,0
  313. mov ecx,1
  314. mov eax,DWORD [4+esi]
  315. jmp NEAR L$0093rdmadd
  316. align 16
  317. L$007common_tail:
  318. mov ebp,DWORD [16+esp]
  319. mov edi,DWORD [4+esp]
  320. lea esi,[32+esp]
  321. mov eax,DWORD [esi]
  322. mov ecx,ebx
  323. xor edx,edx
  324. align 16
  325. L$012sub:
  326. sbb eax,DWORD [edx*4+ebp]
  327. mov DWORD [edx*4+edi],eax
  328. dec ecx
  329. mov eax,DWORD [4+edx*4+esi]
  330. lea edx,[1+edx]
  331. jge NEAR L$012sub
  332. sbb eax,0
  333. mov edx,-1
  334. xor edx,eax
  335. jmp NEAR L$013copy
  336. align 16
  337. L$013copy:
  338. mov esi,DWORD [32+ebx*4+esp]
  339. mov ebp,DWORD [ebx*4+edi]
  340. mov DWORD [32+ebx*4+esp],ecx
  341. and esi,eax
  342. and ebp,edx
  343. or ebp,esi
  344. mov DWORD [ebx*4+edi],ebp
  345. dec ebx
  346. jge NEAR L$013copy
  347. mov esp,DWORD [24+esp]
  348. mov eax,1
  349. L$000just_leave:
  350. pop edi
  351. pop esi
  352. pop ebx
  353. pop ebp
  354. ret
  355. db 77,111,110,116,103,111,109,101,114,121,32,77,117,108,116,105
  356. db 112,108,105,99,97,116,105,111,110,32,102,111,114,32,120,56
  357. db 54,44,32,67,82,89,80,84,79,71,65,77,83,32,98,121
  358. db 32,60,97,112,112,114,111,64,111,112,101,110,115,115,108,46
  359. db 111,114,103,62,0