3 push(@INC,"perlasm","../../perlasm");
6 &asm_init($ARGV[0],$0);
9 for (@ARGV) { $sse2=1 if (/-DOPENSSL_IA32_SSE2/); }
11 &bn_mul_add_words("bn_mul_add_words");
12 &bn_mul_words("bn_mul_words");
13 &bn_sqr_words("bn_sqr_words");
14 &bn_div_words("bn_div_words");
15 &bn_add_words("bn_add_words");
16 &bn_sub_words("bn_sub_words");
17 &bn_sub_part_words("bn_sub_part_words");
25 &function_begin($name,"");
35 &xor($c,$c); # clear carry
36 &mov($r,&wparam(0)); #
38 &mov("ecx",&wparam(2)); #
39 &mov($a,&wparam(1)); #
41 &and("ecx",0xfffffff8); # num / 8
42 &mov($w,&wparam(3)); #
44 &push("ecx"); # Up the stack for a tmp variable
46 &jz(&label("maw_finish"));
49 &picmeup("eax","OPENSSL_ia32cap");
50 &bt(&DWP(0,"eax"),26);
51 &jnc(&label("maw_loop"));
53 &movd("mm0",$w); # mm0 = w
54 &pxor("mm1","mm1"); # mm1 = carry_in
56 &set_label("maw_sse2_loop",0);
57 &movd("mm3",&DWP(0,$r,"",0)); # mm3 = r[0]
58 &paddq("mm1","mm3"); # mm1 = carry_in + r[0]
59 &movd("mm2",&DWP(0,$a,"",0)); # mm2 = a[0]
60 &pmuludq("mm2","mm0"); # mm2 = w*a[0]
61 &movd("mm4",&DWP(4,$a,"",0)); # mm4 = a[1]
62 &pmuludq("mm4","mm0"); # mm4 = w*a[1]
63 &movd("mm6",&DWP(8,$a,"",0)); # mm6 = a[2]
64 &pmuludq("mm6","mm0"); # mm6 = w*a[2]
65 &movd("mm7",&DWP(12,$a,"",0)); # mm7 = a[3]
66 &pmuludq("mm7","mm0"); # mm7 = w*a[3]
67 &paddq("mm1","mm2"); # mm1 = carry_in + r[0] + w*a[0]
68 &movd("mm3",&DWP(4,$r,"",0)); # mm3 = r[1]
69 &paddq("mm3","mm4"); # mm3 = r[1] + w*a[1]
70 &movd("mm5",&DWP(8,$r,"",0)); # mm5 = r[2]
71 &paddq("mm5","mm6"); # mm5 = r[2] + w*a[2]
72 &movd("mm4",&DWP(12,$r,"",0)); # mm4 = r[3]
73 &paddq("mm7","mm4"); # mm7 = r[3] + w*a[3]
74 &movd(&DWP(0,$r,"",0),"mm1");
75 &movd("mm2",&DWP(16,$a,"",0)); # mm2 = a[4]
76 &pmuludq("mm2","mm0"); # mm2 = w*a[4]
77 &psrlq("mm1",32); # mm1 = carry0
78 &movd("mm4",&DWP(20,$a,"",0)); # mm4 = a[5]
79 &pmuludq("mm4","mm0"); # mm4 = w*a[5]
80 &paddq("mm1","mm3"); # mm1 = carry0 + r[1] + w*a[1]
81 &movd("mm6",&DWP(24,$a,"",0)); # mm6 = a[6]
82 &pmuludq("mm6","mm0"); # mm6 = w*a[6]
83 &movd(&DWP(4,$r,"",0),"mm1");
84 &psrlq("mm1",32); # mm1 = carry1
85 &movd("mm3",&DWP(28,$a,"",0)); # mm3 = a[7]
87 &pmuludq("mm3","mm0"); # mm3 = w*a[7]
88 &paddq("mm1","mm5"); # mm1 = carry1 + r[2] + w*a[2]
89 &movd("mm5",&DWP(16,$r,"",0)); # mm5 = r[4]
90 &paddq("mm2","mm5"); # mm2 = r[4] + w*a[4]
91 &movd(&DWP(8,$r,"",0),"mm1");
92 &psrlq("mm1",32); # mm1 = carry2
93 &paddq("mm1","mm7"); # mm1 = carry2 + r[3] + w*a[3]
94 &movd("mm5",&DWP(20,$r,"",0)); # mm5 = r[5]
95 &paddq("mm4","mm5"); # mm4 = r[5] + w*a[5]
96 &movd(&DWP(12,$r,"",0),"mm1");
97 &psrlq("mm1",32); # mm1 = carry3
98 &paddq("mm1","mm2"); # mm1 = carry3 + r[4] + w*a[4]
99 &movd("mm5",&DWP(24,$r,"",0)); # mm5 = r[6]
100 &paddq("mm6","mm5"); # mm6 = r[6] + w*a[6]
101 &movd(&DWP(16,$r,"",0),"mm1");
102 &psrlq("mm1",32); # mm1 = carry4
103 &paddq("mm1","mm4"); # mm1 = carry4 + r[5] + w*a[5]
104 &movd("mm5",&DWP(28,$r,"",0)); # mm5 = r[7]
105 &paddq("mm3","mm5"); # mm3 = r[7] + w*a[7]
106 &movd(&DWP(20,$r,"",0),"mm1");
107 &psrlq("mm1",32); # mm1 = carry5
108 &paddq("mm1","mm6"); # mm1 = carry5 + r[6] + w*a[6]
109 &movd(&DWP(24,$r,"",0),"mm1");
110 &psrlq("mm1",32); # mm1 = carry6
111 &paddq("mm1","mm3"); # mm1 = carry6 + r[7] + w*a[7]
112 &movd(&DWP(28,$r,"",0),"mm1");
114 &psrlq("mm1",32); # mm1 = carry_out
117 &jnz(&label("maw_sse2_loop"));
119 &movd($c,"mm1"); # c = carry_out
122 &jmp(&label("maw_finish"));
125 &set_label("maw_loop",0);
127 &mov(&swtmp(0),"ecx"); #
129 for ($i=0; $i<32; $i+=4)
131 &comment("Round $i");
133 &mov("eax",&DWP($i,$a,"",0)); # *a
135 &add("eax",$c); # L(t)+= *r
136 &mov($c,&DWP($i,$r,"",0)); # L(t)+= *r
137 &adc("edx",0); # H(t)+=carry
138 &add("eax",$c); # L(t)+=c
139 &adc("edx",0); # H(t)+=carry
140 &mov(&DWP($i,$r,"",0),"eax"); # *r= L(t);
141 &mov($c,"edx"); # c= H(t);
145 &mov("ecx",&swtmp(0)); #
149 &jnz(&label("maw_loop"));
151 &set_label("maw_finish",0);
152 &mov("ecx",&wparam(2)); # get num
154 &jnz(&label("maw_finish2")); # helps branch prediction
155 &jmp(&label("maw_end"));
157 &set_label("maw_finish2",1);
158 for ($i=0; $i<7; $i++)
160 &comment("Tail Round $i");
161 &mov("eax",&DWP($i*4,$a,"",0));# *a
163 &add("eax",$c); # L(t)+=c
164 &mov($c,&DWP($i*4,$r,"",0)); # L(t)+= *r
165 &adc("edx",0); # H(t)+=carry
167 &adc("edx",0); # H(t)+=carry
168 &dec("ecx") if ($i != 7-1);
169 &mov(&DWP($i*4,$r,"",0),"eax"); # *r= L(t);
170 &mov($c,"edx"); # c= H(t);
171 &jz(&label("maw_end")) if ($i != 7-1);
173 &set_label("maw_end",0);
176 &pop("ecx"); # clear variable from
178 &function_end($name);
185 &function_begin($name,"");
196 &xor($c,$c); # clear carry
197 &mov($r,&wparam(0)); #
198 &mov($a,&wparam(1)); #
199 &mov($num,&wparam(2)); #
200 &mov($w,&wparam(3)); #
202 &and($num,0xfffffff8); # num / 8
203 &jz(&label("mw_finish"));
205 &set_label("mw_loop",0);
206 for ($i=0; $i<32; $i+=4)
208 &comment("Round $i");
210 &mov("eax",&DWP($i,$a,"",0)); # *a
212 &add("eax",$c); # L(t)+=c
215 &adc("edx",0); # H(t)+=carry
216 &mov(&DWP($i,$r,"",0),"eax"); # *r= L(t);
218 &mov($c,"edx"); # c= H(t);
225 &jz(&label("mw_finish"));
226 &jmp(&label("mw_loop"));
228 &set_label("mw_finish",0);
229 &mov($num,&wparam(2)); # get num
231 &jnz(&label("mw_finish2"));
232 &jmp(&label("mw_end"));
234 &set_label("mw_finish2",1);
235 for ($i=0; $i<7; $i++)
237 &comment("Tail Round $i");
238 &mov("eax",&DWP($i*4,$a,"",0));# *a
240 &add("eax",$c); # L(t)+=c
242 &adc("edx",0); # H(t)+=carry
243 &mov(&DWP($i*4,$r,"",0),"eax");# *r= L(t);
244 &mov($c,"edx"); # c= H(t);
245 &dec($num) if ($i != 7-1);
246 &jz(&label("mw_end")) if ($i != 7-1);
248 &set_label("mw_end",0);
251 &function_end($name);
258 &function_begin($name,"");
265 &mov($r,&wparam(0)); #
266 &mov($a,&wparam(1)); #
267 &mov($num,&wparam(2)); #
269 &and($num,0xfffffff8); # num / 8
270 &jz(&label("sw_finish"));
272 &set_label("sw_loop",0);
273 for ($i=0; $i<32; $i+=4)
275 &comment("Round $i");
276 &mov("eax",&DWP($i,$a,"",0)); # *a
278 &mul("eax"); # *a * *a
279 &mov(&DWP($i*2,$r,"",0),"eax"); #
280 &mov(&DWP($i*2+4,$r,"",0),"edx");#
287 &jnz(&label("sw_loop"));
289 &set_label("sw_finish",0);
290 &mov($num,&wparam(2)); # get num
292 &jz(&label("sw_end"));
294 for ($i=0; $i<7; $i++)
296 &comment("Tail Round $i");
297 &mov("eax",&DWP($i*4,$a,"",0)); # *a
299 &mul("eax"); # *a * *a
300 &mov(&DWP($i*8,$r,"",0),"eax"); #
301 &dec($num) if ($i != 7-1);
302 &mov(&DWP($i*8+4,$r,"",0),"edx");
303 &jz(&label("sw_end")) if ($i != 7-1);
305 &set_label("sw_end",0);
307 &function_end($name);
314 &function_begin($name,"");
315 &mov("edx",&wparam(0)); #
316 &mov("eax",&wparam(1)); #
317 &mov("ebx",&wparam(2)); #
319 &function_end($name);
326 &function_begin($name,"");
337 &mov($r,&wparam(0)); # get r
338 &mov($a,&wparam(1)); # get a
339 &mov($b,&wparam(2)); # get b
340 &mov($num,&wparam(3)); # get num
341 &xor($c,$c); # clear carry
342 &and($num,0xfffffff8); # num / 8
344 &jz(&label("aw_finish"));
346 &set_label("aw_loop",0);
347 for ($i=0; $i<8; $i++)
349 &comment("Round $i");
351 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
352 &mov($tmp2,&DWP($i*4,$b,"",0)); # *b
358 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
366 &jnz(&label("aw_loop"));
368 &set_label("aw_finish",0);
369 &mov($num,&wparam(3)); # get num
371 &jz(&label("aw_end"));
373 for ($i=0; $i<7; $i++)
375 &comment("Tail Round $i");
376 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
377 &mov($tmp2,&DWP($i*4,$b,"",0));# *b
383 &dec($num) if ($i != 6);
384 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
385 &jz(&label("aw_end")) if ($i != 6);
387 &set_label("aw_end",0);
389 # &mov("eax",$c); # $c is "eax"
391 &function_end($name);
398 &function_begin($name,"");
409 &mov($r,&wparam(0)); # get r
410 &mov($a,&wparam(1)); # get a
411 &mov($b,&wparam(2)); # get b
412 &mov($num,&wparam(3)); # get num
413 &xor($c,$c); # clear carry
414 &and($num,0xfffffff8); # num / 8
416 &jz(&label("aw_finish"));
418 &set_label("aw_loop",0);
419 for ($i=0; $i<8; $i++)
421 &comment("Round $i");
423 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
424 &mov($tmp2,&DWP($i*4,$b,"",0)); # *b
430 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
438 &jnz(&label("aw_loop"));
440 &set_label("aw_finish",0);
441 &mov($num,&wparam(3)); # get num
443 &jz(&label("aw_end"));
445 for ($i=0; $i<7; $i++)
447 &comment("Tail Round $i");
448 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
449 &mov($tmp2,&DWP($i*4,$b,"",0));# *b
455 &dec($num) if ($i != 6);
456 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
457 &jz(&label("aw_end")) if ($i != 6);
459 &set_label("aw_end",0);
461 # &mov("eax",$c); # $c is "eax"
463 &function_end($name);
466 sub bn_sub_part_words
470 &function_begin($name,"");
481 &mov($r,&wparam(0)); # get r
482 &mov($a,&wparam(1)); # get a
483 &mov($b,&wparam(2)); # get b
484 &mov($num,&wparam(3)); # get num
485 &xor($c,$c); # clear carry
486 &and($num,0xfffffff8); # num / 8
488 &jz(&label("aw_finish"));
490 &set_label("aw_loop",0);
491 for ($i=0; $i<8; $i++)
493 &comment("Round $i");
495 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
496 &mov($tmp2,&DWP($i*4,$b,"",0)); # *b
502 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
510 &jnz(&label("aw_loop"));
512 &set_label("aw_finish",0);
513 &mov($num,&wparam(3)); # get num
515 &jz(&label("aw_end"));
517 for ($i=0; $i<7; $i++)
519 &comment("Tail Round $i");
520 &mov($tmp1,&DWP(0,$a,"",0)); # *a
521 &mov($tmp2,&DWP(0,$b,"",0));# *b
527 &mov(&DWP(0,$r,"",0),$tmp1); # *r
531 &dec($num) if ($i != 6);
532 &jz(&label("aw_end")) if ($i != 6);
534 &set_label("aw_end",0);
537 &je(&label("pw_end"));
539 &mov($num,&wparam(4)); # get dl
541 &je(&label("pw_end"));
542 &jge(&label("pw_pos"));
548 &and($num,0xfffffff8); # num / 8
549 &jz(&label("pw_neg_finish"));
551 &set_label("pw_neg_loop",0);
552 for ($i=0; $i<8; $i++)
554 &comment("dl<0 Round $i");
557 &mov($tmp2,&DWP($i*4,$b,"",0)); # *b
563 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
570 &jnz(&label("pw_neg_loop"));
572 &set_label("pw_neg_finish",0);
573 &mov($tmp2,&wparam(4)); # get dl
577 &jz(&label("pw_end"));
579 for ($i=0; $i<7; $i++)
581 &comment("dl<0 Tail Round $i");
583 &mov($tmp2,&DWP($i*4,$b,"",0));# *b
589 &dec($num) if ($i != 6);
590 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
591 &jz(&label("pw_end")) if ($i != 6);
594 &jmp(&label("pw_end"));
596 &set_label("pw_pos",0);
598 &and($num,0xfffffff8); # num / 8
599 &jz(&label("pw_pos_finish"));
601 &set_label("pw_pos_loop",0);
603 for ($i=0; $i<8; $i++)
605 &comment("dl>0 Round $i");
607 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
609 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
610 &jnc(&label("pw_nc".$i));
617 &jnz(&label("pw_pos_loop"));
619 &set_label("pw_pos_finish",0);
620 &mov($num,&wparam(4)); # get dl
622 &jz(&label("pw_end"));
624 for ($i=0; $i<7; $i++)
626 &comment("dl>0 Tail Round $i");
627 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
629 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
630 &jnc(&label("pw_tail_nc".$i));
631 &dec($num) if ($i != 6);
632 &jz(&label("pw_end")) if ($i != 6);
635 &jmp(&label("pw_end"));
637 &set_label("pw_nc_loop",0);
638 for ($i=0; $i<8; $i++)
640 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
641 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
642 &set_label("pw_nc".$i,0);
649 &jnz(&label("pw_nc_loop"));
651 &mov($num,&wparam(4)); # get dl
653 &jz(&label("pw_nc_end"));
655 for ($i=0; $i<7; $i++)
657 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
658 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
659 &set_label("pw_tail_nc".$i,0);
660 &dec($num) if ($i != 6);
661 &jz(&label("pw_nc_end")) if ($i != 6);
664 &set_label("pw_nc_end",0);
667 &set_label("pw_end",0);
669 # &mov("eax",$c); # $c is "eax"
671 &function_end($name);