3 push(@INC,"perlasm","../../perlasm");
6 &asm_init($ARGV[0],$0);
9 for (@ARGV) { $sse2=1 if (/-DOPENSSL_IA32_SSE2/); }
11 &external_label("OPENSSL_ia32cap_P") if ($sse2);
13 &bn_mul_add_words("bn_mul_add_words");
14 &bn_mul_words("bn_mul_words");
15 &bn_sqr_words("bn_sqr_words");
16 &bn_div_words("bn_div_words");
17 &bn_add_words("bn_add_words");
18 &bn_sub_words("bn_sub_words");
19 &bn_sub_part_words("bn_sub_part_words");
27 &function_begin($name,$sse2?"EXTRN\t_OPENSSL_ia32cap_P:DWORD":"");
37 &xor($c,$c); # clear carry
38 &mov($r,&wparam(0)); #
40 &mov("ecx",&wparam(2)); #
41 &mov($a,&wparam(1)); #
43 &and("ecx",0xfffffff8); # num / 8
44 &mov($w,&wparam(3)); #
46 &push("ecx"); # Up the stack for a tmp variable
48 &jz(&label("maw_finish"));
51 &picmeup("eax","OPENSSL_ia32cap_P");
52 &bt(&DWP(0,"eax"),26);
53 &jnc(&label("maw_loop"));
55 &movd("mm0",$w); # mm0 = w
56 &pxor("mm1","mm1"); # mm1 = carry_in
58 &set_label("maw_sse2_loop",0);
59 &movd("mm3",&DWP(0,$r,"",0)); # mm3 = r[0]
60 &paddq("mm1","mm3"); # mm1 = carry_in + r[0]
61 &movd("mm2",&DWP(0,$a,"",0)); # mm2 = a[0]
62 &pmuludq("mm2","mm0"); # mm2 = w*a[0]
63 &movd("mm4",&DWP(4,$a,"",0)); # mm4 = a[1]
64 &pmuludq("mm4","mm0"); # mm4 = w*a[1]
65 &movd("mm6",&DWP(8,$a,"",0)); # mm6 = a[2]
66 &pmuludq("mm6","mm0"); # mm6 = w*a[2]
67 &movd("mm7",&DWP(12,$a,"",0)); # mm7 = a[3]
68 &pmuludq("mm7","mm0"); # mm7 = w*a[3]
69 &paddq("mm1","mm2"); # mm1 = carry_in + r[0] + w*a[0]
70 &movd("mm3",&DWP(4,$r,"",0)); # mm3 = r[1]
71 &paddq("mm3","mm4"); # mm3 = r[1] + w*a[1]
72 &movd("mm5",&DWP(8,$r,"",0)); # mm5 = r[2]
73 &paddq("mm5","mm6"); # mm5 = r[2] + w*a[2]
74 &movd("mm4",&DWP(12,$r,"",0)); # mm4 = r[3]
75 &paddq("mm7","mm4"); # mm7 = r[3] + w*a[3]
76 &movd(&DWP(0,$r,"",0),"mm1");
77 &movd("mm2",&DWP(16,$a,"",0)); # mm2 = a[4]
78 &pmuludq("mm2","mm0"); # mm2 = w*a[4]
79 &psrlq("mm1",32); # mm1 = carry0
80 &movd("mm4",&DWP(20,$a,"",0)); # mm4 = a[5]
81 &pmuludq("mm4","mm0"); # mm4 = w*a[5]
82 &paddq("mm1","mm3"); # mm1 = carry0 + r[1] + w*a[1]
83 &movd("mm6",&DWP(24,$a,"",0)); # mm6 = a[6]
84 &pmuludq("mm6","mm0"); # mm6 = w*a[6]
85 &movd(&DWP(4,$r,"",0),"mm1");
86 &psrlq("mm1",32); # mm1 = carry1
87 &movd("mm3",&DWP(28,$a,"",0)); # mm3 = a[7]
89 &pmuludq("mm3","mm0"); # mm3 = w*a[7]
90 &paddq("mm1","mm5"); # mm1 = carry1 + r[2] + w*a[2]
91 &movd("mm5",&DWP(16,$r,"",0)); # mm5 = r[4]
92 &paddq("mm2","mm5"); # mm2 = r[4] + w*a[4]
93 &movd(&DWP(8,$r,"",0),"mm1");
94 &psrlq("mm1",32); # mm1 = carry2
95 &paddq("mm1","mm7"); # mm1 = carry2 + r[3] + w*a[3]
96 &movd("mm5",&DWP(20,$r,"",0)); # mm5 = r[5]
97 &paddq("mm4","mm5"); # mm4 = r[5] + w*a[5]
98 &movd(&DWP(12,$r,"",0),"mm1");
99 &psrlq("mm1",32); # mm1 = carry3
100 &paddq("mm1","mm2"); # mm1 = carry3 + r[4] + w*a[4]
101 &movd("mm5",&DWP(24,$r,"",0)); # mm5 = r[6]
102 &paddq("mm6","mm5"); # mm6 = r[6] + w*a[6]
103 &movd(&DWP(16,$r,"",0),"mm1");
104 &psrlq("mm1",32); # mm1 = carry4
105 &paddq("mm1","mm4"); # mm1 = carry4 + r[5] + w*a[5]
106 &movd("mm5",&DWP(28,$r,"",0)); # mm5 = r[7]
107 &paddq("mm3","mm5"); # mm3 = r[7] + w*a[7]
108 &movd(&DWP(20,$r,"",0),"mm1");
109 &psrlq("mm1",32); # mm1 = carry5
110 &paddq("mm1","mm6"); # mm1 = carry5 + r[6] + w*a[6]
111 &movd(&DWP(24,$r,"",0),"mm1");
112 &psrlq("mm1",32); # mm1 = carry6
113 &paddq("mm1","mm3"); # mm1 = carry6 + r[7] + w*a[7]
114 &movd(&DWP(28,$r,"",0),"mm1");
116 &psrlq("mm1",32); # mm1 = carry_out
119 &jnz(&label("maw_sse2_loop"));
121 &movd($c,"mm1"); # c = carry_out
124 &jmp(&label("maw_finish"));
127 &set_label("maw_loop",0);
129 &mov(&swtmp(0),"ecx"); #
131 for ($i=0; $i<32; $i+=4)
133 &comment("Round $i");
135 &mov("eax",&DWP($i,$a,"",0)); # *a
137 &add("eax",$c); # L(t)+= *r
138 &mov($c,&DWP($i,$r,"",0)); # L(t)+= *r
139 &adc("edx",0); # H(t)+=carry
140 &add("eax",$c); # L(t)+=c
141 &adc("edx",0); # H(t)+=carry
142 &mov(&DWP($i,$r,"",0),"eax"); # *r= L(t);
143 &mov($c,"edx"); # c= H(t);
147 &mov("ecx",&swtmp(0)); #
151 &jnz(&label("maw_loop"));
153 &set_label("maw_finish",0);
154 &mov("ecx",&wparam(2)); # get num
156 &jnz(&label("maw_finish2")); # helps branch prediction
157 &jmp(&label("maw_end"));
159 &set_label("maw_finish2",1);
160 for ($i=0; $i<7; $i++)
162 &comment("Tail Round $i");
163 &mov("eax",&DWP($i*4,$a,"",0));# *a
165 &add("eax",$c); # L(t)+=c
166 &mov($c,&DWP($i*4,$r,"",0)); # L(t)+= *r
167 &adc("edx",0); # H(t)+=carry
169 &adc("edx",0); # H(t)+=carry
170 &dec("ecx") if ($i != 7-1);
171 &mov(&DWP($i*4,$r,"",0),"eax"); # *r= L(t);
172 &mov($c,"edx"); # c= H(t);
173 &jz(&label("maw_end")) if ($i != 7-1);
175 &set_label("maw_end",0);
178 &pop("ecx"); # clear variable from
180 &function_end($name);
187 &function_begin($name,"");
198 &xor($c,$c); # clear carry
199 &mov($r,&wparam(0)); #
200 &mov($a,&wparam(1)); #
201 &mov($num,&wparam(2)); #
202 &mov($w,&wparam(3)); #
204 &and($num,0xfffffff8); # num / 8
205 &jz(&label("mw_finish"));
207 &set_label("mw_loop",0);
208 for ($i=0; $i<32; $i+=4)
210 &comment("Round $i");
212 &mov("eax",&DWP($i,$a,"",0)); # *a
214 &add("eax",$c); # L(t)+=c
217 &adc("edx",0); # H(t)+=carry
218 &mov(&DWP($i,$r,"",0),"eax"); # *r= L(t);
220 &mov($c,"edx"); # c= H(t);
227 &jz(&label("mw_finish"));
228 &jmp(&label("mw_loop"));
230 &set_label("mw_finish",0);
231 &mov($num,&wparam(2)); # get num
233 &jnz(&label("mw_finish2"));
234 &jmp(&label("mw_end"));
236 &set_label("mw_finish2",1);
237 for ($i=0; $i<7; $i++)
239 &comment("Tail Round $i");
240 &mov("eax",&DWP($i*4,$a,"",0));# *a
242 &add("eax",$c); # L(t)+=c
244 &adc("edx",0); # H(t)+=carry
245 &mov(&DWP($i*4,$r,"",0),"eax");# *r= L(t);
246 &mov($c,"edx"); # c= H(t);
247 &dec($num) if ($i != 7-1);
248 &jz(&label("mw_end")) if ($i != 7-1);
250 &set_label("mw_end",0);
253 &function_end($name);
260 &function_begin($name,"");
267 &mov($r,&wparam(0)); #
268 &mov($a,&wparam(1)); #
269 &mov($num,&wparam(2)); #
271 &and($num,0xfffffff8); # num / 8
272 &jz(&label("sw_finish"));
274 &set_label("sw_loop",0);
275 for ($i=0; $i<32; $i+=4)
277 &comment("Round $i");
278 &mov("eax",&DWP($i,$a,"",0)); # *a
280 &mul("eax"); # *a * *a
281 &mov(&DWP($i*2,$r,"",0),"eax"); #
282 &mov(&DWP($i*2+4,$r,"",0),"edx");#
289 &jnz(&label("sw_loop"));
291 &set_label("sw_finish",0);
292 &mov($num,&wparam(2)); # get num
294 &jz(&label("sw_end"));
296 for ($i=0; $i<7; $i++)
298 &comment("Tail Round $i");
299 &mov("eax",&DWP($i*4,$a,"",0)); # *a
301 &mul("eax"); # *a * *a
302 &mov(&DWP($i*8,$r,"",0),"eax"); #
303 &dec($num) if ($i != 7-1);
304 &mov(&DWP($i*8+4,$r,"",0),"edx");
305 &jz(&label("sw_end")) if ($i != 7-1);
307 &set_label("sw_end",0);
309 &function_end($name);
316 &function_begin($name,"");
317 &mov("edx",&wparam(0)); #
318 &mov("eax",&wparam(1)); #
319 &mov("ebx",&wparam(2)); #
321 &function_end($name);
328 &function_begin($name,"");
339 &mov($r,&wparam(0)); # get r
340 &mov($a,&wparam(1)); # get a
341 &mov($b,&wparam(2)); # get b
342 &mov($num,&wparam(3)); # get num
343 &xor($c,$c); # clear carry
344 &and($num,0xfffffff8); # num / 8
346 &jz(&label("aw_finish"));
348 &set_label("aw_loop",0);
349 for ($i=0; $i<8; $i++)
351 &comment("Round $i");
353 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
354 &mov($tmp2,&DWP($i*4,$b,"",0)); # *b
360 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
368 &jnz(&label("aw_loop"));
370 &set_label("aw_finish",0);
371 &mov($num,&wparam(3)); # get num
373 &jz(&label("aw_end"));
375 for ($i=0; $i<7; $i++)
377 &comment("Tail Round $i");
378 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
379 &mov($tmp2,&DWP($i*4,$b,"",0));# *b
385 &dec($num) if ($i != 6);
386 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
387 &jz(&label("aw_end")) if ($i != 6);
389 &set_label("aw_end",0);
391 # &mov("eax",$c); # $c is "eax"
393 &function_end($name);
400 &function_begin($name,"");
411 &mov($r,&wparam(0)); # get r
412 &mov($a,&wparam(1)); # get a
413 &mov($b,&wparam(2)); # get b
414 &mov($num,&wparam(3)); # get num
415 &xor($c,$c); # clear carry
416 &and($num,0xfffffff8); # num / 8
418 &jz(&label("aw_finish"));
420 &set_label("aw_loop",0);
421 for ($i=0; $i<8; $i++)
423 &comment("Round $i");
425 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
426 &mov($tmp2,&DWP($i*4,$b,"",0)); # *b
432 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
440 &jnz(&label("aw_loop"));
442 &set_label("aw_finish",0);
443 &mov($num,&wparam(3)); # get num
445 &jz(&label("aw_end"));
447 for ($i=0; $i<7; $i++)
449 &comment("Tail Round $i");
450 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
451 &mov($tmp2,&DWP($i*4,$b,"",0));# *b
457 &dec($num) if ($i != 6);
458 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
459 &jz(&label("aw_end")) if ($i != 6);
461 &set_label("aw_end",0);
463 # &mov("eax",$c); # $c is "eax"
465 &function_end($name);
468 sub bn_sub_part_words
472 &function_begin($name,"");
483 &mov($r,&wparam(0)); # get r
484 &mov($a,&wparam(1)); # get a
485 &mov($b,&wparam(2)); # get b
486 &mov($num,&wparam(3)); # get num
487 &xor($c,$c); # clear carry
488 &and($num,0xfffffff8); # num / 8
490 &jz(&label("aw_finish"));
492 &set_label("aw_loop",0);
493 for ($i=0; $i<8; $i++)
495 &comment("Round $i");
497 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
498 &mov($tmp2,&DWP($i*4,$b,"",0)); # *b
504 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
512 &jnz(&label("aw_loop"));
514 &set_label("aw_finish",0);
515 &mov($num,&wparam(3)); # get num
517 &jz(&label("aw_end"));
519 for ($i=0; $i<7; $i++)
521 &comment("Tail Round $i");
522 &mov($tmp1,&DWP(0,$a,"",0)); # *a
523 &mov($tmp2,&DWP(0,$b,"",0));# *b
529 &mov(&DWP(0,$r,"",0),$tmp1); # *r
533 &dec($num) if ($i != 6);
534 &jz(&label("aw_end")) if ($i != 6);
536 &set_label("aw_end",0);
539 &je(&label("pw_end"));
541 &mov($num,&wparam(4)); # get dl
543 &je(&label("pw_end"));
544 &jge(&label("pw_pos"));
550 &and($num,0xfffffff8); # num / 8
551 &jz(&label("pw_neg_finish"));
553 &set_label("pw_neg_loop",0);
554 for ($i=0; $i<8; $i++)
556 &comment("dl<0 Round $i");
559 &mov($tmp2,&DWP($i*4,$b,"",0)); # *b
565 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
572 &jnz(&label("pw_neg_loop"));
574 &set_label("pw_neg_finish",0);
575 &mov($tmp2,&wparam(4)); # get dl
579 &jz(&label("pw_end"));
581 for ($i=0; $i<7; $i++)
583 &comment("dl<0 Tail Round $i");
585 &mov($tmp2,&DWP($i*4,$b,"",0));# *b
591 &dec($num) if ($i != 6);
592 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
593 &jz(&label("pw_end")) if ($i != 6);
596 &jmp(&label("pw_end"));
598 &set_label("pw_pos",0);
600 &and($num,0xfffffff8); # num / 8
601 &jz(&label("pw_pos_finish"));
603 &set_label("pw_pos_loop",0);
605 for ($i=0; $i<8; $i++)
607 &comment("dl>0 Round $i");
609 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
611 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
612 &jnc(&label("pw_nc".$i));
619 &jnz(&label("pw_pos_loop"));
621 &set_label("pw_pos_finish",0);
622 &mov($num,&wparam(4)); # get dl
624 &jz(&label("pw_end"));
626 for ($i=0; $i<7; $i++)
628 &comment("dl>0 Tail Round $i");
629 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
631 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
632 &jnc(&label("pw_tail_nc".$i));
633 &dec($num) if ($i != 6);
634 &jz(&label("pw_end")) if ($i != 6);
637 &jmp(&label("pw_end"));
639 &set_label("pw_nc_loop",0);
640 for ($i=0; $i<8; $i++)
642 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
643 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
644 &set_label("pw_nc".$i,0);
651 &jnz(&label("pw_nc_loop"));
653 &mov($num,&wparam(4)); # get dl
655 &jz(&label("pw_nc_end"));
657 for ($i=0; $i<7; $i++)
659 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
660 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
661 &set_label("pw_tail_nc".$i,0);
662 &dec($num) if ($i != 6);
663 &jz(&label("pw_nc_end")) if ($i != 6);
666 &set_label("pw_nc_end",0);
669 &set_label("pw_end",0);
671 # &mov("eax",$c); # $c is "eax"
673 &function_end($name);