1/* 2 * SHA1 routine optimized to do word accesses rather than byte accesses, 3 * and to avoid unnecessary copies into the context array. 4 * 5 * This was initially based on the Mozilla SHA1 implementation, although 6 * none of the original Mozilla code remains. 7 */ 8 9/* this is only to get definitions for memcpy(), ntohl() and htonl() */ 10#include"../git-compat-util.h" 11 12#include"sha1.h" 13 14#if defined(__GNUC__) && (defined(__i386__) || defined(__x86_64__)) 15 16/* 17 * Force usage of rol or ror by selecting the one with the smaller constant. 18 * It _can_ generate slightly smaller code (a constant of 1 is special), but 19 * perhaps more importantly it's possibly faster on any uarch that does a 20 * rotate with a loop. 21 */ 22 23#define SHA_ASM(op, x, n) ({ unsigned int __res; __asm__(op"%1,%0":"=r" (__res):"i" (n),"0" (x)); __res; }) 24#define SHA_ROL(x,n) SHA_ASM("rol", x, n) 25#define SHA_ROR(x,n) SHA_ASM("ror", x, n) 26 27#else 28 29#define SHA_ROT(X,l,r) (((X) << (l)) | ((X) >> (r))) 30#define SHA_ROL(X,n) SHA_ROT(X,n,32-(n)) 31#define SHA_ROR(X,n) SHA_ROT(X,32-(n),n) 32 33#endif 34 35/* 36 * If you have 32 registers or more, the compiler can (and should) 37 * try to change the array[] accesses into registers. However, on 38 * machines with less than ~25 registers, that won't really work, 39 * and at least gcc will make an unholy mess of it. 40 * 41 * So to avoid that mess which just slows things down, we force 42 * the stores to memory to actually happen (we might be better off 43 * with a 'W(t)=(val);asm("":"+m" (W(t))' there instead, as 44 * suggested by Artur Skawina - that will also make gcc unable to 45 * try to do the silly "optimize away loads" part because it won't 46 * see what the value will be). 47 * 48 * Ben Herrenschmidt reports that on PPC, the C version comes close 49 * to the optimized asm with this (ie on PPC you don't want that 50 * 'volatile', since there are lots of registers). 51 * 52 * On ARM we get the best code generation by forcing a full memory barrier 53 * between each SHA_ROUND, otherwise gcc happily get wild with spilling and 54 * the stack frame size simply explode and performance goes down the drain. 55 */ 56 57#if defined(__i386__) || defined(__x86_64__) 58#define setW(x, val) (*(volatile unsigned int *)&W(x) = (val)) 59#elif defined(__GNUC__) && defined(__arm__) 60#define setW(x, val) do { W(x) = (val); __asm__("":::"memory"); } while (0) 61#else 62#define setW(x, val) (W(x) = (val)) 63#endif 64 65/* 66 * Performance might be improved if the CPU architecture is OK with 67 * unaligned 32-bit loads and a fast ntohl() is available. 68 * Otherwise fall back to byte loads and shifts which is portable, 69 * and is faster on architectures with memory alignment issues. 70 */ 71 72#if defined(__i386__) || defined(__x86_64__) || \ 73 defined(_M_IX86) || defined(_M_X64) || \ 74 defined(__ppc__) || defined(__ppc64__) || \ 75 defined(__powerpc__) || defined(__powerpc64__) || \ 76 defined(__s390__) || defined(__s390x__) 77 78#define get_be32(p) ntohl(*(unsigned int *)(p)) 79#define put_be32(p, v) do { *(unsigned int *)(p) = htonl(v); } while (0) 80 81#else 82 83#define get_be32(p) ( \ 84 (*((unsigned char *)(p) + 0) << 24) | \ 85 (*((unsigned char *)(p) + 1) << 16) | \ 86 (*((unsigned char *)(p) + 2) << 8) | \ 87 (*((unsigned char *)(p) + 3) << 0) ) 88#define put_be32(p, v) do { \ 89 unsigned int __v = (v); \ 90 *((unsigned char *)(p) + 0) = __v >> 24; \ 91 *((unsigned char *)(p) + 1) = __v >> 16; \ 92 *((unsigned char *)(p) + 2) = __v >> 8; \ 93 *((unsigned char *)(p) + 3) = __v >> 0; } while (0) 94 95#endif 96 97/* This "rolls" over the 512-bit array */ 98#define W(x) (array[(x)&15]) 99 100/* 101 * Where do we get the source from? The first 16 iterations get it from 102 * the input data, the next mix it from the 512-bit array. 103 */ 104#define SHA_SRC(t) get_be32(data + t) 105#define SHA_MIX(t) SHA_ROL(W(t+13) ^ W(t+8) ^ W(t+2) ^ W(t), 1) 106 107#define SHA_ROUND(t, input, fn, constant, A, B, C, D, E) do { \ 108 unsigned int TEMP = input(t); setW(t, TEMP); \ 109 E += TEMP + SHA_ROL(A,5) + (fn) + (constant); \ 110 B = SHA_ROR(B, 2); } while (0) 111 112#define T_0_15(t, A, B, C, D, E) SHA_ROUND(t, SHA_SRC, (((C^D)&B)^D) , 0x5a827999, A, B, C, D, E ) 113#define T_16_19(t, A, B, C, D, E) SHA_ROUND(t, SHA_MIX, (((C^D)&B)^D) , 0x5a827999, A, B, C, D, E ) 114#define T_20_39(t, A, B, C, D, E) SHA_ROUND(t, SHA_MIX, (B^C^D) , 0x6ed9eba1, A, B, C, D, E ) 115#define T_40_59(t, A, B, C, D, E) SHA_ROUND(t, SHA_MIX, ((B&C)+(D&(B^C))) , 0x8f1bbcdc, A, B, C, D, E ) 116#define T_60_79(t, A, B, C, D, E) SHA_ROUND(t, SHA_MIX, (B^C^D) , 0xca62c1d6, A, B, C, D, E ) 117 118static voidblk_SHA1_Block(blk_SHA_CTX *ctx,const unsigned int*data) 119{ 120unsigned int A,B,C,D,E; 121unsigned int array[16]; 122 123 A = ctx->H[0]; 124 B = ctx->H[1]; 125 C = ctx->H[2]; 126 D = ctx->H[3]; 127 E = ctx->H[4]; 128 129/* Round 1 - iterations 0-16 take their input from 'data' */ 130T_0_15(0, A, B, C, D, E); 131T_0_15(1, E, A, B, C, D); 132T_0_15(2, D, E, A, B, C); 133T_0_15(3, C, D, E, A, B); 134T_0_15(4, B, C, D, E, A); 135T_0_15(5, A, B, C, D, E); 136T_0_15(6, E, A, B, C, D); 137T_0_15(7, D, E, A, B, C); 138T_0_15(8, C, D, E, A, B); 139T_0_15(9, B, C, D, E, A); 140T_0_15(10, A, B, C, D, E); 141T_0_15(11, E, A, B, C, D); 142T_0_15(12, D, E, A, B, C); 143T_0_15(13, C, D, E, A, B); 144T_0_15(14, B, C, D, E, A); 145T_0_15(15, A, B, C, D, E); 146 147/* Round 1 - tail. Input from 512-bit mixing array */ 148T_16_19(16, E, A, B, C, D); 149T_16_19(17, D, E, A, B, C); 150T_16_19(18, C, D, E, A, B); 151T_16_19(19, B, C, D, E, A); 152 153/* Round 2 */ 154T_20_39(20, A, B, C, D, E); 155T_20_39(21, E, A, B, C, D); 156T_20_39(22, D, E, A, B, C); 157T_20_39(23, C, D, E, A, B); 158T_20_39(24, B, C, D, E, A); 159T_20_39(25, A, B, C, D, E); 160T_20_39(26, E, A, B, C, D); 161T_20_39(27, D, E, A, B, C); 162T_20_39(28, C, D, E, A, B); 163T_20_39(29, B, C, D, E, A); 164T_20_39(30, A, B, C, D, E); 165T_20_39(31, E, A, B, C, D); 166T_20_39(32, D, E, A, B, C); 167T_20_39(33, C, D, E, A, B); 168T_20_39(34, B, C, D, E, A); 169T_20_39(35, A, B, C, D, E); 170T_20_39(36, E, A, B, C, D); 171T_20_39(37, D, E, A, B, C); 172T_20_39(38, C, D, E, A, B); 173T_20_39(39, B, C, D, E, A); 174 175/* Round 3 */ 176T_40_59(40, A, B, C, D, E); 177T_40_59(41, E, A, B, C, D); 178T_40_59(42, D, E, A, B, C); 179T_40_59(43, C, D, E, A, B); 180T_40_59(44, B, C, D, E, A); 181T_40_59(45, A, B, C, D, E); 182T_40_59(46, E, A, B, C, D); 183T_40_59(47, D, E, A, B, C); 184T_40_59(48, C, D, E, A, B); 185T_40_59(49, B, C, D, E, A); 186T_40_59(50, A, B, C, D, E); 187T_40_59(51, E, A, B, C, D); 188T_40_59(52, D, E, A, B, C); 189T_40_59(53, C, D, E, A, B); 190T_40_59(54, B, C, D, E, A); 191T_40_59(55, A, B, C, D, E); 192T_40_59(56, E, A, B, C, D); 193T_40_59(57, D, E, A, B, C); 194T_40_59(58, C, D, E, A, B); 195T_40_59(59, B, C, D, E, A); 196 197/* Round 4 */ 198T_60_79(60, A, B, C, D, E); 199T_60_79(61, E, A, B, C, D); 200T_60_79(62, D, E, A, B, C); 201T_60_79(63, C, D, E, A, B); 202T_60_79(64, B, C, D, E, A); 203T_60_79(65, A, B, C, D, E); 204T_60_79(66, E, A, B, C, D); 205T_60_79(67, D, E, A, B, C); 206T_60_79(68, C, D, E, A, B); 207T_60_79(69, B, C, D, E, A); 208T_60_79(70, A, B, C, D, E); 209T_60_79(71, E, A, B, C, D); 210T_60_79(72, D, E, A, B, C); 211T_60_79(73, C, D, E, A, B); 212T_60_79(74, B, C, D, E, A); 213T_60_79(75, A, B, C, D, E); 214T_60_79(76, E, A, B, C, D); 215T_60_79(77, D, E, A, B, C); 216T_60_79(78, C, D, E, A, B); 217T_60_79(79, B, C, D, E, A); 218 219 ctx->H[0] += A; 220 ctx->H[1] += B; 221 ctx->H[2] += C; 222 ctx->H[3] += D; 223 ctx->H[4] += E; 224} 225 226voidblk_SHA1_Init(blk_SHA_CTX *ctx) 227{ 228 ctx->size =0; 229 230/* Initialize H with the magic constants (see FIPS180 for constants) */ 231 ctx->H[0] =0x67452301; 232 ctx->H[1] =0xefcdab89; 233 ctx->H[2] =0x98badcfe; 234 ctx->H[3] =0x10325476; 235 ctx->H[4] =0xc3d2e1f0; 236} 237 238voidblk_SHA1_Update(blk_SHA_CTX *ctx,const void*data,unsigned long len) 239{ 240unsigned int lenW = ctx->size &63; 241 242 ctx->size += len; 243 244/* Read the data into W and process blocks as they get full */ 245if(lenW) { 246unsigned int left =64- lenW; 247if(len < left) 248 left = len; 249memcpy(lenW + (char*)ctx->W, data, left); 250 lenW = (lenW + left) &63; 251 len -= left; 252 data = ((const char*)data + left); 253if(lenW) 254return; 255blk_SHA1_Block(ctx, ctx->W); 256} 257while(len >=64) { 258blk_SHA1_Block(ctx, data); 259 data = ((const char*)data +64); 260 len -=64; 261} 262if(len) 263memcpy(ctx->W, data, len); 264} 265 266voidblk_SHA1_Final(unsigned char hashout[20], blk_SHA_CTX *ctx) 267{ 268static const unsigned char pad[64] = {0x80}; 269unsigned int padlen[2]; 270int i; 271 272/* Pad with a binary 1 (ie 0x80), then zeroes, then length */ 273 padlen[0] =htonl((uint32_t)(ctx->size >>29)); 274 padlen[1] =htonl((uint32_t)(ctx->size <<3)); 275 276 i = ctx->size &63; 277blk_SHA1_Update(ctx, pad,1+ (63& (55- i))); 278blk_SHA1_Update(ctx, padlen,8); 279 280/* Output hash */ 281for(i =0; i <5; i++) 282put_be32(hashout + i*4, ctx->H[i]); 283}