dont memcpy() simpler and same speed
[ffmpeg.git] / libavutil / sha1.c
1 // SHA-1 code Copyright 2007 Michael Nidermayer <michaelni@gmx.at>
2 // license LGPL
3 // based on public domain SHA-1 code by Steve Reid <steve@edmweb.com>
4
5 #include "common.h"
6 #include "sha1.h"
7
8 typedef struct AVSHA1 {
9     uint32_t state[5];
10     uint64_t count;
11     uint8_t buffer[64];
12 } AVSHA1;
13
14 #define rol(value, bits) (((value) << (bits)) | ((value) >> (32 - (bits))))
15
16 /* (R0+R1), R2, R3, R4 are the different operations used in SHA1 */
17 #define blk0(i) (block[i] = be2me_32(((uint32_t*)buffer)[i]))
18 #define blk(i) (block[i] = rol(block[i-3]^block[i-8]^block[i-14]^block[i-16],1))
19
20 #define R0(v,w,x,y,z,i) z+=((w&(x^y))^y)    +blk0(i)+0x5A827999+rol(v,5);w=rol(w,30);
21 #define R1(v,w,x,y,z,i) z+=((w&(x^y))^y)    +blk (i)+0x5A827999+rol(v,5);w=rol(w,30);
22 #define R2(v,w,x,y,z,i) z+=( w^x     ^y)    +blk (i)+0x6ED9EBA1+rol(v,5);w=rol(w,30);
23 #define R3(v,w,x,y,z,i) z+=(((w|x)&y)|(w&x))+blk (i)+0x8F1BBCDC+rol(v,5);w=rol(w,30);
24 #define R4(v,w,x,y,z,i) z+=( w^x     ^y)    +blk (i)+0xCA62C1D6+rol(v,5);w=rol(w,30);
25
26 /* Hash a single 512-bit block. This is the core of the algorithm. */
27
28 //#define VARIANT2
29
30 static void transform(uint32_t state[5], uint8_t buffer[64]){
31     uint32_t block[80];
32     unsigned int i;
33     unsigned int a, b, c, d, e;
34
35     a = state[0];
36     b = state[1];
37     c = state[2];
38     d = state[3];
39     e = state[4];
40 #ifdef VARIANT2
41     for(i=0; i<80; i++){
42         int t;
43         if(i<16) t= be2me_32(((uint32_t*)buffer)[i]);
44         else     t= rol(block[i-3]^block[i-8]^block[i-14]^block[i-16],1);
45         block[i]= t;
46         t+= e+rol(a,5);
47         if(i<40){
48             if(i<20)    t+= ((b&(c^d))^d)    +0x5A827999;
49             else        t+= ( b^c     ^d)    +0x6ED9EBA1;
50         }else{
51             if(i<60)    t+= (((b|c)&d)|(b&c))+0x8F1BBCDC;
52             else        t+= ( b^c     ^d)    +0xCA62C1D6;
53         }
54         e= d;
55         d= c;
56         c= rol(b,30);
57         b= a;
58         a= t;
59     }
60 #else
61     for(i=0; i<15; i+=5){
62         R0(a,b,c,d,e,0+i); R0(e,a,b,c,d,1+i); R0(d,e,a,b,c,2+i); R0(c,d,e,a,b,3+i); R0(b,c,d,e,a,4+i);
63     }
64     R0(a,b,c,d,e,15); R1(e,a,b,c,d,16); R1(d,e,a,b,c,17); R1(c,d,e,a,b,18); R1(b,c,d,e,a,19);
65     for(i=20; i<40; i+=5){
66         R2(a,b,c,d,e,0+i); R2(e,a,b,c,d,1+i); R2(d,e,a,b,c,2+i); R2(c,d,e,a,b,3+i); R2(b,c,d,e,a,4+i);
67     }
68     for(; i<60; i+=5){
69         R3(a,b,c,d,e,0+i); R3(e,a,b,c,d,1+i); R3(d,e,a,b,c,2+i); R3(c,d,e,a,b,3+i); R3(b,c,d,e,a,4+i);
70     }
71     for(; i<80; i+=5){
72         R4(a,b,c,d,e,0+i); R4(e,a,b,c,d,1+i); R4(d,e,a,b,c,2+i); R4(c,d,e,a,b,3+i); R4(b,c,d,e,a,4+i);
73     }
74 #endif
75     state[0] += a;
76     state[1] += b;
77     state[2] += c;
78     state[3] += d;
79     state[4] += e;
80 }
81
82 void av_sha1_init(AVSHA1* context){
83     context->state[0] = 0x67452301;
84     context->state[1] = 0xEFCDAB89;
85     context->state[2] = 0x98BADCFE;
86     context->state[3] = 0x10325476;
87     context->state[4] = 0xC3D2E1F0;
88     context->count    = 0;
89 }
90
91 void av_sha1_update(AVSHA1* context, uint8_t* data, unsigned int len){
92     unsigned int i, j;
93
94     j = context->count & 63;
95     context->count += len;
96     if ((j + len) > 63) {
97         memcpy(&context->buffer[j], data, (i = 64-j));
98         transform(context->state, context->buffer);
99         for ( ; i + 63 < len; i += 64) {
100             transform(context->state, &data[i]);
101         }
102         j=0;
103     }
104     else i = 0;
105     memcpy(&context->buffer[j], &data[i], len - i);
106 }
107
108 void av_sha1_final(AVSHA1* context, uint8_t digest[20]){
109     int i;
110     uint64_t finalcount= be2me_64(context->count<<3);
111
112     av_sha1_update(context, "\200", 1);
113     while ((context->count & 63) != 56) {
114         av_sha1_update(context, "", 1);
115     }
116     av_sha1_update(context, &finalcount, 8);  /* Should cause a transform() */
117     for(i=0; i<5; i++)
118         ((uint32_t*)digest)[i]= be2me_32(context->state[i]);
119 }
120
121 // use the following to test
122 // gcc -DTEST -DHAVE_AV_CONFIG_H -I.. sha1.c -O2 -W -Wall -o sha1 && time ./sha1
123 #ifdef TEST
124 #include <stdio.h>
125 #undef printf
126
127 int main(){
128     int i, k;
129     AVSHA1 context;
130     unsigned char digest[20];
131
132     for(k=0; k<3; k++){
133         av_sha1_init(&context);
134         if(k==0)
135             av_sha1_update(&context, "abc", 3);
136         else if(k==1)
137             av_sha1_update(&context, "abcdbcdecdefdefgefghfghighijhijkijkljklmklmnlmnomnopnopq", 56);
138         else
139             for(i=0; i<1000*1000; i++)
140                 av_sha1_update(&context, "a", 1);
141         av_sha1_final(&context, digest);
142         for (i = 0; i < 20; i++)
143             printf("%02X", digest[i]);
144         putchar('\n');
145     }
146     //Test Vectors (from FIPS PUB 180-1)
147     printf("A9993E36 4706816A BA3E2571 7850C26C 9CD0D89D\n"
148            "84983E44 1C3BD26E BAAE4AA1 F95129E5 E54670F1\n"
149            "34AA973C D4C4DAA4 F61EEB2B DBAD2731 6534016F\n");
150
151     return 0;
152 }
153 #endif