]> granicus.if.org Git - openssl/commitdiff
Add SHA256/512 for PowerISA 2.07.
authorAndy Polyakov <appro@openssl.org>
Tue, 1 Jul 2014 17:06:03 +0000 (19:06 +0200)
committerAndy Polyakov <appro@openssl.org>
Tue, 1 Jul 2014 17:09:11 +0000 (19:09 +0200)
crypto/sha/asm/sha512p8-ppc.pl [new file with mode: 0755]

diff --git a/crypto/sha/asm/sha512p8-ppc.pl b/crypto/sha/asm/sha512p8-ppc.pl
new file mode 100755 (executable)
index 0000000..cb0268c
--- /dev/null
@@ -0,0 +1,423 @@
+#!/usr/bin/env perl
+
+# ====================================================================
+# Written by Andy Polyakov <appro@openssl.org> for the OpenSSL
+# project. The module is, however, dual licensed under OpenSSL and
+# CRYPTOGAMS licenses depending on where you obtain it. For further
+# details see http://www.openssl.org/~appro/cryptogams/.
+# ====================================================================
+
+# SHA256/512 for PowerISA v2.07.
+#
+# Accurate performance measurements are problematic, because it's
+# always virtualized setup with possibly throttled processor.
+# Relative comparison is therefore more informative. This module is
+# ~60% faster than integer-only sha512-ppc.pl. To anchor to something
+# else, SHA256 is 16% slower than sha1-ppc.pl and 2.5x slower than
+# hardware-assisted aes-128-cbc encrypt. SHA512 is 33% faster than
+# sha1-ppc.pl and 1.6x slower than aes-128-cbc. Another interesting
+# result is degree of computational resources' utilization. POWER8 is
+# "massively multi-threaded chip" and difference between single- and
+# maximum multi-process benchmark results tells that utlization is
+# whooping 94%. For sha512-ppc.pl we get [not unimpressive] 84% and
+# for sha1-ppc.pl - 73%. 100% means that multi-process result equals
+# to single-process one, given that all threads end up on the same
+# physical core.
+
+$flavour=shift;
+$output =shift;
+
+if ($flavour =~ /64/) {
+       $SIZE_T=8;
+       $LRSAVE=2*$SIZE_T;
+       $STU="stdu";
+       $POP="ld";
+       $PUSH="std";
+} elsif ($flavour =~ /32/) {
+       $SIZE_T=4;
+       $LRSAVE=$SIZE_T;
+       $STU="stwu";
+       $POP="lwz";
+       $PUSH="stw";
+} else { die "nonsense $flavour"; }
+
+$LENDIAN=($flavour=~/le/);
+
+$0 =~ m/(.*[\/\\])[^\/\\]+$/; $dir=$1;
+( $xlate="${dir}ppc-xlate.pl" and -f $xlate ) or
+( $xlate="${dir}../../perlasm/ppc-xlate.pl" and -f $xlate) or
+die "can't locate ppc-xlate.pl";
+
+open STDOUT,"| $^X $xlate $flavour $output" || die "can't call $xlate: $!";
+
+if ($output =~ /512/) {
+       $bits=512;
+       $SZ=8;
+       $sz="d";
+       $rounds=80;
+} else {
+       $bits=256;
+       $SZ=4;
+       $sz="w";
+       $rounds=64;
+}
+
+$func="sha${bits}_block_p8";
+$FRAME=8*$SIZE_T;
+
+$sp ="r1";
+$toc="r2";
+$ctx="r3";
+$inp="r4";
+$num="r5";
+$Tbl="r6";
+$idx="r7";
+$lrsave="r8";
+$offload="r11";
+$vrsave="r12";
+($x00,$x10,$x20,$x30,$x40,$x50,$x60,$x70)=map("r$_",(0,10,26..31));
+
+@V=($A,$B,$C,$D,$E,$F,$G,$H)=map("v$_",(0..7));
+@X=map("v$_",(8..23));
+($Ki,$Func,$S0,$S1,$s0,$s1,$lemask)=map("v$_",(24..31));
+
+sub ROUND {
+my ($i,$a,$b,$c,$d,$e,$f,$g,$h)=@_;
+my $j=($i+1)%16;
+
+$code.=<<___           if ($i<15 && ($i%(16/$SZ))==(16/$SZ-1));
+       lvx_u           @X[$i+1],0,$inp         ; load X[i] in advance
+       addi            $inp,$inp,16
+___
+$code.=<<___           if ($i<16 && ($i%(16/$SZ)));
+       vsldoi          @X[$i],@X[$i-1],@X[$i-1],$SZ
+___
+$code.=<<___           if ($LENDIAN && $i<16 && ($i%(16/$SZ))==0);
+       vperm           @X[$i],@X[$i],@X[$i],$lemask
+___
+$code.=<<___;
+       `"vshasigma${sz}        $s0,@X[($j+1)%16],0,0"          if ($i>=15)`
+       vsel            $Func,$g,$f,$e          ; Ch(e,f,g)
+       vshasigma${sz}  $S1,$e,1,15             ; Sigma1(e)
+       vaddu${sz}m     $h,$h,@X[$i%16]         ; h+=X[i]
+       vshasigma${sz}  $S0,$a,1,0              ; Sigma0(a)
+       `"vshasigma${sz}        $s1,@X[($j+14)%16],0,15"        if ($i>=15)`
+       vaddu${sz}m     $h,$h,$Func             ; h+=Ch(e,f,g)
+       vxor            $Func,$a,$b
+       `"vaddu${sz}m           @X[$j],@X[$j],@X[($j+9)%16]"    if ($i>=15)`
+       vaddu${sz}m     $h,$h,$S1               ; h+=Sigma1(e)
+       vsel            $Func,$b,$c,$Func       ; Maj(a,b,c)
+       vaddu${sz}m     $g,$g,$Ki               ; future h+=K[i]
+       vaddu${sz}m     $d,$d,$h                ; d+=h
+       vaddu${sz}m     $S0,$S0,$Func           ; Sigma0(a)+Maj(a,b,c)
+       `"vaddu${sz}m           @X[$j],@X[$j],$s0"              if ($i>=15)`
+       lvx             $Ki,$idx,$Tbl           ; load next K[i]
+       addi            $idx,$idx,16
+       vaddu${sz}m     $h,$h,$S0               ; h+=Sigma0(a)+Maj(a,b,c)
+       `"vaddu${sz}m           @X[$j],@X[$j],$s1"              if ($i>=15)`
+___
+}
+
+$code=<<___;
+.machine       "any"
+.text
+
+.globl $func
+.align 6
+$func:
+       $STU            $sp,-`($FRAME+21*16+6*$SIZE_T)`($sp)
+       mflr            $lrsave
+       li              r10,`$FRAME+8*16+15`
+       li              r11,`$FRAME+8*16+31`
+       stvx            v20,r10,$sp             # ABI says so
+       addi            r10,r10,32
+       mfspr           $vrsave,256
+       stvx            v21,r11,$sp
+       addi            r11,r11,32
+       stvx            v22,r10,$sp
+       addi            r10,r10,32
+       stvx            v23,r11,$sp
+       addi            r11,r11,32
+       stvx            v24,r10,$sp
+       addi            r10,r10,32
+       stvx            v25,r11,$sp
+       addi            r11,r11,32
+       stvx            v26,r10,$sp
+       addi            r10,r10,32
+       stvx            v27,r11,$sp
+       addi            r11,r11,32
+       stvx            v28,r10,$sp
+       addi            r10,r10,32
+       stvx            v29,r11,$sp
+       addi            r11,r11,32
+       stvx            v30,r10,$sp
+       stvx            v31,r11,$sp
+       li              r11,-1
+       stw             $vrsave,`$FRAME+21*16-4`($sp)   # save vrsave
+       li              $x10,0x10
+       $PUSH           r26,`$FRAME+21*16+0*$SIZE_T`($sp)
+       li              $x20,0x20
+       $PUSH           r27,`$FRAME+21*16+1*$SIZE_T`($sp)
+       li              $x30,0x30
+       $PUSH           r28,`$FRAME+21*16+2*$SIZE_T`($sp)
+       li              $x40,0x40
+       $PUSH           r29,`$FRAME+21*16+3*$SIZE_T`($sp)
+       li              $x50,0x50
+       $PUSH           r30,`$FRAME+21*16+4*$SIZE_T`($sp)
+       li              $x60,0x60
+       $PUSH           r31,`$FRAME+21*16+5*$SIZE_T`($sp)
+       li              $x70,0x70
+       $PUSH           $lrsave,`$FRAME+21*16+6*$SIZE_T+$LRSAVE`($sp)
+       mtspr           256,r11
+
+       bl              LPICmeup
+       addi            $offload,$sp,$FRAME+15
+___
+$code.=<<___           if ($LENDIAN);
+       li              $idx,8
+       lvsl            $lemask,0,$idx
+       vspltisb        $Ki,0x0f
+       vxor            $lemask,$lemask,$Ki
+___
+$code.=<<___           if ($SZ==4);
+       lvx_4w          $A,$x00,$ctx
+       lvx_4w          $E,$x10,$ctx
+       vsldoi          $B,$A,$A,4              # unpack
+       vsldoi          $C,$A,$A,8
+       vsldoi          $D,$A,$A,12
+       vsldoi          $F,$E,$E,4
+       vsldoi          $G,$E,$E,8
+       vsldoi          $H,$E,$E,12
+___
+$code.=<<___           if ($SZ==8);
+       lvx_u           $A,$x00,$ctx
+       lvx_u           $C,$x10,$ctx
+       lvx_u           $E,$x20,$ctx
+       vsldoi          $B,$A,$A,8              # unpack
+       lvx_u           $G,$x30,$ctx
+       vsldoi          $D,$C,$C,8
+       vsldoi          $F,$E,$E,8
+       vsldoi          $H,$G,$G,8
+___
+$code.=<<___;
+       li              r0,`($rounds-16)/16`    # inner loop counter
+       b               Loop
+.align 5
+Loop:
+       lvx             $Ki,$x00,$Tbl
+       li              $idx,16
+       lvx_u           @X[0],0,$inp
+       addi            $inp,$inp,16
+       stvx            $A,$x00,$offload        # offload $A-$H
+       stvx            $B,$x10,$offload
+       stvx            $C,$x20,$offload
+       stvx            $D,$x30,$offload
+       stvx            $E,$x40,$offload
+       stvx            $F,$x50,$offload
+       stvx            $G,$x60,$offload
+       stvx            $H,$x70,$offload
+       vaddu${sz}m     $H,$H,$Ki               # h+K[i]
+       lvx             $Ki,$idx,$Tbl
+       addi            $idx,$idx,16
+___
+for ($i=0;$i<16;$i++)  { &ROUND($i,@V); unshift(@V,pop(@V)); }
+$code.=<<___;
+       mtctr           r0
+       b               L16_xx
+.align 5
+L16_xx:
+___
+for (;$i<32;$i++)      { &ROUND($i,@V); unshift(@V,pop(@V)); }
+$code.=<<___;
+       bdnz            L16_xx
+
+       lvx             @X[2],$x00,$offload
+       subic.          $num,$num,1
+       lvx             @X[3],$x10,$offload
+       vaddu${sz}m     $A,$A,@X[2]
+       lvx             @X[4],$x20,$offload
+       vaddu${sz}m     $B,$B,@X[3]
+       lvx             @X[5],$x30,$offload
+       vaddu${sz}m     $C,$C,@X[4]
+       lvx             @X[6],$x40,$offload
+       vaddu${sz}m     $D,$D,@X[5]
+       lvx             @X[7],$x50,$offload
+       vaddu${sz}m     $E,$E,@X[6]
+       lvx             @X[8],$x60,$offload
+       vaddu${sz}m     $F,$F,@X[7]
+       lvx             @X[9],$x70,$offload
+       vaddu${sz}m     $G,$G,@X[8]
+       vaddu${sz}m     $H,$H,@X[9]
+       bne             Loop
+___
+$code.=<<___           if ($SZ==4);
+       lvx             @X[0],$idx,$Tbl
+       addi            $idx,$idx,16
+       vperm           $A,$A,$B,$Ki            # pack the answer
+       lvx             @X[1],$idx,$Tbl
+       vperm           $E,$E,$F,$Ki
+       vperm           $A,$A,$C,@X[0]
+       vperm           $E,$E,$G,@X[0]
+       vperm           $A,$A,$D,@X[1]
+       vperm           $E,$E,$H,@X[1]
+       stvx_4w         $A,$x00,$ctx
+       stvx_4w         $E,$x10,$ctx
+___
+$code.=<<___           if ($SZ==8);
+       vperm           $A,$A,$B,$Ki            # pack the answer
+       vperm           $C,$C,$D,$Ki
+       vperm           $E,$E,$F,$Ki
+       vperm           $G,$G,$H,$Ki
+       stvx_u          $A,$x00,$ctx
+       stvx_u          $C,$x10,$ctx
+       stvx_u          $E,$x20,$ctx
+       stvx_u          $G,$x30,$ctx
+___
+$code.=<<___;
+       li              r10,`$FRAME+8*16+15`
+       mtlr            $lrsave
+       li              r11,`$FRAME+8*16+31`
+       mtspr           256,$vrsave
+       lvx             v20,r10,$sp             # ABI says so
+       addi            r10,r10,32
+       lvx             v21,r11,$sp
+       addi            r11,r11,32
+       lvx             v22,r10,$sp
+       addi            r10,r10,32
+       lvx             v23,r11,$sp
+       addi            r11,r11,32
+       lvx             v24,r10,$sp
+       addi            r10,r10,32
+       lvx             v25,r11,$sp
+       addi            r11,r11,32
+       lvx             v26,r10,$sp
+       addi            r10,r10,32
+       lvx             v27,r11,$sp
+       addi            r11,r11,32
+       lvx             v28,r10,$sp
+       addi            r10,r10,32
+       lvx             v29,r11,$sp
+       addi            r11,r11,32
+       lvx             v30,r10,$sp
+       lvx             v31,r11,$sp
+       $POP            r26,`$FRAME+21*16+0*$SIZE_T`($sp)
+       $POP            r27,`$FRAME+21*16+1*$SIZE_T`($sp)
+       $POP            r28,`$FRAME+21*16+2*$SIZE_T`($sp)
+       $POP            r29,`$FRAME+21*16+3*$SIZE_T`($sp)
+       $POP            r30,`$FRAME+21*16+4*$SIZE_T`($sp)
+       $POP            r31,`$FRAME+21*16+5*$SIZE_T`($sp)
+       addi            $sp,$sp,`$FRAME+21*16+6*$SIZE_T`
+       blr
+       .long           0
+       .byte           0,12,4,1,0x80,6,3,0
+       .long           0
+.size  $func,.-$func
+___
+
+# Ugly hack here, because PPC assembler syntax seem to vary too
+# much from platforms to platform...
+$code.=<<___;
+.align 6
+LPICmeup:
+       mflr    r0
+       bcl     20,31,\$+4
+       mflr    $Tbl    ; vvvvvv "distance" between . and 1st data entry
+       addi    $Tbl,$Tbl,`64-8`
+       mtlr    r0
+       blr
+       .long   0
+       .byte   0,12,0x14,0,0,0,0,0
+       .space  `64-9*4`
+___
+
+if ($SZ==8) {
+    local *table = sub {
+       foreach(@_) { $code.=".quad     $_,$_\n"; }
+    };
+    table(
+       "0x428a2f98d728ae22","0x7137449123ef65cd",
+       "0xb5c0fbcfec4d3b2f","0xe9b5dba58189dbbc",
+       "0x3956c25bf348b538","0x59f111f1b605d019",
+       "0x923f82a4af194f9b","0xab1c5ed5da6d8118",
+       "0xd807aa98a3030242","0x12835b0145706fbe",
+       "0x243185be4ee4b28c","0x550c7dc3d5ffb4e2",
+       "0x72be5d74f27b896f","0x80deb1fe3b1696b1",
+       "0x9bdc06a725c71235","0xc19bf174cf692694",
+       "0xe49b69c19ef14ad2","0xefbe4786384f25e3",
+       "0x0fc19dc68b8cd5b5","0x240ca1cc77ac9c65",
+       "0x2de92c6f592b0275","0x4a7484aa6ea6e483",
+       "0x5cb0a9dcbd41fbd4","0x76f988da831153b5",
+       "0x983e5152ee66dfab","0xa831c66d2db43210",
+       "0xb00327c898fb213f","0xbf597fc7beef0ee4",
+       "0xc6e00bf33da88fc2","0xd5a79147930aa725",
+       "0x06ca6351e003826f","0x142929670a0e6e70",
+       "0x27b70a8546d22ffc","0x2e1b21385c26c926",
+       "0x4d2c6dfc5ac42aed","0x53380d139d95b3df",
+       "0x650a73548baf63de","0x766a0abb3c77b2a8",
+       "0x81c2c92e47edaee6","0x92722c851482353b",
+       "0xa2bfe8a14cf10364","0xa81a664bbc423001",
+       "0xc24b8b70d0f89791","0xc76c51a30654be30",
+       "0xd192e819d6ef5218","0xd69906245565a910",
+       "0xf40e35855771202a","0x106aa07032bbd1b8",
+       "0x19a4c116b8d2d0c8","0x1e376c085141ab53",
+       "0x2748774cdf8eeb99","0x34b0bcb5e19b48a8",
+       "0x391c0cb3c5c95a63","0x4ed8aa4ae3418acb",
+       "0x5b9cca4f7763e373","0x682e6ff3d6b2b8a3",
+       "0x748f82ee5defb2fc","0x78a5636f43172f60",
+       "0x84c87814a1f0ab72","0x8cc702081a6439ec",
+       "0x90befffa23631e28","0xa4506cebde82bde9",
+       "0xbef9a3f7b2c67915","0xc67178f2e372532b",
+       "0xca273eceea26619c","0xd186b8c721c0c207",
+       "0xeada7dd6cde0eb1e","0xf57d4f7fee6ed178",
+       "0x06f067aa72176fba","0x0a637dc5a2c898a6",
+       "0x113f9804bef90dae","0x1b710b35131c471b",
+       "0x28db77f523047d84","0x32caab7b40c72493",
+       "0x3c9ebe0a15c9bebc","0x431d67c49c100d4c",
+       "0x4cc5d4becb3e42b6","0x597f299cfc657e2a",
+       "0x5fcb6fab3ad6faec","0x6c44198c4a475817","0");
+$code.=<<___   if (!$LENDIAN);
+.quad  0x0001020304050607,0x1011121314151617
+___
+$code.=<<___   if ($LENDIAN);  # quad-swapped
+.quad  0x1011121314151617,0x0001020304050607
+___
+} else {
+    local *table = sub {
+       foreach(@_) { $code.=".long     $_,$_,$_,$_\n"; }
+    };
+    table(
+       "0x428a2f98","0x71374491","0xb5c0fbcf","0xe9b5dba5",
+       "0x3956c25b","0x59f111f1","0x923f82a4","0xab1c5ed5",
+       "0xd807aa98","0x12835b01","0x243185be","0x550c7dc3",
+       "0x72be5d74","0x80deb1fe","0x9bdc06a7","0xc19bf174",
+       "0xe49b69c1","0xefbe4786","0x0fc19dc6","0x240ca1cc",
+       "0x2de92c6f","0x4a7484aa","0x5cb0a9dc","0x76f988da",
+       "0x983e5152","0xa831c66d","0xb00327c8","0xbf597fc7",
+       "0xc6e00bf3","0xd5a79147","0x06ca6351","0x14292967",
+       "0x27b70a85","0x2e1b2138","0x4d2c6dfc","0x53380d13",
+       "0x650a7354","0x766a0abb","0x81c2c92e","0x92722c85",
+       "0xa2bfe8a1","0xa81a664b","0xc24b8b70","0xc76c51a3",
+       "0xd192e819","0xd6990624","0xf40e3585","0x106aa070",
+       "0x19a4c116","0x1e376c08","0x2748774c","0x34b0bcb5",
+       "0x391c0cb3","0x4ed8aa4a","0x5b9cca4f","0x682e6ff3",
+       "0x748f82ee","0x78a5636f","0x84c87814","0x8cc70208",
+       "0x90befffa","0xa4506ceb","0xbef9a3f7","0xc67178f2","0");
+$code.=<<___   if (!$LENDIAN);
+.long  0x00010203,0x10111213,0x10111213,0x10111213
+.long  0x00010203,0x04050607,0x10111213,0x10111213
+.long  0x00010203,0x04050607,0x08090a0b,0x10111213
+___
+$code.=<<___   if ($LENDIAN);  # word-swapped
+.long  0x10111213,0x10111213,0x10111213,0x00010203
+.long  0x10111213,0x10111213,0x04050607,0x00010203
+.long  0x10111213,0x08090a0b,0x04050607,0x00010203
+___
+}
+$code.=<<___;
+.asciz "SHA${bits} for PowerISA 2.07, CRYPTOGAMS by <appro\@openssl.org>"
+.align 2
+___
+
+$code =~ s/\`([^\`]*)\`/eval $1/gem;
+print $code;
+close STDOUT;