sha512: reduce stack use in transform function by 512 bytes
[libgcrypt.git] / cipher / sha512.c
1 /* sha512.c - SHA384 and SHA512 hash functions
2  * Copyright (C) 2003, 2008, 2009 Free Software Foundation, Inc.
3  *
4  * This file is part of Libgcrypt.
5  *
6  * Libgcrypt is free software; you can redistribute it and/or modify
7  * it under the terms of the GNU Lesser general Public License as
8  * published by the Free Software Foundation; either version 2.1 of
9  * the License, or (at your option) any later version.
10  *
11  * Libgcrypt is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14  * GNU Lesser General Public License for more details.
15  *
16  * You should have received a copy of the GNU Lesser General Public
17  * License along with this program; if not, see <http://www.gnu.org/licenses/>.
18  */
19
20
21 /*  Test vectors from FIPS-180-2:
22  *
23  *  "abc"
24  * 384:
25  *  CB00753F 45A35E8B B5A03D69 9AC65007 272C32AB 0EDED163
26  *  1A8B605A 43FF5BED 8086072B A1E7CC23 58BAECA1 34C825A7
27  * 512:
28  *  DDAF35A1 93617ABA CC417349 AE204131 12E6FA4E 89A97EA2 0A9EEEE6 4B55D39A
29  *  2192992A 274FC1A8 36BA3C23 A3FEEBBD 454D4423 643CE80E 2A9AC94F A54CA49F
30  *
31  *  "abcdefghbcdefghicdefghijdefghijkefghijklfghijklmghijklmnhijklmnoijklmnopjklmnopqklmnopqrlmnopqrsmnopqrstnopqrstu"
32  * 384:
33  *  09330C33 F71147E8 3D192FC7 82CD1B47 53111B17 3B3B05D2
34  *  2FA08086 E3B0F712 FCC7C71A 557E2DB9 66C3E9FA 91746039
35  * 512:
36  *  8E959B75 DAE313DA 8CF4F728 14FC143F 8F7779C6 EB9F7FA1 7299AEAD B6889018
37  *  501D289E 4900F7E4 331B99DE C4B5433A C7D329EE B6DD2654 5E96E55B 874BE909
38  *
39  *  "a" x 1000000
40  * 384:
41  *  9D0E1809 716474CB 086E834E 310A4A1C ED149E9C 00F24852
42  *  7972CEC5 704C2A5B 07B8B3DC 38ECC4EB AE97DDD8 7F3D8985
43  * 512:
44  *  E718483D 0CE76964 4E2E42C7 BC15B463 8E1F98B1 3B204428 5632A803 AFA973EB
45  *  DE0FF244 877EA60A 4CB0432C E577C31B EB009C5C 2C49AA2E 4EADB217 AD8CC09B
46  */
47
48
49 #include <config.h>
50 #include <string.h>
51 #include "g10lib.h"
52 #include "bithelp.h"
53 #include "cipher.h"
54 #include "hash-common.h"
55
56 typedef struct
57 {
58   u64 h0, h1, h2, h3, h4, h5, h6, h7;
59   u64 nblocks;
60   byte buf[128];
61   int count;
62 } SHA512_CONTEXT;
63
64 static void
65 sha512_init (void *context)
66 {
67   SHA512_CONTEXT *hd = context;
68
69   hd->h0 = U64_C(0x6a09e667f3bcc908);
70   hd->h1 = U64_C(0xbb67ae8584caa73b);
71   hd->h2 = U64_C(0x3c6ef372fe94f82b);
72   hd->h3 = U64_C(0xa54ff53a5f1d36f1);
73   hd->h4 = U64_C(0x510e527fade682d1);
74   hd->h5 = U64_C(0x9b05688c2b3e6c1f);
75   hd->h6 = U64_C(0x1f83d9abfb41bd6b);
76   hd->h7 = U64_C(0x5be0cd19137e2179);
77
78   hd->nblocks = 0;
79   hd->count = 0;
80 }
81
82 static void
83 sha384_init (void *context)
84 {
85   SHA512_CONTEXT *hd = context;
86
87   hd->h0 = U64_C(0xcbbb9d5dc1059ed8);
88   hd->h1 = U64_C(0x629a292a367cd507);
89   hd->h2 = U64_C(0x9159015a3070dd17);
90   hd->h3 = U64_C(0x152fecd8f70e5939);
91   hd->h4 = U64_C(0x67332667ffc00b31);
92   hd->h5 = U64_C(0x8eb44a8768581511);
93   hd->h6 = U64_C(0xdb0c2e0d64f98fa7);
94   hd->h7 = U64_C(0x47b5481dbefa4fa4);
95
96   hd->nblocks = 0;
97   hd->count = 0;
98 }
99
100
101 static inline u64
102 ROTR (u64 x, u64 n)
103 {
104   return ((x >> n) | (x << (64 - n)));
105 }
106
107 static inline u64
108 Ch (u64 x, u64 y, u64 z)
109 {
110   return ((x & y) ^ ( ~x & z));
111 }
112
113 static inline u64
114 Maj (u64 x, u64 y, u64 z)
115 {
116   return ((x & y) ^ (x & z) ^ (y & z));
117 }
118
119 static inline u64
120 Sum0 (u64 x)
121 {
122   return (ROTR (x, 28) ^ ROTR (x, 34) ^ ROTR (x, 39));
123 }
124
125 static inline u64
126 Sum1 (u64 x)
127 {
128   return (ROTR (x, 14) ^ ROTR (x, 18) ^ ROTR (x, 41));
129 }
130
131 /****************
132  * Transform the message W which consists of 16 64-bit-words
133  */
134 static void
135 transform (SHA512_CONTEXT *hd, const unsigned char *data)
136 {
137   u64 a, b, c, d, e, f, g, h;
138   u64 w[16];
139   int t;
140   static const u64 k[] =
141     {
142       U64_C(0x428a2f98d728ae22), U64_C(0x7137449123ef65cd),
143       U64_C(0xb5c0fbcfec4d3b2f), U64_C(0xe9b5dba58189dbbc),
144       U64_C(0x3956c25bf348b538), U64_C(0x59f111f1b605d019),
145       U64_C(0x923f82a4af194f9b), U64_C(0xab1c5ed5da6d8118),
146       U64_C(0xd807aa98a3030242), U64_C(0x12835b0145706fbe),
147       U64_C(0x243185be4ee4b28c), U64_C(0x550c7dc3d5ffb4e2),
148       U64_C(0x72be5d74f27b896f), U64_C(0x80deb1fe3b1696b1),
149       U64_C(0x9bdc06a725c71235), U64_C(0xc19bf174cf692694),
150       U64_C(0xe49b69c19ef14ad2), U64_C(0xefbe4786384f25e3),
151       U64_C(0x0fc19dc68b8cd5b5), U64_C(0x240ca1cc77ac9c65),
152       U64_C(0x2de92c6f592b0275), U64_C(0x4a7484aa6ea6e483),
153       U64_C(0x5cb0a9dcbd41fbd4), U64_C(0x76f988da831153b5),
154       U64_C(0x983e5152ee66dfab), U64_C(0xa831c66d2db43210),
155       U64_C(0xb00327c898fb213f), U64_C(0xbf597fc7beef0ee4),
156       U64_C(0xc6e00bf33da88fc2), U64_C(0xd5a79147930aa725),
157       U64_C(0x06ca6351e003826f), U64_C(0x142929670a0e6e70),
158       U64_C(0x27b70a8546d22ffc), U64_C(0x2e1b21385c26c926),
159       U64_C(0x4d2c6dfc5ac42aed), U64_C(0x53380d139d95b3df),
160       U64_C(0x650a73548baf63de), U64_C(0x766a0abb3c77b2a8),
161       U64_C(0x81c2c92e47edaee6), U64_C(0x92722c851482353b),
162       U64_C(0xa2bfe8a14cf10364), U64_C(0xa81a664bbc423001),
163       U64_C(0xc24b8b70d0f89791), U64_C(0xc76c51a30654be30),
164       U64_C(0xd192e819d6ef5218), U64_C(0xd69906245565a910),
165       U64_C(0xf40e35855771202a), U64_C(0x106aa07032bbd1b8),
166       U64_C(0x19a4c116b8d2d0c8), U64_C(0x1e376c085141ab53),
167       U64_C(0x2748774cdf8eeb99), U64_C(0x34b0bcb5e19b48a8),
168       U64_C(0x391c0cb3c5c95a63), U64_C(0x4ed8aa4ae3418acb),
169       U64_C(0x5b9cca4f7763e373), U64_C(0x682e6ff3d6b2b8a3),
170       U64_C(0x748f82ee5defb2fc), U64_C(0x78a5636f43172f60),
171       U64_C(0x84c87814a1f0ab72), U64_C(0x8cc702081a6439ec),
172       U64_C(0x90befffa23631e28), U64_C(0xa4506cebde82bde9),
173       U64_C(0xbef9a3f7b2c67915), U64_C(0xc67178f2e372532b),
174       U64_C(0xca273eceea26619c), U64_C(0xd186b8c721c0c207),
175       U64_C(0xeada7dd6cde0eb1e), U64_C(0xf57d4f7fee6ed178),
176       U64_C(0x06f067aa72176fba), U64_C(0x0a637dc5a2c898a6),
177       U64_C(0x113f9804bef90dae), U64_C(0x1b710b35131c471b),
178       U64_C(0x28db77f523047d84), U64_C(0x32caab7b40c72493),
179       U64_C(0x3c9ebe0a15c9bebc), U64_C(0x431d67c49c100d4c),
180       U64_C(0x4cc5d4becb3e42b6), U64_C(0x597f299cfc657e2a),
181       U64_C(0x5fcb6fab3ad6faec), U64_C(0x6c44198c4a475817)
182     };
183
184   /* get values from the chaining vars */
185   a = hd->h0;
186   b = hd->h1;
187   c = hd->h2;
188   d = hd->h3;
189   e = hd->h4;
190   f = hd->h5;
191   g = hd->h6;
192   h = hd->h7;
193
194 #ifdef WORDS_BIGENDIAN
195   memcpy (w, data, 128);
196 #else
197   {
198     int i;
199     byte *p2;
200
201     for (i = 0, p2 = (byte *) w; i < 16; i++, p2 += 8)
202       {
203         p2[7] = *data++;
204         p2[6] = *data++;
205         p2[5] = *data++;
206         p2[4] = *data++;
207         p2[3] = *data++;
208         p2[2] = *data++;
209         p2[1] = *data++;
210         p2[0] = *data++;
211       }
212   }
213 #endif
214
215 #define S0(x) (ROTR((x),1) ^ ROTR((x),8) ^ ((x)>>7))
216 #define S1(x) (ROTR((x),19) ^ ROTR((x),61) ^ ((x)>>6))
217
218
219   for (t = 0; t < 80 - 16; )
220     {
221       u64 t1, t2;
222
223       /* Performance on a AMD Athlon(tm) Dual Core Processor 4050e
224          with gcc 4.3.3 using gcry_md_hash_buffer of each 10000 bytes
225          initialized to 0,1,2,3...255,0,... and 1000 iterations:
226
227          Not unrolled with macros:  440ms
228          Unrolled with macros:      350ms
229          Unrolled with inline:      330ms
230       */
231 #if 0 /* Not unrolled.  */
232       t1 = h + Sum1 (e) + Ch (e, f, g) + k[t] + w[t%16];
233       w[t%16] += S1 (w[(t - 2)%16]) + w[(t - 7)%16] + S0 (w[(t - 15)%16]);
234       t2 = Sum0 (a) + Maj (a, b, c);
235       h = g;
236       g = f;
237       f = e;
238       e = d + t1;
239       d = c;
240       c = b;
241       b = a;
242       a = t1 + t2;
243       t++;
244 #else /* Unrolled to interweave the chain variables.  */
245       t1 = h + Sum1 (e) + Ch (e, f, g) + k[t] + w[0];
246       w[0] += S1 (w[14]) + w[9] + S0 (w[1]);
247       t2 = Sum0 (a) + Maj (a, b, c);
248       d += t1;
249       h = t1 + t2;
250
251       t1 = g + Sum1 (d) + Ch (d, e, f) + k[t+1] + w[1];
252       w[1] += S1 (w[15]) + w[10] + S0 (w[2]);
253       t2 = Sum0 (h) + Maj (h, a, b);
254       c += t1;
255       g  = t1 + t2;
256
257       t1 = f + Sum1 (c) + Ch (c, d, e) + k[t+2] + w[2];
258       w[2] += S1 (w[0]) + w[11] + S0 (w[3]);
259       t2 = Sum0 (g) + Maj (g, h, a);
260       b += t1;
261       f  = t1 + t2;
262
263       t1 = e + Sum1 (b) + Ch (b, c, d) + k[t+3] + w[3];
264       w[3] += S1 (w[1]) + w[12] + S0 (w[4]);
265       t2 = Sum0 (f) + Maj (f, g, h);
266       a += t1;
267       e  = t1 + t2;
268
269       t1 = d + Sum1 (a) + Ch (a, b, c) + k[t+4] + w[4];
270       w[4] += S1 (w[2]) + w[13] + S0 (w[5]);
271       t2 = Sum0 (e) + Maj (e, f, g);
272       h += t1;
273       d  = t1 + t2;
274
275       t1 = c + Sum1 (h) + Ch (h, a, b) + k[t+5] + w[5];
276       w[5] += S1 (w[3]) + w[14] + S0 (w[6]);
277       t2 = Sum0 (d) + Maj (d, e, f);
278       g += t1;
279       c  = t1 + t2;
280
281       t1 = b + Sum1 (g) + Ch (g, h, a) + k[t+6] + w[6];
282       w[6] += S1 (w[4]) + w[15] + S0 (w[7]);
283       t2 = Sum0 (c) + Maj (c, d, e);
284       f += t1;
285       b  = t1 + t2;
286
287       t1 = a + Sum1 (f) + Ch (f, g, h) + k[t+7] + w[7];
288       w[7] += S1 (w[5]) + w[0] + S0 (w[8]);
289       t2 = Sum0 (b) + Maj (b, c, d);
290       e += t1;
291       a  = t1 + t2;
292
293       t1 = h + Sum1 (e) + Ch (e, f, g) + k[t+8] + w[8];
294       w[8] += S1 (w[6]) + w[1] + S0 (w[9]);
295       t2 = Sum0 (a) + Maj (a, b, c);
296       d += t1;
297       h  = t1 + t2;
298
299       t1 = g + Sum1 (d) + Ch (d, e, f) + k[t+9] + w[9];
300       w[9] += S1 (w[7]) + w[2] + S0 (w[10]);
301       t2 = Sum0 (h) + Maj (h, a, b);
302       c += t1;
303       g  = t1 + t2;
304
305       t1 = f + Sum1 (c) + Ch (c, d, e) + k[t+10] + w[10];
306       w[10] += S1 (w[8]) + w[3] + S0 (w[11]);
307       t2 = Sum0 (g) + Maj (g, h, a);
308       b += t1;
309       f  = t1 + t2;
310
311       t1 = e + Sum1 (b) + Ch (b, c, d) + k[t+11] + w[11];
312       w[11] += S1 (w[9]) + w[4] + S0 (w[12]);
313       t2 = Sum0 (f) + Maj (f, g, h);
314       a += t1;
315       e  = t1 + t2;
316
317       t1 = d + Sum1 (a) + Ch (a, b, c) + k[t+12] + w[12];
318       w[12] += S1 (w[10]) + w[5] + S0 (w[13]);
319       t2 = Sum0 (e) + Maj (e, f, g);
320       h += t1;
321       d  = t1 + t2;
322
323       t1 = c + Sum1 (h) + Ch (h, a, b) + k[t+13] + w[13];
324       w[13] += S1 (w[11]) + w[6] + S0 (w[14]);
325       t2 = Sum0 (d) + Maj (d, e, f);
326       g += t1;
327       c  = t1 + t2;
328
329       t1 = b + Sum1 (g) + Ch (g, h, a) + k[t+14] + w[14];
330       w[14] += S1 (w[12]) + w[7] + S0 (w[15]);
331       t2 = Sum0 (c) + Maj (c, d, e);
332       f += t1;
333       b  = t1 + t2;
334
335       t1 = a + Sum1 (f) + Ch (f, g, h) + k[t+15] + w[15];
336       w[15] += S1 (w[13]) + w[8] + S0 (w[0]);
337       t2 = Sum0 (b) + Maj (b, c, d);
338       e += t1;
339       a  = t1 + t2;
340
341       t += 16;
342 #endif
343     }
344
345   for (; t < 80; )
346     {
347       u64 t1, t2;
348
349 #if 0 /* Not unrolled.  */
350       t1 = h + Sum1 (e) + Ch (e, f, g) + k[t] + w[t%16];
351       t2 = Sum0 (a) + Maj (a, b, c);
352       h = g;
353       g = f;
354       f = e;
355       e = d + t1;
356       d = c;
357       c = b;
358       b = a;
359       a = t1 + t2;
360       t++;
361 #else /* Unrolled to interweave the chain variables.  */
362       t1 = h + Sum1 (e) + Ch (e, f, g) + k[t] + w[0];
363       t2 = Sum0 (a) + Maj (a, b, c);
364       d += t1;
365       h  = t1 + t2;
366
367       t1 = g + Sum1 (d) + Ch (d, e, f) + k[t+1] + w[1];
368       t2 = Sum0 (h) + Maj (h, a, b);
369       c += t1;
370       g  = t1 + t2;
371
372       t1 = f + Sum1 (c) + Ch (c, d, e) + k[t+2] + w[2];
373       t2 = Sum0 (g) + Maj (g, h, a);
374       b += t1;
375       f  = t1 + t2;
376
377       t1 = e + Sum1 (b) + Ch (b, c, d) + k[t+3] + w[3];
378       t2 = Sum0 (f) + Maj (f, g, h);
379       a += t1;
380       e  = t1 + t2;
381
382       t1 = d + Sum1 (a) + Ch (a, b, c) + k[t+4] + w[4];
383       t2 = Sum0 (e) + Maj (e, f, g);
384       h += t1;
385       d  = t1 + t2;
386
387       t1 = c + Sum1 (h) + Ch (h, a, b) + k[t+5] + w[5];
388       t2 = Sum0 (d) + Maj (d, e, f);
389       g += t1;
390       c  = t1 + t2;
391
392       t1 = b + Sum1 (g) + Ch (g, h, a) + k[t+6] + w[6];
393       t2 = Sum0 (c) + Maj (c, d, e);
394       f += t1;
395       b  = t1 + t2;
396
397       t1 = a + Sum1 (f) + Ch (f, g, h) + k[t+7] + w[7];
398       t2 = Sum0 (b) + Maj (b, c, d);
399       e += t1;
400       a  = t1 + t2;
401
402       t1 = h + Sum1 (e) + Ch (e, f, g) + k[t+8] + w[8];
403       t2 = Sum0 (a) + Maj (a, b, c);
404       d += t1;
405       h  = t1 + t2;
406
407       t1 = g + Sum1 (d) + Ch (d, e, f) + k[t+9] + w[9];
408       t2 = Sum0 (h) + Maj (h, a, b);
409       c += t1;
410       g  = t1 + t2;
411
412       t1 = f + Sum1 (c) + Ch (c, d, e) + k[t+10] + w[10];
413       t2 = Sum0 (g) + Maj (g, h, a);
414       b += t1;
415       f  = t1 + t2;
416
417       t1 = e + Sum1 (b) + Ch (b, c, d) + k[t+11] + w[11];
418       t2 = Sum0 (f) + Maj (f, g, h);
419       a += t1;
420       e  = t1 + t2;
421
422       t1 = d + Sum1 (a) + Ch (a, b, c) + k[t+12] + w[12];
423       t2 = Sum0 (e) + Maj (e, f, g);
424       h += t1;
425       d  = t1 + t2;
426
427       t1 = c + Sum1 (h) + Ch (h, a, b) + k[t+13] + w[13];
428       t2 = Sum0 (d) + Maj (d, e, f);
429       g += t1;
430       c  = t1 + t2;
431
432       t1 = b + Sum1 (g) + Ch (g, h, a) + k[t+14] + w[14];
433       t2 = Sum0 (c) + Maj (c, d, e);
434       f += t1;
435       b  = t1 + t2;
436
437       t1 = a + Sum1 (f) + Ch (f, g, h) + k[t+15] + w[15];
438       t2 = Sum0 (b) + Maj (b, c, d);
439       e += t1;
440       a  = t1 + t2;
441
442       t += 16;
443 #endif
444     }
445
446   /* Update chaining vars.  */
447   hd->h0 += a;
448   hd->h1 += b;
449   hd->h2 += c;
450   hd->h3 += d;
451   hd->h4 += e;
452   hd->h5 += f;
453   hd->h6 += g;
454   hd->h7 += h;
455 }
456
457
458 /* Update the message digest with the contents
459  * of INBUF with length INLEN.
460  */
461 static void
462 sha512_write (void *context, const void *inbuf_arg, size_t inlen)
463 {
464   const unsigned char *inbuf = inbuf_arg;
465   SHA512_CONTEXT *hd = context;
466
467   if (hd->count == 128)
468     {                           /* flush the buffer */
469       transform (hd, hd->buf);
470       _gcry_burn_stack (256);
471       hd->count = 0;
472       hd->nblocks++;
473     }
474   if (!inbuf)
475     return;
476   if (hd->count)
477     {
478       for (; inlen && hd->count < 128; inlen--)
479         hd->buf[hd->count++] = *inbuf++;
480       sha512_write (context, NULL, 0);
481       if (!inlen)
482         return;
483     }
484
485   while (inlen >= 128)
486     {
487       transform (hd, inbuf);
488       hd->count = 0;
489       hd->nblocks++;
490       inlen -= 128;
491       inbuf += 128;
492     }
493   _gcry_burn_stack (256);
494   for (; inlen && hd->count < 128; inlen--)
495     hd->buf[hd->count++] = *inbuf++;
496 }
497
498
499 /* The routine final terminates the computation and
500  * returns the digest.
501  * The handle is prepared for a new cycle, but adding bytes to the
502  * handle will the destroy the returned buffer.
503  * Returns: 64 bytes representing the digest.  When used for sha384,
504  * we take the leftmost 48 of those bytes.
505  */
506
507 static void
508 sha512_final (void *context)
509 {
510   SHA512_CONTEXT *hd = context;
511   u64 t, msb, lsb;
512   byte *p;
513
514   sha512_write (context, NULL, 0); /* flush */ ;
515
516   t = hd->nblocks;
517   /* multiply by 128 to make a byte count */
518   lsb = t << 7;
519   msb = t >> 57;
520   /* add the count */
521   t = lsb;
522   if ((lsb += hd->count) < t)
523     msb++;
524   /* multiply by 8 to make a bit count */
525   t = lsb;
526   lsb <<= 3;
527   msb <<= 3;
528   msb |= t >> 61;
529
530   if (hd->count < 112)
531     {                           /* enough room */
532       hd->buf[hd->count++] = 0x80;      /* pad */
533       while (hd->count < 112)
534         hd->buf[hd->count++] = 0;       /* pad */
535     }
536   else
537     {                           /* need one extra block */
538       hd->buf[hd->count++] = 0x80;      /* pad character */
539       while (hd->count < 128)
540         hd->buf[hd->count++] = 0;
541       sha512_write (context, NULL, 0); /* flush */ ;
542       memset (hd->buf, 0, 112); /* fill next block with zeroes */
543     }
544   /* append the 128 bit count */
545   hd->buf[112] = msb >> 56;
546   hd->buf[113] = msb >> 48;
547   hd->buf[114] = msb >> 40;
548   hd->buf[115] = msb >> 32;
549   hd->buf[116] = msb >> 24;
550   hd->buf[117] = msb >> 16;
551   hd->buf[118] = msb >> 8;
552   hd->buf[119] = msb;
553
554   hd->buf[120] = lsb >> 56;
555   hd->buf[121] = lsb >> 48;
556   hd->buf[122] = lsb >> 40;
557   hd->buf[123] = lsb >> 32;
558   hd->buf[124] = lsb >> 24;
559   hd->buf[125] = lsb >> 16;
560   hd->buf[126] = lsb >> 8;
561   hd->buf[127] = lsb;
562   transform (hd, hd->buf);
563   _gcry_burn_stack (256);
564
565   p = hd->buf;
566 #ifdef WORDS_BIGENDIAN
567 #define X(a) do { *(u64*)p = hd->h##a ; p += 8; } while (0)
568 #else /* little endian */
569 #define X(a) do { *p++ = hd->h##a >> 56; *p++ = hd->h##a >> 48;       \
570                   *p++ = hd->h##a >> 40; *p++ = hd->h##a >> 32;       \
571                   *p++ = hd->h##a >> 24; *p++ = hd->h##a >> 16;       \
572                   *p++ = hd->h##a >> 8;  *p++ = hd->h##a; } while (0)
573 #endif
574   X (0);
575   X (1);
576   X (2);
577   X (3);
578   X (4);
579   X (5);
580   /* Note that these last two chunks are included even for SHA384.
581      We just ignore them. */
582   X (6);
583   X (7);
584 #undef X
585 }
586
587 static byte *
588 sha512_read (void *context)
589 {
590   SHA512_CONTEXT *hd = (SHA512_CONTEXT *) context;
591   return hd->buf;
592 }
593
594
595 \f
596 /*
597      Self-test section.
598  */
599
600
601 static gpg_err_code_t
602 selftests_sha384 (int extended, selftest_report_func_t report)
603 {
604   const char *what;
605   const char *errtxt;
606
607   what = "short string";
608   errtxt = _gcry_hash_selftest_check_one
609     (GCRY_MD_SHA384, 0,
610      "abc", 3,
611      "\xcb\x00\x75\x3f\x45\xa3\x5e\x8b\xb5\xa0\x3d\x69\x9a\xc6\x50\x07"
612      "\x27\x2c\x32\xab\x0e\xde\xd1\x63\x1a\x8b\x60\x5a\x43\xff\x5b\xed"
613      "\x80\x86\x07\x2b\xa1\xe7\xcc\x23\x58\xba\xec\xa1\x34\xc8\x25\xa7", 48);
614   if (errtxt)
615     goto failed;
616
617   if (extended)
618     {
619       what = "long string";
620       errtxt = _gcry_hash_selftest_check_one
621         (GCRY_MD_SHA384, 0,
622          "abcdefghbcdefghicdefghijdefghijkefghijklfghijklmghijklmn"
623          "hijklmnoijklmnopjklmnopqklmnopqrlmnopqrsmnopqrstnopqrstu", 112,
624          "\x09\x33\x0C\x33\xF7\x11\x47\xE8\x3D\x19\x2F\xC7\x82\xCD\x1B\x47"
625          "\x53\x11\x1B\x17\x3B\x3B\x05\xD2\x2F\xA0\x80\x86\xE3\xB0\xF7\x12"
626          "\xFC\xC7\xC7\x1A\x55\x7E\x2D\xB9\x66\xC3\xE9\xFA\x91\x74\x60\x39",
627          48);
628       if (errtxt)
629         goto failed;
630
631       what = "one million \"a\"";
632       errtxt = _gcry_hash_selftest_check_one
633         (GCRY_MD_SHA384, 1,
634          NULL, 0,
635          "\x9D\x0E\x18\x09\x71\x64\x74\xCB\x08\x6E\x83\x4E\x31\x0A\x4A\x1C"
636          "\xED\x14\x9E\x9C\x00\xF2\x48\x52\x79\x72\xCE\xC5\x70\x4C\x2A\x5B"
637          "\x07\xB8\xB3\xDC\x38\xEC\xC4\xEB\xAE\x97\xDD\xD8\x7F\x3D\x89\x85",
638          48);
639       if (errtxt)
640         goto failed;
641     }
642
643   return 0; /* Succeeded. */
644
645  failed:
646   if (report)
647     report ("digest", GCRY_MD_SHA384, what, errtxt);
648   return GPG_ERR_SELFTEST_FAILED;
649 }
650
651 static gpg_err_code_t
652 selftests_sha512 (int extended, selftest_report_func_t report)
653 {
654   const char *what;
655   const char *errtxt;
656
657   what = "short string";
658   errtxt = _gcry_hash_selftest_check_one
659     (GCRY_MD_SHA512, 0,
660      "abc", 3,
661      "\xDD\xAF\x35\xA1\x93\x61\x7A\xBA\xCC\x41\x73\x49\xAE\x20\x41\x31"
662      "\x12\xE6\xFA\x4E\x89\xA9\x7E\xA2\x0A\x9E\xEE\xE6\x4B\x55\xD3\x9A"
663      "\x21\x92\x99\x2A\x27\x4F\xC1\xA8\x36\xBA\x3C\x23\xA3\xFE\xEB\xBD"
664      "\x45\x4D\x44\x23\x64\x3C\xE8\x0E\x2A\x9A\xC9\x4F\xA5\x4C\xA4\x9F", 64);
665   if (errtxt)
666     goto failed;
667
668   if (extended)
669     {
670       what = "long string";
671       errtxt = _gcry_hash_selftest_check_one
672         (GCRY_MD_SHA512, 0,
673          "abcdefghbcdefghicdefghijdefghijkefghijklfghijklmghijklmn"
674          "hijklmnoijklmnopjklmnopqklmnopqrlmnopqrsmnopqrstnopqrstu", 112,
675          "\x8E\x95\x9B\x75\xDA\xE3\x13\xDA\x8C\xF4\xF7\x28\x14\xFC\x14\x3F"
676          "\x8F\x77\x79\xC6\xEB\x9F\x7F\xA1\x72\x99\xAE\xAD\xB6\x88\x90\x18"
677          "\x50\x1D\x28\x9E\x49\x00\xF7\xE4\x33\x1B\x99\xDE\xC4\xB5\x43\x3A"
678          "\xC7\xD3\x29\xEE\xB6\xDD\x26\x54\x5E\x96\xE5\x5B\x87\x4B\xE9\x09",
679          64);
680       if (errtxt)
681         goto failed;
682
683       what = "one million \"a\"";
684       errtxt = _gcry_hash_selftest_check_one
685         (GCRY_MD_SHA512, 1,
686          NULL, 0,
687          "\xE7\x18\x48\x3D\x0C\xE7\x69\x64\x4E\x2E\x42\xC7\xBC\x15\xB4\x63"
688          "\x8E\x1F\x98\xB1\x3B\x20\x44\x28\x56\x32\xA8\x03\xAF\xA9\x73\xEB"
689          "\xDE\x0F\xF2\x44\x87\x7E\xA6\x0A\x4C\xB0\x43\x2C\xE5\x77\xC3\x1B"
690          "\xEB\x00\x9C\x5C\x2C\x49\xAA\x2E\x4E\xAD\xB2\x17\xAD\x8C\xC0\x9B",
691          64);
692       if (errtxt)
693         goto failed;
694     }
695
696   return 0; /* Succeeded. */
697
698  failed:
699   if (report)
700     report ("digest", GCRY_MD_SHA512, what, errtxt);
701   return GPG_ERR_SELFTEST_FAILED;
702 }
703
704
705 /* Run a full self-test for ALGO and return 0 on success.  */
706 static gpg_err_code_t
707 run_selftests (int algo, int extended, selftest_report_func_t report)
708 {
709   gpg_err_code_t ec;
710
711   switch (algo)
712     {
713     case GCRY_MD_SHA384:
714       ec = selftests_sha384 (extended, report);
715       break;
716     case GCRY_MD_SHA512:
717       ec = selftests_sha512 (extended, report);
718       break;
719     default:
720       ec = GPG_ERR_DIGEST_ALGO;
721       break;
722
723     }
724   return ec;
725 }
726
727
728
729 \f
730 static byte sha512_asn[] =      /* Object ID is 2.16.840.1.101.3.4.2.3 */
731   {
732     0x30, 0x51, 0x30, 0x0d, 0x06, 0x09, 0x60, 0x86,
733     0x48, 0x01, 0x65, 0x03, 0x04, 0x02, 0x03, 0x05,
734     0x00, 0x04, 0x40
735   };
736
737 static gcry_md_oid_spec_t oid_spec_sha512[] =
738   {
739     { "2.16.840.1.101.3.4.2.3" },
740
741     /* PKCS#1 sha512WithRSAEncryption */
742     { "1.2.840.113549.1.1.13" },
743
744     { NULL }
745   };
746
747 gcry_md_spec_t _gcry_digest_spec_sha512 =
748   {
749     "SHA512", sha512_asn, DIM (sha512_asn), oid_spec_sha512, 64,
750     sha512_init, sha512_write, sha512_final, sha512_read,
751     sizeof (SHA512_CONTEXT),
752   };
753 md_extra_spec_t _gcry_digest_extraspec_sha512 =
754   {
755     run_selftests
756   };
757
758 static byte sha384_asn[] =      /* Object ID is 2.16.840.1.101.3.4.2.2 */
759   {
760     0x30, 0x41, 0x30, 0x0d, 0x06, 0x09, 0x60, 0x86,
761     0x48, 0x01, 0x65, 0x03, 0x04, 0x02, 0x02, 0x05,
762     0x00, 0x04, 0x30
763   };
764
765 static gcry_md_oid_spec_t oid_spec_sha384[] =
766   {
767     { "2.16.840.1.101.3.4.2.2" },
768
769     /* PKCS#1 sha384WithRSAEncryption */
770     { "1.2.840.113549.1.1.12" },
771
772     { NULL },
773   };
774
775 gcry_md_spec_t _gcry_digest_spec_sha384 =
776   {
777     "SHA384", sha384_asn, DIM (sha384_asn), oid_spec_sha384, 48,
778     sha384_init, sha512_write, sha512_final, sha512_read,
779     sizeof (SHA512_CONTEXT),
780   };
781 md_extra_spec_t _gcry_digest_extraspec_sha384 =
782   {
783     run_selftests
784   };