1238384Sjkim#!/usr/bin/env perl 2238384Sjkim 3238384Sjkim# ==================================================================== 4290207Sjkim# Written by Andy Polyakov <appro@openssl.org> for the OpenSSL 5238384Sjkim# project. The module is, however, dual licensed under OpenSSL and 6238384Sjkim# CRYPTOGAMS licenses depending on where you obtain it. For further 7238384Sjkim# details see http://www.openssl.org/~appro/cryptogams/. 8238384Sjkim# ==================================================================== 9238384Sjkim 10238384Sjkim# SHA512 block procedure for ARMv4. September 2007. 11238384Sjkim 12238384Sjkim# This code is ~4.5 (four and a half) times faster than code generated 13238384Sjkim# by gcc 3.4 and it spends ~72 clock cycles per byte [on single-issue 14238384Sjkim# Xscale PXA250 core]. 15238384Sjkim# 16238384Sjkim# July 2010. 17238384Sjkim# 18238384Sjkim# Rescheduling for dual-issue pipeline resulted in 6% improvement on 19238384Sjkim# Cortex A8 core and ~40 cycles per processed byte. 20238384Sjkim 21238384Sjkim# February 2011. 22238384Sjkim# 23238384Sjkim# Profiler-assisted and platform-specific optimization resulted in 7% 24238384Sjkim# improvement on Coxtex A8 core and ~38 cycles per byte. 25238384Sjkim 26238384Sjkim# March 2011. 27238384Sjkim# 28238384Sjkim# Add NEON implementation. On Cortex A8 it was measured to process 29290207Sjkim# one byte in 23.3 cycles or ~60% faster than integer-only code. 30238384Sjkim 31290207Sjkim# August 2012. 32290207Sjkim# 33290207Sjkim# Improve NEON performance by 12% on Snapdragon S4. In absolute 34290207Sjkim# terms it's 22.6 cycles per byte, which is disappointing result. 35290207Sjkim# Technical writers asserted that 3-way S4 pipeline can sustain 36290207Sjkim# multiple NEON instructions per cycle, but dual NEON issue could 37290207Sjkim# not be observed, and for NEON-only sequences IPC(*) was found to 38290207Sjkim# be limited by 1:-( 0.33 and 0.66 were measured for sequences with 39290207Sjkim# ILPs(*) of 1 and 2 respectively. This in turn means that you can 40290207Sjkim# even find yourself striving, as I did here, for achieving IPC 41290207Sjkim# adequate to one delivered by Cortex A8 [for reference, it's 42290207Sjkim# 0.5 for ILP of 1, and 1 for higher ILPs]. 43290207Sjkim# 44290207Sjkim# (*) ILP, instruction-level parallelism, how many instructions 45290207Sjkim# *can* execute at the same time. IPC, instructions per cycle, 46290207Sjkim# indicates how many instructions actually execute. 47290207Sjkim 48238384Sjkim# Byte order [in]dependence. ========================================= 49238384Sjkim# 50238384Sjkim# Originally caller was expected to maintain specific *dword* order in 51238384Sjkim# h[0-7], namely with most significant dword at *lower* address, which 52238384Sjkim# was reflected in below two parameters as 0 and 4. Now caller is 53238384Sjkim# expected to maintain native byte order for whole 64-bit values. 54238384Sjkim$hi="HI"; 55238384Sjkim$lo="LO"; 56238384Sjkim# ==================================================================== 57238384Sjkim 58238384Sjkimwhile (($output=shift) && ($output!~/^\w[\w\-]*\.\w+$/)) {} 59238384Sjkimopen STDOUT,">$output"; 60238384Sjkim 61238384Sjkim$ctx="r0"; # parameter block 62238384Sjkim$inp="r1"; 63238384Sjkim$len="r2"; 64238384Sjkim 65238384Sjkim$Tlo="r3"; 66238384Sjkim$Thi="r4"; 67238384Sjkim$Alo="r5"; 68238384Sjkim$Ahi="r6"; 69238384Sjkim$Elo="r7"; 70238384Sjkim$Ehi="r8"; 71238384Sjkim$t0="r9"; 72238384Sjkim$t1="r10"; 73238384Sjkim$t2="r11"; 74238384Sjkim$t3="r12"; 75238384Sjkim############ r13 is stack pointer 76238384Sjkim$Ktbl="r14"; 77238384Sjkim############ r15 is program counter 78238384Sjkim 79238384Sjkim$Aoff=8*0; 80238384Sjkim$Boff=8*1; 81238384Sjkim$Coff=8*2; 82238384Sjkim$Doff=8*3; 83238384Sjkim$Eoff=8*4; 84238384Sjkim$Foff=8*5; 85238384Sjkim$Goff=8*6; 86238384Sjkim$Hoff=8*7; 87238384Sjkim$Xoff=8*8; 88238384Sjkim 89238384Sjkimsub BODY_00_15() { 90238384Sjkimmy $magic = shift; 91238384Sjkim$code.=<<___; 92238384Sjkim @ Sigma1(x) (ROTR((x),14) ^ ROTR((x),18) ^ ROTR((x),41)) 93238384Sjkim @ LO lo>>14^hi<<18 ^ lo>>18^hi<<14 ^ hi>>9^lo<<23 94238384Sjkim @ HI hi>>14^lo<<18 ^ hi>>18^lo<<14 ^ lo>>9^hi<<23 95238384Sjkim mov $t0,$Elo,lsr#14 96238384Sjkim str $Tlo,[sp,#$Xoff+0] 97238384Sjkim mov $t1,$Ehi,lsr#14 98238384Sjkim str $Thi,[sp,#$Xoff+4] 99238384Sjkim eor $t0,$t0,$Ehi,lsl#18 100238384Sjkim ldr $t2,[sp,#$Hoff+0] @ h.lo 101238384Sjkim eor $t1,$t1,$Elo,lsl#18 102238384Sjkim ldr $t3,[sp,#$Hoff+4] @ h.hi 103238384Sjkim eor $t0,$t0,$Elo,lsr#18 104238384Sjkim eor $t1,$t1,$Ehi,lsr#18 105238384Sjkim eor $t0,$t0,$Ehi,lsl#14 106238384Sjkim eor $t1,$t1,$Elo,lsl#14 107238384Sjkim eor $t0,$t0,$Ehi,lsr#9 108238384Sjkim eor $t1,$t1,$Elo,lsr#9 109238384Sjkim eor $t0,$t0,$Elo,lsl#23 110238384Sjkim eor $t1,$t1,$Ehi,lsl#23 @ Sigma1(e) 111238384Sjkim adds $Tlo,$Tlo,$t0 112238384Sjkim ldr $t0,[sp,#$Foff+0] @ f.lo 113238384Sjkim adc $Thi,$Thi,$t1 @ T += Sigma1(e) 114238384Sjkim ldr $t1,[sp,#$Foff+4] @ f.hi 115238384Sjkim adds $Tlo,$Tlo,$t2 116238384Sjkim ldr $t2,[sp,#$Goff+0] @ g.lo 117238384Sjkim adc $Thi,$Thi,$t3 @ T += h 118238384Sjkim ldr $t3,[sp,#$Goff+4] @ g.hi 119238384Sjkim 120238384Sjkim eor $t0,$t0,$t2 121238384Sjkim str $Elo,[sp,#$Eoff+0] 122238384Sjkim eor $t1,$t1,$t3 123238384Sjkim str $Ehi,[sp,#$Eoff+4] 124238384Sjkim and $t0,$t0,$Elo 125238384Sjkim str $Alo,[sp,#$Aoff+0] 126238384Sjkim and $t1,$t1,$Ehi 127238384Sjkim str $Ahi,[sp,#$Aoff+4] 128238384Sjkim eor $t0,$t0,$t2 129238384Sjkim ldr $t2,[$Ktbl,#$lo] @ K[i].lo 130238384Sjkim eor $t1,$t1,$t3 @ Ch(e,f,g) 131238384Sjkim ldr $t3,[$Ktbl,#$hi] @ K[i].hi 132238384Sjkim 133238384Sjkim adds $Tlo,$Tlo,$t0 134238384Sjkim ldr $Elo,[sp,#$Doff+0] @ d.lo 135238384Sjkim adc $Thi,$Thi,$t1 @ T += Ch(e,f,g) 136238384Sjkim ldr $Ehi,[sp,#$Doff+4] @ d.hi 137238384Sjkim adds $Tlo,$Tlo,$t2 138238384Sjkim and $t0,$t2,#0xff 139238384Sjkim adc $Thi,$Thi,$t3 @ T += K[i] 140238384Sjkim adds $Elo,$Elo,$Tlo 141238384Sjkim ldr $t2,[sp,#$Boff+0] @ b.lo 142238384Sjkim adc $Ehi,$Ehi,$Thi @ d += T 143238384Sjkim teq $t0,#$magic 144238384Sjkim 145238384Sjkim ldr $t3,[sp,#$Coff+0] @ c.lo 146238384Sjkim orreq $Ktbl,$Ktbl,#1 147238384Sjkim @ Sigma0(x) (ROTR((x),28) ^ ROTR((x),34) ^ ROTR((x),39)) 148238384Sjkim @ LO lo>>28^hi<<4 ^ hi>>2^lo<<30 ^ hi>>7^lo<<25 149238384Sjkim @ HI hi>>28^lo<<4 ^ lo>>2^hi<<30 ^ lo>>7^hi<<25 150238384Sjkim mov $t0,$Alo,lsr#28 151238384Sjkim mov $t1,$Ahi,lsr#28 152238384Sjkim eor $t0,$t0,$Ahi,lsl#4 153238384Sjkim eor $t1,$t1,$Alo,lsl#4 154238384Sjkim eor $t0,$t0,$Ahi,lsr#2 155238384Sjkim eor $t1,$t1,$Alo,lsr#2 156238384Sjkim eor $t0,$t0,$Alo,lsl#30 157238384Sjkim eor $t1,$t1,$Ahi,lsl#30 158238384Sjkim eor $t0,$t0,$Ahi,lsr#7 159238384Sjkim eor $t1,$t1,$Alo,lsr#7 160238384Sjkim eor $t0,$t0,$Alo,lsl#25 161238384Sjkim eor $t1,$t1,$Ahi,lsl#25 @ Sigma0(a) 162238384Sjkim adds $Tlo,$Tlo,$t0 163238384Sjkim and $t0,$Alo,$t2 164238384Sjkim adc $Thi,$Thi,$t1 @ T += Sigma0(a) 165238384Sjkim 166238384Sjkim ldr $t1,[sp,#$Boff+4] @ b.hi 167238384Sjkim orr $Alo,$Alo,$t2 168238384Sjkim ldr $t2,[sp,#$Coff+4] @ c.hi 169238384Sjkim and $Alo,$Alo,$t3 170238384Sjkim and $t3,$Ahi,$t1 171238384Sjkim orr $Ahi,$Ahi,$t1 172238384Sjkim orr $Alo,$Alo,$t0 @ Maj(a,b,c).lo 173238384Sjkim and $Ahi,$Ahi,$t2 174238384Sjkim adds $Alo,$Alo,$Tlo 175238384Sjkim orr $Ahi,$Ahi,$t3 @ Maj(a,b,c).hi 176238384Sjkim sub sp,sp,#8 177238384Sjkim adc $Ahi,$Ahi,$Thi @ h += T 178238384Sjkim tst $Ktbl,#1 179238384Sjkim add $Ktbl,$Ktbl,#8 180238384Sjkim___ 181238384Sjkim} 182238384Sjkim$code=<<___; 183238384Sjkim#include "arm_arch.h" 184238384Sjkim#ifdef __ARMEL__ 185238384Sjkim# define LO 0 186238384Sjkim# define HI 4 187238384Sjkim# define WORD64(hi0,lo0,hi1,lo1) .word lo0,hi0, lo1,hi1 188238384Sjkim#else 189238384Sjkim# define HI 0 190238384Sjkim# define LO 4 191238384Sjkim# define WORD64(hi0,lo0,hi1,lo1) .word hi0,lo0, hi1,lo1 192238384Sjkim#endif 193238384Sjkim 194238384Sjkim.text 195238384Sjkim.code 32 196238384Sjkim.type K512,%object 197238384Sjkim.align 5 198238384SjkimK512: 199238384SjkimWORD64(0x428a2f98,0xd728ae22, 0x71374491,0x23ef65cd) 200238384SjkimWORD64(0xb5c0fbcf,0xec4d3b2f, 0xe9b5dba5,0x8189dbbc) 201238384SjkimWORD64(0x3956c25b,0xf348b538, 0x59f111f1,0xb605d019) 202238384SjkimWORD64(0x923f82a4,0xaf194f9b, 0xab1c5ed5,0xda6d8118) 203238384SjkimWORD64(0xd807aa98,0xa3030242, 0x12835b01,0x45706fbe) 204238384SjkimWORD64(0x243185be,0x4ee4b28c, 0x550c7dc3,0xd5ffb4e2) 205238384SjkimWORD64(0x72be5d74,0xf27b896f, 0x80deb1fe,0x3b1696b1) 206238384SjkimWORD64(0x9bdc06a7,0x25c71235, 0xc19bf174,0xcf692694) 207238384SjkimWORD64(0xe49b69c1,0x9ef14ad2, 0xefbe4786,0x384f25e3) 208238384SjkimWORD64(0x0fc19dc6,0x8b8cd5b5, 0x240ca1cc,0x77ac9c65) 209238384SjkimWORD64(0x2de92c6f,0x592b0275, 0x4a7484aa,0x6ea6e483) 210238384SjkimWORD64(0x5cb0a9dc,0xbd41fbd4, 0x76f988da,0x831153b5) 211238384SjkimWORD64(0x983e5152,0xee66dfab, 0xa831c66d,0x2db43210) 212238384SjkimWORD64(0xb00327c8,0x98fb213f, 0xbf597fc7,0xbeef0ee4) 213238384SjkimWORD64(0xc6e00bf3,0x3da88fc2, 0xd5a79147,0x930aa725) 214238384SjkimWORD64(0x06ca6351,0xe003826f, 0x14292967,0x0a0e6e70) 215238384SjkimWORD64(0x27b70a85,0x46d22ffc, 0x2e1b2138,0x5c26c926) 216238384SjkimWORD64(0x4d2c6dfc,0x5ac42aed, 0x53380d13,0x9d95b3df) 217238384SjkimWORD64(0x650a7354,0x8baf63de, 0x766a0abb,0x3c77b2a8) 218238384SjkimWORD64(0x81c2c92e,0x47edaee6, 0x92722c85,0x1482353b) 219238384SjkimWORD64(0xa2bfe8a1,0x4cf10364, 0xa81a664b,0xbc423001) 220238384SjkimWORD64(0xc24b8b70,0xd0f89791, 0xc76c51a3,0x0654be30) 221238384SjkimWORD64(0xd192e819,0xd6ef5218, 0xd6990624,0x5565a910) 222238384SjkimWORD64(0xf40e3585,0x5771202a, 0x106aa070,0x32bbd1b8) 223238384SjkimWORD64(0x19a4c116,0xb8d2d0c8, 0x1e376c08,0x5141ab53) 224238384SjkimWORD64(0x2748774c,0xdf8eeb99, 0x34b0bcb5,0xe19b48a8) 225238384SjkimWORD64(0x391c0cb3,0xc5c95a63, 0x4ed8aa4a,0xe3418acb) 226238384SjkimWORD64(0x5b9cca4f,0x7763e373, 0x682e6ff3,0xd6b2b8a3) 227238384SjkimWORD64(0x748f82ee,0x5defb2fc, 0x78a5636f,0x43172f60) 228238384SjkimWORD64(0x84c87814,0xa1f0ab72, 0x8cc70208,0x1a6439ec) 229238384SjkimWORD64(0x90befffa,0x23631e28, 0xa4506ceb,0xde82bde9) 230238384SjkimWORD64(0xbef9a3f7,0xb2c67915, 0xc67178f2,0xe372532b) 231238384SjkimWORD64(0xca273ece,0xea26619c, 0xd186b8c7,0x21c0c207) 232238384SjkimWORD64(0xeada7dd6,0xcde0eb1e, 0xf57d4f7f,0xee6ed178) 233238384SjkimWORD64(0x06f067aa,0x72176fba, 0x0a637dc5,0xa2c898a6) 234238384SjkimWORD64(0x113f9804,0xbef90dae, 0x1b710b35,0x131c471b) 235238384SjkimWORD64(0x28db77f5,0x23047d84, 0x32caab7b,0x40c72493) 236238384SjkimWORD64(0x3c9ebe0a,0x15c9bebc, 0x431d67c4,0x9c100d4c) 237238384SjkimWORD64(0x4cc5d4be,0xcb3e42b6, 0x597f299c,0xfc657e2a) 238238384SjkimWORD64(0x5fcb6fab,0x3ad6faec, 0x6c44198c,0x4a475817) 239238384Sjkim.size K512,.-K512 240290207Sjkim#if __ARM_MAX_ARCH__>=7 241238384Sjkim.LOPENSSL_armcap: 242238384Sjkim.word OPENSSL_armcap_P-sha512_block_data_order 243238384Sjkim.skip 32-4 244290207Sjkim#else 245290207Sjkim.skip 32 246290207Sjkim#endif 247238384Sjkim 248238384Sjkim.global sha512_block_data_order 249238384Sjkim.type sha512_block_data_order,%function 250238384Sjkimsha512_block_data_order: 251238384Sjkim sub r3,pc,#8 @ sha512_block_data_order 252238384Sjkim add $len,$inp,$len,lsl#7 @ len to point at the end of inp 253290207Sjkim#if __ARM_MAX_ARCH__>=7 254238384Sjkim ldr r12,.LOPENSSL_armcap 255238384Sjkim ldr r12,[r3,r12] @ OPENSSL_armcap_P 256238384Sjkim tst r12,#1 257238384Sjkim bne .LNEON 258238384Sjkim#endif 259238384Sjkim stmdb sp!,{r4-r12,lr} 260238384Sjkim sub $Ktbl,r3,#672 @ K512 261238384Sjkim sub sp,sp,#9*8 262238384Sjkim 263238384Sjkim ldr $Elo,[$ctx,#$Eoff+$lo] 264238384Sjkim ldr $Ehi,[$ctx,#$Eoff+$hi] 265238384Sjkim ldr $t0, [$ctx,#$Goff+$lo] 266238384Sjkim ldr $t1, [$ctx,#$Goff+$hi] 267238384Sjkim ldr $t2, [$ctx,#$Hoff+$lo] 268238384Sjkim ldr $t3, [$ctx,#$Hoff+$hi] 269238384Sjkim.Loop: 270238384Sjkim str $t0, [sp,#$Goff+0] 271238384Sjkim str $t1, [sp,#$Goff+4] 272238384Sjkim str $t2, [sp,#$Hoff+0] 273238384Sjkim str $t3, [sp,#$Hoff+4] 274238384Sjkim ldr $Alo,[$ctx,#$Aoff+$lo] 275238384Sjkim ldr $Ahi,[$ctx,#$Aoff+$hi] 276238384Sjkim ldr $Tlo,[$ctx,#$Boff+$lo] 277238384Sjkim ldr $Thi,[$ctx,#$Boff+$hi] 278238384Sjkim ldr $t0, [$ctx,#$Coff+$lo] 279238384Sjkim ldr $t1, [$ctx,#$Coff+$hi] 280238384Sjkim ldr $t2, [$ctx,#$Doff+$lo] 281238384Sjkim ldr $t3, [$ctx,#$Doff+$hi] 282238384Sjkim str $Tlo,[sp,#$Boff+0] 283238384Sjkim str $Thi,[sp,#$Boff+4] 284238384Sjkim str $t0, [sp,#$Coff+0] 285238384Sjkim str $t1, [sp,#$Coff+4] 286238384Sjkim str $t2, [sp,#$Doff+0] 287238384Sjkim str $t3, [sp,#$Doff+4] 288238384Sjkim ldr $Tlo,[$ctx,#$Foff+$lo] 289238384Sjkim ldr $Thi,[$ctx,#$Foff+$hi] 290238384Sjkim str $Tlo,[sp,#$Foff+0] 291238384Sjkim str $Thi,[sp,#$Foff+4] 292238384Sjkim 293238384Sjkim.L00_15: 294238384Sjkim#if __ARM_ARCH__<7 295238384Sjkim ldrb $Tlo,[$inp,#7] 296238384Sjkim ldrb $t0, [$inp,#6] 297238384Sjkim ldrb $t1, [$inp,#5] 298238384Sjkim ldrb $t2, [$inp,#4] 299238384Sjkim ldrb $Thi,[$inp,#3] 300238384Sjkim ldrb $t3, [$inp,#2] 301238384Sjkim orr $Tlo,$Tlo,$t0,lsl#8 302238384Sjkim ldrb $t0, [$inp,#1] 303238384Sjkim orr $Tlo,$Tlo,$t1,lsl#16 304238384Sjkim ldrb $t1, [$inp],#8 305238384Sjkim orr $Tlo,$Tlo,$t2,lsl#24 306238384Sjkim orr $Thi,$Thi,$t3,lsl#8 307238384Sjkim orr $Thi,$Thi,$t0,lsl#16 308238384Sjkim orr $Thi,$Thi,$t1,lsl#24 309238384Sjkim#else 310238384Sjkim ldr $Tlo,[$inp,#4] 311238384Sjkim ldr $Thi,[$inp],#8 312238384Sjkim#ifdef __ARMEL__ 313238384Sjkim rev $Tlo,$Tlo 314238384Sjkim rev $Thi,$Thi 315238384Sjkim#endif 316238384Sjkim#endif 317238384Sjkim___ 318238384Sjkim &BODY_00_15(0x94); 319238384Sjkim$code.=<<___; 320238384Sjkim tst $Ktbl,#1 321238384Sjkim beq .L00_15 322238384Sjkim ldr $t0,[sp,#`$Xoff+8*(16-1)`+0] 323238384Sjkim ldr $t1,[sp,#`$Xoff+8*(16-1)`+4] 324238384Sjkim bic $Ktbl,$Ktbl,#1 325238384Sjkim.L16_79: 326238384Sjkim @ sigma0(x) (ROTR((x),1) ^ ROTR((x),8) ^ ((x)>>7)) 327238384Sjkim @ LO lo>>1^hi<<31 ^ lo>>8^hi<<24 ^ lo>>7^hi<<25 328238384Sjkim @ HI hi>>1^lo<<31 ^ hi>>8^lo<<24 ^ hi>>7 329238384Sjkim mov $Tlo,$t0,lsr#1 330238384Sjkim ldr $t2,[sp,#`$Xoff+8*(16-14)`+0] 331238384Sjkim mov $Thi,$t1,lsr#1 332238384Sjkim ldr $t3,[sp,#`$Xoff+8*(16-14)`+4] 333238384Sjkim eor $Tlo,$Tlo,$t1,lsl#31 334238384Sjkim eor $Thi,$Thi,$t0,lsl#31 335238384Sjkim eor $Tlo,$Tlo,$t0,lsr#8 336238384Sjkim eor $Thi,$Thi,$t1,lsr#8 337238384Sjkim eor $Tlo,$Tlo,$t1,lsl#24 338238384Sjkim eor $Thi,$Thi,$t0,lsl#24 339238384Sjkim eor $Tlo,$Tlo,$t0,lsr#7 340238384Sjkim eor $Thi,$Thi,$t1,lsr#7 341238384Sjkim eor $Tlo,$Tlo,$t1,lsl#25 342238384Sjkim 343238384Sjkim @ sigma1(x) (ROTR((x),19) ^ ROTR((x),61) ^ ((x)>>6)) 344238384Sjkim @ LO lo>>19^hi<<13 ^ hi>>29^lo<<3 ^ lo>>6^hi<<26 345238384Sjkim @ HI hi>>19^lo<<13 ^ lo>>29^hi<<3 ^ hi>>6 346238384Sjkim mov $t0,$t2,lsr#19 347238384Sjkim mov $t1,$t3,lsr#19 348238384Sjkim eor $t0,$t0,$t3,lsl#13 349238384Sjkim eor $t1,$t1,$t2,lsl#13 350238384Sjkim eor $t0,$t0,$t3,lsr#29 351238384Sjkim eor $t1,$t1,$t2,lsr#29 352238384Sjkim eor $t0,$t0,$t2,lsl#3 353238384Sjkim eor $t1,$t1,$t3,lsl#3 354238384Sjkim eor $t0,$t0,$t2,lsr#6 355238384Sjkim eor $t1,$t1,$t3,lsr#6 356238384Sjkim ldr $t2,[sp,#`$Xoff+8*(16-9)`+0] 357238384Sjkim eor $t0,$t0,$t3,lsl#26 358238384Sjkim 359238384Sjkim ldr $t3,[sp,#`$Xoff+8*(16-9)`+4] 360238384Sjkim adds $Tlo,$Tlo,$t0 361238384Sjkim ldr $t0,[sp,#`$Xoff+8*16`+0] 362238384Sjkim adc $Thi,$Thi,$t1 363238384Sjkim 364238384Sjkim ldr $t1,[sp,#`$Xoff+8*16`+4] 365238384Sjkim adds $Tlo,$Tlo,$t2 366238384Sjkim adc $Thi,$Thi,$t3 367238384Sjkim adds $Tlo,$Tlo,$t0 368238384Sjkim adc $Thi,$Thi,$t1 369238384Sjkim___ 370238384Sjkim &BODY_00_15(0x17); 371238384Sjkim$code.=<<___; 372238384Sjkim ldreq $t0,[sp,#`$Xoff+8*(16-1)`+0] 373238384Sjkim ldreq $t1,[sp,#`$Xoff+8*(16-1)`+4] 374238384Sjkim beq .L16_79 375238384Sjkim bic $Ktbl,$Ktbl,#1 376238384Sjkim 377238384Sjkim ldr $Tlo,[sp,#$Boff+0] 378238384Sjkim ldr $Thi,[sp,#$Boff+4] 379238384Sjkim ldr $t0, [$ctx,#$Aoff+$lo] 380238384Sjkim ldr $t1, [$ctx,#$Aoff+$hi] 381238384Sjkim ldr $t2, [$ctx,#$Boff+$lo] 382238384Sjkim ldr $t3, [$ctx,#$Boff+$hi] 383238384Sjkim adds $t0,$Alo,$t0 384238384Sjkim str $t0, [$ctx,#$Aoff+$lo] 385238384Sjkim adc $t1,$Ahi,$t1 386238384Sjkim str $t1, [$ctx,#$Aoff+$hi] 387238384Sjkim adds $t2,$Tlo,$t2 388238384Sjkim str $t2, [$ctx,#$Boff+$lo] 389238384Sjkim adc $t3,$Thi,$t3 390238384Sjkim str $t3, [$ctx,#$Boff+$hi] 391238384Sjkim 392238384Sjkim ldr $Alo,[sp,#$Coff+0] 393238384Sjkim ldr $Ahi,[sp,#$Coff+4] 394238384Sjkim ldr $Tlo,[sp,#$Doff+0] 395238384Sjkim ldr $Thi,[sp,#$Doff+4] 396238384Sjkim ldr $t0, [$ctx,#$Coff+$lo] 397238384Sjkim ldr $t1, [$ctx,#$Coff+$hi] 398238384Sjkim ldr $t2, [$ctx,#$Doff+$lo] 399238384Sjkim ldr $t3, [$ctx,#$Doff+$hi] 400238384Sjkim adds $t0,$Alo,$t0 401238384Sjkim str $t0, [$ctx,#$Coff+$lo] 402238384Sjkim adc $t1,$Ahi,$t1 403238384Sjkim str $t1, [$ctx,#$Coff+$hi] 404238384Sjkim adds $t2,$Tlo,$t2 405238384Sjkim str $t2, [$ctx,#$Doff+$lo] 406238384Sjkim adc $t3,$Thi,$t3 407238384Sjkim str $t3, [$ctx,#$Doff+$hi] 408238384Sjkim 409238384Sjkim ldr $Tlo,[sp,#$Foff+0] 410238384Sjkim ldr $Thi,[sp,#$Foff+4] 411238384Sjkim ldr $t0, [$ctx,#$Eoff+$lo] 412238384Sjkim ldr $t1, [$ctx,#$Eoff+$hi] 413238384Sjkim ldr $t2, [$ctx,#$Foff+$lo] 414238384Sjkim ldr $t3, [$ctx,#$Foff+$hi] 415238384Sjkim adds $Elo,$Elo,$t0 416238384Sjkim str $Elo,[$ctx,#$Eoff+$lo] 417238384Sjkim adc $Ehi,$Ehi,$t1 418238384Sjkim str $Ehi,[$ctx,#$Eoff+$hi] 419238384Sjkim adds $t2,$Tlo,$t2 420238384Sjkim str $t2, [$ctx,#$Foff+$lo] 421238384Sjkim adc $t3,$Thi,$t3 422238384Sjkim str $t3, [$ctx,#$Foff+$hi] 423238384Sjkim 424238384Sjkim ldr $Alo,[sp,#$Goff+0] 425238384Sjkim ldr $Ahi,[sp,#$Goff+4] 426238384Sjkim ldr $Tlo,[sp,#$Hoff+0] 427238384Sjkim ldr $Thi,[sp,#$Hoff+4] 428238384Sjkim ldr $t0, [$ctx,#$Goff+$lo] 429238384Sjkim ldr $t1, [$ctx,#$Goff+$hi] 430238384Sjkim ldr $t2, [$ctx,#$Hoff+$lo] 431238384Sjkim ldr $t3, [$ctx,#$Hoff+$hi] 432238384Sjkim adds $t0,$Alo,$t0 433238384Sjkim str $t0, [$ctx,#$Goff+$lo] 434238384Sjkim adc $t1,$Ahi,$t1 435238384Sjkim str $t1, [$ctx,#$Goff+$hi] 436238384Sjkim adds $t2,$Tlo,$t2 437238384Sjkim str $t2, [$ctx,#$Hoff+$lo] 438238384Sjkim adc $t3,$Thi,$t3 439238384Sjkim str $t3, [$ctx,#$Hoff+$hi] 440238384Sjkim 441238384Sjkim add sp,sp,#640 442238384Sjkim sub $Ktbl,$Ktbl,#640 443238384Sjkim 444238384Sjkim teq $inp,$len 445238384Sjkim bne .Loop 446238384Sjkim 447238384Sjkim add sp,sp,#8*9 @ destroy frame 448238384Sjkim#if __ARM_ARCH__>=5 449238384Sjkim ldmia sp!,{r4-r12,pc} 450238384Sjkim#else 451238384Sjkim ldmia sp!,{r4-r12,lr} 452238384Sjkim tst lr,#1 453238384Sjkim moveq pc,lr @ be binary compatible with V4, yet 454238384Sjkim bx lr @ interoperable with Thumb ISA:-) 455238384Sjkim#endif 456238384Sjkim___ 457238384Sjkim 458238384Sjkim{ 459238384Sjkimmy @Sigma0=(28,34,39); 460238384Sjkimmy @Sigma1=(14,18,41); 461238384Sjkimmy @sigma0=(1, 8, 7); 462238384Sjkimmy @sigma1=(19,61,6); 463238384Sjkim 464238384Sjkimmy $Ktbl="r3"; 465238384Sjkimmy $cnt="r12"; # volatile register known as ip, intra-procedure-call scratch 466238384Sjkim 467238384Sjkimmy @X=map("d$_",(0..15)); 468238384Sjkimmy @V=($A,$B,$C,$D,$E,$F,$G,$H)=map("d$_",(16..23)); 469238384Sjkim 470238384Sjkimsub NEON_00_15() { 471238384Sjkimmy $i=shift; 472238384Sjkimmy ($a,$b,$c,$d,$e,$f,$g,$h)=@_; 473238384Sjkimmy ($t0,$t1,$t2,$T1,$K,$Ch,$Maj)=map("d$_",(24..31)); # temps 474238384Sjkim 475238384Sjkim$code.=<<___ if ($i<16 || $i&1); 476238384Sjkim vshr.u64 $t0,$e,#@Sigma1[0] @ $i 477238384Sjkim#if $i<16 478238384Sjkim vld1.64 {@X[$i%16]},[$inp]! @ handles unaligned 479238384Sjkim#endif 480238384Sjkim vshr.u64 $t1,$e,#@Sigma1[1] 481290207Sjkim#if $i>0 482290207Sjkim vadd.i64 $a,$Maj @ h+=Maj from the past 483290207Sjkim#endif 484238384Sjkim vshr.u64 $t2,$e,#@Sigma1[2] 485238384Sjkim___ 486238384Sjkim$code.=<<___; 487238384Sjkim vld1.64 {$K},[$Ktbl,:64]! @ K[i++] 488238384Sjkim vsli.64 $t0,$e,#`64-@Sigma1[0]` 489238384Sjkim vsli.64 $t1,$e,#`64-@Sigma1[1]` 490290207Sjkim vmov $Ch,$e 491238384Sjkim vsli.64 $t2,$e,#`64-@Sigma1[2]` 492238384Sjkim#if $i<16 && defined(__ARMEL__) 493238384Sjkim vrev64.8 @X[$i],@X[$i] 494238384Sjkim#endif 495290207Sjkim veor $t1,$t0 496290207Sjkim vbsl $Ch,$f,$g @ Ch(e,f,g) 497238384Sjkim vshr.u64 $t0,$a,#@Sigma0[0] 498290207Sjkim veor $t2,$t1 @ Sigma1(e) 499290207Sjkim vadd.i64 $T1,$Ch,$h 500238384Sjkim vshr.u64 $t1,$a,#@Sigma0[1] 501290207Sjkim vsli.64 $t0,$a,#`64-@Sigma0[0]` 502290207Sjkim vadd.i64 $T1,$t2 503238384Sjkim vshr.u64 $t2,$a,#@Sigma0[2] 504290207Sjkim vadd.i64 $K,@X[$i%16] 505238384Sjkim vsli.64 $t1,$a,#`64-@Sigma0[1]` 506290207Sjkim veor $Maj,$a,$b 507238384Sjkim vsli.64 $t2,$a,#`64-@Sigma0[2]` 508238384Sjkim veor $h,$t0,$t1 509290207Sjkim vadd.i64 $T1,$K 510290207Sjkim vbsl $Maj,$c,$b @ Maj(a,b,c) 511238384Sjkim veor $h,$t2 @ Sigma0(a) 512238384Sjkim vadd.i64 $d,$T1 513290207Sjkim vadd.i64 $Maj,$T1 514290207Sjkim @ vadd.i64 $h,$Maj 515238384Sjkim___ 516238384Sjkim} 517238384Sjkim 518238384Sjkimsub NEON_16_79() { 519238384Sjkimmy $i=shift; 520238384Sjkim 521238384Sjkimif ($i&1) { &NEON_00_15($i,@_); return; } 522238384Sjkim 523238384Sjkim# 2x-vectorized, therefore runs every 2nd round 524238384Sjkimmy @X=map("q$_",(0..7)); # view @X as 128-bit vector 525238384Sjkimmy ($t0,$t1,$s0,$s1) = map("q$_",(12..15)); # temps 526238384Sjkimmy ($d0,$d1,$d2) = map("d$_",(24..26)); # temps from NEON_00_15 527238384Sjkimmy $e=@_[4]; # $e from NEON_00_15 528238384Sjkim$i /= 2; 529238384Sjkim$code.=<<___; 530238384Sjkim vshr.u64 $t0,@X[($i+7)%8],#@sigma1[0] 531238384Sjkim vshr.u64 $t1,@X[($i+7)%8],#@sigma1[1] 532290207Sjkim vadd.i64 @_[0],d30 @ h+=Maj from the past 533238384Sjkim vshr.u64 $s1,@X[($i+7)%8],#@sigma1[2] 534238384Sjkim vsli.64 $t0,@X[($i+7)%8],#`64-@sigma1[0]` 535238384Sjkim vext.8 $s0,@X[$i%8],@X[($i+1)%8],#8 @ X[i+1] 536238384Sjkim vsli.64 $t1,@X[($i+7)%8],#`64-@sigma1[1]` 537238384Sjkim veor $s1,$t0 538238384Sjkim vshr.u64 $t0,$s0,#@sigma0[0] 539238384Sjkim veor $s1,$t1 @ sigma1(X[i+14]) 540238384Sjkim vshr.u64 $t1,$s0,#@sigma0[1] 541238384Sjkim vadd.i64 @X[$i%8],$s1 542238384Sjkim vshr.u64 $s1,$s0,#@sigma0[2] 543238384Sjkim vsli.64 $t0,$s0,#`64-@sigma0[0]` 544238384Sjkim vsli.64 $t1,$s0,#`64-@sigma0[1]` 545238384Sjkim vext.8 $s0,@X[($i+4)%8],@X[($i+5)%8],#8 @ X[i+9] 546238384Sjkim veor $s1,$t0 547238384Sjkim vshr.u64 $d0,$e,#@Sigma1[0] @ from NEON_00_15 548238384Sjkim vadd.i64 @X[$i%8],$s0 549238384Sjkim vshr.u64 $d1,$e,#@Sigma1[1] @ from NEON_00_15 550238384Sjkim veor $s1,$t1 @ sigma0(X[i+1]) 551238384Sjkim vshr.u64 $d2,$e,#@Sigma1[2] @ from NEON_00_15 552238384Sjkim vadd.i64 @X[$i%8],$s1 553238384Sjkim___ 554238384Sjkim &NEON_00_15(2*$i,@_); 555238384Sjkim} 556238384Sjkim 557238384Sjkim$code.=<<___; 558290207Sjkim#if __ARM_MAX_ARCH__>=7 559290207Sjkim.arch armv7-a 560238384Sjkim.fpu neon 561238384Sjkim 562238384Sjkim.align 4 563238384Sjkim.LNEON: 564238384Sjkim dmb @ errata #451034 on early Cortex A8 565238384Sjkim vstmdb sp!,{d8-d15} @ ABI specification says so 566238384Sjkim sub $Ktbl,r3,#672 @ K512 567238384Sjkim vldmia $ctx,{$A-$H} @ load context 568238384Sjkim.Loop_neon: 569238384Sjkim___ 570238384Sjkimfor($i=0;$i<16;$i++) { &NEON_00_15($i,@V); unshift(@V,pop(@V)); } 571238384Sjkim$code.=<<___; 572238384Sjkim mov $cnt,#4 573238384Sjkim.L16_79_neon: 574238384Sjkim subs $cnt,#1 575238384Sjkim___ 576238384Sjkimfor(;$i<32;$i++) { &NEON_16_79($i,@V); unshift(@V,pop(@V)); } 577238384Sjkim$code.=<<___; 578238384Sjkim bne .L16_79_neon 579238384Sjkim 580290207Sjkim vadd.i64 $A,d30 @ h+=Maj from the past 581238384Sjkim vldmia $ctx,{d24-d31} @ load context to temp 582238384Sjkim vadd.i64 q8,q12 @ vectorized accumulate 583238384Sjkim vadd.i64 q9,q13 584238384Sjkim vadd.i64 q10,q14 585238384Sjkim vadd.i64 q11,q15 586238384Sjkim vstmia $ctx,{$A-$H} @ save context 587238384Sjkim teq $inp,$len 588238384Sjkim sub $Ktbl,#640 @ rewind K512 589238384Sjkim bne .Loop_neon 590238384Sjkim 591238384Sjkim vldmia sp!,{d8-d15} @ epilogue 592290207Sjkim ret @ bx lr 593238384Sjkim#endif 594238384Sjkim___ 595238384Sjkim} 596238384Sjkim$code.=<<___; 597238384Sjkim.size sha512_block_data_order,.-sha512_block_data_order 598238384Sjkim.asciz "SHA512 block transform for ARMv4/NEON, CRYPTOGAMS by <appro\@openssl.org>" 599238384Sjkim.align 2 600290207Sjkim#if __ARM_MAX_ARCH__>=7 601238384Sjkim.comm OPENSSL_armcap_P,4,4 602290207Sjkim#endif 603238384Sjkim___ 604238384Sjkim 605238384Sjkim$code =~ s/\`([^\`]*)\`/eval $1/gem; 606238384Sjkim$code =~ s/\bbx\s+lr\b/.word\t0xe12fff1e/gm; # make it possible to compile with -march=armv4 607290207Sjkim$code =~ s/\bret\b/bx lr/gm; 608238384Sjkimprint $code; 609238384Sjkimclose STDOUT; # enforce flush 610