crypto/modes/asm/ghashp8-ppc.pl

   1 #!/usr/bin/env perl
   2 #
   3 # ====================================================================
   4 # Written by Andy Polyakov <appro@openssl.org> for the OpenSSL
   5 # project. The module is, however, dual licensed under OpenSSL and
   6 # CRYPTOGAMS licenses depending on where you obtain it. For further
   7 # details see http://www.openssl.org/~appro/cryptogams/.
   8 # ====================================================================
   9 #
  10 # GHASH for for PowerISA v2.07.
  11 #
  12 # July 2014
  13 #
  14 # Accurate performance measurements are problematic, because it's
  15 # always virtualized setup with possibly throttled processor.
  16 # Relative comparison is therefore more informative. This initial
  17 # version is ~2.1x slower than hardware-assisted AES-128-CTR, ~12x
  18 # faster than "4-bit" integer-only compiler-generated 64-bit code.
  19 # "Initial version" means that there is room for futher improvement.
  20
  21 $flavour=shift;
  22 $output =shift;
  23
  24 if ($flavour =~ /64/) {
  25         $SIZE_T=8;
  26         $LRSAVE=2*$SIZE_T;
  27         $STU="stdu";
  28         $POP="ld";
  29         $PUSH="std";
  30 } elsif ($flavour =~ /32/) {
  31         $SIZE_T=4;
  32         $LRSAVE=$SIZE_T;
  33         $STU="stwu";
  34         $POP="lwz";
  35         $PUSH="stw";
  36 } else { die "nonsense $flavour"; }
  37
  38 $0 =~ m/(.*[\/\\])[^\/\\]+$/; $dir=$1;
  39 ( $xlate="${dir}ppc-xlate.pl" and -f $xlate ) or
  40 ( $xlate="${dir}../../perlasm/ppc-xlate.pl" and -f $xlate) or
  41 die "can't locate ppc-xlate.pl";
  42
  43 open STDOUT,"| $^X $xlate $flavour $output" || die "can't call $xlate: $!";
  44
  45 my ($Xip,$Htbl,$inp,$len)=map("r$_",(3..6));    # argument block
  46
  47 my ($Xl,$Xm,$Xh,$IN)=map("v$_",(0..3));
  48 my ($zero,$t0,$t1,$t2,$xC2,$H,$Hh,$Hl,$lemask)=map("v$_",(4..12));
  49 my $vrsave="r12";
  50
  51 $code=<<___;
  52 .machine        "any"
  53
  54 .text
  55
  56 .globl  .gcm_init_p8
  57 .align  5
  58 .gcm_init_p8:
  59         lis             r0,0xfff0
  60         li              r8,0x10
  61         mfspr           $vrsave,256
  62         li              r9,0x20
  63         mtspr           256,r0
  64         li              r10,0x30
  65         lvx_u           $H,0,r4                 # load H
  66
  67         vspltisb        $xC2,-16                # 0xf0
  68         vspltisb        $t0,1                   # one
  69         vaddubm         $xC2,$xC2,$xC2          # 0xe0
  70         vxor            $zero,$zero,$zero
  71         vor             $xC2,$xC2,$t0           # 0xe1
  72         vsldoi          $xC2,$xC2,$zero,15      # 0xe1...
  73         vsldoi          $t1,$zero,$t0,1         # ...1
  74         vaddubm         $xC2,$xC2,$xC2          # 0xc2...
  75         vspltisb        $t2,7
  76         vor             $xC2,$xC2,$t1           # 0xc2....01
  77         vspltb          $t1,$H,0                # most significant byte
  78         vsl             $H,$H,$t0               # H<<=1
  79         vsrab           $t1,$t1,$t2             # broadcast carry bit
  80         vand            $t1,$t1,$xC2
  81         vxor            $H,$H,$t1               # twisted H
  82
  83         vsldoi          $H,$H,$H,8              # twist even more ...
  84         vsldoi          $xC2,$zero,$xC2,8       # 0xc2.0
  85         vsldoi          $Hl,$zero,$H,8          # ... and split
  86         vsldoi          $Hh,$H,$zero,8
  87
  88         stvx_u          $xC2,0,r3               # save pre-computed table
  89         stvx_u          $Hl,r8,r3
  90         stvx_u          $H, r9,r3
  91         stvx_u          $Hh,r10,r3
  92
  93         mtspr           256,$vrsave
  94         blr
  95         .long           0
  96         .byte           0,12,0x14,0,0,0,2,0
  97         .long           0
  98 .size   .gcm_init_p8,.-.gcm_init_p8
  99
 100 .globl  .gcm_gmult_p8
 101 .align  5
 102 .gcm_gmult_p8:
 103         lis             r0,0xfff8
 104         li              r8,0x10
 105         mfspr           $vrsave,256
 106         li              r9,0x20
 107         mtspr           256,r0
 108         li              r10,0x30
 109         lvx_u           $IN,0,$Xip              # load Xi
 110
 111         lvx_u           $Hl,r8,$Htbl            # load pre-computed table
 112          le?lvsl        $lemask,r0,r0
 113         lvx_u           $H, r9,$Htbl
 114          le?vspltisb    $t0,0x07
 115         lvx_u           $Hh,r10,$Htbl
 116          le?vxor        $lemask,$lemask,$t0
 117         lvx_u           $xC2,0,$Htbl
 118          le?vperm       $IN,$IN,$IN,$lemask
 119         vxor            $zero,$zero,$zero
 120
 121         vpmsumd         $Xl,$IN,$Hl             # H.lo·Xi.lo
 122         vpmsumd         $Xm,$IN,$H              # H.hi·Xi.lo+H.lo·Xi.hi
 123         vpmsumd         $Xh,$IN,$Hh             # H.hi·Xi.hi
 124
 125         vpmsumd         $t2,$Xl,$xC2            # 1st phase
 126
 127         vsldoi          $t0,$Xm,$zero,8
 128         vsldoi          $t1,$zero,$Xm,8
 129         vxor            $Xl,$Xl,$t0
 130         vxor            $Xh,$Xh,$t1
 131
 132         vsldoi          $Xl,$Xl,$Xl,8
 133         vxor            $Xl,$Xl,$t2
 134
 135         vsldoi          $t1,$Xl,$Xl,8           # 2nd phase
 136         vpmsumd         $Xl,$Xl,$xC2
 137         vxor            $t1,$t1,$Xh
 138         vxor            $Xl,$Xl,$t1
 139
 140         le?vperm        $Xl,$Xl,$Xl,$lemask
 141         stvx_u          $Xl,0,$Xip              # write out Xi
 142
 143         mtspr           256,$vrsave
 144         blr
 145         .long           0
 146         .byte           0,12,0x14,0,0,0,2,0
 147         .long           0
 148 .size   .gcm_gmult_p8,.-.gcm_gmult_p8
 149
 150 .globl  .gcm_ghash_p8
 151 .align  5
 152 .gcm_ghash_p8:
 153         lis             r0,0xfff8
 154         li              r8,0x10
 155         mfspr           $vrsave,256
 156         li              r9,0x20
 157         mtspr           256,r0
 158         li              r10,0x30
 159         lvx_u           $Xl,0,$Xip              # load Xi
 160
 161         lvx_u           $Hl,r8,$Htbl            # load pre-computed table
 162          le?lvsl        $lemask,r0,r0
 163         lvx_u           $H, r9,$Htbl
 164          le?vspltisb    $t0,0x07
 165         lvx_u           $Hh,r10,$Htbl
 166          le?vxor        $lemask,$lemask,$t0
 167         lvx_u           $xC2,0,$Htbl
 168          le?vperm       $Xl,$Xl,$Xl,$lemask
 169         vxor            $zero,$zero,$zero
 170
 171         lvx_u           $IN,0,$inp
 172         addi            $inp,$inp,16
 173         subi            $len,$len,16
 174          le?vperm       $IN,$IN,$IN,$lemask
 175         vxor            $IN,$IN,$Xl
 176         b               Loop
 177
 178 .align  5
 179 Loop:
 180          subic          $len,$len,16
 181         vpmsumd         $Xl,$IN,$Hl             # H.lo·Xi.lo
 182          subfe.         r0,r0,r0                # borrow?-1:0
 183         vpmsumd         $Xm,$IN,$H              # H.hi·Xi.lo+H.lo·Xi.hi
 184          and            r0,r0,$len
 185         vpmsumd         $Xh,$IN,$Hh             # H.hi·Xi.hi
 186          add            $inp,$inp,r0
 187
 188         vpmsumd         $t2,$Xl,$xC2            # 1st phase
 189
 190         vsldoi          $t0,$Xm,$zero,8
 191         vsldoi          $t1,$zero,$Xm,8
 192         vxor            $Xl,$Xl,$t0
 193         vxor            $Xh,$Xh,$t1
 194
 195         vsldoi          $Xl,$Xl,$Xl,8
 196         vxor            $Xl,$Xl,$t2
 197          lvx_u          $IN,0,$inp
 198          addi           $inp,$inp,16
 199
 200         vsldoi          $t1,$Xl,$Xl,8           # 2nd phase
 201         vpmsumd         $Xl,$Xl,$xC2
 202          le?vperm       $IN,$IN,$IN,$lemask
 203         vxor            $t1,$t1,$Xh
 204         vxor            $IN,$IN,$t1
 205         vxor            $IN,$IN,$Xl
 206         beq             Loop                    # did $len-=16 borrow?
 207
 208         vxor            $Xl,$Xl,$t1
 209         le?vperm        $Xl,$Xl,$Xl,$lemask
 210         stvx_u          $Xl,0,$Xip              # write out Xi
 211
 212         mtspr           256,$vrsave
 213         blr
 214         .long           0
 215         .byte           0,12,0x14,0,0,0,4,0
 216         .long           0
 217 .size   .gcm_ghash_p8,.-.gcm_ghash_p8
 218
 219 .asciz  "GHASH for PowerISA 2.07, CRYPTOGAMS by <appro\@openssl.org>"
 220 .align  2
 221 ___
 222
 223 foreach (split("\n",$code)) {
 224         if ($flavour =~ /le$/o) {       # little-endian
 225             s/le\?//o           or
 226             s/be\?/#be#/o;
 227         } else {
 228             s/le\?/#le#/o       or
 229             s/be\?//o;
 230         }
 231         print $_,"\n";
 232 }
 233
 234 close STDOUT; # enforce flush