2 # Copyright 2016 The OpenSSL Project Authors. All Rights Reserved.
4 # Licensed under the OpenSSL license (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
10 # ====================================================================
11 # Written by Andy Polyakov <appro@openssl.org> 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 # ====================================================================
17 # This module implements Poly1305 hash for x86_64.
25 # Add AVX512F+VL+BW code path.
27 # Numbers are cycles per processed byte with poly1305_blocks alone,
28 # measured with rdtsc at fixed clock frequency.
30 # IALU/gcc-4.8(*) AVX(**) AVX2 AVX-512
33 # Westmere 1.88/+120% -
34 # Sandy Bridge 1.39/+140% 1.10
35 # Haswell 1.14/+175% 1.11 0.65
36 # Skylake[-X] 1.13/+120% 0.96 0.51 [0.35]
37 # Silvermont 2.83/+95% -
38 # Knights L 3.60/- 1.65 1.10 (***)
39 # Goldmont 1.70/+180% -
40 # VIA Nano 1.82/+150% -
41 # Sledgehammer 1.38/+160% -
42 # Bulldozer 2.30/+130% 0.97
43 # Ryzen 1.15/+200% 1.08 1.18
45 # (*) improvement coefficients relative to clang are more modest and
46 # are ~50% on most processors, in both cases we are comparing to
48 # (**) SSE2 implementation was attempted, but among non-AVX processors
49 # it was faster than integer-only code only on older Intel P4 and
50 # Core processors, 50-30%, less newer processor is, but slower on
51 # contemporary ones, for example almost 2x slower on Atom, and as
52 # former are naturally disappearing, SSE2 is deemed unnecessary;
53 # (***) Current AVX-512 code requires BW and VL extensions and can not
54 # execute on Knights Landing;
58 if ($flavour =~ /\./) { $output = $flavour; undef $flavour; }
60 $win64=0; $win64=1 if ($flavour =~ /[nm]asm|mingw64/ || $output =~ /\.asm$/);
62 $0 =~ m/(.*[\/\\])[^\/\\]+$/; $dir=$1;
63 ( $xlate="${dir}x86_64-xlate.pl" and -f $xlate ) or
64 ( $xlate="${dir}../../perlasm/x86_64-xlate.pl" and -f $xlate) or
65 die "can't locate x86_64-xlate.pl";
67 if (`$ENV{CC} -Wa,-v -c -o /dev/null -x assembler /dev/null 2>&1`
68 =~ /GNU assembler version ([2-9]\.[0-9]+)/) {
69 $avx = ($1>=2.19) + ($1>=2.22) + ($1>=2.25) + ($1>=2.26);
72 if (!$avx && $win64 && ($flavour =~ /nasm/ || $ENV{ASM} =~ /nasm/) &&
73 `nasm -v 2>&1` =~ /NASM version ([2-9]\.[0-9]+)(?:\.([0-9]+))?/) {
74 $avx = ($1>=2.09) + ($1>=2.10) + 2 * ($1>=2.12);
75 $avx += 2 if ($1==2.11 && $2>=8);
78 if (!$avx && $win64 && ($flavour =~ /masm/ || $ENV{ASM} =~ /ml64/) &&
79 `ml64 2>&1` =~ /Version ([0-9]+)\./) {
80 $avx = ($1>=10) + ($1>=12);
83 if (!$avx && `$ENV{CC} -v 2>&1` =~ /((?:^clang|LLVM) version|.*based on LLVM) ([3-9]\.[0-9]+)/) {
84 $avx = ($2>=3.0) + ($2>3.0);
87 open OUT,"| \"$^X\" \"$xlate\" $flavour \"$output\"";
90 my ($ctx,$inp,$len,$padbit)=("%rdi","%rsi","%rdx","%rcx");
91 my ($mac,$nonce)=($inp,$len); # *_emit arguments
92 my ($d1,$d2,$d3, $r0,$r1,$s1)=map("%r$_",(8..13));
93 my ($h0,$h1,$h2)=("%r14","%rbx","%rbp");
95 sub poly1305_iteration {
96 # input: copy of $r1 in %rax, $h0-$h2, $r0-$r1
97 # output: $h0-$h2 *= $r0-$r1
105 mov %rax,$h0 # future $h0
115 mov $h2,$h1 # borrow $h1
119 imulq $s1,$h1 # h2*s1
124 imulq $r0,$h2 # h2*r0
126 mov \$-4,%rax # mask value
129 and $d3,%rax # last reduction step
140 ########################################################################
141 # Layout of opaque area is following.
143 # unsigned __int64 h[3]; # current hash value base 2^64
144 # unsigned __int64 r[2]; # key value base 2^64
149 .extern OPENSSL_ia32cap_P
152 .hidden poly1305_init
153 .globl poly1305_blocks
154 .hidden poly1305_blocks
156 .hidden poly1305_emit
158 .type poly1305_init,\@function,3
162 mov %rax,0($ctx) # initialize hash value
169 lea poly1305_blocks(%rip),%r10
170 lea poly1305_emit(%rip),%r11
172 $code.=<<___ if ($avx);
173 mov OPENSSL_ia32cap_P+4(%rip),%r9
174 lea poly1305_blocks_avx(%rip),%rax
175 lea poly1305_emit_avx(%rip),%rcx
176 bt \$`60-32`,%r9 # AVX?
180 $code.=<<___ if ($avx>1);
181 lea poly1305_blocks_avx2(%rip),%rax
182 bt \$`5+32`,%r9 # AVX2?
185 $code.=<<___ if ($avx>3);
186 mov \$`(1<<31|1<<21|1<<16)`,%rax
193 mov \$0x0ffffffc0fffffff,%rax
194 mov \$0x0ffffffc0ffffffc,%rcx
200 $code.=<<___ if ($flavour !~ /elf32/);
204 $code.=<<___ if ($flavour =~ /elf32/);
212 .size poly1305_init,.-poly1305_init
214 .type poly1305_blocks,\@function,4
220 jz .Lno_data # too short
236 mov $len,%r15 # reassign $len
238 mov 24($ctx),$r0 # load r
241 mov 0($ctx),$h0 # load hash value
248 add $r1,$s1 # s1 = r1 + (r1 >> 2)
253 add 0($inp),$h0 # accumulate input
258 &poly1305_iteration();
264 mov $h0,0($ctx) # store hash value
281 .cfi_adjust_cfa_offset -48
286 .size poly1305_blocks,.-poly1305_blocks
288 .type poly1305_emit,\@function,3
292 mov 0($ctx),%r8 # load hash value
297 add \$5,%r8 # compare to modulus
301 shr \$2,%r10 # did 130-bit value overfow?
305 add 0($nonce),%rax # accumulate nonce
307 mov %rax,0($mac) # write result
311 .size poly1305_emit,.-poly1305_emit
315 ########################################################################
316 # Layout of opaque area is following.
318 # unsigned __int32 h[5]; # current hash value base 2^26
319 # unsigned __int32 is_base2_26;
320 # unsigned __int64 r[2]; # key value base 2^64
321 # unsigned __int64 pad;
322 # struct { unsigned __int32 r^2, r^1, r^4, r^3; } r[9];
324 # where r^n are base 2^26 digits of degrees of multiplier key. There are
325 # 5 digits, but last four are interleaved with multiples of 5, totalling
326 # in 9 elements: r0, r1, 5*r1, r2, 5*r2, r3, 5*r3, r4, 5*r4.
328 my ($H0,$H1,$H2,$H3,$H4, $T0,$T1,$T2,$T3,$T4, $D0,$D1,$D2,$D3,$D4, $MASK) =
329 map("%xmm$_",(0..15));
332 .type __poly1305_block,\@abi-omnipotent
336 &poly1305_iteration();
339 .size __poly1305_block,.-__poly1305_block
341 .type __poly1305_init_avx,\@abi-omnipotent
348 lea 48+64($ctx),$ctx # size optimization
351 call __poly1305_block # r^2
353 mov \$0x3ffffff,%eax # save interleaved r^2 and r base 2^26
359 mov %eax,`16*0+0-64`($ctx)
361 mov %edx,`16*0+4-64`($ctx)
368 mov %eax,`16*1+0-64`($ctx)
369 lea (%rax,%rax,4),%eax # *5
370 mov %edx,`16*1+4-64`($ctx)
371 lea (%rdx,%rdx,4),%edx # *5
372 mov %eax,`16*2+0-64`($ctx)
374 mov %edx,`16*2+4-64`($ctx)
385 mov %eax,`16*3+0-64`($ctx)
386 lea (%rax,%rax,4),%eax # *5
387 mov %edx,`16*3+4-64`($ctx)
388 lea (%rdx,%rdx,4),%edx # *5
389 mov %eax,`16*4+0-64`($ctx)
391 mov %edx,`16*4+4-64`($ctx)
400 mov %eax,`16*5+0-64`($ctx)
401 lea (%rax,%rax,4),%eax # *5
402 mov %edx,`16*5+4-64`($ctx)
403 lea (%rdx,%rdx,4),%edx # *5
404 mov %eax,`16*6+0-64`($ctx)
406 mov %edx,`16*6+4-64`($ctx)
412 mov $d1#d,`16*7+0-64`($ctx)
413 lea ($d1,$d1,4),$d1 # *5
414 mov $d2#d,`16*7+4-64`($ctx)
415 lea ($d2,$d2,4),$d2 # *5
416 mov $d1#d,`16*8+0-64`($ctx)
417 mov $d2#d,`16*8+4-64`($ctx)
420 call __poly1305_block # r^3
422 mov \$0x3ffffff,%eax # save r^3 base 2^26
426 mov %eax,`16*0+12-64`($ctx)
430 mov %edx,`16*1+12-64`($ctx)
431 lea (%rdx,%rdx,4),%edx # *5
433 mov %edx,`16*2+12-64`($ctx)
439 mov %eax,`16*3+12-64`($ctx)
440 lea (%rax,%rax,4),%eax # *5
442 mov %eax,`16*4+12-64`($ctx)
447 mov %edx,`16*5+12-64`($ctx)
448 lea (%rdx,%rdx,4),%edx # *5
450 mov %edx,`16*6+12-64`($ctx)
455 mov $d1#d,`16*7+12-64`($ctx)
456 lea ($d1,$d1,4),$d1 # *5
457 mov $d1#d,`16*8+12-64`($ctx)
460 call __poly1305_block # r^4
462 mov \$0x3ffffff,%eax # save r^4 base 2^26
466 mov %eax,`16*0+8-64`($ctx)
470 mov %edx,`16*1+8-64`($ctx)
471 lea (%rdx,%rdx,4),%edx # *5
473 mov %edx,`16*2+8-64`($ctx)
479 mov %eax,`16*3+8-64`($ctx)
480 lea (%rax,%rax,4),%eax # *5
482 mov %eax,`16*4+8-64`($ctx)
487 mov %edx,`16*5+8-64`($ctx)
488 lea (%rdx,%rdx,4),%edx # *5
490 mov %edx,`16*6+8-64`($ctx)
495 mov $d1#d,`16*7+8-64`($ctx)
496 lea ($d1,$d1,4),$d1 # *5
497 mov $d1#d,`16*8+8-64`($ctx)
499 lea -48-64($ctx),$ctx # size [de-]optimization
501 .size __poly1305_init_avx,.-__poly1305_init_avx
503 .type poly1305_blocks_avx,\@function,4
507 mov 20($ctx),%r8d # is_base2_26
539 mov $len,%r15 # reassign $len
541 mov 0($ctx),$d1 # load hash value
545 mov 24($ctx),$r0 # load r
548 ################################# base 2^26 -> base 2^64
550 and \$`-1*(1<<31)`,$d1
551 mov $d2,$r1 # borrow $r1
553 and \$`-1*(1<<31)`,$d2
567 adc \$0,$h2 # can be partially reduced...
569 mov \$-4,$d2 # ... so reduce
582 add $r1,$s1 # s1 = r1 + (r1 >> 2)
584 add 0($inp),$h0 # accumulate input
589 call __poly1305_block
591 test $padbit,$padbit # if $padbit is zero,
592 jz .Lstore_base2_64_avx # store hash in base 2^64 format
594 ################################# base 2^64 -> base 2^26
601 and \$0x3ffffff,%rax # h[0]
603 and \$0x3ffffff,%rdx # h[1]
607 and \$0x3ffffff,$h0 # h[2]
609 and \$0x3ffffff,$h1 # h[3]
613 jz .Lstore_base2_26_avx
623 .Lstore_base2_64_avx:
626 mov $h2,16($ctx) # note that is_base2_26 is zeroed
630 .Lstore_base2_26_avx:
631 mov %rax#d,0($ctx) # store hash value base 2^26
651 .cfi_adjust_cfa_offset -48
653 .Lblocks_avx_epilogue:
674 mov $len,%r15 # reassign $len
676 mov 24($ctx),$r0 # load r
679 mov 0($ctx),$h0 # load hash value
686 add $r1,$s1 # s1 = r1 + (r1 >> 2)
691 add 0($inp),$h0 # accumulate input
697 call __poly1305_block
700 ################################# base 2^64 -> base 2^26
707 and \$0x3ffffff,%rax # h[0]
709 and \$0x3ffffff,%rdx # h[1]
713 and \$0x3ffffff,$h0 # h[2]
715 and \$0x3ffffff,$h1 # h[3]
723 movl \$1,20($ctx) # set is_base2_26
725 call __poly1305_init_avx
744 .cfi_adjust_cfa_offset -48
745 .Lbase2_64_avx_epilogue:
752 vmovd 4*0($ctx),$H0 # load hash value
760 $code.=<<___ if (!$win64);
762 .cfi_def_cfa %r11,0x60
765 $code.=<<___ if ($win64);
768 vmovdqa %xmm6,0x50(%r11)
769 vmovdqa %xmm7,0x60(%r11)
770 vmovdqa %xmm8,0x70(%r11)
771 vmovdqa %xmm9,0x80(%r11)
772 vmovdqa %xmm10,0x90(%r11)
773 vmovdqa %xmm11,0xa0(%r11)
774 vmovdqa %xmm12,0xb0(%r11)
775 vmovdqa %xmm13,0xc0(%r11)
776 vmovdqa %xmm14,0xd0(%r11)
777 vmovdqa %xmm15,0xe0(%r11)
785 vmovdqu `16*3`($ctx),$D4 # preload r0^2
786 lea `16*3+64`($ctx),$ctx # size optimization
787 lea .Lconst(%rip),%rcx
789 ################################################################
791 vmovdqu 16*2($inp),$T0
792 vmovdqu 16*3($inp),$T1
793 vmovdqa 64(%rcx),$MASK # .Lmask26
795 vpsrldq \$6,$T0,$T2 # splat input
797 vpunpckhqdq $T1,$T0,$T4 # 4
798 vpunpcklqdq $T1,$T0,$T0 # 0:1
799 vpunpcklqdq $T3,$T2,$T3 # 2:3
801 vpsrlq \$40,$T4,$T4 # 4
803 vpand $MASK,$T0,$T0 # 0
805 vpand $MASK,$T1,$T1 # 1
807 vpand $MASK,$T2,$T2 # 2
808 vpand $MASK,$T3,$T3 # 3
809 vpor 32(%rcx),$T4,$T4 # padbit, yes, always
813 # expand and copy pre-calculated table to stack
814 vmovdqu `16*1-64`($ctx),$D1
815 vmovdqu `16*2-64`($ctx),$D2
816 vpshufd \$0xEE,$D4,$D3 # 34xx -> 3434
817 vpshufd \$0x44,$D4,$D0 # xx12 -> 1212
818 vmovdqa $D3,-0x90(%r11)
819 vmovdqa $D0,0x00(%rsp)
820 vpshufd \$0xEE,$D1,$D4
821 vmovdqu `16*3-64`($ctx),$D0
822 vpshufd \$0x44,$D1,$D1
823 vmovdqa $D4,-0x80(%r11)
824 vmovdqa $D1,0x10(%rsp)
825 vpshufd \$0xEE,$D2,$D3
826 vmovdqu `16*4-64`($ctx),$D1
827 vpshufd \$0x44,$D2,$D2
828 vmovdqa $D3,-0x70(%r11)
829 vmovdqa $D2,0x20(%rsp)
830 vpshufd \$0xEE,$D0,$D4
831 vmovdqu `16*5-64`($ctx),$D2
832 vpshufd \$0x44,$D0,$D0
833 vmovdqa $D4,-0x60(%r11)
834 vmovdqa $D0,0x30(%rsp)
835 vpshufd \$0xEE,$D1,$D3
836 vmovdqu `16*6-64`($ctx),$D0
837 vpshufd \$0x44,$D1,$D1
838 vmovdqa $D3,-0x50(%r11)
839 vmovdqa $D1,0x40(%rsp)
840 vpshufd \$0xEE,$D2,$D4
841 vmovdqu `16*7-64`($ctx),$D1
842 vpshufd \$0x44,$D2,$D2
843 vmovdqa $D4,-0x40(%r11)
844 vmovdqa $D2,0x50(%rsp)
845 vpshufd \$0xEE,$D0,$D3
846 vmovdqu `16*8-64`($ctx),$D2
847 vpshufd \$0x44,$D0,$D0
848 vmovdqa $D3,-0x30(%r11)
849 vmovdqa $D0,0x60(%rsp)
850 vpshufd \$0xEE,$D1,$D4
851 vpshufd \$0x44,$D1,$D1
852 vmovdqa $D4,-0x20(%r11)
853 vmovdqa $D1,0x70(%rsp)
854 vpshufd \$0xEE,$D2,$D3
855 vmovdqa 0x00(%rsp),$D4 # preload r0^2
856 vpshufd \$0x44,$D2,$D2
857 vmovdqa $D3,-0x10(%r11)
858 vmovdqa $D2,0x80(%rsp)
864 ################################################################
865 # ((inp[0]*r^4+inp[2]*r^2+inp[4])*r^4+inp[6]*r^2
866 # ((inp[1]*r^4+inp[3]*r^2+inp[5])*r^3+inp[7]*r
867 # \___________________/
868 # ((inp[0]*r^4+inp[2]*r^2+inp[4])*r^4+inp[6]*r^2+inp[8])*r^2
869 # ((inp[1]*r^4+inp[3]*r^2+inp[5])*r^4+inp[7]*r^2+inp[9])*r
870 # \___________________/ \____________________/
872 # Note that we start with inp[2:3]*r^2. This is because it
873 # doesn't depend on reduction in previous iteration.
874 ################################################################
875 # d4 = h4*r0 + h3*r1 + h2*r2 + h1*r3 + h0*r4
876 # d3 = h3*r0 + h2*r1 + h1*r2 + h0*r3 + h4*5*r4
877 # d2 = h2*r0 + h1*r1 + h0*r2 + h4*5*r3 + h3*5*r4
878 # d1 = h1*r0 + h0*r1 + h4*5*r2 + h3*5*r3 + h2*5*r4
879 # d0 = h0*r0 + h4*5*r1 + h3*5*r2 + h2*5*r3 + h1*5*r4
881 # though note that $Tx and $Hx are "reversed" in this section,
882 # and $D4 is preloaded with r0^2...
884 vpmuludq $T0,$D4,$D0 # d0 = h0*r0
885 vpmuludq $T1,$D4,$D1 # d1 = h1*r0
886 vmovdqa $H2,0x20(%r11) # offload hash
887 vpmuludq $T2,$D4,$D2 # d3 = h2*r0
888 vmovdqa 0x10(%rsp),$H2 # r1^2
889 vpmuludq $T3,$D4,$D3 # d3 = h3*r0
890 vpmuludq $T4,$D4,$D4 # d4 = h4*r0
892 vmovdqa $H0,0x00(%r11) #
893 vpmuludq 0x20(%rsp),$T4,$H0 # h4*s1
894 vmovdqa $H1,0x10(%r11) #
895 vpmuludq $T3,$H2,$H1 # h3*r1
896 vpaddq $H0,$D0,$D0 # d0 += h4*s1
897 vpaddq $H1,$D4,$D4 # d4 += h3*r1
898 vmovdqa $H3,0x30(%r11) #
899 vpmuludq $T2,$H2,$H0 # h2*r1
900 vpmuludq $T1,$H2,$H1 # h1*r1
901 vpaddq $H0,$D3,$D3 # d3 += h2*r1
902 vmovdqa 0x30(%rsp),$H3 # r2^2
903 vpaddq $H1,$D2,$D2 # d2 += h1*r1
904 vmovdqa $H4,0x40(%r11) #
905 vpmuludq $T0,$H2,$H2 # h0*r1
906 vpmuludq $T2,$H3,$H0 # h2*r2
907 vpaddq $H2,$D1,$D1 # d1 += h0*r1
909 vmovdqa 0x40(%rsp),$H4 # s2^2
910 vpaddq $H0,$D4,$D4 # d4 += h2*r2
911 vpmuludq $T1,$H3,$H1 # h1*r2
912 vpmuludq $T0,$H3,$H3 # h0*r2
913 vpaddq $H1,$D3,$D3 # d3 += h1*r2
914 vmovdqa 0x50(%rsp),$H2 # r3^2
915 vpaddq $H3,$D2,$D2 # d2 += h0*r2
916 vpmuludq $T4,$H4,$H0 # h4*s2
917 vpmuludq $T3,$H4,$H4 # h3*s2
918 vpaddq $H0,$D1,$D1 # d1 += h4*s2
919 vmovdqa 0x60(%rsp),$H3 # s3^2
920 vpaddq $H4,$D0,$D0 # d0 += h3*s2
922 vmovdqa 0x80(%rsp),$H4 # s4^2
923 vpmuludq $T1,$H2,$H1 # h1*r3
924 vpmuludq $T0,$H2,$H2 # h0*r3
925 vpaddq $H1,$D4,$D4 # d4 += h1*r3
926 vpaddq $H2,$D3,$D3 # d3 += h0*r3
927 vpmuludq $T4,$H3,$H0 # h4*s3
928 vpmuludq $T3,$H3,$H1 # h3*s3
929 vpaddq $H0,$D2,$D2 # d2 += h4*s3
930 vmovdqu 16*0($inp),$H0 # load input
931 vpaddq $H1,$D1,$D1 # d1 += h3*s3
932 vpmuludq $T2,$H3,$H3 # h2*s3
933 vpmuludq $T2,$H4,$T2 # h2*s4
934 vpaddq $H3,$D0,$D0 # d0 += h2*s3
936 vmovdqu 16*1($inp),$H1 #
937 vpaddq $T2,$D1,$D1 # d1 += h2*s4
938 vpmuludq $T3,$H4,$T3 # h3*s4
939 vpmuludq $T4,$H4,$T4 # h4*s4
940 vpsrldq \$6,$H0,$H2 # splat input
941 vpaddq $T3,$D2,$D2 # d2 += h3*s4
942 vpaddq $T4,$D3,$D3 # d3 += h4*s4
943 vpsrldq \$6,$H1,$H3 #
944 vpmuludq 0x70(%rsp),$T0,$T4 # h0*r4
945 vpmuludq $T1,$H4,$T0 # h1*s4
946 vpunpckhqdq $H1,$H0,$H4 # 4
947 vpaddq $T4,$D4,$D4 # d4 += h0*r4
948 vmovdqa -0x90(%r11),$T4 # r0^4
949 vpaddq $T0,$D0,$D0 # d0 += h1*s4
951 vpunpcklqdq $H1,$H0,$H0 # 0:1
952 vpunpcklqdq $H3,$H2,$H3 # 2:3
954 #vpsrlq \$40,$H4,$H4 # 4
955 vpsrldq \$`40/8`,$H4,$H4 # 4
957 vpand $MASK,$H0,$H0 # 0
959 vpand $MASK,$H1,$H1 # 1
960 vpand 0(%rcx),$H4,$H4 # .Lmask24
962 vpand $MASK,$H2,$H2 # 2
963 vpand $MASK,$H3,$H3 # 3
964 vpor 32(%rcx),$H4,$H4 # padbit, yes, always
966 vpaddq 0x00(%r11),$H0,$H0 # add hash value
967 vpaddq 0x10(%r11),$H1,$H1
968 vpaddq 0x20(%r11),$H2,$H2
969 vpaddq 0x30(%r11),$H3,$H3
970 vpaddq 0x40(%r11),$H4,$H4
977 ################################################################
978 # Now we accumulate (inp[0:1]+hash)*r^4
979 ################################################################
980 # d4 = h4*r0 + h3*r1 + h2*r2 + h1*r3 + h0*r4
981 # d3 = h3*r0 + h2*r1 + h1*r2 + h0*r3 + h4*5*r4
982 # d2 = h2*r0 + h1*r1 + h0*r2 + h4*5*r3 + h3*5*r4
983 # d1 = h1*r0 + h0*r1 + h4*5*r2 + h3*5*r3 + h2*5*r4
984 # d0 = h0*r0 + h4*5*r1 + h3*5*r2 + h2*5*r3 + h1*5*r4
986 vpmuludq $H0,$T4,$T0 # h0*r0
987 vpmuludq $H1,$T4,$T1 # h1*r0
990 vmovdqa -0x80(%r11),$T2 # r1^4
991 vpmuludq $H2,$T4,$T0 # h2*r0
992 vpmuludq $H3,$T4,$T1 # h3*r0
995 vpmuludq $H4,$T4,$T4 # h4*r0
996 vpmuludq -0x70(%r11),$H4,$T0 # h4*s1
999 vpaddq $T0,$D0,$D0 # d0 += h4*s1
1000 vpmuludq $H2,$T2,$T1 # h2*r1
1001 vpmuludq $H3,$T2,$T0 # h3*r1
1002 vpaddq $T1,$D3,$D3 # d3 += h2*r1
1003 vmovdqa -0x60(%r11),$T3 # r2^4
1004 vpaddq $T0,$D4,$D4 # d4 += h3*r1
1005 vpmuludq $H1,$T2,$T1 # h1*r1
1006 vpmuludq $H0,$T2,$T2 # h0*r1
1007 vpaddq $T1,$D2,$D2 # d2 += h1*r1
1008 vpaddq $T2,$D1,$D1 # d1 += h0*r1
1010 vmovdqa -0x50(%r11),$T4 # s2^4
1011 vpmuludq $H2,$T3,$T0 # h2*r2
1012 vpmuludq $H1,$T3,$T1 # h1*r2
1013 vpaddq $T0,$D4,$D4 # d4 += h2*r2
1014 vpaddq $T1,$D3,$D3 # d3 += h1*r2
1015 vmovdqa -0x40(%r11),$T2 # r3^4
1016 vpmuludq $H0,$T3,$T3 # h0*r2
1017 vpmuludq $H4,$T4,$T0 # h4*s2
1018 vpaddq $T3,$D2,$D2 # d2 += h0*r2
1019 vpaddq $T0,$D1,$D1 # d1 += h4*s2
1020 vmovdqa -0x30(%r11),$T3 # s3^4
1021 vpmuludq $H3,$T4,$T4 # h3*s2
1022 vpmuludq $H1,$T2,$T1 # h1*r3
1023 vpaddq $T4,$D0,$D0 # d0 += h3*s2
1025 vmovdqa -0x10(%r11),$T4 # s4^4
1026 vpaddq $T1,$D4,$D4 # d4 += h1*r3
1027 vpmuludq $H0,$T2,$T2 # h0*r3
1028 vpmuludq $H4,$T3,$T0 # h4*s3
1029 vpaddq $T2,$D3,$D3 # d3 += h0*r3
1030 vpaddq $T0,$D2,$D2 # d2 += h4*s3
1031 vmovdqu 16*2($inp),$T0 # load input
1032 vpmuludq $H3,$T3,$T2 # h3*s3
1033 vpmuludq $H2,$T3,$T3 # h2*s3
1034 vpaddq $T2,$D1,$D1 # d1 += h3*s3
1035 vmovdqu 16*3($inp),$T1 #
1036 vpaddq $T3,$D0,$D0 # d0 += h2*s3
1038 vpmuludq $H2,$T4,$H2 # h2*s4
1039 vpmuludq $H3,$T4,$H3 # h3*s4
1040 vpsrldq \$6,$T0,$T2 # splat input
1041 vpaddq $H2,$D1,$D1 # d1 += h2*s4
1042 vpmuludq $H4,$T4,$H4 # h4*s4
1043 vpsrldq \$6,$T1,$T3 #
1044 vpaddq $H3,$D2,$H2 # h2 = d2 + h3*s4
1045 vpaddq $H4,$D3,$H3 # h3 = d3 + h4*s4
1046 vpmuludq -0x20(%r11),$H0,$H4 # h0*r4
1047 vpmuludq $H1,$T4,$H0
1048 vpunpckhqdq $T1,$T0,$T4 # 4
1049 vpaddq $H4,$D4,$H4 # h4 = d4 + h0*r4
1050 vpaddq $H0,$D0,$H0 # h0 = d0 + h1*s4
1052 vpunpcklqdq $T1,$T0,$T0 # 0:1
1053 vpunpcklqdq $T3,$T2,$T3 # 2:3
1055 #vpsrlq \$40,$T4,$T4 # 4
1056 vpsrldq \$`40/8`,$T4,$T4 # 4
1058 vmovdqa 0x00(%rsp),$D4 # preload r0^2
1059 vpand $MASK,$T0,$T0 # 0
1061 vpand $MASK,$T1,$T1 # 1
1062 vpand 0(%rcx),$T4,$T4 # .Lmask24
1064 vpand $MASK,$T2,$T2 # 2
1065 vpand $MASK,$T3,$T3 # 3
1066 vpor 32(%rcx),$T4,$T4 # padbit, yes, always
1068 ################################################################
1069 # lazy reduction as discussed in "NEON crypto" by D.J. Bernstein
1074 vpaddq $D3,$H4,$H4 # h3 -> h4
1078 vpaddq $D0,$D1,$H1 # h0 -> h1
1085 vpaddq $D1,$H2,$H2 # h1 -> h2
1089 vpaddq $D0,$H0,$H0 # h4 -> h0
1093 vpaddq $D2,$H3,$H3 # h2 -> h3
1097 vpaddq $D0,$H1,$H1 # h0 -> h1
1101 vpaddq $D3,$H4,$H4 # h3 -> h4
1106 ################################################################
1107 # multiply (inp[0:1]+hash) or inp[2:3] by r^2:r^1
1109 vpshufd \$0x10,$D4,$D4 # r0^n, xx12 -> x1x2
1120 vmovdqa $H2,0x20(%r11)
1121 vmovdqa $H0,0x00(%r11)
1122 vmovdqa $H1,0x10(%r11)
1123 vmovdqa $H3,0x30(%r11)
1124 vmovdqa $H4,0x40(%r11)
1126 # d4 = h4*r0 + h3*r1 + h2*r2 + h1*r3 + h0*r4
1127 # d3 = h3*r0 + h2*r1 + h1*r2 + h0*r3 + h4*5*r4
1128 # d2 = h2*r0 + h1*r1 + h0*r2 + h4*5*r3 + h3*5*r4
1129 # d1 = h1*r0 + h0*r1 + h4*5*r2 + h3*5*r3 + h2*5*r4
1130 # d0 = h0*r0 + h4*5*r1 + h3*5*r2 + h2*5*r3 + h1*5*r4
1132 vpmuludq $T2,$D4,$D2 # d2 = h2*r0
1133 vpmuludq $T0,$D4,$D0 # d0 = h0*r0
1134 vpshufd \$0x10,`16*1-64`($ctx),$H2 # r1^n
1135 vpmuludq $T1,$D4,$D1 # d1 = h1*r0
1136 vpmuludq $T3,$D4,$D3 # d3 = h3*r0
1137 vpmuludq $T4,$D4,$D4 # d4 = h4*r0
1139 vpmuludq $T3,$H2,$H0 # h3*r1
1140 vpaddq $H0,$D4,$D4 # d4 += h3*r1
1141 vpshufd \$0x10,`16*2-64`($ctx),$H3 # s1^n
1142 vpmuludq $T2,$H2,$H1 # h2*r1
1143 vpaddq $H1,$D3,$D3 # d3 += h2*r1
1144 vpshufd \$0x10,`16*3-64`($ctx),$H4 # r2^n
1145 vpmuludq $T1,$H2,$H0 # h1*r1
1146 vpaddq $H0,$D2,$D2 # d2 += h1*r1
1147 vpmuludq $T0,$H2,$H2 # h0*r1
1148 vpaddq $H2,$D1,$D1 # d1 += h0*r1
1149 vpmuludq $T4,$H3,$H3 # h4*s1
1150 vpaddq $H3,$D0,$D0 # d0 += h4*s1
1152 vpshufd \$0x10,`16*4-64`($ctx),$H2 # s2^n
1153 vpmuludq $T2,$H4,$H1 # h2*r2
1154 vpaddq $H1,$D4,$D4 # d4 += h2*r2
1155 vpmuludq $T1,$H4,$H0 # h1*r2
1156 vpaddq $H0,$D3,$D3 # d3 += h1*r2
1157 vpshufd \$0x10,`16*5-64`($ctx),$H3 # r3^n
1158 vpmuludq $T0,$H4,$H4 # h0*r2
1159 vpaddq $H4,$D2,$D2 # d2 += h0*r2
1160 vpmuludq $T4,$H2,$H1 # h4*s2
1161 vpaddq $H1,$D1,$D1 # d1 += h4*s2
1162 vpshufd \$0x10,`16*6-64`($ctx),$H4 # s3^n
1163 vpmuludq $T3,$H2,$H2 # h3*s2
1164 vpaddq $H2,$D0,$D0 # d0 += h3*s2
1166 vpmuludq $T1,$H3,$H0 # h1*r3
1167 vpaddq $H0,$D4,$D4 # d4 += h1*r3
1168 vpmuludq $T0,$H3,$H3 # h0*r3
1169 vpaddq $H3,$D3,$D3 # d3 += h0*r3
1170 vpshufd \$0x10,`16*7-64`($ctx),$H2 # r4^n
1171 vpmuludq $T4,$H4,$H1 # h4*s3
1172 vpaddq $H1,$D2,$D2 # d2 += h4*s3
1173 vpshufd \$0x10,`16*8-64`($ctx),$H3 # s4^n
1174 vpmuludq $T3,$H4,$H0 # h3*s3
1175 vpaddq $H0,$D1,$D1 # d1 += h3*s3
1176 vpmuludq $T2,$H4,$H4 # h2*s3
1177 vpaddq $H4,$D0,$D0 # d0 += h2*s3
1179 vpmuludq $T0,$H2,$H2 # h0*r4
1180 vpaddq $H2,$D4,$D4 # h4 = d4 + h0*r4
1181 vpmuludq $T4,$H3,$H1 # h4*s4
1182 vpaddq $H1,$D3,$D3 # h3 = d3 + h4*s4
1183 vpmuludq $T3,$H3,$H0 # h3*s4
1184 vpaddq $H0,$D2,$D2 # h2 = d2 + h3*s4
1185 vpmuludq $T2,$H3,$H1 # h2*s4
1186 vpaddq $H1,$D1,$D1 # h1 = d1 + h2*s4
1187 vpmuludq $T1,$H3,$H3 # h1*s4
1188 vpaddq $H3,$D0,$D0 # h0 = d0 + h1*s4
1192 vmovdqu 16*0($inp),$H0 # load input
1193 vmovdqu 16*1($inp),$H1
1195 vpsrldq \$6,$H0,$H2 # splat input
1197 vpunpckhqdq $H1,$H0,$H4 # 4
1198 vpunpcklqdq $H1,$H0,$H0 # 0:1
1199 vpunpcklqdq $H3,$H2,$H3 # 2:3
1201 vpsrlq \$40,$H4,$H4 # 4
1203 vpand $MASK,$H0,$H0 # 0
1205 vpand $MASK,$H1,$H1 # 1
1207 vpand $MASK,$H2,$H2 # 2
1208 vpand $MASK,$H3,$H3 # 3
1209 vpor 32(%rcx),$H4,$H4 # padbit, yes, always
1211 vpshufd \$0x32,`16*0-64`($ctx),$T4 # r0^n, 34xx -> x3x4
1212 vpaddq 0x00(%r11),$H0,$H0
1213 vpaddq 0x10(%r11),$H1,$H1
1214 vpaddq 0x20(%r11),$H2,$H2
1215 vpaddq 0x30(%r11),$H3,$H3
1216 vpaddq 0x40(%r11),$H4,$H4
1218 ################################################################
1219 # multiply (inp[0:1]+hash) by r^4:r^3 and accumulate
1221 vpmuludq $H0,$T4,$T0 # h0*r0
1222 vpaddq $T0,$D0,$D0 # d0 += h0*r0
1223 vpmuludq $H1,$T4,$T1 # h1*r0
1224 vpaddq $T1,$D1,$D1 # d1 += h1*r0
1225 vpmuludq $H2,$T4,$T0 # h2*r0
1226 vpaddq $T0,$D2,$D2 # d2 += h2*r0
1227 vpshufd \$0x32,`16*1-64`($ctx),$T2 # r1^n
1228 vpmuludq $H3,$T4,$T1 # h3*r0
1229 vpaddq $T1,$D3,$D3 # d3 += h3*r0
1230 vpmuludq $H4,$T4,$T4 # h4*r0
1231 vpaddq $T4,$D4,$D4 # d4 += h4*r0
1233 vpmuludq $H3,$T2,$T0 # h3*r1
1234 vpaddq $T0,$D4,$D4 # d4 += h3*r1
1235 vpshufd \$0x32,`16*2-64`($ctx),$T3 # s1
1236 vpmuludq $H2,$T2,$T1 # h2*r1
1237 vpaddq $T1,$D3,$D3 # d3 += h2*r1
1238 vpshufd \$0x32,`16*3-64`($ctx),$T4 # r2
1239 vpmuludq $H1,$T2,$T0 # h1*r1
1240 vpaddq $T0,$D2,$D2 # d2 += h1*r1
1241 vpmuludq $H0,$T2,$T2 # h0*r1
1242 vpaddq $T2,$D1,$D1 # d1 += h0*r1
1243 vpmuludq $H4,$T3,$T3 # h4*s1
1244 vpaddq $T3,$D0,$D0 # d0 += h4*s1
1246 vpshufd \$0x32,`16*4-64`($ctx),$T2 # s2
1247 vpmuludq $H2,$T4,$T1 # h2*r2
1248 vpaddq $T1,$D4,$D4 # d4 += h2*r2
1249 vpmuludq $H1,$T4,$T0 # h1*r2
1250 vpaddq $T0,$D3,$D3 # d3 += h1*r2
1251 vpshufd \$0x32,`16*5-64`($ctx),$T3 # r3
1252 vpmuludq $H0,$T4,$T4 # h0*r2
1253 vpaddq $T4,$D2,$D2 # d2 += h0*r2
1254 vpmuludq $H4,$T2,$T1 # h4*s2
1255 vpaddq $T1,$D1,$D1 # d1 += h4*s2
1256 vpshufd \$0x32,`16*6-64`($ctx),$T4 # s3
1257 vpmuludq $H3,$T2,$T2 # h3*s2
1258 vpaddq $T2,$D0,$D0 # d0 += h3*s2
1260 vpmuludq $H1,$T3,$T0 # h1*r3
1261 vpaddq $T0,$D4,$D4 # d4 += h1*r3
1262 vpmuludq $H0,$T3,$T3 # h0*r3
1263 vpaddq $T3,$D3,$D3 # d3 += h0*r3
1264 vpshufd \$0x32,`16*7-64`($ctx),$T2 # r4
1265 vpmuludq $H4,$T4,$T1 # h4*s3
1266 vpaddq $T1,$D2,$D2 # d2 += h4*s3
1267 vpshufd \$0x32,`16*8-64`($ctx),$T3 # s4
1268 vpmuludq $H3,$T4,$T0 # h3*s3
1269 vpaddq $T0,$D1,$D1 # d1 += h3*s3
1270 vpmuludq $H2,$T4,$T4 # h2*s3
1271 vpaddq $T4,$D0,$D0 # d0 += h2*s3
1273 vpmuludq $H0,$T2,$T2 # h0*r4
1274 vpaddq $T2,$D4,$D4 # d4 += h0*r4
1275 vpmuludq $H4,$T3,$T1 # h4*s4
1276 vpaddq $T1,$D3,$D3 # d3 += h4*s4
1277 vpmuludq $H3,$T3,$T0 # h3*s4
1278 vpaddq $T0,$D2,$D2 # d2 += h3*s4
1279 vpmuludq $H2,$T3,$T1 # h2*s4
1280 vpaddq $T1,$D1,$D1 # d1 += h2*s4
1281 vpmuludq $H1,$T3,$T3 # h1*s4
1282 vpaddq $T3,$D0,$D0 # d0 += h1*s4
1285 ################################################################
1286 # horizontal addition
1299 ################################################################
1304 vpaddq $H3,$D4,$D4 # h3 -> h4
1308 vpaddq $H0,$D1,$D1 # h0 -> h1
1315 vpaddq $H1,$D2,$D2 # h1 -> h2
1319 vpaddq $H4,$D0,$D0 # h4 -> h0
1323 vpaddq $H2,$D3,$D3 # h2 -> h3
1327 vpaddq $H0,$D1,$D1 # h0 -> h1
1331 vpaddq $H3,$D4,$D4 # h3 -> h4
1333 vmovd $D0,`4*0-48-64`($ctx) # save partially reduced
1334 vmovd $D1,`4*1-48-64`($ctx)
1335 vmovd $D2,`4*2-48-64`($ctx)
1336 vmovd $D3,`4*3-48-64`($ctx)
1337 vmovd $D4,`4*4-48-64`($ctx)
1339 $code.=<<___ if ($win64);
1340 vmovdqa 0x50(%r11),%xmm6
1341 vmovdqa 0x60(%r11),%xmm7
1342 vmovdqa 0x70(%r11),%xmm8
1343 vmovdqa 0x80(%r11),%xmm9
1344 vmovdqa 0x90(%r11),%xmm10
1345 vmovdqa 0xa0(%r11),%xmm11
1346 vmovdqa 0xb0(%r11),%xmm12
1347 vmovdqa 0xc0(%r11),%xmm13
1348 vmovdqa 0xd0(%r11),%xmm14
1349 vmovdqa 0xe0(%r11),%xmm15
1353 $code.=<<___ if (!$win64);
1361 .size poly1305_blocks_avx,.-poly1305_blocks_avx
1363 .type poly1305_emit_avx,\@function,3
1366 cmpl \$0,20($ctx) # is_base2_26?
1369 mov 0($ctx),%eax # load hash value base 2^26
1375 shl \$26,%rcx # base 2^26 -> base 2^64
1391 mov %r10,%rax # could be partially reduced, so reduce
1402 add \$5,%r8 # compare to modulus
1406 shr \$2,%r10 # did 130-bit value overfow?
1410 add 0($nonce),%rax # accumulate nonce
1412 mov %rax,0($mac) # write result
1416 .size poly1305_emit_avx,.-poly1305_emit_avx
1420 my ($H0,$H1,$H2,$H3,$H4, $MASK, $T4,$T0,$T1,$T2,$T3, $D0,$D1,$D2,$D3,$D4) =
1421 map("%ymm$_",(0..15));
1425 .type poly1305_blocks_avx2,\@function,4
1427 poly1305_blocks_avx2:
1429 mov 20($ctx),%r8d # is_base2_26
1461 mov $len,%r15 # reassign $len
1463 mov 0($ctx),$d1 # load hash value
1467 mov 24($ctx),$r0 # load r
1470 ################################# base 2^26 -> base 2^64
1472 and \$`-1*(1<<31)`,$d1
1473 mov $d2,$r1 # borrow $r1
1475 and \$`-1*(1<<31)`,$d2
1489 adc \$0,$h2 # can be partially reduced...
1491 mov \$-4,$d2 # ... so reduce
1504 add $r1,$s1 # s1 = r1 + (r1 >> 2)
1506 .Lbase2_26_pre_avx2:
1507 add 0($inp),$h0 # accumulate input
1513 call __poly1305_block
1517 jnz .Lbase2_26_pre_avx2
1519 test $padbit,$padbit # if $padbit is zero,
1520 jz .Lstore_base2_64_avx2 # store hash in base 2^64 format
1522 ################################# base 2^64 -> base 2^26
1529 and \$0x3ffffff,%rax # h[0]
1531 and \$0x3ffffff,%rdx # h[1]
1535 and \$0x3ffffff,$h0 # h[2]
1537 and \$0x3ffffff,$h1 # h[3]
1541 jz .Lstore_base2_26_avx2
1551 .Lstore_base2_64_avx2:
1554 mov $h2,16($ctx) # note that is_base2_26 is zeroed
1558 .Lstore_base2_26_avx2:
1559 mov %rax#d,0($ctx) # store hash value base 2^26
1579 .cfi_adjust_cfa_offset -48
1581 .Lblocks_avx2_epilogue:
1600 .Lbase2_64_avx2_body:
1602 mov $len,%r15 # reassign $len
1604 mov 24($ctx),$r0 # load r
1607 mov 0($ctx),$h0 # load hash value
1614 add $r1,$s1 # s1 = r1 + (r1 >> 2)
1619 .Lbase2_64_pre_avx2:
1620 add 0($inp),$h0 # accumulate input
1626 call __poly1305_block
1630 jnz .Lbase2_64_pre_avx2
1633 ################################# base 2^64 -> base 2^26
1640 and \$0x3ffffff,%rax # h[0]
1642 and \$0x3ffffff,%rdx # h[1]
1646 and \$0x3ffffff,$h0 # h[2]
1648 and \$0x3ffffff,$h1 # h[3]
1656 movl \$1,20($ctx) # set is_base2_26
1658 call __poly1305_init_avx
1661 mov %r15,$len # restore $len
1662 mov OPENSSL_ia32cap_P+8(%rip),%r10d
1663 mov \$`(1<<31|1<<30|1<<16)`,%r11d
1679 .cfi_adjust_cfa_offset -48
1680 .Lbase2_64_avx2_epilogue:
1687 mov OPENSSL_ia32cap_P+8(%rip),%r10d
1688 mov \$`(1<<31|1<<30|1<<16)`,%r11d
1689 vmovd 4*0($ctx),%x#$H0 # load hash value base 2^26
1690 vmovd 4*1($ctx),%x#$H1
1691 vmovd 4*2($ctx),%x#$H2
1692 vmovd 4*3($ctx),%x#$H3
1693 vmovd 4*4($ctx),%x#$H4
1697 $code.=<<___ if ($avx>2);
1701 cmp %r11d,%r10d # check for AVX512F+BW+VL
1705 $code.=<<___ if (!$win64);
1707 .cfi_def_cfa %r11,16
1710 $code.=<<___ if ($win64);
1711 lea -0xf8(%rsp),%r11
1713 vmovdqa %xmm6,0x50(%r11)
1714 vmovdqa %xmm7,0x60(%r11)
1715 vmovdqa %xmm8,0x70(%r11)
1716 vmovdqa %xmm9,0x80(%r11)
1717 vmovdqa %xmm10,0x90(%r11)
1718 vmovdqa %xmm11,0xa0(%r11)
1719 vmovdqa %xmm12,0xb0(%r11)
1720 vmovdqa %xmm13,0xc0(%r11)
1721 vmovdqa %xmm14,0xd0(%r11)
1722 vmovdqa %xmm15,0xe0(%r11)
1726 lea .Lconst(%rip),%rcx
1727 lea 48+64($ctx),$ctx # size optimization
1728 vmovdqa 96(%rcx),$T0 # .Lpermd_avx2
1730 # expand and copy pre-calculated table to stack
1731 vmovdqu `16*0-64`($ctx),%x#$T2
1733 vmovdqu `16*1-64`($ctx),%x#$T3
1734 vmovdqu `16*2-64`($ctx),%x#$T4
1735 vmovdqu `16*3-64`($ctx),%x#$D0
1736 vmovdqu `16*4-64`($ctx),%x#$D1
1737 vmovdqu `16*5-64`($ctx),%x#$D2
1738 lea 0x90(%rsp),%rax # size optimization
1739 vmovdqu `16*6-64`($ctx),%x#$D3
1740 vpermd $T2,$T0,$T2 # 00003412 -> 14243444
1741 vmovdqu `16*7-64`($ctx),%x#$D4
1743 vmovdqu `16*8-64`($ctx),%x#$MASK
1745 vmovdqa $T2,0x00(%rsp)
1747 vmovdqa $T3,0x20-0x90(%rax)
1749 vmovdqa $T4,0x40-0x90(%rax)
1751 vmovdqa $D0,0x60-0x90(%rax)
1753 vmovdqa $D1,0x80-0x90(%rax)
1755 vmovdqa $D2,0xa0-0x90(%rax)
1756 vpermd $MASK,$T0,$MASK
1757 vmovdqa $D3,0xc0-0x90(%rax)
1758 vmovdqa $D4,0xe0-0x90(%rax)
1759 vmovdqa $MASK,0x100-0x90(%rax)
1760 vmovdqa 64(%rcx),$MASK # .Lmask26
1762 ################################################################
1764 vmovdqu 16*0($inp),%x#$T0
1765 vmovdqu 16*1($inp),%x#$T1
1766 vinserti128 \$1,16*2($inp),$T0,$T0
1767 vinserti128 \$1,16*3($inp),$T1,$T1
1770 vpsrldq \$6,$T0,$T2 # splat input
1772 vpunpckhqdq $T1,$T0,$T4 # 4
1773 vpunpcklqdq $T3,$T2,$T2 # 2:3
1774 vpunpcklqdq $T1,$T0,$T0 # 0:1
1779 vpsrlq \$40,$T4,$T4 # 4
1780 vpand $MASK,$T2,$T2 # 2
1781 vpand $MASK,$T0,$T0 # 0
1782 vpand $MASK,$T1,$T1 # 1
1783 vpand $MASK,$T3,$T3 # 3
1784 vpor 32(%rcx),$T4,$T4 # padbit, yes, always
1786 vpaddq $H2,$T2,$H2 # accumulate input
1793 ################################################################
1794 # ((inp[0]*r^4+inp[4])*r^4+inp[ 8])*r^4
1795 # ((inp[1]*r^4+inp[5])*r^4+inp[ 9])*r^3
1796 # ((inp[2]*r^4+inp[6])*r^4+inp[10])*r^2
1797 # ((inp[3]*r^4+inp[7])*r^4+inp[11])*r^1
1798 # \________/\__________/
1799 ################################################################
1800 #vpaddq $H2,$T2,$H2 # accumulate input
1802 vmovdqa `32*0`(%rsp),$T0 # r0^4
1804 vmovdqa `32*1`(%rsp),$T1 # r1^4
1806 vmovdqa `32*3`(%rsp),$T2 # r2^4
1808 vmovdqa `32*6-0x90`(%rax),$T3 # s3^4
1809 vmovdqa `32*8-0x90`(%rax),$S4 # s4^4
1811 # d4 = h4*r0 + h3*r1 + h2*r2 + h1*r3 + h0*r4
1812 # d3 = h3*r0 + h2*r1 + h1*r2 + h0*r3 + h4*5*r4
1813 # d2 = h2*r0 + h1*r1 + h0*r2 + h4*5*r3 + h3*5*r4
1814 # d1 = h1*r0 + h0*r1 + h4*5*r2 + h3*5*r3 + h2*5*r4
1815 # d0 = h0*r0 + h4*5*r1 + h3*5*r2 + h2*5*r3 + h1*5*r4
1817 # however, as h2 is "chronologically" first one available pull
1818 # corresponding operations up, so it's
1820 # d4 = h2*r2 + h4*r0 + h3*r1 + h1*r3 + h0*r4
1821 # d3 = h2*r1 + h3*r0 + h1*r2 + h0*r3 + h4*5*r4
1822 # d2 = h2*r0 + h1*r1 + h0*r2 + h4*5*r3 + h3*5*r4
1823 # d1 = h2*5*r4 + h1*r0 + h0*r1 + h4*5*r2 + h3*5*r3
1824 # d0 = h2*5*r3 + h0*r0 + h4*5*r1 + h3*5*r2 + h1*5*r4
1826 vpmuludq $H2,$T0,$D2 # d2 = h2*r0
1827 vpmuludq $H2,$T1,$D3 # d3 = h2*r1
1828 vpmuludq $H2,$T2,$D4 # d4 = h2*r2
1829 vpmuludq $H2,$T3,$D0 # d0 = h2*s3
1830 vpmuludq $H2,$S4,$D1 # d1 = h2*s4
1832 vpmuludq $H0,$T1,$T4 # h0*r1
1833 vpmuludq $H1,$T1,$H2 # h1*r1, borrow $H2 as temp
1834 vpaddq $T4,$D1,$D1 # d1 += h0*r1
1835 vpaddq $H2,$D2,$D2 # d2 += h1*r1
1836 vpmuludq $H3,$T1,$T4 # h3*r1
1837 vpmuludq `32*2`(%rsp),$H4,$H2 # h4*s1
1838 vpaddq $T4,$D4,$D4 # d4 += h3*r1
1839 vpaddq $H2,$D0,$D0 # d0 += h4*s1
1840 vmovdqa `32*4-0x90`(%rax),$T1 # s2
1842 vpmuludq $H0,$T0,$T4 # h0*r0
1843 vpmuludq $H1,$T0,$H2 # h1*r0
1844 vpaddq $T4,$D0,$D0 # d0 += h0*r0
1845 vpaddq $H2,$D1,$D1 # d1 += h1*r0
1846 vpmuludq $H3,$T0,$T4 # h3*r0
1847 vpmuludq $H4,$T0,$H2 # h4*r0
1848 vmovdqu 16*0($inp),%x#$T0 # load input
1849 vpaddq $T4,$D3,$D3 # d3 += h3*r0
1850 vpaddq $H2,$D4,$D4 # d4 += h4*r0
1851 vinserti128 \$1,16*2($inp),$T0,$T0
1853 vpmuludq $H3,$T1,$T4 # h3*s2
1854 vpmuludq $H4,$T1,$H2 # h4*s2
1855 vmovdqu 16*1($inp),%x#$T1
1856 vpaddq $T4,$D0,$D0 # d0 += h3*s2
1857 vpaddq $H2,$D1,$D1 # d1 += h4*s2
1858 vmovdqa `32*5-0x90`(%rax),$H2 # r3
1859 vpmuludq $H1,$T2,$T4 # h1*r2
1860 vpmuludq $H0,$T2,$T2 # h0*r2
1861 vpaddq $T4,$D3,$D3 # d3 += h1*r2
1862 vpaddq $T2,$D2,$D2 # d2 += h0*r2
1863 vinserti128 \$1,16*3($inp),$T1,$T1
1866 vpmuludq $H1,$H2,$T4 # h1*r3
1867 vpmuludq $H0,$H2,$H2 # h0*r3
1868 vpsrldq \$6,$T0,$T2 # splat input
1869 vpaddq $T4,$D4,$D4 # d4 += h1*r3
1870 vpaddq $H2,$D3,$D3 # d3 += h0*r3
1871 vpmuludq $H3,$T3,$T4 # h3*s3
1872 vpmuludq $H4,$T3,$H2 # h4*s3
1874 vpaddq $T4,$D1,$D1 # d1 += h3*s3
1875 vpaddq $H2,$D2,$D2 # d2 += h4*s3
1876 vpunpckhqdq $T1,$T0,$T4 # 4
1878 vpmuludq $H3,$S4,$H3 # h3*s4
1879 vpmuludq $H4,$S4,$H4 # h4*s4
1880 vpunpcklqdq $T1,$T0,$T0 # 0:1
1881 vpaddq $H3,$D2,$H2 # h2 = d2 + h3*r4
1882 vpaddq $H4,$D3,$H3 # h3 = d3 + h4*r4
1883 vpunpcklqdq $T3,$T2,$T3 # 2:3
1884 vpmuludq `32*7-0x90`(%rax),$H0,$H4 # h0*r4
1885 vpmuludq $H1,$S4,$H0 # h1*s4
1886 vmovdqa 64(%rcx),$MASK # .Lmask26
1887 vpaddq $H4,$D4,$H4 # h4 = d4 + h0*r4
1888 vpaddq $H0,$D0,$H0 # h0 = d0 + h1*s4
1890 ################################################################
1891 # lazy reduction (interleaved with tail of input splat)
1895 vpaddq $D3,$H4,$H4 # h3 -> h4
1899 vpaddq $D0,$D1,$H1 # h0 -> h1
1908 vpaddq $D1,$H2,$H2 # h1 -> h2
1912 vpaddq $D4,$H0,$H0 # h4 -> h0
1914 vpand $MASK,$T2,$T2 # 2
1919 vpaddq $D2,$H3,$H3 # h2 -> h3
1921 vpaddq $T2,$H2,$H2 # modulo-scheduled
1926 vpaddq $D0,$H1,$H1 # h0 -> h1
1928 vpsrlq \$40,$T4,$T4 # 4
1932 vpaddq $D3,$H4,$H4 # h3 -> h4
1934 vpand $MASK,$T0,$T0 # 0
1935 vpand $MASK,$T1,$T1 # 1
1936 vpand $MASK,$T3,$T3 # 3
1937 vpor 32(%rcx),$T4,$T4 # padbit, yes, always
1944 ################################################################
1945 # while above multiplications were by r^4 in all lanes, in last
1946 # iteration we multiply least significant lane by r^4 and most
1947 # significant one by r, so copy of above except that references
1948 # to the precomputed table are displaced by 4...
1950 #vpaddq $H2,$T2,$H2 # accumulate input
1952 vmovdqu `32*0+4`(%rsp),$T0 # r0^4
1954 vmovdqu `32*1+4`(%rsp),$T1 # r1^4
1956 vmovdqu `32*3+4`(%rsp),$T2 # r2^4
1958 vmovdqu `32*6+4-0x90`(%rax),$T3 # s3^4
1959 vmovdqu `32*8+4-0x90`(%rax),$S4 # s4^4
1961 vpmuludq $H2,$T0,$D2 # d2 = h2*r0
1962 vpmuludq $H2,$T1,$D3 # d3 = h2*r1
1963 vpmuludq $H2,$T2,$D4 # d4 = h2*r2
1964 vpmuludq $H2,$T3,$D0 # d0 = h2*s3
1965 vpmuludq $H2,$S4,$D1 # d1 = h2*s4
1967 vpmuludq $H0,$T1,$T4 # h0*r1
1968 vpmuludq $H1,$T1,$H2 # h1*r1
1969 vpaddq $T4,$D1,$D1 # d1 += h0*r1
1970 vpaddq $H2,$D2,$D2 # d2 += h1*r1
1971 vpmuludq $H3,$T1,$T4 # h3*r1
1972 vpmuludq `32*2+4`(%rsp),$H4,$H2 # h4*s1
1973 vpaddq $T4,$D4,$D4 # d4 += h3*r1
1974 vpaddq $H2,$D0,$D0 # d0 += h4*s1
1976 vpmuludq $H0,$T0,$T4 # h0*r0
1977 vpmuludq $H1,$T0,$H2 # h1*r0
1978 vpaddq $T4,$D0,$D0 # d0 += h0*r0
1979 vmovdqu `32*4+4-0x90`(%rax),$T1 # s2
1980 vpaddq $H2,$D1,$D1 # d1 += h1*r0
1981 vpmuludq $H3,$T0,$T4 # h3*r0
1982 vpmuludq $H4,$T0,$H2 # h4*r0
1983 vpaddq $T4,$D3,$D3 # d3 += h3*r0
1984 vpaddq $H2,$D4,$D4 # d4 += h4*r0
1986 vpmuludq $H3,$T1,$T4 # h3*s2
1987 vpmuludq $H4,$T1,$H2 # h4*s2
1988 vpaddq $T4,$D0,$D0 # d0 += h3*s2
1989 vpaddq $H2,$D1,$D1 # d1 += h4*s2
1990 vmovdqu `32*5+4-0x90`(%rax),$H2 # r3
1991 vpmuludq $H1,$T2,$T4 # h1*r2
1992 vpmuludq $H0,$T2,$T2 # h0*r2
1993 vpaddq $T4,$D3,$D3 # d3 += h1*r2
1994 vpaddq $T2,$D2,$D2 # d2 += h0*r2
1996 vpmuludq $H1,$H2,$T4 # h1*r3
1997 vpmuludq $H0,$H2,$H2 # h0*r3
1998 vpaddq $T4,$D4,$D4 # d4 += h1*r3
1999 vpaddq $H2,$D3,$D3 # d3 += h0*r3
2000 vpmuludq $H3,$T3,$T4 # h3*s3
2001 vpmuludq $H4,$T3,$H2 # h4*s3
2002 vpaddq $T4,$D1,$D1 # d1 += h3*s3
2003 vpaddq $H2,$D2,$D2 # d2 += h4*s3
2005 vpmuludq $H3,$S4,$H3 # h3*s4
2006 vpmuludq $H4,$S4,$H4 # h4*s4
2007 vpaddq $H3,$D2,$H2 # h2 = d2 + h3*r4
2008 vpaddq $H4,$D3,$H3 # h3 = d3 + h4*r4
2009 vpmuludq `32*7+4-0x90`(%rax),$H0,$H4 # h0*r4
2010 vpmuludq $H1,$S4,$H0 # h1*s4
2011 vmovdqa 64(%rcx),$MASK # .Lmask26
2012 vpaddq $H4,$D4,$H4 # h4 = d4 + h0*r4
2013 vpaddq $H0,$D0,$H0 # h0 = d0 + h1*s4
2015 ################################################################
2016 # horizontal addition
2029 vpermq \$0x2,$H3,$T3
2030 vpermq \$0x2,$H4,$T4
2031 vpermq \$0x2,$H0,$T0
2032 vpermq \$0x2,$D1,$T1
2033 vpermq \$0x2,$H2,$T2
2040 ################################################################
2045 vpaddq $D3,$H4,$H4 # h3 -> h4
2049 vpaddq $D0,$D1,$H1 # h0 -> h1
2056 vpaddq $D1,$H2,$H2 # h1 -> h2
2060 vpaddq $D4,$H0,$H0 # h4 -> h0
2064 vpaddq $D2,$H3,$H3 # h2 -> h3
2068 vpaddq $D0,$H1,$H1 # h0 -> h1
2072 vpaddq $D3,$H4,$H4 # h3 -> h4
2074 vmovd %x#$H0,`4*0-48-64`($ctx)# save partially reduced
2075 vmovd %x#$H1,`4*1-48-64`($ctx)
2076 vmovd %x#$H2,`4*2-48-64`($ctx)
2077 vmovd %x#$H3,`4*3-48-64`($ctx)
2078 vmovd %x#$H4,`4*4-48-64`($ctx)
2080 $code.=<<___ if ($win64);
2081 vmovdqa 0x50(%r11),%xmm6
2082 vmovdqa 0x60(%r11),%xmm7
2083 vmovdqa 0x70(%r11),%xmm8
2084 vmovdqa 0x80(%r11),%xmm9
2085 vmovdqa 0x90(%r11),%xmm10
2086 vmovdqa 0xa0(%r11),%xmm11
2087 vmovdqa 0xb0(%r11),%xmm12
2088 vmovdqa 0xc0(%r11),%xmm13
2089 vmovdqa 0xd0(%r11),%xmm14
2090 vmovdqa 0xe0(%r11),%xmm15
2094 $code.=<<___ if (!$win64);
2102 .size poly1305_blocks_avx2,.-poly1305_blocks_avx2
2104 #######################################################################
2106 # On entry we have input length divisible by 64. But since inner loop
2107 # processes 128 bytes per iteration, cases when length is not divisible
2108 # by 128 are handled by passing tail 64 bytes to .Ltail_avx2. For this
2109 # reason stack layout is kept identical to poly1305_blocks_avx2. If not
2110 # for this tail, we wouldn't have to even allocate stack frame...
2112 my ($R0,$R1,$R2,$R3,$R4, $S1,$S2,$S3,$S4) = map("%ymm$_",(16..24));
2113 my ($M0,$M1,$M2,$M3,$M4) = map("%ymm$_",(25..29));
2114 my $PADBIT="%zmm30";
2115 my $GATHER="%ymm31";
2118 .type poly1305_blocks_avx512,\@function,4
2120 poly1305_blocks_avx512:
2125 $code.=<<___ if (!$win64);
2127 .cfi_def_cfa %r11,16
2130 $code.=<<___ if ($win64);
2131 lea -0xf8(%rsp),%r11
2133 vmovdqa %xmm6,0x50(%r11)
2134 vmovdqa %xmm7,0x60(%r11)
2135 vmovdqa %xmm8,0x70(%r11)
2136 vmovdqa32 %xmm9,0x80(%r11)
2137 vmovdqa32 %xmm10,0x90(%r11)
2138 vmovdqa32 %xmm11,0xa0(%r11)
2139 vmovdqa32 %xmm12,0xb0(%r11)
2140 vmovdqa32 %xmm13,0xc0(%r11)
2141 vmovdqa32 %xmm14,0xd0(%r11)
2142 vmovdqa32 %xmm15,0xe0(%r11)
2146 lea .Lconst(%rip),%rcx
2147 lea 48+64($ctx),$ctx # size optimization
2148 vmovdqa 96(%rcx),$T2 # .Lpermd_avx2
2150 # expand pre-calculated table
2151 vmovdqu32 `16*0-64`($ctx),%x#$R0
2153 vmovdqu32 `16*1-64`($ctx),%x#$R1
2154 vmovdqu32 `16*2-64`($ctx),%x#$S1
2155 vmovdqu32 `16*3-64`($ctx),%x#$R2
2156 vmovdqu32 `16*4-64`($ctx),%x#$S2
2157 vmovdqu32 `16*5-64`($ctx),%x#$R3
2158 vmovdqu32 `16*6-64`($ctx),%x#$S3
2159 vmovdqu32 `16*7-64`($ctx),%x#$R4
2160 vmovdqu32 `16*8-64`($ctx),%x#$S4
2161 vpermd $R0,$T2,$R0 # 00003412 -> 14243444
2162 vmovdqa64 64(%rcx),$MASK # .Lmask26
2166 vmovdqa32 $R0,0x00(%rsp) # save in case $len%128 != 0
2167 vpsrlq \$32,$R0,$T0 # 14243444 -> 01020304
2169 vmovdqa32 $R1,0x20(%rsp)
2172 vmovdqa32 $S1,0x40(%rsp)
2175 vmovdqa32 $R2,0x60(%rsp)
2177 vmovdqa32 $S2,0x80(%rsp)
2178 vmovdqa32 $R3,0xa0(%rsp)
2179 vmovdqa32 $S3,0xc0(%rsp)
2180 vmovdqa32 $R4,0xe0(%rsp)
2181 vmovdqa32 $S4,0x100(%rsp)
2183 ################################################################
2184 # calculate 5th through 8th powers of the key
2186 # d0 = r0'*r0 + r1'*5*r4 + r2'*5*r3 + r3'*5*r2 + r4'*5*r1
2187 # d1 = r0'*r1 + r1'*r0 + r2'*5*r4 + r3'*5*r3 + r4'*5*r2
2188 # d2 = r0'*r2 + r1'*r1 + r2'*r0 + r3'*5*r4 + r4'*5*r3
2189 # d3 = r0'*r3 + r1'*r2 + r2'*r1 + r3'*r0 + r4'*5*r4
2190 # d4 = r0'*r4 + r1'*r3 + r2'*r2 + r3'*r1 + r4'*r0
2192 vpmuludq $T0,$R0,$D0 # d0 = r0'*r0
2193 vpmuludq $T0,$R1,$D1 # d1 = r0'*r1
2194 vpmuludq $T0,$R2,$D2 # d2 = r0'*r2
2195 vpmuludq $T0,$R3,$D3 # d3 = r0'*r3
2196 vpmuludq $T0,$R4,$D4 # d4 = r0'*r4
2199 vpmuludq $T1,$S4,$M0
2200 vpmuludq $T1,$R0,$M1
2201 vpmuludq $T1,$R1,$M2
2202 vpmuludq $T1,$R2,$M3
2203 vpmuludq $T1,$R3,$M4
2205 vpaddq $M0,$D0,$D0 # d0 += r1'*5*r4
2206 vpaddq $M1,$D1,$D1 # d1 += r1'*r0
2207 vpaddq $M2,$D2,$D2 # d2 += r1'*r1
2208 vpaddq $M3,$D3,$D3 # d3 += r1'*r2
2209 vpaddq $M4,$D4,$D4 # d4 += r1'*r3
2211 vpmuludq $T2,$S3,$M0
2212 vpmuludq $T2,$S4,$M1
2213 vpmuludq $T2,$R1,$M3
2214 vpmuludq $T2,$R2,$M4
2215 vpmuludq $T2,$R0,$M2
2217 vpaddq $M0,$D0,$D0 # d0 += r2'*5*r3
2218 vpaddq $M1,$D1,$D1 # d1 += r2'*5*r4
2219 vpaddq $M3,$D3,$D3 # d3 += r2'*r1
2220 vpaddq $M4,$D4,$D4 # d4 += r2'*r2
2221 vpaddq $M2,$D2,$D2 # d2 += r2'*r0
2223 vpmuludq $T3,$S2,$M0
2224 vpmuludq $T3,$R0,$M3
2225 vpmuludq $T3,$R1,$M4
2226 vpmuludq $T3,$S3,$M1
2227 vpmuludq $T3,$S4,$M2
2228 vpaddq $M0,$D0,$D0 # d0 += r3'*5*r2
2229 vpaddq $M3,$D3,$D3 # d3 += r3'*r0
2230 vpaddq $M4,$D4,$D4 # d4 += r3'*r1
2231 vpaddq $M1,$D1,$D1 # d1 += r3'*5*r3
2232 vpaddq $M2,$D2,$D2 # d2 += r3'*5*r4
2234 vpmuludq $T4,$S4,$M3
2235 vpmuludq $T4,$R0,$M4
2236 vpmuludq $T4,$S1,$M0
2237 vpmuludq $T4,$S2,$M1
2238 vpmuludq $T4,$S3,$M2
2239 vpaddq $M3,$D3,$D3 # d3 += r2'*5*r4
2240 vpaddq $M4,$D4,$D4 # d4 += r2'*r0
2241 vpaddq $M0,$D0,$D0 # d0 += r2'*5*r1
2242 vpaddq $M1,$D1,$D1 # d1 += r2'*5*r2
2243 vpaddq $M2,$D2,$D2 # d2 += r2'*5*r3
2245 ################################################################
2247 vmovdqu64 16*0($inp),%z#$T3
2248 vmovdqu64 16*4($inp),%z#$T4
2251 ################################################################
2255 vpandq $MASK,$D3,$D3
2256 vpaddq $M3,$D4,$D4 # d3 -> d4
2259 vpandq $MASK,$D0,$D0
2260 vpaddq $M0,$D1,$D1 # d0 -> d1
2263 vpandq $MASK,$D4,$D4
2266 vpandq $MASK,$D1,$D1
2267 vpaddq $M1,$D2,$D2 # d1 -> d2
2271 vpaddq $M4,$D0,$D0 # d4 -> d0
2274 vpandq $MASK,$D2,$D2
2275 vpaddq $M2,$D3,$D3 # d2 -> d3
2278 vpandq $MASK,$D0,$D0
2279 vpaddq $M0,$D1,$D1 # d0 -> d1
2282 vpandq $MASK,$D3,$D3
2283 vpaddq $M3,$D4,$D4 # d3 -> d4
2286 map(s/%y/%z/,($T4,$T0,$T1,$T2,$T3)); # switch to %zmm domain
2287 map(s/%y/%z/,($M4,$M0,$M1,$M2,$M3));
2288 map(s/%y/%z/,($D0,$D1,$D2,$D3,$D4));
2289 map(s/%y/%z/,($R0,$R1,$R2,$R3,$R4, $S1,$S2,$S3,$S4));
2290 map(s/%y/%z/,($H0,$H1,$H2,$H3,$H4));
2291 map(s/%y/%z/,($MASK));
2293 ################################################################
2294 # at this point we have 14243444 in $R0-$S4 and 05060708 in
2297 vpunpcklqdq $T4,$T3,$T0 # transpose input
2298 vpunpckhqdq $T4,$T3,$T4
2300 # ... since input 64-bit lanes are ordered as 73625140, we could
2301 # "vperm" it to 76543210 (here and in each loop iteration), *or*
2302 # we could just flow along, hence the goal for $R0-$S4 is
2303 # 1858286838784888 ...
2305 vmovdqa32 128(%rcx),$M0 # .Lpermd_avx512:
2309 vpermd $R0,$M0,$R0 # 14243444 -> 1---2---3---4---
2315 vpermd $D0,$M0,${R0}{%k1} # 05060708 -> 1858286838784888
2316 vpermd $D1,$M0,${R1}{%k1}
2317 vpermd $D2,$M0,${R2}{%k1}
2318 vpermd $D3,$M0,${R3}{%k1}
2319 vpermd $D4,$M0,${R4}{%k1}
2321 vpslld \$2,$R1,$S1 # *5
2330 vpbroadcastq %x#$MASK,$MASK
2331 vpbroadcastq 32(%rcx),$PADBIT # .L129
2333 vpsrlq \$52,$T0,$T2 # splat input
2338 vpsrlq \$40,$T4,$T4 # 4
2339 vpandq $MASK,$T2,$T2 # 2
2340 vpandq $MASK,$T0,$T0 # 0
2341 #vpandq $MASK,$T1,$T1 # 1
2342 #vpandq $MASK,$T3,$T3 # 3
2343 #vporq $PADBIT,$T4,$T4 # padbit, yes, always
2345 vpaddq $H2,$T2,$H2 # accumulate input
2352 ################################################################
2353 # ((inp[0]*r^8+inp[ 8])*r^8+inp[16])*r^8
2354 # ((inp[1]*r^8+inp[ 9])*r^8+inp[17])*r^7
2355 # ((inp[2]*r^8+inp[10])*r^8+inp[18])*r^6
2356 # ((inp[3]*r^8+inp[11])*r^8+inp[19])*r^5
2357 # ((inp[4]*r^8+inp[12])*r^8+inp[20])*r^4
2358 # ((inp[5]*r^8+inp[13])*r^8+inp[21])*r^3
2359 # ((inp[6]*r^8+inp[14])*r^8+inp[22])*r^2
2360 # ((inp[7]*r^8+inp[15])*r^8+inp[23])*r^1
2361 # \________/\___________/
2362 ################################################################
2363 #vpaddq $H2,$T2,$H2 # accumulate input
2365 # d4 = h4*r0 + h3*r1 + h2*r2 + h1*r3 + h0*r4
2366 # d3 = h3*r0 + h2*r1 + h1*r2 + h0*r3 + h4*5*r4
2367 # d2 = h2*r0 + h1*r1 + h0*r2 + h4*5*r3 + h3*5*r4
2368 # d1 = h1*r0 + h0*r1 + h4*5*r2 + h3*5*r3 + h2*5*r4
2369 # d0 = h0*r0 + h4*5*r1 + h3*5*r2 + h2*5*r3 + h1*5*r4
2371 # however, as h2 is "chronologically" first one available pull
2372 # corresponding operations up, so it's
2374 # d3 = h2*r1 + h0*r3 + h1*r2 + h3*r0 + h4*5*r4
2375 # d4 = h2*r2 + h0*r4 + h1*r3 + h3*r1 + h4*r0
2376 # d0 = h2*5*r3 + h0*r0 + h1*5*r4 + h3*5*r2 + h4*5*r1
2377 # d1 = h2*5*r4 + h0*r1 + h1*r0 + h3*5*r3 + h4*5*r2
2378 # d2 = h2*r0 + h0*r2 + h1*r1 + h3*5*r4 + h4*5*r3
2380 vpmuludq $H2,$R1,$D3 # d3 = h2*r1
2382 vpmuludq $H2,$R2,$D4 # d4 = h2*r2
2383 vpandq $MASK,$T1,$T1 # 1
2384 vpmuludq $H2,$S3,$D0 # d0 = h2*s3
2385 vpandq $MASK,$T3,$T3 # 3
2386 vpmuludq $H2,$S4,$D1 # d1 = h2*s4
2387 vporq $PADBIT,$T4,$T4 # padbit, yes, always
2388 vpmuludq $H2,$R0,$D2 # d2 = h2*r0
2389 vpaddq $H1,$T1,$H1 # accumulate input
2393 vmovdqu64 16*0($inp),$T3 # load input
2394 vmovdqu64 16*4($inp),$T4
2396 vpmuludq $H0,$R3,$M3
2397 vpmuludq $H0,$R4,$M4
2398 vpmuludq $H0,$R0,$M0
2399 vpmuludq $H0,$R1,$M1
2400 vpaddq $M3,$D3,$D3 # d3 += h0*r3
2401 vpaddq $M4,$D4,$D4 # d4 += h0*r4
2402 vpaddq $M0,$D0,$D0 # d0 += h0*r0
2403 vpaddq $M1,$D1,$D1 # d1 += h0*r1
2405 vpmuludq $H1,$R2,$M3
2406 vpmuludq $H1,$R3,$M4
2407 vpmuludq $H1,$S4,$M0
2408 vpmuludq $H0,$R2,$M2
2409 vpaddq $M3,$D3,$D3 # d3 += h1*r2
2410 vpaddq $M4,$D4,$D4 # d4 += h1*r3
2411 vpaddq $M0,$D0,$D0 # d0 += h1*s4
2412 vpaddq $M2,$D2,$D2 # d2 += h0*r2
2414 vpunpcklqdq $T4,$T3,$T0 # transpose input
2415 vpunpckhqdq $T4,$T3,$T4
2417 vpmuludq $H3,$R0,$M3
2418 vpmuludq $H3,$R1,$M4
2419 vpmuludq $H1,$R0,$M1
2420 vpmuludq $H1,$R1,$M2
2421 vpaddq $M3,$D3,$D3 # d3 += h3*r0
2422 vpaddq $M4,$D4,$D4 # d4 += h3*r1
2423 vpaddq $M1,$D1,$D1 # d1 += h1*r0
2424 vpaddq $M2,$D2,$D2 # d2 += h1*r1
2426 vpmuludq $H4,$S4,$M3
2427 vpmuludq $H4,$R0,$M4
2428 vpmuludq $H3,$S2,$M0
2429 vpmuludq $H3,$S3,$M1
2430 vpaddq $M3,$D3,$D3 # d3 += h4*s4
2431 vpmuludq $H3,$S4,$M2
2432 vpaddq $M4,$D4,$D4 # d4 += h4*r0
2433 vpaddq $M0,$D0,$D0 # d0 += h3*s2
2434 vpaddq $M1,$D1,$D1 # d1 += h3*s3
2435 vpaddq $M2,$D2,$D2 # d2 += h3*s4
2437 vpmuludq $H4,$S1,$M0
2438 vpmuludq $H4,$S2,$M1
2439 vpmuludq $H4,$S3,$M2
2440 vpaddq $M0,$D0,$H0 # h0 = d0 + h4*s1
2441 vpaddq $M1,$D1,$H1 # h1 = d2 + h4*s2
2442 vpaddq $M2,$D2,$H2 # h2 = d3 + h4*s3
2444 ################################################################
2445 # lazy reduction (interleaved with input splat)
2447 vpsrlq \$52,$T0,$T2 # splat input
2451 vpandq $MASK,$D3,$D3
2452 vpaddq $H3,$D4,$H4 # h3 -> h4
2457 vpandq $MASK,$H0,$H0
2458 vpaddq $D0,$H1,$H1 # h0 -> h1
2460 vpandq $MASK,$T2,$T2 # 2
2463 vpandq $MASK,$H4,$H4
2466 vpandq $MASK,$H1,$H1
2467 vpaddq $D1,$H2,$H2 # h1 -> h2
2471 vpaddq $D4,$H0,$H0 # h4 -> h0
2473 vpaddq $T2,$H2,$H2 # modulo-scheduled
2477 vpandq $MASK,$H2,$H2
2478 vpaddq $D2,$D3,$H3 # h2 -> h3
2483 vpandq $MASK,$H0,$H0
2484 vpaddq $D0,$H1,$H1 # h0 -> h1
2486 vpsrlq \$40,$T4,$T4 # 4
2489 vpandq $MASK,$H3,$H3
2490 vpaddq $D3,$H4,$H4 # h3 -> h4
2492 vpandq $MASK,$T0,$T0 # 0
2493 #vpandq $MASK,$T1,$T1 # 1
2494 #vpandq $MASK,$T3,$T3 # 3
2495 #vporq $PADBIT,$T4,$T4 # padbit, yes, always
2501 ################################################################
2502 # while above multiplications were by r^8 in all lanes, in last
2503 # iteration we multiply least significant lane by r^8 and most
2504 # significant one by r, that's why table gets shifted...
2506 vpsrlq \$32,$R0,$R0 # 0105020603070408
2516 ################################################################
2517 # load either next or last 64 byte of input
2518 lea ($inp,$len),$inp
2520 #vpaddq $H2,$T2,$H2 # accumulate input
2523 vpmuludq $H2,$R1,$D3 # d3 = h2*r1
2524 vpmuludq $H2,$R2,$D4 # d4 = h2*r2
2525 vpmuludq $H2,$S3,$D0 # d0 = h2*s3
2526 vpandq $MASK,$T1,$T1 # 1
2527 vpmuludq $H2,$S4,$D1 # d1 = h2*s4
2528 vpandq $MASK,$T3,$T3 # 3
2529 vpmuludq $H2,$R0,$D2 # d2 = h2*r0
2530 vporq $PADBIT,$T4,$T4 # padbit, yes, always
2531 vpaddq $H1,$T1,$H1 # accumulate input
2535 vmovdqu64 16*0($inp),%x#$T0
2536 vpmuludq $H0,$R3,$M3
2537 vpmuludq $H0,$R4,$M4
2538 vpmuludq $H0,$R0,$M0
2539 vpmuludq $H0,$R1,$M1
2540 vpaddq $M3,$D3,$D3 # d3 += h0*r3
2541 vpaddq $M4,$D4,$D4 # d4 += h0*r4
2542 vpaddq $M0,$D0,$D0 # d0 += h0*r0
2543 vpaddq $M1,$D1,$D1 # d1 += h0*r1
2545 vmovdqu64 16*1($inp),%x#$T1
2546 vpmuludq $H1,$R2,$M3
2547 vpmuludq $H1,$R3,$M4
2548 vpmuludq $H1,$S4,$M0
2549 vpmuludq $H0,$R2,$M2
2550 vpaddq $M3,$D3,$D3 # d3 += h1*r2
2551 vpaddq $M4,$D4,$D4 # d4 += h1*r3
2552 vpaddq $M0,$D0,$D0 # d0 += h1*s4
2553 vpaddq $M2,$D2,$D2 # d2 += h0*r2
2555 vinserti64x2 \$1,16*2($inp),$T0,$T0
2556 vpmuludq $H3,$R0,$M3
2557 vpmuludq $H3,$R1,$M4
2558 vpmuludq $H1,$R0,$M1
2559 vpmuludq $H1,$R1,$M2
2560 vpaddq $M3,$D3,$D3 # d3 += h3*r0
2561 vpaddq $M4,$D4,$D4 # d4 += h3*r1
2562 vpaddq $M1,$D1,$D1 # d1 += h1*r0
2563 vpaddq $M2,$D2,$D2 # d2 += h1*r1
2565 vinserti64x2 \$1,16*3($inp),$T1,$T1
2566 vpmuludq $H4,$S4,$M3
2567 vpmuludq $H4,$R0,$M4
2568 vpmuludq $H3,$S2,$M0
2569 vpmuludq $H3,$S3,$M1
2570 vpmuludq $H3,$S4,$M2
2571 vpaddq $M3,$D3,$H3 # h3 = d3 + h4*s4
2572 vpaddq $M4,$D4,$D4 # d4 += h4*r0
2573 vpaddq $M0,$D0,$D0 # d0 += h3*s2
2574 vpaddq $M1,$D1,$D1 # d1 += h3*s3
2575 vpaddq $M2,$D2,$D2 # d2 += h3*s4
2577 vpmuludq $H4,$S1,$M0
2578 vpmuludq $H4,$S2,$M1
2579 vpmuludq $H4,$S3,$M2
2580 vpaddq $M0,$D0,$H0 # h0 = d0 + h4*s1
2581 vpaddq $M1,$D1,$H1 # h1 = d2 + h4*s2
2582 vpaddq $M2,$D2,$H2 # h2 = d3 + h4*s3
2584 ################################################################
2585 # horizontal addition
2600 vpermq \$0x2,$H3,$D3
2601 vpermq \$0x2,$H4,$D4
2602 vpermq \$0x2,$H0,$D0
2603 vpermq \$0x2,$H1,$D1
2604 vpermq \$0x2,$H2,$D2
2611 vextracti64x4 \$0x1,$H3,%y#$D3
2612 vextracti64x4 \$0x1,$H4,%y#$D4
2613 vextracti64x4 \$0x1,$H0,%y#$D0
2614 vextracti64x4 \$0x1,$H1,%y#$D1
2615 vextracti64x4 \$0x1,$H2,%y#$D2
2616 vpaddq $D3,$H3,${H3}{%k3}{z} # keep single qword in case
2617 vpaddq $D4,$H4,${H4}{%k3}{z} # it's passed to .Ltail_avx2
2618 vpaddq $D0,$H0,${H0}{%k3}{z}
2619 vpaddq $D1,$H1,${H1}{%k3}{z}
2620 vpaddq $D2,$H2,${H2}{%k3}{z}
2622 map(s/%z/%y/,($T0,$T1,$T2,$T3,$T4, $PADBIT));
2623 map(s/%z/%y/,($H0,$H1,$H2,$H3,$H4, $D0,$D1,$D2,$D3,$D4, $MASK));
2625 ################################################################
2626 # lazy reduction (interleaved with input splat)
2629 vpandq $MASK,$H3,$H3
2630 vpsrldq \$6,$T0,$T2 # splat input
2632 vpunpckhqdq $T1,$T0,$T4 # 4
2633 vpaddq $D3,$H4,$H4 # h3 -> h4
2636 vpandq $MASK,$H0,$H0
2637 vpunpcklqdq $T3,$T2,$T2 # 2:3
2638 vpunpcklqdq $T1,$T0,$T0 # 0:1
2639 vpaddq $D0,$H1,$H1 # h0 -> h1
2642 vpandq $MASK,$H4,$H4
2645 vpandq $MASK,$H1,$H1
2648 vpaddq $D1,$H2,$H2 # h1 -> h2
2653 vpsrlq \$40,$T4,$T4 # 4
2654 vpaddq $D4,$H0,$H0 # h4 -> h0
2657 vpandq $MASK,$H2,$H2
2658 vpandq $MASK,$T2,$T2 # 2
2659 vpandq $MASK,$T0,$T0 # 0
2660 vpaddq $D2,$H3,$H3 # h2 -> h3
2663 vpandq $MASK,$H0,$H0
2664 vpaddq $H2,$T2,$H2 # accumulate input for .Ltail_avx2
2665 vpandq $MASK,$T1,$T1 # 1
2666 vpaddq $D0,$H1,$H1 # h0 -> h1
2669 vpandq $MASK,$H3,$H3
2670 vpandq $MASK,$T3,$T3 # 3
2671 vporq $PADBIT,$T4,$T4 # padbit, yes, always
2672 vpaddq $D3,$H4,$H4 # h3 -> h4
2674 lea 0x90(%rsp),%rax # size optimization for .Ltail_avx2
2678 vpsubq $T2,$H2,$H2 # undo input accumulation
2679 vmovd %x#$H0,`4*0-48-64`($ctx)# save partially reduced
2680 vmovd %x#$H1,`4*1-48-64`($ctx)
2681 vmovd %x#$H2,`4*2-48-64`($ctx)
2682 vmovd %x#$H3,`4*3-48-64`($ctx)
2683 vmovd %x#$H4,`4*4-48-64`($ctx)
2686 $code.=<<___ if ($win64);
2687 movdqa 0x50(%r11),%xmm6
2688 movdqa 0x60(%r11),%xmm7
2689 movdqa 0x70(%r11),%xmm8
2690 movdqa 0x80(%r11),%xmm9
2691 movdqa 0x90(%r11),%xmm10
2692 movdqa 0xa0(%r11),%xmm11
2693 movdqa 0xb0(%r11),%xmm12
2694 movdqa 0xc0(%r11),%xmm13
2695 movdqa 0xd0(%r11),%xmm14
2696 movdqa 0xe0(%r11),%xmm15
2698 .Ldo_avx512_epilogue:
2700 $code.=<<___ if (!$win64);
2707 .size poly1305_blocks_avx512,.-poly1305_blocks_avx512
2710 ########################################################################
2711 # VPMADD52 version using 2^44 radix.
2713 # One can argue that base 2^52 would be more natural. Well, even though
2714 # some operations would be more natural, one has to recognize couple of
2715 # things. Base 2^52 doesn't provide advantage over base 2^44 if you look
2716 # at amount of multiply-n-accumulate operations. Secondly, it makes it
2717 # impossible to pre-compute multiples of 5 [referred to as s[]/sN in
2718 # reference implementations], which means that more such operations
2719 # would have to be performed in inner loop, which in turn makes critical
2720 # path longer. In other words, even though base 2^44 reduction might
2721 # look less elegant, overall critical path is actually shorter...
2723 ########################################################################
2724 # Layout of opaque area is following.
2726 # unsigned __int64 h[3]; # current hash value base 2^44
2727 # unsigned __int64 s[2]; # key value*20 base 2^44
2728 # unsigned __int64 r[3]; # key value base 2^44
2729 # struct { unsigned __int64 r^1, r^3, r^2, r^4; } R[4];
2730 # # r^n positions reflect
2731 # # placement in register, not
2732 # # memory, R[3] is R[1]*20
2735 .type poly1305_init_base2_44,\@function,3
2737 poly1305_init_base2_44:
2739 mov %rax,0($ctx) # initialize hash value
2744 lea poly1305_blocks_vpmadd52(%rip),%r10
2745 lea poly1305_emit_base2_44(%rip),%r11
2747 mov \$0x0ffffffc0fffffff,%rax
2748 mov \$0x0ffffffc0ffffffc,%rcx
2750 mov \$0x00000fffffffffff,%r8
2752 mov \$0x00000fffffffffff,%r9
2755 mov %r8,40($ctx) # r0
2758 mov %rax,48($ctx) # r1
2759 lea (%rax,%rax,4),%rax # *5
2760 mov %rcx,56($ctx) # r2
2761 shl \$2,%rax # magic <<2
2762 lea (%rcx,%rcx,4),%rcx # *5
2763 shl \$2,%rcx # magic <<2
2764 mov %rax,24($ctx) # s1
2765 mov %rcx,32($ctx) # s2
2766 movq \$-1,64($ctx) # write impossible value
2768 $code.=<<___ if ($flavour !~ /elf32/);
2772 $code.=<<___ if ($flavour =~ /elf32/);
2779 .size poly1305_init_base2_44,.-poly1305_init_base2_44
2782 my ($H0,$H1,$H2,$r2r1r0,$r1r0s2,$r0s2s1,$Dlo,$Dhi) = map("%ymm$_",(0..5,16,17));
2783 my ($T0,$inp_permd,$inp_shift,$PAD) = map("%ymm$_",(18..21));
2784 my ($reduc_mask,$reduc_rght,$reduc_left) = map("%ymm$_",(22..25));
2787 .type poly1305_blocks_vpmadd52,\@function,4
2789 poly1305_blocks_vpmadd52:
2791 jz .Lno_data_vpmadd52 # too short
2794 mov 64($ctx),%r8 # peek on power of the key
2796 # if powers of the key are not calculated yet, process up to 3
2797 # blocks with this single-block subroutine, otherwise ensure that
2798 # length is divisible by 2 blocks and pass the rest down to next
2803 cmp \$4,$len # is input long
2805 test %r8,%r8 # is power value impossible?
2808 and $len,%rax # is input of favourable length?
2809 jz .Lblocks_vpmadd52_4x
2815 lea .L2_44_inp_permd(%rip),%r10
2818 vmovq $padbit,%x#$PAD
2819 vmovdqa64 0(%r10),$inp_permd # .L2_44_inp_permd
2820 vmovdqa64 32(%r10),$inp_shift # .L2_44_inp_shift
2821 vpermq \$0xcf,$PAD,$PAD
2822 vmovdqa64 64(%r10),$reduc_mask # .L2_44_mask
2824 vmovdqu64 0($ctx),${Dlo}{%k7}{z} # load hash value
2825 vmovdqu64 40($ctx),${r2r1r0}{%k7}{z} # load keys
2826 vmovdqu64 32($ctx),${r1r0s2}{%k7}{z}
2827 vmovdqu64 24($ctx),${r0s2s1}{%k7}{z}
2829 vmovdqa64 96(%r10),$reduc_rght # .L2_44_shift_rgt
2830 vmovdqa64 128(%r10),$reduc_left # .L2_44_shift_lft
2836 vmovdqu32 0($inp),%x#$T0 # load input as ----3210
2839 vpermd $T0,$inp_permd,$T0 # ----3210 -> --322110
2840 vpsrlvq $inp_shift,$T0,$T0
2841 vpandq $reduc_mask,$T0,$T0
2844 vpaddq $T0,$Dlo,$Dlo # accumulate input
2846 vpermq \$0,$Dlo,${H0}{%k7}{z} # smash hash value
2847 vpermq \$0b01010101,$Dlo,${H1}{%k7}{z}
2848 vpermq \$0b10101010,$Dlo,${H2}{%k7}{z}
2850 vpxord $Dlo,$Dlo,$Dlo
2851 vpxord $Dhi,$Dhi,$Dhi
2853 vpmadd52luq $r2r1r0,$H0,$Dlo
2854 vpmadd52huq $r2r1r0,$H0,$Dhi
2856 vpmadd52luq $r1r0s2,$H1,$Dlo
2857 vpmadd52huq $r1r0s2,$H1,$Dhi
2859 vpmadd52luq $r0s2s1,$H2,$Dlo
2860 vpmadd52huq $r0s2s1,$H2,$Dhi
2862 vpsrlvq $reduc_rght,$Dlo,$T0 # 0 in topmost qword
2863 vpsllvq $reduc_left,$Dhi,$Dhi # 0 in topmost qword
2864 vpandq $reduc_mask,$Dlo,$Dlo
2866 vpaddq $T0,$Dhi,$Dhi
2868 vpermq \$0b10010011,$Dhi,$Dhi # 0 in lowest qword
2870 vpaddq $Dhi,$Dlo,$Dlo # note topmost qword :-)
2872 vpsrlvq $reduc_rght,$Dlo,$T0 # 0 in topmost word
2873 vpandq $reduc_mask,$Dlo,$Dlo
2875 vpermq \$0b10010011,$T0,$T0
2877 vpaddq $T0,$Dlo,$Dlo
2879 vpermq \$0b10010011,$Dlo,${T0}{%k1}{z}
2881 vpaddq $T0,$Dlo,$Dlo
2884 vpaddq $T0,$Dlo,$Dlo
2889 vmovdqu64 $Dlo,0($ctx){%k7} # store hash value
2892 jnz .Lblocks_vpmadd52_4x
2896 .size poly1305_blocks_vpmadd52,.-poly1305_blocks_vpmadd52
2900 ########################################################################
2901 # As implied by its name 4x subroutine processes 4 blocks in parallel
2902 # (but handles even 4*n+2 blocks lengths). It takes up to 4th key power
2903 # and is handled in 256-bit %ymm registers.
2905 my ($H0,$H1,$H2,$R0,$R1,$R2,$S1,$S2) = map("%ymm$_",(0..5,16,17));
2906 my ($D0lo,$D0hi,$D1lo,$D1hi,$D2lo,$D2hi) = map("%ymm$_",(18..23));
2907 my ($T0,$T1,$T2,$T3,$mask44,$mask42,$tmp,$PAD) = map("%ymm$_",(24..31));
2910 .type poly1305_blocks_vpmadd52_4x,\@function,4
2912 poly1305_blocks_vpmadd52_4x:
2914 jz .Lno_data_vpmadd52_4x # too short
2917 mov 64($ctx),%r8 # peek on power of the key
2919 .Lblocks_vpmadd52_4x:
2920 vpbroadcastq $padbit,$PAD
2922 vmovdqa64 .Lx_mask44(%rip),$mask44
2924 vmovdqa64 .Lx_mask42(%rip),$mask42
2925 kmovw %eax,%k1 # used in 2x path
2927 test %r8,%r8 # is power value impossible?
2928 js .Linit_vpmadd52 # if it is, then init R[4]
2930 vmovq 0($ctx),%x#$H0 # load current hash value
2931 vmovq 8($ctx),%x#$H1
2932 vmovq 16($ctx),%x#$H2
2934 test \$3,$len # is length 4*n+2?
2935 jnz .Lblocks_vpmadd52_2x_do
2937 .Lblocks_vpmadd52_4x_do:
2938 vpbroadcastq 64($ctx),$R0 # load 4th power of the key
2939 vpbroadcastq 96($ctx),$R1
2940 vpbroadcastq 128($ctx),$R2
2941 vpbroadcastq 160($ctx),$S1
2943 .Lblocks_vpmadd52_4x_key_loaded:
2944 vpsllq \$2,$R2,$S2 # S2 = R2*5*4
2948 test \$7,$len # is len 8*n?
2949 jz .Lblocks_vpmadd52_8x
2951 vmovdqu64 16*0($inp),$T2 # load data
2952 vmovdqu64 16*2($inp),$T3
2955 vpunpcklqdq $T3,$T2,$T1 # transpose data
2956 vpunpckhqdq $T3,$T2,$T3
2958 # at this point 64-bit lanes are ordered as 3-1-2-0
2960 vpsrlq \$24,$T3,$T2 # splat the data
2962 vpaddq $T2,$H2,$H2 # accumulate input
2963 vpandq $mask44,$T1,$T0
2967 vpandq $mask44,$T1,$T1
2970 jz .Ltail_vpmadd52_4x
2971 jmp .Loop_vpmadd52_4x
2976 vmovq 24($ctx),%x#$S1 # load key
2977 vmovq 56($ctx),%x#$H2
2978 vmovq 32($ctx),%x#$S2
2979 vmovq 40($ctx),%x#$R0
2980 vmovq 48($ctx),%x#$R1
2988 .Lmul_init_vpmadd52:
2989 vpxorq $D0lo,$D0lo,$D0lo
2990 vpmadd52luq $H2,$S1,$D0lo
2991 vpxorq $D0hi,$D0hi,$D0hi
2992 vpmadd52huq $H2,$S1,$D0hi
2993 vpxorq $D1lo,$D1lo,$D1lo
2994 vpmadd52luq $H2,$S2,$D1lo
2995 vpxorq $D1hi,$D1hi,$D1hi
2996 vpmadd52huq $H2,$S2,$D1hi
2997 vpxorq $D2lo,$D2lo,$D2lo
2998 vpmadd52luq $H2,$R0,$D2lo
2999 vpxorq $D2hi,$D2hi,$D2hi
3000 vpmadd52huq $H2,$R0,$D2hi
3002 vpmadd52luq $H0,$R0,$D0lo
3003 vpmadd52huq $H0,$R0,$D0hi
3004 vpmadd52luq $H0,$R1,$D1lo
3005 vpmadd52huq $H0,$R1,$D1hi
3006 vpmadd52luq $H0,$R2,$D2lo
3007 vpmadd52huq $H0,$R2,$D2hi
3009 vpmadd52luq $H1,$S2,$D0lo
3010 vpmadd52huq $H1,$S2,$D0hi
3011 vpmadd52luq $H1,$R0,$D1lo
3012 vpmadd52huq $H1,$R0,$D1hi
3013 vpmadd52luq $H1,$R1,$D2lo
3014 vpmadd52huq $H1,$R1,$D2hi
3016 ################################################################
3018 vpsrlq \$44,$D0lo,$tmp
3019 vpsllq \$8,$D0hi,$D0hi
3020 vpandq $mask44,$D0lo,$H0
3021 vpaddq $tmp,$D0hi,$D0hi
3023 vpaddq $D0hi,$D1lo,$D1lo
3025 vpsrlq \$44,$D1lo,$tmp
3026 vpsllq \$8,$D1hi,$D1hi
3027 vpandq $mask44,$D1lo,$H1
3028 vpaddq $tmp,$D1hi,$D1hi
3030 vpaddq $D1hi,$D2lo,$D2lo
3032 vpsrlq \$42,$D2lo,$tmp
3033 vpsllq \$10,$D2hi,$D2hi
3034 vpandq $mask42,$D2lo,$H2
3035 vpaddq $tmp,$D2hi,$D2hi
3037 vpaddq $D2hi,$H0,$H0
3038 vpsllq \$2,$D2hi,$D2hi
3040 vpaddq $D2hi,$H0,$H0
3042 vpsrlq \$44,$H0,$tmp # additional step
3043 vpandq $mask44,$H0,$H0
3048 jz .Ldone_init_vpmadd52
3050 vpunpcklqdq $R1,$H1,$R1 # 1,2
3051 vpbroadcastq %x#$H1,%x#$H1 # 2,2
3052 vpunpcklqdq $R2,$H2,$R2
3053 vpbroadcastq %x#$H2,%x#$H2
3054 vpunpcklqdq $R0,$H0,$R0
3055 vpbroadcastq %x#$H0,%x#$H0
3057 vpsllq \$2,$R1,$S1 # S1 = R1*5*4
3058 vpsllq \$2,$R2,$S2 # S2 = R2*5*4
3064 jmp .Lmul_init_vpmadd52
3068 .Ldone_init_vpmadd52:
3069 vinserti128 \$1,%x#$R1,$H1,$R1 # 1,2,3,4
3070 vinserti128 \$1,%x#$R2,$H2,$R2
3071 vinserti128 \$1,%x#$R0,$H0,$R0
3073 vpermq \$0b11011000,$R1,$R1 # 1,3,2,4
3074 vpermq \$0b11011000,$R2,$R2
3075 vpermq \$0b11011000,$R0,$R0
3077 vpsllq \$2,$R1,$S1 # S1 = R1*5*4
3081 vmovq 0($ctx),%x#$H0 # load current hash value
3082 vmovq 8($ctx),%x#$H1
3083 vmovq 16($ctx),%x#$H2
3085 test \$3,$len # is length 4*n+2?
3086 jnz .Ldone_init_vpmadd52_2x
3088 vmovdqu64 $R0,64($ctx) # save key powers
3089 vpbroadcastq %x#$R0,$R0 # broadcast 4th power
3090 vmovdqu64 $R1,96($ctx)
3091 vpbroadcastq %x#$R1,$R1
3092 vmovdqu64 $R2,128($ctx)
3093 vpbroadcastq %x#$R2,$R2
3094 vmovdqu64 $S1,160($ctx)
3095 vpbroadcastq %x#$S1,$S1
3097 jmp .Lblocks_vpmadd52_4x_key_loaded
3101 .Ldone_init_vpmadd52_2x:
3102 vmovdqu64 $R0,64($ctx) # save key powers
3103 vpsrldq \$8,$R0,$R0 # 0-1-0-2
3104 vmovdqu64 $R1,96($ctx)
3106 vmovdqu64 $R2,128($ctx)
3108 vmovdqu64 $S1,160($ctx)
3110 jmp .Lblocks_vpmadd52_2x_key_loaded
3114 .Lblocks_vpmadd52_2x_do:
3115 vmovdqu64 128+8($ctx),${R2}{%k1}{z}# load 2nd and 1st key powers
3116 vmovdqu64 160+8($ctx),${S1}{%k1}{z}
3117 vmovdqu64 64+8($ctx),${R0}{%k1}{z}
3118 vmovdqu64 96+8($ctx),${R1}{%k1}{z}
3120 .Lblocks_vpmadd52_2x_key_loaded:
3121 vmovdqu64 16*0($inp),$T2 # load data
3125 vpunpcklqdq $T3,$T2,$T1 # transpose data
3126 vpunpckhqdq $T3,$T2,$T3
3128 # at this point 64-bit lanes are ordered as x-1-x-0
3130 vpsrlq \$24,$T3,$T2 # splat the data
3132 vpaddq $T2,$H2,$H2 # accumulate input
3133 vpandq $mask44,$T1,$T0
3137 vpandq $mask44,$T1,$T1
3139 jmp .Ltail_vpmadd52_2x
3144 #vpaddq $T2,$H2,$H2 # accumulate input
3148 vpxorq $D0lo,$D0lo,$D0lo
3149 vpmadd52luq $H2,$S1,$D0lo
3150 vpxorq $D0hi,$D0hi,$D0hi
3151 vpmadd52huq $H2,$S1,$D0hi
3152 vpxorq $D1lo,$D1lo,$D1lo
3153 vpmadd52luq $H2,$S2,$D1lo
3154 vpxorq $D1hi,$D1hi,$D1hi
3155 vpmadd52huq $H2,$S2,$D1hi
3156 vpxorq $D2lo,$D2lo,$D2lo
3157 vpmadd52luq $H2,$R0,$D2lo
3158 vpxorq $D2hi,$D2hi,$D2hi
3159 vpmadd52huq $H2,$R0,$D2hi
3161 vmovdqu64 16*0($inp),$T2 # load data
3162 vmovdqu64 16*2($inp),$T3
3164 vpmadd52luq $H0,$R0,$D0lo
3165 vpmadd52huq $H0,$R0,$D0hi
3166 vpmadd52luq $H0,$R1,$D1lo
3167 vpmadd52huq $H0,$R1,$D1hi
3168 vpmadd52luq $H0,$R2,$D2lo
3169 vpmadd52huq $H0,$R2,$D2hi
3171 vpunpcklqdq $T3,$T2,$T1 # transpose data
3172 vpunpckhqdq $T3,$T2,$T3
3173 vpmadd52luq $H1,$S2,$D0lo
3174 vpmadd52huq $H1,$S2,$D0hi
3175 vpmadd52luq $H1,$R0,$D1lo
3176 vpmadd52huq $H1,$R0,$D1hi
3177 vpmadd52luq $H1,$R1,$D2lo
3178 vpmadd52huq $H1,$R1,$D2hi
3180 ################################################################
3181 # partial reduction (interleaved with data splat)
3182 vpsrlq \$44,$D0lo,$tmp
3183 vpsllq \$8,$D0hi,$D0hi
3184 vpandq $mask44,$D0lo,$H0
3185 vpaddq $tmp,$D0hi,$D0hi
3189 vpaddq $D0hi,$D1lo,$D1lo
3191 vpsrlq \$44,$D1lo,$tmp
3192 vpsllq \$8,$D1hi,$D1hi
3193 vpandq $mask44,$D1lo,$H1
3194 vpaddq $tmp,$D1hi,$D1hi
3196 vpandq $mask44,$T1,$T0
3199 vpaddq $D1hi,$D2lo,$D2lo
3201 vpsrlq \$42,$D2lo,$tmp
3202 vpsllq \$10,$D2hi,$D2hi
3203 vpandq $mask42,$D2lo,$H2
3204 vpaddq $tmp,$D2hi,$D2hi
3206 vpaddq $T2,$H2,$H2 # accumulate input
3207 vpaddq $D2hi,$H0,$H0
3208 vpsllq \$2,$D2hi,$D2hi
3210 vpaddq $D2hi,$H0,$H0
3212 vpandq $mask44,$T1,$T1
3214 vpsrlq \$44,$H0,$tmp # additional step
3215 vpandq $mask44,$H0,$H0
3219 sub \$4,$len # len-=64
3220 jnz .Loop_vpmadd52_4x
3223 vmovdqu64 128($ctx),$R2 # load all key powers
3224 vmovdqu64 160($ctx),$S1
3225 vmovdqu64 64($ctx),$R0
3226 vmovdqu64 96($ctx),$R1
3229 vpsllq \$2,$R2,$S2 # S2 = R2*5*4
3233 #vpaddq $T2,$H2,$H2 # accumulate input
3237 vpxorq $D0lo,$D0lo,$D0lo
3238 vpmadd52luq $H2,$S1,$D0lo
3239 vpxorq $D0hi,$D0hi,$D0hi
3240 vpmadd52huq $H2,$S1,$D0hi
3241 vpxorq $D1lo,$D1lo,$D1lo
3242 vpmadd52luq $H2,$S2,$D1lo
3243 vpxorq $D1hi,$D1hi,$D1hi
3244 vpmadd52huq $H2,$S2,$D1hi
3245 vpxorq $D2lo,$D2lo,$D2lo
3246 vpmadd52luq $H2,$R0,$D2lo
3247 vpxorq $D2hi,$D2hi,$D2hi
3248 vpmadd52huq $H2,$R0,$D2hi
3250 vpmadd52luq $H0,$R0,$D0lo
3251 vpmadd52huq $H0,$R0,$D0hi
3252 vpmadd52luq $H0,$R1,$D1lo
3253 vpmadd52huq $H0,$R1,$D1hi
3254 vpmadd52luq $H0,$R2,$D2lo
3255 vpmadd52huq $H0,$R2,$D2hi
3257 vpmadd52luq $H1,$S2,$D0lo
3258 vpmadd52huq $H1,$S2,$D0hi
3259 vpmadd52luq $H1,$R0,$D1lo
3260 vpmadd52huq $H1,$R0,$D1hi
3261 vpmadd52luq $H1,$R1,$D2lo
3262 vpmadd52huq $H1,$R1,$D2hi
3264 ################################################################
3265 # horizontal addition
3269 vpsrldq \$8,$D0lo,$T0
3270 vpsrldq \$8,$D0hi,$H0
3271 vpsrldq \$8,$D1lo,$T1
3272 vpsrldq \$8,$D1hi,$H1
3273 vpaddq $T0,$D0lo,$D0lo
3274 vpaddq $H0,$D0hi,$D0hi
3275 vpsrldq \$8,$D2lo,$T2
3276 vpsrldq \$8,$D2hi,$H2
3277 vpaddq $T1,$D1lo,$D1lo
3278 vpaddq $H1,$D1hi,$D1hi
3279 vpermq \$0x2,$D0lo,$T0
3280 vpermq \$0x2,$D0hi,$H0
3281 vpaddq $T2,$D2lo,$D2lo
3282 vpaddq $H2,$D2hi,$D2hi
3284 vpermq \$0x2,$D1lo,$T1
3285 vpermq \$0x2,$D1hi,$H1
3286 vpaddq $T0,$D0lo,${D0lo}{%k1}{z}
3287 vpaddq $H0,$D0hi,${D0hi}{%k1}{z}
3288 vpermq \$0x2,$D2lo,$T2
3289 vpermq \$0x2,$D2hi,$H2
3290 vpaddq $T1,$D1lo,${D1lo}{%k1}{z}
3291 vpaddq $H1,$D1hi,${D1hi}{%k1}{z}
3292 vpaddq $T2,$D2lo,${D2lo}{%k1}{z}
3293 vpaddq $H2,$D2hi,${D2hi}{%k1}{z}
3295 ################################################################
3297 vpsrlq \$44,$D0lo,$tmp
3298 vpsllq \$8,$D0hi,$D0hi
3299 vpandq $mask44,$D0lo,$H0
3300 vpaddq $tmp,$D0hi,$D0hi
3302 vpaddq $D0hi,$D1lo,$D1lo
3304 vpsrlq \$44,$D1lo,$tmp
3305 vpsllq \$8,$D1hi,$D1hi
3306 vpandq $mask44,$D1lo,$H1
3307 vpaddq $tmp,$D1hi,$D1hi
3309 vpaddq $D1hi,$D2lo,$D2lo
3311 vpsrlq \$42,$D2lo,$tmp
3312 vpsllq \$10,$D2hi,$D2hi
3313 vpandq $mask42,$D2lo,$H2
3314 vpaddq $tmp,$D2hi,$D2hi
3316 vpaddq $D2hi,$H0,$H0
3317 vpsllq \$2,$D2hi,$D2hi
3319 vpaddq $D2hi,$H0,$H0
3321 vpsrlq \$44,$H0,$tmp # additional step
3322 vpandq $mask44,$H0,$H0
3325 # at this point $len is
3326 # either 4*n+2 or 0...
3327 sub \$2,$len # len-=32
3328 ja .Lblocks_vpmadd52_4x_do
3330 vmovq %x#$H0,0($ctx)
3331 vmovq %x#$H1,8($ctx)
3332 vmovq %x#$H2,16($ctx)
3335 .Lno_data_vpmadd52_4x:
3337 .size poly1305_blocks_vpmadd52_4x,.-poly1305_blocks_vpmadd52_4x
3341 ########################################################################
3342 # As implied by its name 8x subroutine processes 8 blocks in parallel...
3343 # This is intermediate version, as it's used only in cases when input
3344 # length is either 8*n, 8*n+1 or 8*n+2...
3346 my ($H0,$H1,$H2,$R0,$R1,$R2,$S1,$S2) = map("%ymm$_",(0..5,16,17));
3347 my ($D0lo,$D0hi,$D1lo,$D1hi,$D2lo,$D2hi) = map("%ymm$_",(18..23));
3348 my ($T0,$T1,$T2,$T3,$mask44,$mask42,$tmp,$PAD) = map("%ymm$_",(24..31));
3349 my ($RR0,$RR1,$RR2,$SS1,$SS2) = map("%ymm$_",(6..10));
3352 .type poly1305_blocks_vpmadd52_8x,\@function,4
3354 poly1305_blocks_vpmadd52_8x:
3356 jz .Lno_data_vpmadd52_8x # too short
3359 mov 64($ctx),%r8 # peek on power of the key
3361 vmovdqa64 .Lx_mask44(%rip),$mask44
3362 vmovdqa64 .Lx_mask42(%rip),$mask42
3364 test %r8,%r8 # is power value impossible?
3365 js .Linit_vpmadd52 # if it is, then init R[4]
3367 vmovq 0($ctx),%x#$H0 # load current hash value
3368 vmovq 8($ctx),%x#$H1
3369 vmovq 16($ctx),%x#$H2
3371 .Lblocks_vpmadd52_8x:
3372 ################################################################
3373 # fist we calculate more key powers
3375 vmovdqu64 128($ctx),$R2 # load 1-3-2-4 powers
3376 vmovdqu64 160($ctx),$S1
3377 vmovdqu64 64($ctx),$R0
3378 vmovdqu64 96($ctx),$R1
3380 vpsllq \$2,$R2,$S2 # S2 = R2*5*4
3384 vpbroadcastq %x#$R2,$RR2 # broadcast 4th power
3385 vpbroadcastq %x#$R0,$RR0
3386 vpbroadcastq %x#$R1,$RR1
3388 vpxorq $D0lo,$D0lo,$D0lo
3389 vpmadd52luq $RR2,$S1,$D0lo
3390 vpxorq $D0hi,$D0hi,$D0hi
3391 vpmadd52huq $RR2,$S1,$D0hi
3392 vpxorq $D1lo,$D1lo,$D1lo
3393 vpmadd52luq $RR2,$S2,$D1lo
3394 vpxorq $D1hi,$D1hi,$D1hi
3395 vpmadd52huq $RR2,$S2,$D1hi
3396 vpxorq $D2lo,$D2lo,$D2lo
3397 vpmadd52luq $RR2,$R0,$D2lo
3398 vpxorq $D2hi,$D2hi,$D2hi
3399 vpmadd52huq $RR2,$R0,$D2hi
3401 vpmadd52luq $RR0,$R0,$D0lo
3402 vpmadd52huq $RR0,$R0,$D0hi
3403 vpmadd52luq $RR0,$R1,$D1lo
3404 vpmadd52huq $RR0,$R1,$D1hi
3405 vpmadd52luq $RR0,$R2,$D2lo
3406 vpmadd52huq $RR0,$R2,$D2hi
3408 vpmadd52luq $RR1,$S2,$D0lo
3409 vpmadd52huq $RR1,$S2,$D0hi
3410 vpmadd52luq $RR1,$R0,$D1lo
3411 vpmadd52huq $RR1,$R0,$D1hi
3412 vpmadd52luq $RR1,$R1,$D2lo
3413 vpmadd52huq $RR1,$R1,$D2hi
3415 ################################################################
3417 vpsrlq \$44,$D0lo,$tmp
3418 vpsllq \$8,$D0hi,$D0hi
3419 vpandq $mask44,$D0lo,$RR0
3420 vpaddq $tmp,$D0hi,$D0hi
3422 vpaddq $D0hi,$D1lo,$D1lo
3424 vpsrlq \$44,$D1lo,$tmp
3425 vpsllq \$8,$D1hi,$D1hi
3426 vpandq $mask44,$D1lo,$RR1
3427 vpaddq $tmp,$D1hi,$D1hi
3429 vpaddq $D1hi,$D2lo,$D2lo
3431 vpsrlq \$42,$D2lo,$tmp
3432 vpsllq \$10,$D2hi,$D2hi
3433 vpandq $mask42,$D2lo,$RR2
3434 vpaddq $tmp,$D2hi,$D2hi
3436 vpaddq $D2hi,$RR0,$RR0
3437 vpsllq \$2,$D2hi,$D2hi
3439 vpaddq $D2hi,$RR0,$RR0
3441 vpsrlq \$44,$RR0,$tmp # additional step
3442 vpandq $mask44,$RR0,$RR0
3444 vpaddq $tmp,$RR1,$RR1
3446 ################################################################
3447 # At this point Rx holds 1324 powers, RRx - 5768, and the goal
3448 # is 15263748, which reflects how data is loaded...
3450 vpunpcklqdq $R2,$RR2,$T2 # 3748
3451 vpunpckhqdq $R2,$RR2,$R2 # 1526
3452 vpunpcklqdq $R0,$RR0,$T0
3453 vpunpckhqdq $R0,$RR0,$R0
3454 vpunpcklqdq $R1,$RR1,$T1
3455 vpunpckhqdq $R1,$RR1,$R1
3457 ######## switch to %zmm
3458 map(s/%y/%z/, $H0,$H1,$H2,$R0,$R1,$R2,$S1,$S2);
3459 map(s/%y/%z/, $D0lo,$D0hi,$D1lo,$D1hi,$D2lo,$D2hi);
3460 map(s/%y/%z/, $T0,$T1,$T2,$T3,$mask44,$mask42,$tmp,$PAD);
3461 map(s/%y/%z/, $RR0,$RR1,$RR2,$SS1,$SS2);
3464 vshufi64x2 \$0x44,$R2,$T2,$RR2 # 15263748
3465 vshufi64x2 \$0x44,$R0,$T0,$RR0
3466 vshufi64x2 \$0x44,$R1,$T1,$RR1
3468 vmovdqu64 16*0($inp),$T2 # load data
3469 vmovdqu64 16*4($inp),$T3
3472 vpsllq \$2,$RR2,$SS2 # S2 = R2*5*4
3473 vpsllq \$2,$RR1,$SS1 # S1 = R1*5*4
3474 vpaddq $RR2,$SS2,$SS2
3475 vpaddq $RR1,$SS1,$SS1
3476 vpsllq \$2,$SS2,$SS2
3477 vpsllq \$2,$SS1,$SS1
3479 vpbroadcastq $padbit,$PAD
3480 vpbroadcastq %x#$mask44,$mask44
3481 vpbroadcastq %x#$mask42,$mask42
3483 vpbroadcastq %x#$SS1,$S1 # broadcast 8th power
3484 vpbroadcastq %x#$SS2,$S2
3485 vpbroadcastq %x#$RR0,$R0
3486 vpbroadcastq %x#$RR1,$R1
3487 vpbroadcastq %x#$RR2,$R2
3489 vpunpcklqdq $T3,$T2,$T1 # transpose data
3490 vpunpckhqdq $T3,$T2,$T3
3492 # at this point 64-bit lanes are ordered as 73625140
3494 vpsrlq \$24,$T3,$T2 # splat the data
3496 vpaddq $T2,$H2,$H2 # accumulate input
3497 vpandq $mask44,$T1,$T0
3501 vpandq $mask44,$T1,$T1
3504 jz .Ltail_vpmadd52_8x
3505 jmp .Loop_vpmadd52_8x
3509 #vpaddq $T2,$H2,$H2 # accumulate input
3513 vpxorq $D0lo,$D0lo,$D0lo
3514 vpmadd52luq $H2,$S1,$D0lo
3515 vpxorq $D0hi,$D0hi,$D0hi
3516 vpmadd52huq $H2,$S1,$D0hi
3517 vpxorq $D1lo,$D1lo,$D1lo
3518 vpmadd52luq $H2,$S2,$D1lo
3519 vpxorq $D1hi,$D1hi,$D1hi
3520 vpmadd52huq $H2,$S2,$D1hi
3521 vpxorq $D2lo,$D2lo,$D2lo
3522 vpmadd52luq $H2,$R0,$D2lo
3523 vpxorq $D2hi,$D2hi,$D2hi
3524 vpmadd52huq $H2,$R0,$D2hi
3526 vmovdqu64 16*0($inp),$T2 # load data
3527 vmovdqu64 16*4($inp),$T3
3529 vpmadd52luq $H0,$R0,$D0lo
3530 vpmadd52huq $H0,$R0,$D0hi
3531 vpmadd52luq $H0,$R1,$D1lo
3532 vpmadd52huq $H0,$R1,$D1hi
3533 vpmadd52luq $H0,$R2,$D2lo
3534 vpmadd52huq $H0,$R2,$D2hi
3536 vpunpcklqdq $T3,$T2,$T1 # transpose data
3537 vpunpckhqdq $T3,$T2,$T3
3538 vpmadd52luq $H1,$S2,$D0lo
3539 vpmadd52huq $H1,$S2,$D0hi
3540 vpmadd52luq $H1,$R0,$D1lo
3541 vpmadd52huq $H1,$R0,$D1hi
3542 vpmadd52luq $H1,$R1,$D2lo
3543 vpmadd52huq $H1,$R1,$D2hi
3545 ################################################################
3546 # partial reduction (interleaved with data splat)
3547 vpsrlq \$44,$D0lo,$tmp
3548 vpsllq \$8,$D0hi,$D0hi
3549 vpandq $mask44,$D0lo,$H0
3550 vpaddq $tmp,$D0hi,$D0hi
3554 vpaddq $D0hi,$D1lo,$D1lo
3556 vpsrlq \$44,$D1lo,$tmp
3557 vpsllq \$8,$D1hi,$D1hi
3558 vpandq $mask44,$D1lo,$H1
3559 vpaddq $tmp,$D1hi,$D1hi
3561 vpandq $mask44,$T1,$T0
3564 vpaddq $D1hi,$D2lo,$D2lo
3566 vpsrlq \$42,$D2lo,$tmp
3567 vpsllq \$10,$D2hi,$D2hi
3568 vpandq $mask42,$D2lo,$H2
3569 vpaddq $tmp,$D2hi,$D2hi
3571 vpaddq $T2,$H2,$H2 # accumulate input
3572 vpaddq $D2hi,$H0,$H0
3573 vpsllq \$2,$D2hi,$D2hi
3575 vpaddq $D2hi,$H0,$H0
3577 vpandq $mask44,$T1,$T1
3579 vpsrlq \$44,$H0,$tmp # additional step
3580 vpandq $mask44,$H0,$H0
3584 sub \$8,$len # len-=128
3585 jnz .Loop_vpmadd52_8x
3588 #vpaddq $T2,$H2,$H2 # accumulate input
3592 vpxorq $D0lo,$D0lo,$D0lo
3593 vpmadd52luq $H2,$SS1,$D0lo
3594 vpxorq $D0hi,$D0hi,$D0hi
3595 vpmadd52huq $H2,$SS1,$D0hi
3596 vpxorq $D1lo,$D1lo,$D1lo
3597 vpmadd52luq $H2,$SS2,$D1lo
3598 vpxorq $D1hi,$D1hi,$D1hi
3599 vpmadd52huq $H2,$SS2,$D1hi
3600 vpxorq $D2lo,$D2lo,$D2lo
3601 vpmadd52luq $H2,$RR0,$D2lo
3602 vpxorq $D2hi,$D2hi,$D2hi
3603 vpmadd52huq $H2,$RR0,$D2hi
3605 vpmadd52luq $H0,$RR0,$D0lo
3606 vpmadd52huq $H0,$RR0,$D0hi
3607 vpmadd52luq $H0,$RR1,$D1lo
3608 vpmadd52huq $H0,$RR1,$D1hi
3609 vpmadd52luq $H0,$RR2,$D2lo
3610 vpmadd52huq $H0,$RR2,$D2hi
3612 vpmadd52luq $H1,$SS2,$D0lo
3613 vpmadd52huq $H1,$SS2,$D0hi
3614 vpmadd52luq $H1,$RR0,$D1lo
3615 vpmadd52huq $H1,$RR0,$D1hi
3616 vpmadd52luq $H1,$RR1,$D2lo
3617 vpmadd52huq $H1,$RR1,$D2hi
3619 ################################################################
3620 # horizontal addition
3624 vpsrldq \$8,$D0lo,$T0
3625 vpsrldq \$8,$D0hi,$H0
3626 vpsrldq \$8,$D1lo,$T1
3627 vpsrldq \$8,$D1hi,$H1
3628 vpaddq $T0,$D0lo,$D0lo
3629 vpaddq $H0,$D0hi,$D0hi
3630 vpsrldq \$8,$D2lo,$T2
3631 vpsrldq \$8,$D2hi,$H2
3632 vpaddq $T1,$D1lo,$D1lo
3633 vpaddq $H1,$D1hi,$D1hi
3634 vpermq \$0x2,$D0lo,$T0
3635 vpermq \$0x2,$D0hi,$H0
3636 vpaddq $T2,$D2lo,$D2lo
3637 vpaddq $H2,$D2hi,$D2hi
3639 vpermq \$0x2,$D1lo,$T1
3640 vpermq \$0x2,$D1hi,$H1
3641 vpaddq $T0,$D0lo,$D0lo
3642 vpaddq $H0,$D0hi,$D0hi
3643 vpermq \$0x2,$D2lo,$T2
3644 vpermq \$0x2,$D2hi,$H2
3645 vpaddq $T1,$D1lo,$D1lo
3646 vpaddq $H1,$D1hi,$D1hi
3647 vextracti64x4 \$1,$D0lo,%y#$T0
3648 vextracti64x4 \$1,$D0hi,%y#$H0
3649 vpaddq $T2,$D2lo,$D2lo
3650 vpaddq $H2,$D2hi,$D2hi
3652 vextracti64x4 \$1,$D1lo,%y#$T1
3653 vextracti64x4 \$1,$D1hi,%y#$H1
3654 vextracti64x4 \$1,$D2lo,%y#$T2
3655 vextracti64x4 \$1,$D2hi,%y#$H2
3657 ######## switch back to %ymm
3658 map(s/%z/%y/, $H0,$H1,$H2,$R0,$R1,$R2,$S1,$S2);
3659 map(s/%z/%y/, $D0lo,$D0hi,$D1lo,$D1hi,$D2lo,$D2hi);
3660 map(s/%z/%y/, $T0,$T1,$T2,$T3,$mask44,$mask42,$tmp,$PAD);
3663 vpaddq $T0,$D0lo,${D0lo}{%k1}{z}
3664 vpaddq $H0,$D0hi,${D0hi}{%k1}{z}
3665 vpaddq $T1,$D1lo,${D1lo}{%k1}{z}
3666 vpaddq $H1,$D1hi,${D1hi}{%k1}{z}
3667 vpaddq $T2,$D2lo,${D2lo}{%k1}{z}
3668 vpaddq $H2,$D2hi,${D2hi}{%k1}{z}
3670 ################################################################
3672 vpsrlq \$44,$D0lo,$tmp
3673 vpsllq \$8,$D0hi,$D0hi
3674 vpandq $mask44,$D0lo,$H0
3675 vpaddq $tmp,$D0hi,$D0hi
3677 vpaddq $D0hi,$D1lo,$D1lo
3679 vpsrlq \$44,$D1lo,$tmp
3680 vpsllq \$8,$D1hi,$D1hi
3681 vpandq $mask44,$D1lo,$H1
3682 vpaddq $tmp,$D1hi,$D1hi
3684 vpaddq $D1hi,$D2lo,$D2lo
3686 vpsrlq \$42,$D2lo,$tmp
3687 vpsllq \$10,$D2hi,$D2hi
3688 vpandq $mask42,$D2lo,$H2
3689 vpaddq $tmp,$D2hi,$D2hi
3691 vpaddq $D2hi,$H0,$H0
3692 vpsllq \$2,$D2hi,$D2hi
3694 vpaddq $D2hi,$H0,$H0
3696 vpsrlq \$44,$H0,$tmp # additional step
3697 vpandq $mask44,$H0,$H0
3701 ################################################################
3703 vmovq %x#$H0,0($ctx)
3704 vmovq %x#$H1,8($ctx)
3705 vmovq %x#$H2,16($ctx)
3708 .Lno_data_vpmadd52_8x:
3710 .size poly1305_blocks_vpmadd52_8x,.-poly1305_blocks_vpmadd52_8x
3714 .type poly1305_emit_base2_44,\@function,3
3716 poly1305_emit_base2_44:
3717 mov 0($ctx),%r8 # load hash value
3733 add \$5,%r8 # compare to modulus
3737 shr \$2,%r10 # did 130-bit value overfow?
3741 add 0($nonce),%rax # accumulate nonce
3743 mov %rax,0($mac) # write result
3747 .size poly1305_emit_base2_44,.-poly1305_emit_base2_44
3754 .long 0x0ffffff,0,0x0ffffff,0,0x0ffffff,0,0x0ffffff,0
3756 .long `1<<24`,0,`1<<24`,0,`1<<24`,0,`1<<24`,0
3758 .long 0x3ffffff,0,0x3ffffff,0,0x3ffffff,0,0x3ffffff,0
3760 .long 2,2,2,3,2,0,2,1
3762 .long 0,0,0,1, 0,2,0,3, 0,4,0,5, 0,6,0,7
3765 .long 0,1,1,2,2,3,7,7
3769 .quad 0xfffffffffff,0xfffffffffff,0x3ffffffffff,0xffffffffffffffff
3777 .quad 0xfffffffffff,0xfffffffffff,0xfffffffffff,0xfffffffffff
3778 .quad 0xfffffffffff,0xfffffffffff,0xfffffffffff,0xfffffffffff
3780 .quad 0x3ffffffffff,0x3ffffffffff,0x3ffffffffff,0x3ffffffffff
3781 .quad 0x3ffffffffff,0x3ffffffffff,0x3ffffffffff,0x3ffffffffff
3786 .asciz "Poly1305 for x86_64, CRYPTOGAMS by <appro\@openssl.org>"
3790 # EXCEPTION_DISPOSITION handler (EXCEPTION_RECORD *rec,ULONG64 frame,
3791 # CONTEXT *context,DISPATCHER_CONTEXT *disp)
3799 .extern __imp_RtlVirtualUnwind
3800 .type se_handler,\@abi-omnipotent
3814 mov 120($context),%rax # pull context->Rax
3815 mov 248($context),%rbx # pull context->Rip
3817 mov 8($disp),%rsi # disp->ImageBase
3818 mov 56($disp),%r11 # disp->HandlerData
3820 mov 0(%r11),%r10d # HandlerData[0]
3821 lea (%rsi,%r10),%r10 # prologue label
3822 cmp %r10,%rbx # context->Rip<.Lprologue
3823 jb .Lcommon_seh_tail
3825 mov 152($context),%rax # pull context->Rsp
3827 mov 4(%r11),%r10d # HandlerData[1]
3828 lea (%rsi,%r10),%r10 # epilogue label
3829 cmp %r10,%rbx # context->Rip>=.Lepilogue
3830 jae .Lcommon_seh_tail
3840 mov %rbx,144($context) # restore context->Rbx
3841 mov %rbp,160($context) # restore context->Rbp
3842 mov %r12,216($context) # restore context->R12
3843 mov %r13,224($context) # restore context->R13
3844 mov %r14,232($context) # restore context->R14
3845 mov %r15,240($context) # restore context->R14
3847 jmp .Lcommon_seh_tail
3848 .size se_handler,.-se_handler
3850 .type avx_handler,\@abi-omnipotent
3864 mov 120($context),%rax # pull context->Rax
3865 mov 248($context),%rbx # pull context->Rip
3867 mov 8($disp),%rsi # disp->ImageBase
3868 mov 56($disp),%r11 # disp->HandlerData
3870 mov 0(%r11),%r10d # HandlerData[0]
3871 lea (%rsi,%r10),%r10 # prologue label
3872 cmp %r10,%rbx # context->Rip<prologue label
3873 jb .Lcommon_seh_tail
3875 mov 152($context),%rax # pull context->Rsp
3877 mov 4(%r11),%r10d # HandlerData[1]
3878 lea (%rsi,%r10),%r10 # epilogue label
3879 cmp %r10,%rbx # context->Rip>=epilogue label
3880 jae .Lcommon_seh_tail
3882 mov 208($context),%rax # pull context->R11
3886 lea 512($context),%rdi # &context.Xmm6
3888 .long 0xa548f3fc # cld; rep movsq
3893 mov %rax,152($context) # restore context->Rsp
3894 mov %rsi,168($context) # restore context->Rsi
3895 mov %rdi,176($context) # restore context->Rdi
3897 mov 40($disp),%rdi # disp->ContextRecord
3898 mov $context,%rsi # context
3899 mov \$154,%ecx # sizeof(CONTEXT)
3900 .long 0xa548f3fc # cld; rep movsq
3903 xor %rcx,%rcx # arg1, UNW_FLAG_NHANDLER
3904 mov 8(%rsi),%rdx # arg2, disp->ImageBase
3905 mov 0(%rsi),%r8 # arg3, disp->ControlPc
3906 mov 16(%rsi),%r9 # arg4, disp->FunctionEntry
3907 mov 40(%rsi),%r10 # disp->ContextRecord
3908 lea 56(%rsi),%r11 # &disp->HandlerData
3909 lea 24(%rsi),%r12 # &disp->EstablisherFrame
3910 mov %r10,32(%rsp) # arg5
3911 mov %r11,40(%rsp) # arg6
3912 mov %r12,48(%rsp) # arg7
3913 mov %rcx,56(%rsp) # arg8, (NULL)
3914 call *__imp_RtlVirtualUnwind(%rip)
3916 mov \$1,%eax # ExceptionContinueSearch
3928 .size avx_handler,.-avx_handler
3932 .rva .LSEH_begin_poly1305_init
3933 .rva .LSEH_end_poly1305_init
3934 .rva .LSEH_info_poly1305_init
3936 .rva .LSEH_begin_poly1305_blocks
3937 .rva .LSEH_end_poly1305_blocks
3938 .rva .LSEH_info_poly1305_blocks
3940 .rva .LSEH_begin_poly1305_emit
3941 .rva .LSEH_end_poly1305_emit
3942 .rva .LSEH_info_poly1305_emit
3944 $code.=<<___ if ($avx);
3945 .rva .LSEH_begin_poly1305_blocks_avx
3947 .rva .LSEH_info_poly1305_blocks_avx_1
3951 .rva .LSEH_info_poly1305_blocks_avx_2
3954 .rva .LSEH_end_poly1305_blocks_avx
3955 .rva .LSEH_info_poly1305_blocks_avx_3
3957 .rva .LSEH_begin_poly1305_emit_avx
3958 .rva .LSEH_end_poly1305_emit_avx
3959 .rva .LSEH_info_poly1305_emit_avx
3961 $code.=<<___ if ($avx>1);
3962 .rva .LSEH_begin_poly1305_blocks_avx2
3963 .rva .Lbase2_64_avx2
3964 .rva .LSEH_info_poly1305_blocks_avx2_1
3966 .rva .Lbase2_64_avx2
3968 .rva .LSEH_info_poly1305_blocks_avx2_2
3971 .rva .LSEH_end_poly1305_blocks_avx2
3972 .rva .LSEH_info_poly1305_blocks_avx2_3
3974 $code.=<<___ if ($avx>2);
3975 .rva .LSEH_begin_poly1305_blocks_avx512
3976 .rva .LSEH_end_poly1305_blocks_avx512
3977 .rva .LSEH_info_poly1305_blocks_avx512
3982 .LSEH_info_poly1305_init:
3985 .rva .LSEH_begin_poly1305_init,.LSEH_begin_poly1305_init
3987 .LSEH_info_poly1305_blocks:
3990 .rva .Lblocks_body,.Lblocks_epilogue
3992 .LSEH_info_poly1305_emit:
3995 .rva .LSEH_begin_poly1305_emit,.LSEH_begin_poly1305_emit
3997 $code.=<<___ if ($avx);
3998 .LSEH_info_poly1305_blocks_avx_1:
4001 .rva .Lblocks_avx_body,.Lblocks_avx_epilogue # HandlerData[]
4003 .LSEH_info_poly1305_blocks_avx_2:
4006 .rva .Lbase2_64_avx_body,.Lbase2_64_avx_epilogue # HandlerData[]
4008 .LSEH_info_poly1305_blocks_avx_3:
4011 .rva .Ldo_avx_body,.Ldo_avx_epilogue # HandlerData[]
4013 .LSEH_info_poly1305_emit_avx:
4016 .rva .LSEH_begin_poly1305_emit_avx,.LSEH_begin_poly1305_emit_avx
4018 $code.=<<___ if ($avx>1);
4019 .LSEH_info_poly1305_blocks_avx2_1:
4022 .rva .Lblocks_avx2_body,.Lblocks_avx2_epilogue # HandlerData[]
4024 .LSEH_info_poly1305_blocks_avx2_2:
4027 .rva .Lbase2_64_avx2_body,.Lbase2_64_avx2_epilogue # HandlerData[]
4029 .LSEH_info_poly1305_blocks_avx2_3:
4032 .rva .Ldo_avx2_body,.Ldo_avx2_epilogue # HandlerData[]
4034 $code.=<<___ if ($avx>2);
4035 .LSEH_info_poly1305_blocks_avx512:
4038 .rva .Ldo_avx512_body,.Ldo_avx512_epilogue # HandlerData[]
4042 foreach (split('\n',$code)) {
4043 s/\`([^\`]*)\`/eval($1)/ge;
4044 s/%r([a-z]+)#d/%e$1/g;
4045 s/%r([0-9]+)#d/%r$1d/g;
4046 s/%x#%[yz]/%x/g or s/%y#%z/%y/g or s/%z#%[yz]/%z/g;