VirtualBox

source: vbox/trunk/src/libs/openssl-3.0.1/crypto/sha/asm/sha256-586.pl@ 94082

Last change on this file since 94082 was 94082, checked in by vboxsync, 3 years ago

libs/openssl-3.0.1: started applying and adjusting our OpenSSL changes to 3.0.1. bugref:10128
File size: 36.0 KB
Line 
1#! /usr/bin/env perl
2# Copyright 2007-2020 The OpenSSL Project Authors. All Rights Reserved.
3#
4# Licensed under the Apache License 2.0 (the "License"). You may not use
5# this file except in compliance with the License. You can obtain a copy
6# in the file LICENSE in the source distribution or at
7# https://www.openssl.org/source/license.html
8
9#
10# ====================================================================
11# Written by Andy Polyakov <[email protected]> for the OpenSSL
12# project. The module is, however, dual licensed under OpenSSL and
13# CRYPTOGAMS licenses depending on where you obtain it. For further
14# details see http://www.openssl.org/~appro/cryptogams/.
15# ====================================================================
16#
17# SHA256 block transform for x86. September 2007.
18#
19# Performance improvement over compiler generated code varies from
20# 10% to 40% [see below]. Not very impressive on some µ-archs, but
21# it's 5 times smaller and optimizes amount of writes.
22#
23# May 2012.
24#
25# Optimization including two of Pavel Semjanov's ideas, alternative
26# Maj and full unroll, resulted in ~20-25% improvement on most CPUs,
27# ~7% on Pentium, ~40% on Atom. As fully unrolled loop body is almost
28# 15x larger, 8KB vs. 560B, it's fired only for longer inputs. But not
29# on P4, where it kills performance, nor Sandy Bridge, where folded
30# loop is approximately as fast...
31#
32# June 2012.
33#
34# Add AMD XOP-specific code path, >30% improvement on Bulldozer over
35# May version, >60% over original. Add AVX+shrd code path, >25%
36# improvement on Sandy Bridge over May version, 60% over original.
37#
38# May 2013.
39#
40# Replace AMD XOP code path with SSSE3 to cover more processors.
41# (Biggest improvement coefficient is on upcoming Atom Silvermont,
42# not shown.) Add AVX+BMI code path.
43#
44# March 2014.
45#
46# Add support for Intel SHA Extensions.
47#
48# Performance in clock cycles per processed byte (less is better):
49#
50# gcc icc x86 asm(*) SIMD x86_64 asm(**)
51# Pentium 46 57 40/38 - -
52# PIII 36 33 27/24 - -
53# P4 41 38 28 - 17.3
54# AMD K8 27 25 19/15.5 - 14.9
55# Core2 26 23 18/15.6 14.3 13.8
56# Westmere 27 - 19/15.7 13.4 12.3
57# Sandy Bridge 25 - 15.9 12.4 11.6
58# Ivy Bridge 24 - 15.0 11.4 10.3
59# Haswell 22 - 13.9 9.46 7.80
60# Skylake 20 - 14.9 9.50 7.70
61# Bulldozer 36 - 27/22 17.0 13.6
62# VIA Nano 36 - 25/22 16.8 16.5
63# Atom 50 - 30/25 21.9 18.9
64# Silvermont 40 - 34/31 22.9 20.6
65# Goldmont 29 - 20 16.3(***)
66#
67# (*) numbers after slash are for unrolled loop, where applicable;
68# (**) x86_64 assembly performance is presented for reference
69# purposes, results are best-available;
70# (***) SHAEXT result is 4.1, strangely enough better than 64-bit one;
71
72$0 =~ m/(.*[\/\\])[^\/\\]+$/; $dir=$1;
73push(@INC,"${dir}","${dir}../../perlasm");
74require "x86asm.pl";
75
76$output=pop and open STDOUT,">$output";
77
78&asm_init($ARGV[0],$ARGV[$#ARGV] eq "386");
79
80$xmm=$avx=0;
81for (@ARGV) { $xmm=1 if (/-DOPENSSL_IA32_SSE2/); }
82
83if ($xmm && `$ENV{CC} -Wa,-v -c -o /dev/null -x assembler /dev/null 2>&1`
84 =~ /GNU assembler version ([2-9]\.[0-9]+)/) {
85 $avx = ($1>=2.19) + ($1>=2.22);
86}
87
88if ($xmm && !$avx && $ARGV[0] eq "win32n" &&
89 `nasm -v 2>&1` =~ /NASM version ([2-9]\.[0-9]+)/) {
90 $avx = ($1>=2.03) + ($1>=2.10);
91}
92
93if ($xmm && !$avx && $ARGV[0] eq "win32" &&
94 `ml 2>&1` =~ /Version ([0-9]+)\./) {
95 $avx = ($1>=10) + ($1>=11);
96}
97
98if ($xmm && !$avx && `$ENV{CC} -v 2>&1` =~ /((?:clang|LLVM) version|based on LLVM) ([0-9]+\.[0-9]+)/) {
99 $avx = ($2>=3.0) + ($2>3.0);
100}
101
102$shaext=$xmm; ### set to zero if compiling for 1.0.1
103
104$unroll_after = 64*4; # If pre-evicted from L1P cache first spin of
105 # fully unrolled loop was measured to run about
106 # 3-4x slower. If slowdown coefficient is N and
107 # unrolled loop is m times faster, then you break
108 # even at (N-1)/(m-1) blocks. Then it needs to be
109 # adjusted for probability of code being evicted,
110 # code size/cache size=1/4. Typical m is 1.15...
111
112$A="eax";
113$E="edx";
114$T="ebx";
115$Aoff=&DWP(4,"esp");
116$Boff=&DWP(8,"esp");
117$Coff=&DWP(12,"esp");
118$Doff=&DWP(16,"esp");
119$Eoff=&DWP(20,"esp");
120$Foff=&DWP(24,"esp");
121$Goff=&DWP(28,"esp");
122$Hoff=&DWP(32,"esp");
123$Xoff=&DWP(36,"esp");
124$K256="ebp";
125
126sub BODY_16_63() {
127 &mov ($T,"ecx"); # "ecx" is preloaded
128 &mov ("esi",&DWP(4*(9+15+16-14),"esp"));
129 &ror ("ecx",18-7);
130 &mov ("edi","esi");
131 &ror ("esi",19-17);
132 &xor ("ecx",$T);
133 &shr ($T,3);
134 &ror ("ecx",7);
135 &xor ("esi","edi");
136 &xor ($T,"ecx"); # T = sigma0(X[-15])
137 &ror ("esi",17);
138 &add ($T,&DWP(4*(9+15+16),"esp")); # T += X[-16]
139 &shr ("edi",10);
140 &add ($T,&DWP(4*(9+15+16-9),"esp")); # T += X[-7]
141 #&xor ("edi","esi") # sigma1(X[-2])
142 # &add ($T,"edi"); # T += sigma1(X[-2])
143 # &mov (&DWP(4*(9+15),"esp"),$T); # save X[0]
144
145 &BODY_00_15(1);
146}
147sub BODY_00_15() {
148 my $in_16_63=shift;
149
150 &mov ("ecx",$E);
151 &xor ("edi","esi") if ($in_16_63); # sigma1(X[-2])
152 &mov ("esi",$Foff);
153 &ror ("ecx",25-11);
154 &add ($T,"edi") if ($in_16_63); # T += sigma1(X[-2])
155 &mov ("edi",$Goff);
156 &xor ("ecx",$E);
157 &xor ("esi","edi");
158 &mov ($T,&DWP(4*(9+15),"esp")) if (!$in_16_63);
159 &mov (&DWP(4*(9+15),"esp"),$T) if ($in_16_63); # save X[0]
160 &ror ("ecx",11-6);
161 &and ("esi",$E);
162 &mov ($Eoff,$E); # modulo-scheduled
163 &xor ($E,"ecx");
164 &add ($T,$Hoff); # T += h
165 &xor ("esi","edi"); # Ch(e,f,g)
166 &ror ($E,6); # Sigma1(e)
167 &mov ("ecx",$A);
168 &add ($T,"esi"); # T += Ch(e,f,g)
169
170 &ror ("ecx",22-13);
171 &add ($T,$E); # T += Sigma1(e)
172 &mov ("edi",$Boff);
173 &xor ("ecx",$A);
174 &mov ($Aoff,$A); # modulo-scheduled
175 &lea ("esp",&DWP(-4,"esp"));
176 &ror ("ecx",13-2);
177 &mov ("esi",&DWP(0,$K256));
178 &xor ("ecx",$A);
179 &mov ($E,$Eoff); # e in next iteration, d in this one
180 &xor ($A,"edi"); # a ^= b
181 &ror ("ecx",2); # Sigma0(a)
182
183 &add ($T,"esi"); # T+= K[i]
184 &mov (&DWP(0,"esp"),$A); # (b^c) in next round
185 &add ($E,$T); # d += T
186 &and ($A,&DWP(4,"esp")); # a &= (b^c)
187 &add ($T,"ecx"); # T += Sigma0(a)
188 &xor ($A,"edi"); # h = Maj(a,b,c) = Ch(a^b,c,b)
189 &mov ("ecx",&DWP(4*(9+15+16-1),"esp")) if ($in_16_63); # preload T
190 &add ($K256,4);
191 &add ($A,$T); # h += T
192}
193
194&external_label("OPENSSL_ia32cap_P") if (!$i386);
195
196&function_begin("sha256_block_data_order");
197 &mov ("esi",wparam(0)); # ctx
198 &mov ("edi",wparam(1)); # inp
199 &mov ("eax",wparam(2)); # num
200 &mov ("ebx","esp"); # saved sp
201
202 &call (&label("pic_point")); # make it PIC!
203&set_label("pic_point");
204 &blindpop($K256);
205 &lea ($K256,&DWP(&label("K256")."-".&label("pic_point"),$K256));
206
207 &sub ("esp",16);
208 &and ("esp",-64);
209
210 &shl ("eax",6);
211 &add ("eax","edi");
212 &mov (&DWP(0,"esp"),"esi"); # ctx
213 &mov (&DWP(4,"esp"),"edi"); # inp
214 &mov (&DWP(8,"esp"),"eax"); # inp+num*128
215 &mov (&DWP(12,"esp"),"ebx"); # saved sp
216 if (!$i386 && $xmm) {
217 &picmeup("edx","OPENSSL_ia32cap_P",$K256,&label("K256"));
218 &mov ("ecx",&DWP(0,"edx"));
219 &mov ("ebx",&DWP(4,"edx"));
220 &test ("ecx",1<<20); # check for P4
221 &jnz (&label("loop"));
222 &mov ("edx",&DWP(8,"edx")) if ($xmm);
223 &test ("ecx",1<<24); # check for FXSR
224 &jz ($unroll_after?&label("no_xmm"):&label("loop"));
225 &and ("ecx",1<<30); # mask "Intel CPU" bit
226 &and ("ebx",1<<28|1<<9); # mask AVX and SSSE3 bits
227 &test ("edx",1<<29) if ($shaext); # check for SHA
228 &jnz (&label("shaext")) if ($shaext);
229 &or ("ecx","ebx");
230 &and ("ecx",1<<28|1<<30);
231 &cmp ("ecx",1<<28|1<<30);
232 if ($xmm) {
233 &je (&label("AVX")) if ($avx);
234 &test ("ebx",1<<9); # check for SSSE3
235 &jnz (&label("SSSE3"));
236 } else {
237 &je (&label("loop_shrd"));
238 }
239 if ($unroll_after) {
240&set_label("no_xmm");
241 &sub ("eax","edi");
242 &cmp ("eax",$unroll_after);
243 &jae (&label("unrolled"));
244 } }
245 &jmp (&label("loop"));
246
247sub COMPACT_LOOP() {
248my $suffix=shift;
249
250&set_label("loop$suffix",$suffix?32:16);
251 # copy input block to stack reversing byte and dword order
252 for($i=0;$i<4;$i++) {
253 &mov ("eax",&DWP($i*16+0,"edi"));
254 &mov ("ebx",&DWP($i*16+4,"edi"));
255 &mov ("ecx",&DWP($i*16+8,"edi"));
256 &bswap ("eax");
257 &mov ("edx",&DWP($i*16+12,"edi"));
258 &bswap ("ebx");
259 &push ("eax");
260 &bswap ("ecx");
261 &push ("ebx");
262 &bswap ("edx");
263 &push ("ecx");
264 &push ("edx");
265 }
266 &add ("edi",64);
267 &lea ("esp",&DWP(-4*9,"esp"));# place for A,B,C,D,E,F,G,H
268 &mov (&DWP(4*(9+16)+4,"esp"),"edi");
269
270 # copy ctx->h[0-7] to A,B,C,D,E,F,G,H on stack
271 &mov ($A,&DWP(0,"esi"));
272 &mov ("ebx",&DWP(4,"esi"));
273 &mov ("ecx",&DWP(8,"esi"));
274 &mov ("edi",&DWP(12,"esi"));
275 # &mov ($Aoff,$A);
276 &mov ($Boff,"ebx");
277 &xor ("ebx","ecx");
278 &mov ($Coff,"ecx");
279 &mov ($Doff,"edi");
280 &mov (&DWP(0,"esp"),"ebx"); # magic
281 &mov ($E,&DWP(16,"esi"));
282 &mov ("ebx",&DWP(20,"esi"));
283 &mov ("ecx",&DWP(24,"esi"));
284 &mov ("edi",&DWP(28,"esi"));
285 # &mov ($Eoff,$E);
286 &mov ($Foff,"ebx");
287 &mov ($Goff,"ecx");
288 &mov ($Hoff,"edi");
289
290&set_label("00_15$suffix",16);
291
292 &BODY_00_15();
293
294 &cmp ("esi",0xc19bf174);
295 &jne (&label("00_15$suffix"));
296
297 &mov ("ecx",&DWP(4*(9+15+16-1),"esp")); # preloaded in BODY_00_15(1)
298 &jmp (&label("16_63$suffix"));
299
300&set_label("16_63$suffix",16);
301
302 &BODY_16_63();
303
304 &cmp ("esi",0xc67178f2);
305 &jne (&label("16_63$suffix"));
306
307 &mov ("esi",&DWP(4*(9+16+64)+0,"esp"));#ctx
308 # &mov ($A,$Aoff);
309 &mov ("ebx",$Boff);
310 # &mov ("edi",$Coff);
311 &mov ("ecx",$Doff);
312 &add ($A,&DWP(0,"esi"));
313 &add ("ebx",&DWP(4,"esi"));
314 &add ("edi",&DWP(8,"esi"));
315 &add ("ecx",&DWP(12,"esi"));
316 &mov (&DWP(0,"esi"),$A);
317 &mov (&DWP(4,"esi"),"ebx");
318 &mov (&DWP(8,"esi"),"edi");
319 &mov (&DWP(12,"esi"),"ecx");
320 # &mov ($E,$Eoff);
321 &mov ("eax",$Foff);
322 &mov ("ebx",$Goff);
323 &mov ("ecx",$Hoff);
324 &mov ("edi",&DWP(4*(9+16+64)+4,"esp"));#inp
325 &add ($E,&DWP(16,"esi"));
326 &add ("eax",&DWP(20,"esi"));
327 &add ("ebx",&DWP(24,"esi"));
328 &add ("ecx",&DWP(28,"esi"));
329 &mov (&DWP(16,"esi"),$E);
330 &mov (&DWP(20,"esi"),"eax");
331 &mov (&DWP(24,"esi"),"ebx");
332 &mov (&DWP(28,"esi"),"ecx");
333
334 &lea ("esp",&DWP(4*(9+16+64),"esp"));# destroy frame
335 &sub ($K256,4*64); # rewind K
336
337 &cmp ("edi",&DWP(8,"esp")); # are we done yet?
338 &jb (&label("loop$suffix"));
339}
340 &COMPACT_LOOP();
341 &mov ("esp",&DWP(12,"esp")); # restore sp
342&function_end_A();
343 if (!$i386 && !$xmm) {
344 # ~20% improvement on Sandy Bridge
345 local *ror = sub { &shrd(@_[0],@_) };
346 &COMPACT_LOOP("_shrd");
347 &mov ("esp",&DWP(12,"esp")); # restore sp
348&function_end_A();
349 }
350
351&set_label("K256",64); # Yes! I keep it in the code segment!
352@K256=( 0x428a2f98,0x71374491,0xb5c0fbcf,0xe9b5dba5,
353 0x3956c25b,0x59f111f1,0x923f82a4,0xab1c5ed5,
354 0xd807aa98,0x12835b01,0x243185be,0x550c7dc3,
355 0x72be5d74,0x80deb1fe,0x9bdc06a7,0xc19bf174,
356 0xe49b69c1,0xefbe4786,0x0fc19dc6,0x240ca1cc,
357 0x2de92c6f,0x4a7484aa,0x5cb0a9dc,0x76f988da,
358 0x983e5152,0xa831c66d,0xb00327c8,0xbf597fc7,
359 0xc6e00bf3,0xd5a79147,0x06ca6351,0x14292967,
360 0x27b70a85,0x2e1b2138,0x4d2c6dfc,0x53380d13,
361 0x650a7354,0x766a0abb,0x81c2c92e,0x92722c85,
362 0xa2bfe8a1,0xa81a664b,0xc24b8b70,0xc76c51a3,
363 0xd192e819,0xd6990624,0xf40e3585,0x106aa070,
364 0x19a4c116,0x1e376c08,0x2748774c,0x34b0bcb5,
365 0x391c0cb3,0x4ed8aa4a,0x5b9cca4f,0x682e6ff3,
366 0x748f82ee,0x78a5636f,0x84c87814,0x8cc70208,
367 0x90befffa,0xa4506ceb,0xbef9a3f7,0xc67178f2 );
368&data_word(@K256);
369&data_word(0x00010203,0x04050607,0x08090a0b,0x0c0d0e0f); # byte swap mask
370&asciz("SHA256 block transform for x86, CRYPTOGAMS by <appro\@openssl.org>");
371
372($a,$b,$c,$d,$e,$f,$g,$h)=(0..7); # offsets
373sub off { &DWP(4*(((shift)-$i)&7),"esp"); }
374
375if (!$i386 && $unroll_after) {
376my @AH=($A,$K256);
377
378&set_label("unrolled",16);
379 &lea ("esp",&DWP(-96,"esp"));
380 # copy ctx->h[0-7] to A,B,C,D,E,F,G,H on stack
381 &mov ($AH[0],&DWP(0,"esi"));
382 &mov ($AH[1],&DWP(4,"esi"));
383 &mov ("ecx",&DWP(8,"esi"));
384 &mov ("ebx",&DWP(12,"esi"));
385 #&mov (&DWP(0,"esp"),$AH[0]);
386 &mov (&DWP(4,"esp"),$AH[1]);
387 &xor ($AH[1],"ecx"); # magic
388 &mov (&DWP(8,"esp"),"ecx");
389 &mov (&DWP(12,"esp"),"ebx");
390 &mov ($E,&DWP(16,"esi"));
391 &mov ("ebx",&DWP(20,"esi"));
392 &mov ("ecx",&DWP(24,"esi"));
393 &mov ("esi",&DWP(28,"esi"));
394 #&mov (&DWP(16,"esp"),$E);
395 &mov (&DWP(20,"esp"),"ebx");
396 &mov (&DWP(24,"esp"),"ecx");
397 &mov (&DWP(28,"esp"),"esi");
398 &jmp (&label("grand_loop"));
399
400&set_label("grand_loop",16);
401 # copy input block to stack reversing byte order
402 for($i=0;$i<5;$i++) {
403 &mov ("ebx",&DWP(12*$i+0,"edi"));
404 &mov ("ecx",&DWP(12*$i+4,"edi"));
405 &bswap ("ebx");
406 &mov ("esi",&DWP(12*$i+8,"edi"));
407 &bswap ("ecx");
408 &mov (&DWP(32+12*$i+0,"esp"),"ebx");
409 &bswap ("esi");
410 &mov (&DWP(32+12*$i+4,"esp"),"ecx");
411 &mov (&DWP(32+12*$i+8,"esp"),"esi");
412 }
413 &mov ("ebx",&DWP($i*12,"edi"));
414 &add ("edi",64);
415 &bswap ("ebx");
416 &mov (&DWP(96+4,"esp"),"edi");
417 &mov (&DWP(32+12*$i,"esp"),"ebx");
418
419 my ($t1,$t2) = ("ecx","esi");
420
421 for ($i=0;$i<64;$i++) {
422
423 if ($i>=16) {
424 &mov ($T,$t1); # $t1 is preloaded
425 # &mov ($t2,&DWP(32+4*(($i+14)&15),"esp"));
426 &ror ($t1,18-7);
427 &mov ("edi",$t2);
428 &ror ($t2,19-17);
429 &xor ($t1,$T);
430 &shr ($T,3);
431 &ror ($t1,7);
432 &xor ($t2,"edi");
433 &xor ($T,$t1); # T = sigma0(X[-15])
434 &ror ($t2,17);
435 &add ($T,&DWP(32+4*($i&15),"esp")); # T += X[-16]
436 &shr ("edi",10);
437 &add ($T,&DWP(32+4*(($i+9)&15),"esp")); # T += X[-7]
438 #&xor ("edi",$t2) # sigma1(X[-2])
439 # &add ($T,"edi"); # T += sigma1(X[-2])
440 # &mov (&DWP(4*(9+15),"esp"),$T); # save X[0]
441 }
442 &mov ($t1,$E);
443 &xor ("edi",$t2) if ($i>=16); # sigma1(X[-2])
444 &mov ($t2,&off($f));
445 &ror ($E,25-11);
446 &add ($T,"edi") if ($i>=16); # T += sigma1(X[-2])
447 &mov ("edi",&off($g));
448 &xor ($E,$t1);
449 &mov ($T,&DWP(32+4*($i&15),"esp")) if ($i<16); # X[i]
450 &mov (&DWP(32+4*($i&15),"esp"),$T) if ($i>=16 && $i<62); # save X[0]
451 &xor ($t2,"edi");
452 &ror ($E,11-6);
453 &and ($t2,$t1);
454 &mov (&off($e),$t1); # save $E, modulo-scheduled
455 &xor ($E,$t1);
456 &add ($T,&off($h)); # T += h
457 &xor ("edi",$t2); # Ch(e,f,g)
458 &ror ($E,6); # Sigma1(e)
459 &mov ($t1,$AH[0]);
460 &add ($T,"edi"); # T += Ch(e,f,g)
461
462 &ror ($t1,22-13);
463 &mov ($t2,$AH[0]);
464 &mov ("edi",&off($b));
465 &xor ($t1,$AH[0]);
466 &mov (&off($a),$AH[0]); # save $A, modulo-scheduled
467 &xor ($AH[0],"edi"); # a ^= b, (b^c) in next round
468 &ror ($t1,13-2);
469 &and ($AH[1],$AH[0]); # (b^c) &= (a^b)
470 &lea ($E,&DWP(@K256[$i],$T,$E)); # T += Sigma1(1)+K[i]
471 &xor ($t1,$t2);
472 &xor ($AH[1],"edi"); # h = Maj(a,b,c) = Ch(a^b,c,b)
473 &mov ($t2,&DWP(32+4*(($i+2)&15),"esp")) if ($i>=15 && $i<63);
474 &ror ($t1,2); # Sigma0(a)
475
476 &add ($AH[1],$E); # h += T
477 &add ($E,&off($d)); # d += T
478 &add ($AH[1],$t1); # h += Sigma0(a)
479 &mov ($t1,&DWP(32+4*(($i+15)&15),"esp")) if ($i>=15 && $i<63);
480
481 @AH = reverse(@AH); # rotate(a,h)
482 ($t1,$t2) = ($t2,$t1); # rotate(t1,t2)
483 }
484 &mov ("esi",&DWP(96,"esp")); #ctx
485 #&mov ($AH[0],&DWP(0,"esp"));
486 &xor ($AH[1],"edi"); #&mov ($AH[1],&DWP(4,"esp"));
487 #&mov ("edi", &DWP(8,"esp"));
488 &mov ("ecx",&DWP(12,"esp"));
489 &add ($AH[0],&DWP(0,"esi"));
490 &add ($AH[1],&DWP(4,"esi"));
491 &add ("edi",&DWP(8,"esi"));
492 &add ("ecx",&DWP(12,"esi"));
493 &mov (&DWP(0,"esi"),$AH[0]);
494 &mov (&DWP(4,"esi"),$AH[1]);
495 &mov (&DWP(8,"esi"),"edi");
496 &mov (&DWP(12,"esi"),"ecx");
497 #&mov (&DWP(0,"esp"),$AH[0]);
498 &mov (&DWP(4,"esp"),$AH[1]);
499 &xor ($AH[1],"edi"); # magic
500 &mov (&DWP(8,"esp"),"edi");
501 &mov (&DWP(12,"esp"),"ecx");
502 #&mov ($E,&DWP(16,"esp"));
503 &mov ("edi",&DWP(20,"esp"));
504 &mov ("ebx",&DWP(24,"esp"));
505 &mov ("ecx",&DWP(28,"esp"));
506 &add ($E,&DWP(16,"esi"));
507 &add ("edi",&DWP(20,"esi"));
508 &add ("ebx",&DWP(24,"esi"));
509 &add ("ecx",&DWP(28,"esi"));
510 &mov (&DWP(16,"esi"),$E);
511 &mov (&DWP(20,"esi"),"edi");
512 &mov (&DWP(24,"esi"),"ebx");
513 &mov (&DWP(28,"esi"),"ecx");
514 #&mov (&DWP(16,"esp"),$E);
515 &mov (&DWP(20,"esp"),"edi");
516 &mov ("edi",&DWP(96+4,"esp")); # inp
517 &mov (&DWP(24,"esp"),"ebx");
518 &mov (&DWP(28,"esp"),"ecx");
519
520 &cmp ("edi",&DWP(96+8,"esp")); # are we done yet?
521 &jb (&label("grand_loop"));
522
523 &mov ("esp",&DWP(96+12,"esp")); # restore sp
524&function_end_A();
525}
526 if (!$i386 && $xmm) {{{
527if ($shaext) {
528######################################################################
529# Intel SHA Extensions implementation of SHA256 update function.
530#
531my ($ctx,$inp,$end)=("esi","edi","eax");
532my ($Wi,$ABEF,$CDGH,$TMP)=map("xmm$_",(0..2,7));
533my @MSG=map("xmm$_",(3..6));
534
535sub sha256op38 {
536 my ($opcodelet,$dst,$src)=@_;
537 if ("$dst:$src" =~ /xmm([0-7]):xmm([0-7])/)
538 { &data_byte(0x0f,0x38,$opcodelet,0xc0|($1<<3)|$2); }
539}
540sub sha256rnds2 { sha256op38(0xcb,@_); }
541sub sha256msg1 { sha256op38(0xcc,@_); }
542sub sha256msg2 { sha256op38(0xcd,@_); }
543
544&set_label("shaext",32);
545 &sub ("esp",32);
546
547 &movdqu ($ABEF,&QWP(0,$ctx)); # DCBA
548 &lea ($K256,&DWP(0x80,$K256));
549 &movdqu ($CDGH,&QWP(16,$ctx)); # HGFE
550 &movdqa ($TMP,&QWP(0x100-0x80,$K256)); # byte swap mask
551
552 &pshufd ($Wi,$ABEF,0x1b); # ABCD
553 &pshufd ($ABEF,$ABEF,0xb1); # CDAB
554 &pshufd ($CDGH,$CDGH,0x1b); # EFGH
555 &palignr ($ABEF,$CDGH,8); # ABEF
556 &punpcklqdq ($CDGH,$Wi); # CDGH
557 &jmp (&label("loop_shaext"));
558
559&set_label("loop_shaext",16);
560 &movdqu (@MSG[0],&QWP(0,$inp));
561 &movdqu (@MSG[1],&QWP(0x10,$inp));
562 &movdqu (@MSG[2],&QWP(0x20,$inp));
563 &pshufb (@MSG[0],$TMP);
564 &movdqu (@MSG[3],&QWP(0x30,$inp));
565 &movdqa (&QWP(16,"esp"),$CDGH); # offload
566
567 &movdqa ($Wi,&QWP(0*16-0x80,$K256));
568 &paddd ($Wi,@MSG[0]);
569 &pshufb (@MSG[1],$TMP);
570 &sha256rnds2 ($CDGH,$ABEF); # 0-3
571 &pshufd ($Wi,$Wi,0x0e);
572 &nop ();
573 &movdqa (&QWP(0,"esp"),$ABEF); # offload
574 &sha256rnds2 ($ABEF,$CDGH);
575
576 &movdqa ($Wi,&QWP(1*16-0x80,$K256));
577 &paddd ($Wi,@MSG[1]);
578 &pshufb (@MSG[2],$TMP);
579 &sha256rnds2 ($CDGH,$ABEF); # 4-7
580 &pshufd ($Wi,$Wi,0x0e);
581 &lea ($inp,&DWP(0x40,$inp));
582 &sha256msg1 (@MSG[0],@MSG[1]);
583 &sha256rnds2 ($ABEF,$CDGH);
584
585 &movdqa ($Wi,&QWP(2*16-0x80,$K256));
586 &paddd ($Wi,@MSG[2]);
587 &pshufb (@MSG[3],$TMP);
588 &sha256rnds2 ($CDGH,$ABEF); # 8-11
589 &pshufd ($Wi,$Wi,0x0e);
590 &movdqa ($TMP,@MSG[3]);
591 &palignr ($TMP,@MSG[2],4);
592 &nop ();
593 &paddd (@MSG[0],$TMP);
594 &sha256msg1 (@MSG[1],@MSG[2]);
595 &sha256rnds2 ($ABEF,$CDGH);
596
597 &movdqa ($Wi,&QWP(3*16-0x80,$K256));
598 &paddd ($Wi,@MSG[3]);
599 &sha256msg2 (@MSG[0],@MSG[3]);
600 &sha256rnds2 ($CDGH,$ABEF); # 12-15
601 &pshufd ($Wi,$Wi,0x0e);
602 &movdqa ($TMP,@MSG[0]);
603 &palignr ($TMP,@MSG[3],4);
604 &nop ();
605 &paddd (@MSG[1],$TMP);
606 &sha256msg1 (@MSG[2],@MSG[3]);
607 &sha256rnds2 ($ABEF,$CDGH);
608
609for($i=4;$i<16-3;$i++) {
610 &movdqa ($Wi,&QWP($i*16-0x80,$K256));
611 &paddd ($Wi,@MSG[0]);
612 &sha256msg2 (@MSG[1],@MSG[0]);
613 &sha256rnds2 ($CDGH,$ABEF); # 16-19...
614 &pshufd ($Wi,$Wi,0x0e);
615 &movdqa ($TMP,@MSG[1]);
616 &palignr ($TMP,@MSG[0],4);
617 &nop ();
618 &paddd (@MSG[2],$TMP);
619 &sha256msg1 (@MSG[3],@MSG[0]);
620 &sha256rnds2 ($ABEF,$CDGH);
621
622 push(@MSG,shift(@MSG));
623}
624 &movdqa ($Wi,&QWP(13*16-0x80,$K256));
625 &paddd ($Wi,@MSG[0]);
626 &sha256msg2 (@MSG[1],@MSG[0]);
627 &sha256rnds2 ($CDGH,$ABEF); # 52-55
628 &pshufd ($Wi,$Wi,0x0e);
629 &movdqa ($TMP,@MSG[1])
630 &palignr ($TMP,@MSG[0],4);
631 &sha256rnds2 ($ABEF,$CDGH);
632 &paddd (@MSG[2],$TMP);
633
634 &movdqa ($Wi,&QWP(14*16-0x80,$K256));
635 &paddd ($Wi,@MSG[1]);
636 &sha256rnds2 ($CDGH,$ABEF); # 56-59
637 &pshufd ($Wi,$Wi,0x0e);
638 &sha256msg2 (@MSG[2],@MSG[1]);
639 &movdqa ($TMP,&QWP(0x100-0x80,$K256)); # byte swap mask
640 &sha256rnds2 ($ABEF,$CDGH);
641
642 &movdqa ($Wi,&QWP(15*16-0x80,$K256));
643 &paddd ($Wi,@MSG[2]);
644 &nop ();
645 &sha256rnds2 ($CDGH,$ABEF); # 60-63
646 &pshufd ($Wi,$Wi,0x0e);
647 &cmp ($end,$inp);
648 &nop ();
649 &sha256rnds2 ($ABEF,$CDGH);
650
651 &paddd ($CDGH,&QWP(16,"esp"));
652 &paddd ($ABEF,&QWP(0,"esp"));
653 &jnz (&label("loop_shaext"));
654
655 &pshufd ($CDGH,$CDGH,0xb1); # DCHG
656 &pshufd ($TMP,$ABEF,0x1b); # FEBA
657 &pshufd ($ABEF,$ABEF,0xb1); # BAFE
658 &punpckhqdq ($ABEF,$CDGH); # DCBA
659 &palignr ($CDGH,$TMP,8); # HGFE
660
661 &mov ("esp",&DWP(32+12,"esp"));
662 &movdqu (&QWP(0,$ctx),$ABEF);
663 &movdqu (&QWP(16,$ctx),$CDGH);
664&function_end_A();
665}
666
667my @X = map("xmm$_",(0..3));
668my ($t0,$t1,$t2,$t3) = map("xmm$_",(4..7));
669my @AH = ($A,$T);
670
671&set_label("SSSE3",32);
672 &lea ("esp",&DWP(-96,"esp"));
673 # copy ctx->h[0-7] to A,B,C,D,E,F,G,H on stack
674 &mov ($AH[0],&DWP(0,"esi"));
675 &mov ($AH[1],&DWP(4,"esi"));
676 &mov ("ecx",&DWP(8,"esi"));
677 &mov ("edi",&DWP(12,"esi"));
678 #&mov (&DWP(0,"esp"),$AH[0]);
679 &mov (&DWP(4,"esp"),$AH[1]);
680 &xor ($AH[1],"ecx"); # magic
681 &mov (&DWP(8,"esp"),"ecx");
682 &mov (&DWP(12,"esp"),"edi");
683 &mov ($E,&DWP(16,"esi"));
684 &mov ("edi",&DWP(20,"esi"));
685 &mov ("ecx",&DWP(24,"esi"));
686 &mov ("esi",&DWP(28,"esi"));
687 #&mov (&DWP(16,"esp"),$E);
688 &mov (&DWP(20,"esp"),"edi");
689 &mov ("edi",&DWP(96+4,"esp")); # inp
690 &mov (&DWP(24,"esp"),"ecx");
691 &mov (&DWP(28,"esp"),"esi");
692 &movdqa ($t3,&QWP(256,$K256));
693 &jmp (&label("grand_ssse3"));
694
695&set_label("grand_ssse3",16);
696 # load input, reverse byte order, add K256[0..15], save to stack
697 &movdqu (@X[0],&QWP(0,"edi"));
698 &movdqu (@X[1],&QWP(16,"edi"));
699 &movdqu (@X[2],&QWP(32,"edi"));
700 &movdqu (@X[3],&QWP(48,"edi"));
701 &add ("edi",64);
702 &pshufb (@X[0],$t3);
703 &mov (&DWP(96+4,"esp"),"edi");
704 &pshufb (@X[1],$t3);
705 &movdqa ($t0,&QWP(0,$K256));
706 &pshufb (@X[2],$t3);
707 &movdqa ($t1,&QWP(16,$K256));
708 &paddd ($t0,@X[0]);
709 &pshufb (@X[3],$t3);
710 &movdqa ($t2,&QWP(32,$K256));
711 &paddd ($t1,@X[1]);
712 &movdqa ($t3,&QWP(48,$K256));
713 &movdqa (&QWP(32+0,"esp"),$t0);
714 &paddd ($t2,@X[2]);
715 &movdqa (&QWP(32+16,"esp"),$t1);
716 &paddd ($t3,@X[3]);
717 &movdqa (&QWP(32+32,"esp"),$t2);
718 &movdqa (&QWP(32+48,"esp"),$t3);
719 &jmp (&label("ssse3_00_47"));
720
721&set_label("ssse3_00_47",16);
722 &add ($K256,64);
723
724sub SSSE3_00_47 () {
725my $j = shift;
726my $body = shift;
727my @X = @_;
728my @insns = (&$body,&$body,&$body,&$body); # 120 instructions
729
730 eval(shift(@insns));
731 &movdqa ($t0,@X[1]);
732 eval(shift(@insns)); # @
733 eval(shift(@insns));
734 &movdqa ($t3,@X[3]);
735 eval(shift(@insns));
736 eval(shift(@insns));
737 &palignr ($t0,@X[0],4); # X[1..4]
738 eval(shift(@insns));
739 eval(shift(@insns)); # @
740 eval(shift(@insns));
741 &palignr ($t3,@X[2],4); # X[9..12]
742 eval(shift(@insns));
743 eval(shift(@insns));
744 eval(shift(@insns));
745 &movdqa ($t1,$t0);
746 eval(shift(@insns)); # @
747 eval(shift(@insns));
748 &movdqa ($t2,$t0);
749 eval(shift(@insns));
750 eval(shift(@insns));
751 &psrld ($t0,3);
752 eval(shift(@insns));
753 eval(shift(@insns)); # @
754 &paddd (@X[0],$t3); # X[0..3] += X[9..12]
755 eval(shift(@insns));
756 eval(shift(@insns));
757 &psrld ($t2,7);
758 eval(shift(@insns));
759 eval(shift(@insns));
760 eval(shift(@insns)); # @
761 eval(shift(@insns));
762 &pshufd ($t3,@X[3],0b11111010); # X[14..15]
763 eval(shift(@insns));
764 eval(shift(@insns));
765 &pslld ($t1,32-18);
766 eval(shift(@insns));
767 eval(shift(@insns)); # @
768 &pxor ($t0,$t2);
769 eval(shift(@insns));
770 eval(shift(@insns));
771 &psrld ($t2,18-7);
772 eval(shift(@insns));
773 eval(shift(@insns));
774 eval(shift(@insns)); # @
775 &pxor ($t0,$t1);
776 eval(shift(@insns));
777 eval(shift(@insns));
778 &pslld ($t1,18-7);
779 eval(shift(@insns));
780 eval(shift(@insns));
781 eval(shift(@insns)); # @
782 &pxor ($t0,$t2);
783 eval(shift(@insns));
784 eval(shift(@insns));
785 &movdqa ($t2,$t3);
786 eval(shift(@insns));
787 eval(shift(@insns));
788 eval(shift(@insns)); # @
789 &pxor ($t0,$t1); # sigma0(X[1..4])
790 eval(shift(@insns));
791 eval(shift(@insns));
792 &psrld ($t3,10);
793 eval(shift(@insns));
794 eval(shift(@insns));
795 eval(shift(@insns)); # @
796 &paddd (@X[0],$t0); # X[0..3] += sigma0(X[1..4])
797 eval(shift(@insns));
798 eval(shift(@insns));
799 &psrlq ($t2,17);
800 eval(shift(@insns));
801 eval(shift(@insns));
802 eval(shift(@insns)); # @
803 &pxor ($t3,$t2);
804 eval(shift(@insns));
805 eval(shift(@insns));
806 &psrlq ($t2,19-17);
807 eval(shift(@insns));
808 eval(shift(@insns));
809 eval(shift(@insns)); # @
810 &pxor ($t3,$t2);
811 eval(shift(@insns));
812 eval(shift(@insns));
813 &pshufd ($t3,$t3,0b10000000);
814 eval(shift(@insns));
815 eval(shift(@insns));
816 eval(shift(@insns)); # @
817 eval(shift(@insns));
818 eval(shift(@insns));
819 eval(shift(@insns));
820 eval(shift(@insns));
821 eval(shift(@insns)); # @
822 eval(shift(@insns));
823 &psrldq ($t3,8);
824 eval(shift(@insns));
825 eval(shift(@insns));
826 eval(shift(@insns));
827 &paddd (@X[0],$t3); # X[0..1] += sigma1(X[14..15])
828 eval(shift(@insns)); # @
829 eval(shift(@insns));
830 eval(shift(@insns));
831 eval(shift(@insns));
832 eval(shift(@insns));
833 eval(shift(@insns)); # @
834 eval(shift(@insns));
835 &pshufd ($t3,@X[0],0b01010000); # X[16..17]
836 eval(shift(@insns));
837 eval(shift(@insns));
838 eval(shift(@insns));
839 &movdqa ($t2,$t3);
840 eval(shift(@insns)); # @
841 &psrld ($t3,10);
842 eval(shift(@insns));
843 &psrlq ($t2,17);
844 eval(shift(@insns));
845 eval(shift(@insns));
846 eval(shift(@insns));
847 eval(shift(@insns)); # @
848 &pxor ($t3,$t2);
849 eval(shift(@insns));
850 eval(shift(@insns));
851 &psrlq ($t2,19-17);
852 eval(shift(@insns));
853 eval(shift(@insns));
854 eval(shift(@insns)); # @
855 &pxor ($t3,$t2);
856 eval(shift(@insns));
857 eval(shift(@insns));
858 eval(shift(@insns));
859 &pshufd ($t3,$t3,0b00001000);
860 eval(shift(@insns));
861 eval(shift(@insns)); # @
862 &movdqa ($t2,&QWP(16*$j,$K256));
863 eval(shift(@insns));
864 eval(shift(@insns));
865 &pslldq ($t3,8);
866 eval(shift(@insns));
867 eval(shift(@insns));
868 eval(shift(@insns)); # @
869 eval(shift(@insns));
870 eval(shift(@insns));
871 eval(shift(@insns));
872 eval(shift(@insns));
873 eval(shift(@insns)); # @
874 &paddd (@X[0],$t3); # X[2..3] += sigma1(X[16..17])
875 eval(shift(@insns));
876 eval(shift(@insns));
877 eval(shift(@insns));
878 eval(shift(@insns));
879 &paddd ($t2,@X[0]);
880 eval(shift(@insns)); # @
881
882 foreach (@insns) { eval; } # remaining instructions
883
884 &movdqa (&QWP(32+16*$j,"esp"),$t2);
885}
886
887sub body_00_15 () {
888 (
889 '&mov ("ecx",$E);',
890 '&ror ($E,25-11);',
891 '&mov ("esi",&off($f));',
892 '&xor ($E,"ecx");',
893 '&mov ("edi",&off($g));',
894 '&xor ("esi","edi");',
895 '&ror ($E,11-6);',
896 '&and ("esi","ecx");',
897 '&mov (&off($e),"ecx");', # save $E, modulo-scheduled
898 '&xor ($E,"ecx");',
899 '&xor ("edi","esi");', # Ch(e,f,g)
900 '&ror ($E,6);', # T = Sigma1(e)
901 '&mov ("ecx",$AH[0]);',
902 '&add ($E,"edi");', # T += Ch(e,f,g)
903 '&mov ("edi",&off($b));',
904 '&mov ("esi",$AH[0]);',
905
906 '&ror ("ecx",22-13);',
907 '&mov (&off($a),$AH[0]);', # save $A, modulo-scheduled
908 '&xor ("ecx",$AH[0]);',
909 '&xor ($AH[0],"edi");', # a ^= b, (b^c) in next round
910 '&add ($E,&off($h));', # T += h
911 '&ror ("ecx",13-2);',
912 '&and ($AH[1],$AH[0]);', # (b^c) &= (a^b)
913 '&xor ("ecx","esi");',
914 '&add ($E,&DWP(32+4*($i&15),"esp"));', # T += K[i]+X[i]
915 '&xor ($AH[1],"edi");', # h = Maj(a,b,c) = Ch(a^b,c,b)
916 '&ror ("ecx",2);', # Sigma0(a)
917
918 '&add ($AH[1],$E);', # h += T
919 '&add ($E,&off($d));', # d += T
920 '&add ($AH[1],"ecx");'. # h += Sigma0(a)
921
922 '@AH = reverse(@AH); $i++;' # rotate(a,h)
923 );
924}
925
926 for ($i=0,$j=0; $j<4; $j++) {
927 &SSSE3_00_47($j,\&body_00_15,@X);
928 push(@X,shift(@X)); # rotate(@X)
929 }
930 &cmp (&DWP(16*$j,$K256),0x00010203);
931 &jne (&label("ssse3_00_47"));
932
933 for ($i=0; $i<16; ) {
934 foreach(body_00_15()) { eval; }
935 }
936
937 &mov ("esi",&DWP(96,"esp")); #ctx
938 #&mov ($AH[0],&DWP(0,"esp"));
939 &xor ($AH[1],"edi"); #&mov ($AH[1],&DWP(4,"esp"));
940 #&mov ("edi", &DWP(8,"esp"));
941 &mov ("ecx",&DWP(12,"esp"));
942 &add ($AH[0],&DWP(0,"esi"));
943 &add ($AH[1],&DWP(4,"esi"));
944 &add ("edi",&DWP(8,"esi"));
945 &add ("ecx",&DWP(12,"esi"));
946 &mov (&DWP(0,"esi"),$AH[0]);
947 &mov (&DWP(4,"esi"),$AH[1]);
948 &mov (&DWP(8,"esi"),"edi");
949 &mov (&DWP(12,"esi"),"ecx");
950 #&mov (&DWP(0,"esp"),$AH[0]);
951 &mov (&DWP(4,"esp"),$AH[1]);
952 &xor ($AH[1],"edi"); # magic
953 &mov (&DWP(8,"esp"),"edi");
954 &mov (&DWP(12,"esp"),"ecx");
955 #&mov ($E,&DWP(16,"esp"));
956 &mov ("edi",&DWP(20,"esp"));
957 &mov ("ecx",&DWP(24,"esp"));
958 &add ($E,&DWP(16,"esi"));
959 &add ("edi",&DWP(20,"esi"));
960 &add ("ecx",&DWP(24,"esi"));
961 &mov (&DWP(16,"esi"),$E);
962 &mov (&DWP(20,"esi"),"edi");
963 &mov (&DWP(20,"esp"),"edi");
964 &mov ("edi",&DWP(28,"esp"));
965 &mov (&DWP(24,"esi"),"ecx");
966 #&mov (&DWP(16,"esp"),$E);
967 &add ("edi",&DWP(28,"esi"));
968 &mov (&DWP(24,"esp"),"ecx");
969 &mov (&DWP(28,"esi"),"edi");
970 &mov (&DWP(28,"esp"),"edi");
971 &mov ("edi",&DWP(96+4,"esp")); # inp
972
973 &movdqa ($t3,&QWP(64,$K256));
974 &sub ($K256,3*64); # rewind K
975 &cmp ("edi",&DWP(96+8,"esp")); # are we done yet?
976 &jb (&label("grand_ssse3"));
977
978 &mov ("esp",&DWP(96+12,"esp")); # restore sp
979&function_end_A();
980 if ($avx) {
981&set_label("AVX",32);
982 if ($avx>1) {
983 &and ("edx",1<<8|1<<3); # check for BMI2+BMI1
984 &cmp ("edx",1<<8|1<<3);
985 &je (&label("AVX_BMI"));
986 }
987 &lea ("esp",&DWP(-96,"esp"));
988 &vzeroall ();
989 # copy ctx->h[0-7] to A,B,C,D,E,F,G,H on stack
990 &mov ($AH[0],&DWP(0,"esi"));
991 &mov ($AH[1],&DWP(4,"esi"));
992 &mov ("ecx",&DWP(8,"esi"));
993 &mov ("edi",&DWP(12,"esi"));
994 #&mov (&DWP(0,"esp"),$AH[0]);
995 &mov (&DWP(4,"esp"),$AH[1]);
996 &xor ($AH[1],"ecx"); # magic
997 &mov (&DWP(8,"esp"),"ecx");
998 &mov (&DWP(12,"esp"),"edi");
999 &mov ($E,&DWP(16,"esi"));
1000 &mov ("edi",&DWP(20,"esi"));
1001 &mov ("ecx",&DWP(24,"esi"));
1002 &mov ("esi",&DWP(28,"esi"));
1003 #&mov (&DWP(16,"esp"),$E);
1004 &mov (&DWP(20,"esp"),"edi");
1005 &mov ("edi",&DWP(96+4,"esp")); # inp
1006 &mov (&DWP(24,"esp"),"ecx");
1007 &mov (&DWP(28,"esp"),"esi");
1008 &vmovdqa ($t3,&QWP(256,$K256));
1009 &jmp (&label("grand_avx"));
1010
1011&set_label("grand_avx",32);
1012 # load input, reverse byte order, add K256[0..15], save to stack
1013 &vmovdqu (@X[0],&QWP(0,"edi"));
1014 &vmovdqu (@X[1],&QWP(16,"edi"));
1015 &vmovdqu (@X[2],&QWP(32,"edi"));
1016 &vmovdqu (@X[3],&QWP(48,"edi"));
1017 &add ("edi",64);
1018 &vpshufb (@X[0],@X[0],$t3);
1019 &mov (&DWP(96+4,"esp"),"edi");
1020 &vpshufb (@X[1],@X[1],$t3);
1021 &vpshufb (@X[2],@X[2],$t3);
1022 &vpaddd ($t0,@X[0],&QWP(0,$K256));
1023 &vpshufb (@X[3],@X[3],$t3);
1024 &vpaddd ($t1,@X[1],&QWP(16,$K256));
1025 &vpaddd ($t2,@X[2],&QWP(32,$K256));
1026 &vpaddd ($t3,@X[3],&QWP(48,$K256));
1027 &vmovdqa (&QWP(32+0,"esp"),$t0);
1028 &vmovdqa (&QWP(32+16,"esp"),$t1);
1029 &vmovdqa (&QWP(32+32,"esp"),$t2);
1030 &vmovdqa (&QWP(32+48,"esp"),$t3);
1031 &jmp (&label("avx_00_47"));
1032
1033&set_label("avx_00_47",16);
1034 &add ($K256,64);
1035
1036sub Xupdate_AVX () {
1037 (
1038 '&vpalignr ($t0,@X[1],@X[0],4);', # X[1..4]
1039 '&vpalignr ($t3,@X[3],@X[2],4);', # X[9..12]
1040 '&vpsrld ($t2,$t0,7);',
1041 '&vpaddd (@X[0],@X[0],$t3);', # X[0..3] += X[9..16]
1042 '&vpsrld ($t3,$t0,3);',
1043 '&vpslld ($t1,$t0,14);',
1044 '&vpxor ($t0,$t3,$t2);',
1045 '&vpshufd ($t3,@X[3],0b11111010)',# X[14..15]
1046 '&vpsrld ($t2,$t2,18-7);',
1047 '&vpxor ($t0,$t0,$t1);',
1048 '&vpslld ($t1,$t1,25-14);',
1049 '&vpxor ($t0,$t0,$t2);',
1050 '&vpsrld ($t2,$t3,10);',
1051 '&vpxor ($t0,$t0,$t1);', # sigma0(X[1..4])
1052 '&vpsrlq ($t1,$t3,17);',
1053 '&vpaddd (@X[0],@X[0],$t0);', # X[0..3] += sigma0(X[1..4])
1054 '&vpxor ($t2,$t2,$t1);',
1055 '&vpsrlq ($t3,$t3,19);',
1056 '&vpxor ($t2,$t2,$t3);', # sigma1(X[14..15]
1057 '&vpshufd ($t3,$t2,0b10000100);',
1058 '&vpsrldq ($t3,$t3,8);',
1059 '&vpaddd (@X[0],@X[0],$t3);', # X[0..1] += sigma1(X[14..15])
1060 '&vpshufd ($t3,@X[0],0b01010000)',# X[16..17]
1061 '&vpsrld ($t2,$t3,10);',
1062 '&vpsrlq ($t1,$t3,17);',
1063 '&vpxor ($t2,$t2,$t1);',
1064 '&vpsrlq ($t3,$t3,19);',
1065 '&vpxor ($t2,$t2,$t3);', # sigma1(X[16..17]
1066 '&vpshufd ($t3,$t2,0b11101000);',
1067 '&vpslldq ($t3,$t3,8);',
1068 '&vpaddd (@X[0],@X[0],$t3);' # X[2..3] += sigma1(X[16..17])
1069 );
1070}
1071
1072local *ror = sub { &shrd(@_[0],@_) };
1073sub AVX_00_47 () {
1074my $j = shift;
1075my $body = shift;
1076my @X = @_;
1077my @insns = (&$body,&$body,&$body,&$body); # 120 instructions
1078my $insn;
1079
1080 foreach (Xupdate_AVX()) { # 31 instructions
1081 eval;
1082 eval(shift(@insns));
1083 eval(shift(@insns));
1084 eval($insn = shift(@insns));
1085 eval(shift(@insns)) if ($insn =~ /rorx/ && @insns[0] =~ /rorx/);
1086 }
1087 &vpaddd ($t2,@X[0],&QWP(16*$j,$K256));
1088 foreach (@insns) { eval; } # remaining instructions
1089 &vmovdqa (&QWP(32+16*$j,"esp"),$t2);
1090}
1091
1092 for ($i=0,$j=0; $j<4; $j++) {
1093 &AVX_00_47($j,\&body_00_15,@X);
1094 push(@X,shift(@X)); # rotate(@X)
1095 }
1096 &cmp (&DWP(16*$j,$K256),0x00010203);
1097 &jne (&label("avx_00_47"));
1098
1099 for ($i=0; $i<16; ) {
1100 foreach(body_00_15()) { eval; }
1101 }
1102
1103 &mov ("esi",&DWP(96,"esp")); #ctx
1104 #&mov ($AH[0],&DWP(0,"esp"));
1105 &xor ($AH[1],"edi"); #&mov ($AH[1],&DWP(4,"esp"));
1106 #&mov ("edi", &DWP(8,"esp"));
1107 &mov ("ecx",&DWP(12,"esp"));
1108 &add ($AH[0],&DWP(0,"esi"));
1109 &add ($AH[1],&DWP(4,"esi"));
1110 &add ("edi",&DWP(8,"esi"));
1111 &add ("ecx",&DWP(12,"esi"));
1112 &mov (&DWP(0,"esi"),$AH[0]);
1113 &mov (&DWP(4,"esi"),$AH[1]);
1114 &mov (&DWP(8,"esi"),"edi");
1115 &mov (&DWP(12,"esi"),"ecx");
1116 #&mov (&DWP(0,"esp"),$AH[0]);
1117 &mov (&DWP(4,"esp"),$AH[1]);
1118 &xor ($AH[1],"edi"); # magic
1119 &mov (&DWP(8,"esp"),"edi");
1120 &mov (&DWP(12,"esp"),"ecx");
1121 #&mov ($E,&DWP(16,"esp"));
1122 &mov ("edi",&DWP(20,"esp"));
1123 &mov ("ecx",&DWP(24,"esp"));
1124 &add ($E,&DWP(16,"esi"));
1125 &add ("edi",&DWP(20,"esi"));
1126 &add ("ecx",&DWP(24,"esi"));
1127 &mov (&DWP(16,"esi"),$E);
1128 &mov (&DWP(20,"esi"),"edi");
1129 &mov (&DWP(20,"esp"),"edi");
1130 &mov ("edi",&DWP(28,"esp"));
1131 &mov (&DWP(24,"esi"),"ecx");
1132 #&mov (&DWP(16,"esp"),$E);
1133 &add ("edi",&DWP(28,"esi"));
1134 &mov (&DWP(24,"esp"),"ecx");
1135 &mov (&DWP(28,"esi"),"edi");
1136 &mov (&DWP(28,"esp"),"edi");
1137 &mov ("edi",&DWP(96+4,"esp")); # inp
1138
1139 &vmovdqa ($t3,&QWP(64,$K256));
1140 &sub ($K256,3*64); # rewind K
1141 &cmp ("edi",&DWP(96+8,"esp")); # are we done yet?
1142 &jb (&label("grand_avx"));
1143
1144 &mov ("esp",&DWP(96+12,"esp")); # restore sp
1145 &vzeroall ();
1146&function_end_A();
1147 if ($avx>1) {
1148sub bodyx_00_15 () { # +10%
1149 (
1150 '&rorx ("ecx",$E,6)',
1151 '&rorx ("esi",$E,11)',
1152 '&mov (&off($e),$E)', # save $E, modulo-scheduled
1153 '&rorx ("edi",$E,25)',
1154 '&xor ("ecx","esi")',
1155 '&andn ("esi",$E,&off($g))',
1156 '&xor ("ecx","edi")', # Sigma1(e)
1157 '&and ($E,&off($f))',
1158 '&mov (&off($a),$AH[0]);', # save $A, modulo-scheduled
1159 '&or ($E,"esi")', # T = Ch(e,f,g)
1160
1161 '&rorx ("edi",$AH[0],2)',
1162 '&rorx ("esi",$AH[0],13)',
1163 '&lea ($E,&DWP(0,$E,"ecx"))', # T += Sigma1(e)
1164 '&rorx ("ecx",$AH[0],22)',
1165 '&xor ("esi","edi")',
1166 '&mov ("edi",&off($b))',
1167 '&xor ("ecx","esi")', # Sigma0(a)
1168
1169 '&xor ($AH[0],"edi")', # a ^= b, (b^c) in next round
1170 '&add ($E,&off($h))', # T += h
1171 '&and ($AH[1],$AH[0])', # (b^c) &= (a^b)
1172 '&add ($E,&DWP(32+4*($i&15),"esp"))', # T += K[i]+X[i]
1173 '&xor ($AH[1],"edi")', # h = Maj(a,b,c) = Ch(a^b,c,b)
1174
1175 '&add ("ecx",$E)', # h += T
1176 '&add ($E,&off($d))', # d += T
1177 '&lea ($AH[1],&DWP(0,$AH[1],"ecx"));'. # h += Sigma0(a)
1178
1179 '@AH = reverse(@AH); $i++;' # rotate(a,h)
1180 );
1181}
1182
1183&set_label("AVX_BMI",32);
1184 &lea ("esp",&DWP(-96,"esp"));
1185 &vzeroall ();
1186 # copy ctx->h[0-7] to A,B,C,D,E,F,G,H on stack
1187 &mov ($AH[0],&DWP(0,"esi"));
1188 &mov ($AH[1],&DWP(4,"esi"));
1189 &mov ("ecx",&DWP(8,"esi"));
1190 &mov ("edi",&DWP(12,"esi"));
1191 #&mov (&DWP(0,"esp"),$AH[0]);
1192 &mov (&DWP(4,"esp"),$AH[1]);
1193 &xor ($AH[1],"ecx"); # magic
1194 &mov (&DWP(8,"esp"),"ecx");
1195 &mov (&DWP(12,"esp"),"edi");
1196 &mov ($E,&DWP(16,"esi"));
1197 &mov ("edi",&DWP(20,"esi"));
1198 &mov ("ecx",&DWP(24,"esi"));
1199 &mov ("esi",&DWP(28,"esi"));
1200 #&mov (&DWP(16,"esp"),$E);
1201 &mov (&DWP(20,"esp"),"edi");
1202 &mov ("edi",&DWP(96+4,"esp")); # inp
1203 &mov (&DWP(24,"esp"),"ecx");
1204 &mov (&DWP(28,"esp"),"esi");
1205 &vmovdqa ($t3,&QWP(256,$K256));
1206 &jmp (&label("grand_avx_bmi"));
1207
1208&set_label("grand_avx_bmi",32);
1209 # load input, reverse byte order, add K256[0..15], save to stack
1210 &vmovdqu (@X[0],&QWP(0,"edi"));
1211 &vmovdqu (@X[1],&QWP(16,"edi"));
1212 &vmovdqu (@X[2],&QWP(32,"edi"));
1213 &vmovdqu (@X[3],&QWP(48,"edi"));
1214 &add ("edi",64);
1215 &vpshufb (@X[0],@X[0],$t3);
1216 &mov (&DWP(96+4,"esp"),"edi");
1217 &vpshufb (@X[1],@X[1],$t3);
1218 &vpshufb (@X[2],@X[2],$t3);
1219 &vpaddd ($t0,@X[0],&QWP(0,$K256));
1220 &vpshufb (@X[3],@X[3],$t3);
1221 &vpaddd ($t1,@X[1],&QWP(16,$K256));
1222 &vpaddd ($t2,@X[2],&QWP(32,$K256));
1223 &vpaddd ($t3,@X[3],&QWP(48,$K256));
1224 &vmovdqa (&QWP(32+0,"esp"),$t0);
1225 &vmovdqa (&QWP(32+16,"esp"),$t1);
1226 &vmovdqa (&QWP(32+32,"esp"),$t2);
1227 &vmovdqa (&QWP(32+48,"esp"),$t3);
1228 &jmp (&label("avx_bmi_00_47"));
1229
1230&set_label("avx_bmi_00_47",16);
1231 &add ($K256,64);
1232
1233 for ($i=0,$j=0; $j<4; $j++) {
1234 &AVX_00_47($j,\&bodyx_00_15,@X);
1235 push(@X,shift(@X)); # rotate(@X)
1236 }
1237 &cmp (&DWP(16*$j,$K256),0x00010203);
1238 &jne (&label("avx_bmi_00_47"));
1239
1240 for ($i=0; $i<16; ) {
1241 foreach(bodyx_00_15()) { eval; }
1242 }
1243
1244 &mov ("esi",&DWP(96,"esp")); #ctx
1245 #&mov ($AH[0],&DWP(0,"esp"));
1246 &xor ($AH[1],"edi"); #&mov ($AH[1],&DWP(4,"esp"));
1247 #&mov ("edi", &DWP(8,"esp"));
1248 &mov ("ecx",&DWP(12,"esp"));
1249 &add ($AH[0],&DWP(0,"esi"));
1250 &add ($AH[1],&DWP(4,"esi"));
1251 &add ("edi",&DWP(8,"esi"));
1252 &add ("ecx",&DWP(12,"esi"));
1253 &mov (&DWP(0,"esi"),$AH[0]);
1254 &mov (&DWP(4,"esi"),$AH[1]);
1255 &mov (&DWP(8,"esi"),"edi");
1256 &mov (&DWP(12,"esi"),"ecx");
1257 #&mov (&DWP(0,"esp"),$AH[0]);
1258 &mov (&DWP(4,"esp"),$AH[1]);
1259 &xor ($AH[1],"edi"); # magic
1260 &mov (&DWP(8,"esp"),"edi");
1261 &mov (&DWP(12,"esp"),"ecx");
1262 #&mov ($E,&DWP(16,"esp"));
1263 &mov ("edi",&DWP(20,"esp"));
1264 &mov ("ecx",&DWP(24,"esp"));
1265 &add ($E,&DWP(16,"esi"));
1266 &add ("edi",&DWP(20,"esi"));
1267 &add ("ecx",&DWP(24,"esi"));
1268 &mov (&DWP(16,"esi"),$E);
1269 &mov (&DWP(20,"esi"),"edi");
1270 &mov (&DWP(20,"esp"),"edi");
1271 &mov ("edi",&DWP(28,"esp"));
1272 &mov (&DWP(24,"esi"),"ecx");
1273 #&mov (&DWP(16,"esp"),$E);
1274 &add ("edi",&DWP(28,"esi"));
1275 &mov (&DWP(24,"esp"),"ecx");
1276 &mov (&DWP(28,"esi"),"edi");
1277 &mov (&DWP(28,"esp"),"edi");
1278 &mov ("edi",&DWP(96+4,"esp")); # inp
1279
1280 &vmovdqa ($t3,&QWP(64,$K256));
1281 &sub ($K256,3*64); # rewind K
1282 &cmp ("edi",&DWP(96+8,"esp")); # are we done yet?
1283 &jb (&label("grand_avx_bmi"));
1284
1285 &mov ("esp",&DWP(96+12,"esp")); # restore sp
1286 &vzeroall ();
1287&function_end_A();
1288 }
1289 }
1290 }}}
1291&function_end_B("sha256_block_data_order");
1292
1293&asm_finish();
1294
1295close STDOUT or die "error closing STDOUT: $!";
Note: See TracBrowser for help on using the repository browser.

© 2024 Oracle Support Privacy / Do Not Sell My Info Terms of Use Trademark Policy Automated Access Etiquette