3 push(@INC,"perlasm","../../perlasm");
6 &asm_init($ARGV[0],$0);
8 $sse2=0; # SSE2 is temporarily disabled...
10 &bn_mul_add_words("bn_mul_add_words");
11 &bn_mul_words("bn_mul_words");
12 &bn_sqr_words("bn_sqr_words");
13 &bn_div_words("bn_div_words");
14 &bn_add_words("bn_add_words");
15 &bn_sub_words("bn_sub_words");
16 &bn_sub_part_words("bn_sub_part_words");
24 &function_begin($name,"");
34 &xor($c,$c); # clear carry
35 &mov($r,&wparam(0)); #
37 &mov("ecx",&wparam(2)); #
38 &mov($a,&wparam(1)); #
40 &and("ecx",0xfffffff8); # num / 8
41 &mov($w,&wparam(3)); #
43 &push("ecx"); # Up the stack for a tmp variable
45 &jz(&label("maw_finish"));
48 &picmeup("eax","OPENSSL_ia32cap");
49 &bt(&DWP(0,"eax"),26);
50 &jnc(&label("maw_loop"));
52 &movd("mm0",$w); # mm0 = w
53 &pxor("mm1","mm1"); # mm1 = carry_in
55 &set_label("maw_sse2_loop",0);
56 &movd("mm3",&DWP(0,$r,"",0)); # mm3 = r[0]
57 &paddq("mm1","mm3"); # mm1 = carry_in + r[0]
58 &movd("mm2",&DWP(0,$a,"",0)); # mm2 = a[0]
59 &pmuludq("mm2","mm0"); # mm2 = w*a[0]
60 &movd("mm4",&DWP(4,$a,"",0)); # mm4 = a[1]
61 &pmuludq("mm4","mm0"); # mm4 = w*a[1]
62 &movd("mm6",&DWP(8,$a,"",0)); # mm6 = a[2]
63 &pmuludq("mm6","mm0"); # mm6 = w*a[2]
64 &movd("mm7",&DWP(12,$a,"",0)); # mm7 = a[3]
65 &pmuludq("mm7","mm0"); # mm7 = w*a[3]
66 &paddq("mm1","mm2"); # mm1 = carry_in + r[0] + w*a[0]
67 &movd("mm3",&DWP(4,$r,"",0)); # mm3 = r[1]
68 &paddq("mm3","mm4"); # mm3 = r[1] + w*a[1]
69 &movd("mm5",&DWP(8,$r,"",0)); # mm5 = r[2]
70 &paddq("mm5","mm6"); # mm5 = r[2] + w*a[2]
71 &movd("mm4",&DWP(12,$r,"",0)); # mm4 = r[3]
72 &paddq("mm7","mm4"); # mm7 = r[3] + w*a[3]
73 &movd(&DWP(0,$r,"",0),"mm1");
74 &movd("mm2",&DWP(16,$a,"",0)); # mm2 = a[4]
75 &pmuludq("mm2","mm0"); # mm2 = w*a[4]
76 &psrlq("mm1",32); # mm1 = carry0
77 &movd("mm4",&DWP(20,$a,"",0)); # mm4 = a[5]
78 &pmuludq("mm4","mm0"); # mm4 = w*a[5]
79 &paddq("mm1","mm3"); # mm1 = carry0 + r[1] + w*a[1]
80 &movd("mm6",&DWP(24,$a,"",0)); # mm6 = a[6]
81 &pmuludq("mm6","mm0"); # mm6 = w*a[6]
82 &movd(&DWP(4,$r,"",0),"mm1");
83 &psrlq("mm1",32); # mm1 = carry1
84 &movd("mm3",&DWP(28,$a,"",0)); # mm3 = a[7]
86 &pmuludq("mm3","mm0"); # mm3 = w*a[7]
87 &paddq("mm1","mm5"); # mm1 = carry1 + r[2] + w*a[2]
88 &movd("mm5",&DWP(16,$r,"",0)); # mm5 = r[4]
89 &paddq("mm2","mm5"); # mm2 = r[4] + w*a[4]
90 &movd(&DWP(8,$r,"",0),"mm1");
91 &psrlq("mm1",32); # mm1 = carry2
92 &paddq("mm1","mm7"); # mm1 = carry2 + r[3] + w*a[3]
93 &movd("mm5",&DWP(20,$r,"",0)); # mm5 = r[5]
94 &paddq("mm4","mm5"); # mm4 = r[5] + w*a[5]
95 &movd(&DWP(12,$r,"",0),"mm1");
96 &psrlq("mm1",32); # mm1 = carry3
97 &paddq("mm1","mm2"); # mm1 = carry3 + r[4] + w*a[4]
98 &movd("mm5",&DWP(24,$r,"",0)); # mm5 = r[6]
99 &paddq("mm6","mm5"); # mm6 = r[6] + w*a[6]
100 &movd(&DWP(16,$r,"",0),"mm1");
101 &psrlq("mm1",32); # mm1 = carry4
102 &paddq("mm1","mm4"); # mm1 = carry4 + r[5] + w*a[5]
103 &movd("mm5",&DWP(28,$r,"",0)); # mm5 = r[7]
104 &paddq("mm3","mm5"); # mm3 = r[7] + w*a[7]
105 &movd(&DWP(20,$r,"",0),"mm1");
106 &psrlq("mm1",32); # mm1 = carry5
107 &paddq("mm1","mm6"); # mm1 = carry5 + r[6] + w*a[6]
108 &movd(&DWP(24,$r,"",0),"mm1");
109 &psrlq("mm1",32); # mm1 = carry6
110 &paddq("mm1","mm3"); # mm1 = carry6 + r[7] + w*a[7]
111 &movd(&DWP(28,$r,"",0),"mm1");
113 &psrlq("mm1",32); # mm1 = carry_out
116 &jnz(&label("maw_sse2_loop"));
118 &movd($c,"mm1"); # c = carry_out
121 &jmp(&label("maw_finish"));
124 &set_label("maw_loop",0);
126 &mov(&swtmp(0),"ecx"); #
128 for ($i=0; $i<32; $i+=4)
130 &comment("Round $i");
132 &mov("eax",&DWP($i,$a,"",0)); # *a
134 &add("eax",$c); # L(t)+= *r
135 &mov($c,&DWP($i,$r,"",0)); # L(t)+= *r
136 &adc("edx",0); # H(t)+=carry
137 &add("eax",$c); # L(t)+=c
138 &adc("edx",0); # H(t)+=carry
139 &mov(&DWP($i,$r,"",0),"eax"); # *r= L(t);
140 &mov($c,"edx"); # c= H(t);
144 &mov("ecx",&swtmp(0)); #
148 &jnz(&label("maw_loop"));
150 &set_label("maw_finish",0);
151 &mov("ecx",&wparam(2)); # get num
153 &jnz(&label("maw_finish2")); # helps branch prediction
154 &jmp(&label("maw_end"));
156 &set_label("maw_finish2",1);
157 for ($i=0; $i<7; $i++)
159 &comment("Tail Round $i");
160 &mov("eax",&DWP($i*4,$a,"",0));# *a
162 &add("eax",$c); # L(t)+=c
163 &mov($c,&DWP($i*4,$r,"",0)); # L(t)+= *r
164 &adc("edx",0); # H(t)+=carry
166 &adc("edx",0); # H(t)+=carry
167 &dec("ecx") if ($i != 7-1);
168 &mov(&DWP($i*4,$r,"",0),"eax"); # *r= L(t);
169 &mov($c,"edx"); # c= H(t);
170 &jz(&label("maw_end")) if ($i != 7-1);
172 &set_label("maw_end",0);
175 &pop("ecx"); # clear variable from
177 &function_end($name);
184 &function_begin($name,"");
195 &xor($c,$c); # clear carry
196 &mov($r,&wparam(0)); #
197 &mov($a,&wparam(1)); #
198 &mov($num,&wparam(2)); #
199 &mov($w,&wparam(3)); #
201 &and($num,0xfffffff8); # num / 8
202 &jz(&label("mw_finish"));
204 &set_label("mw_loop",0);
205 for ($i=0; $i<32; $i+=4)
207 &comment("Round $i");
209 &mov("eax",&DWP($i,$a,"",0)); # *a
211 &add("eax",$c); # L(t)+=c
214 &adc("edx",0); # H(t)+=carry
215 &mov(&DWP($i,$r,"",0),"eax"); # *r= L(t);
217 &mov($c,"edx"); # c= H(t);
224 &jz(&label("mw_finish"));
225 &jmp(&label("mw_loop"));
227 &set_label("mw_finish",0);
228 &mov($num,&wparam(2)); # get num
230 &jnz(&label("mw_finish2"));
231 &jmp(&label("mw_end"));
233 &set_label("mw_finish2",1);
234 for ($i=0; $i<7; $i++)
236 &comment("Tail Round $i");
237 &mov("eax",&DWP($i*4,$a,"",0));# *a
239 &add("eax",$c); # L(t)+=c
241 &adc("edx",0); # H(t)+=carry
242 &mov(&DWP($i*4,$r,"",0),"eax");# *r= L(t);
243 &mov($c,"edx"); # c= H(t);
244 &dec($num) if ($i != 7-1);
245 &jz(&label("mw_end")) if ($i != 7-1);
247 &set_label("mw_end",0);
250 &function_end($name);
257 &function_begin($name,"");
264 &mov($r,&wparam(0)); #
265 &mov($a,&wparam(1)); #
266 &mov($num,&wparam(2)); #
268 &and($num,0xfffffff8); # num / 8
269 &jz(&label("sw_finish"));
271 &set_label("sw_loop",0);
272 for ($i=0; $i<32; $i+=4)
274 &comment("Round $i");
275 &mov("eax",&DWP($i,$a,"",0)); # *a
277 &mul("eax"); # *a * *a
278 &mov(&DWP($i*2,$r,"",0),"eax"); #
279 &mov(&DWP($i*2+4,$r,"",0),"edx");#
286 &jnz(&label("sw_loop"));
288 &set_label("sw_finish",0);
289 &mov($num,&wparam(2)); # get num
291 &jz(&label("sw_end"));
293 for ($i=0; $i<7; $i++)
295 &comment("Tail Round $i");
296 &mov("eax",&DWP($i*4,$a,"",0)); # *a
298 &mul("eax"); # *a * *a
299 &mov(&DWP($i*8,$r,"",0),"eax"); #
300 &dec($num) if ($i != 7-1);
301 &mov(&DWP($i*8+4,$r,"",0),"edx");
302 &jz(&label("sw_end")) if ($i != 7-1);
304 &set_label("sw_end",0);
306 &function_end($name);
313 &function_begin($name,"");
314 &mov("edx",&wparam(0)); #
315 &mov("eax",&wparam(1)); #
316 &mov("ebx",&wparam(2)); #
318 &function_end($name);
325 &function_begin($name,"");
336 &mov($r,&wparam(0)); # get r
337 &mov($a,&wparam(1)); # get a
338 &mov($b,&wparam(2)); # get b
339 &mov($num,&wparam(3)); # get num
340 &xor($c,$c); # clear carry
341 &and($num,0xfffffff8); # num / 8
343 &jz(&label("aw_finish"));
345 &set_label("aw_loop",0);
346 for ($i=0; $i<8; $i++)
348 &comment("Round $i");
350 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
351 &mov($tmp2,&DWP($i*4,$b,"",0)); # *b
357 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
365 &jnz(&label("aw_loop"));
367 &set_label("aw_finish",0);
368 &mov($num,&wparam(3)); # get num
370 &jz(&label("aw_end"));
372 for ($i=0; $i<7; $i++)
374 &comment("Tail Round $i");
375 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
376 &mov($tmp2,&DWP($i*4,$b,"",0));# *b
382 &dec($num) if ($i != 6);
383 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
384 &jz(&label("aw_end")) if ($i != 6);
386 &set_label("aw_end",0);
388 # &mov("eax",$c); # $c is "eax"
390 &function_end($name);
397 &function_begin($name,"");
408 &mov($r,&wparam(0)); # get r
409 &mov($a,&wparam(1)); # get a
410 &mov($b,&wparam(2)); # get b
411 &mov($num,&wparam(3)); # get num
412 &xor($c,$c); # clear carry
413 &and($num,0xfffffff8); # num / 8
415 &jz(&label("aw_finish"));
417 &set_label("aw_loop",0);
418 for ($i=0; $i<8; $i++)
420 &comment("Round $i");
422 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
423 &mov($tmp2,&DWP($i*4,$b,"",0)); # *b
429 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
437 &jnz(&label("aw_loop"));
439 &set_label("aw_finish",0);
440 &mov($num,&wparam(3)); # get num
442 &jz(&label("aw_end"));
444 for ($i=0; $i<7; $i++)
446 &comment("Tail Round $i");
447 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
448 &mov($tmp2,&DWP($i*4,$b,"",0));# *b
454 &dec($num) if ($i != 6);
455 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
456 &jz(&label("aw_end")) if ($i != 6);
458 &set_label("aw_end",0);
460 # &mov("eax",$c); # $c is "eax"
462 &function_end($name);
465 sub bn_sub_part_words
469 &function_begin($name,"");
480 &mov($r,&wparam(0)); # get r
481 &mov($a,&wparam(1)); # get a
482 &mov($b,&wparam(2)); # get b
483 &mov($num,&wparam(3)); # get num
484 &xor($c,$c); # clear carry
485 &and($num,0xfffffff8); # num / 8
487 &jz(&label("aw_finish"));
489 &set_label("aw_loop",0);
490 for ($i=0; $i<8; $i++)
492 &comment("Round $i");
494 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
495 &mov($tmp2,&DWP($i*4,$b,"",0)); # *b
501 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
509 &jnz(&label("aw_loop"));
511 &set_label("aw_finish",0);
512 &mov($num,&wparam(3)); # get num
514 &jz(&label("aw_end"));
516 for ($i=0; $i<7; $i++)
518 &comment("Tail Round $i");
519 &mov($tmp1,&DWP(0,$a,"",0)); # *a
520 &mov($tmp2,&DWP(0,$b,"",0));# *b
526 &mov(&DWP(0,$r,"",0),$tmp1); # *r
530 &dec($num) if ($i != 6);
531 &jz(&label("aw_end")) if ($i != 6);
533 &set_label("aw_end",0);
536 &je(&label("pw_end"));
538 &mov($num,&wparam(4)); # get dl
540 &je(&label("pw_end"));
541 &jge(&label("pw_pos"));
547 &and($num,0xfffffff8); # num / 8
548 &jz(&label("pw_neg_finish"));
550 &set_label("pw_neg_loop",0);
551 for ($i=0; $i<8; $i++)
553 &comment("dl<0 Round $i");
556 &mov($tmp2,&DWP($i*4,$b,"",0)); # *b
562 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
569 &jnz(&label("pw_neg_loop"));
571 &set_label("pw_neg_finish",0);
572 &mov($tmp2,&wparam(4)); # get dl
576 &jz(&label("pw_end"));
578 for ($i=0; $i<7; $i++)
580 &comment("dl<0 Tail Round $i");
582 &mov($tmp2,&DWP($i*4,$b,"",0));# *b
588 &dec($num) if ($i != 6);
589 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
590 &jz(&label("pw_end")) if ($i != 6);
593 &jmp(&label("pw_end"));
595 &set_label("pw_pos",0);
597 &and($num,0xfffffff8); # num / 8
598 &jz(&label("pw_pos_finish"));
600 &set_label("pw_pos_loop",0);
602 for ($i=0; $i<8; $i++)
604 &comment("dl>0 Round $i");
606 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
608 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
609 &jnc(&label("pw_nc".$i));
616 &jnz(&label("pw_pos_loop"));
618 &set_label("pw_pos_finish",0);
619 &mov($num,&wparam(4)); # get dl
621 &jz(&label("pw_end"));
623 for ($i=0; $i<7; $i++)
625 &comment("dl>0 Tail Round $i");
626 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
628 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
629 &jnc(&label("pw_tail_nc".$i));
630 &dec($num) if ($i != 6);
631 &jz(&label("pw_end")) if ($i != 6);
634 &jmp(&label("pw_end"));
636 &set_label("pw_nc_loop",0);
637 for ($i=0; $i<8; $i++)
639 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
640 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
641 &set_label("pw_nc".$i,0);
648 &jnz(&label("pw_nc_loop"));
650 &mov($num,&wparam(4)); # get dl
652 &jz(&label("pw_nc_end"));
654 for ($i=0; $i<7; $i++)
656 &mov($tmp1,&DWP($i*4,$a,"",0)); # *a
657 &mov(&DWP($i*4,$r,"",0),$tmp1); # *r
658 &set_label("pw_tail_nc".$i,0);
659 &dec($num) if ($i != 6);
660 &jz(&label("pw_nc_end")) if ($i != 6);
663 &set_label("pw_nc_end",0);
666 &set_label("pw_end",0);
668 # &mov("eax",$c); # $c is "eax"
670 &function_end($name);