Enable AMD64 SHA512 implementations for WIN64
[libgcrypt.git] / cipher / md4.c
index d649f68..c9b4154 100644 (file)
@@ -1,5 +1,5 @@
 /* md4.c - MD4 Message-Digest Algorithm
- * Copyright (C) 2002 Free Software Foundation, Inc.
+ * Copyright (C) 2002, 2003 Free Software Foundation, Inc.
  *
  * This file is part of Libgcrypt.
  *
 #include <stdio.h>
 #include <stdlib.h>
 #include <string.h>
-#include <assert.h>
+
 #include "g10lib.h"
-#include "memory.h"
-#include "dynload.h"
+#include "cipher.h"
 
 #include "bithelp.h"
+#include "bufhelp.h"
+#include "hash-common.h"
 
 
 typedef struct {
+    gcry_md_block_ctx_t bctx;
     u32 A,B,C,D;         /* chaining variables */
-    u32  nblocks;
-    byte buf[64];
-    int  count;
 } MD4_CONTEXT;
 
+static unsigned int
+transform ( void *c, const unsigned char *data, size_t nblks );
 
 static void
-md4_init( MD4_CONTEXT *ctx )
+md4_init (void *context, unsigned int flags)
 {
-    ctx->A = 0x67452301;
-    ctx->B = 0xefcdab89;
-    ctx->C = 0x98badcfe;
-    ctx->D = 0x10325476;
+  MD4_CONTEXT *ctx = context;
 
-    ctx->nblocks = 0;
-    ctx->count = 0;
-}
+  (void)flags;
 
-static void
-burn_stack (int bytes)
-{
-    char buf[128];
-    
-    memset (buf, 0, sizeof buf);
-    bytes -= sizeof buf;
-    if (bytes > 0)
-        burn_stack (bytes);
+  ctx->A = 0x67452301;
+  ctx->B = 0xefcdab89;
+  ctx->C = 0x98badcfe;
+  ctx->D = 0x10325476;
+
+  ctx->bctx.nblocks = 0;
+  ctx->bctx.nblocks_high = 0;
+  ctx->bctx.count = 0;
+  ctx->bctx.blocksize = 64;
+  ctx->bctx.bwrite = transform;
 }
 
 #define F(x, y, z) ((z) ^ ((x) & ((y) ^ (z))))
@@ -98,140 +95,110 @@ burn_stack (int bytes)
 /****************
  * transform 64 bytes
  */
-static void
-transform( MD4_CONTEXT *ctx, byte *data )
+static unsigned int
+transform_blk ( void *c, const unsigned char *data )
 {
-    u32 in[16];
-    register u32 A = ctx->A;
-    register u32 B = ctx->B;
-    register u32 C = ctx->C;
-    register u32 D = ctx->D;
-
-#ifdef BIG_ENDIAN_HOST
-    { int i;
-      byte *p2, *p1;
-      for(i=0, p1=data, p2=(byte*)in; i < 16; i++, p2 += 4 ) {
-       p2[3] = *p1++;
-       p2[2] = *p1++;
-       p2[1] = *p1++;
-       p2[0] = *p1++;
-      }
-    }
-#else
-    memcpy (in, data, 64);
-#endif
-
-    /* Round 1.  */
+  MD4_CONTEXT *ctx = c;
+  u32 in[16];
+  register u32 A = ctx->A;
+  register u32 B = ctx->B;
+  register u32 C = ctx->C;
+  register u32 D = ctx->D;
+  int i;
+
+  for ( i = 0; i < 16; i++ )
+    in[i] = buf_get_le32(data + i * 4);
+
+  /* Round 1.  */
 #define function(a,b,c,d,k,s) a=rol(a+F(b,c,d)+in[k],s);
-          function(A,B,C,D, 0, 3);
-          function(D,A,B,C, 1, 7);
-          function(C,D,A,B, 2,11);
-          function(B,C,D,A, 3,19);
-          function(A,B,C,D, 4, 3);
-          function(D,A,B,C, 5, 7);
-          function(C,D,A,B, 6,11);
-          function(B,C,D,A, 7,19);
-          function(A,B,C,D, 8, 3);
-          function(D,A,B,C, 9, 7);
-          function(C,D,A,B,10,11);
-          function(B,C,D,A,11,19);
-          function(A,B,C,D,12, 3);
-          function(D,A,B,C,13, 7);
-          function(C,D,A,B,14,11);
-          function(B,C,D,A,15,19);
+  function(A,B,C,D, 0, 3);
+  function(D,A,B,C, 1, 7);
+  function(C,D,A,B, 2,11);
+  function(B,C,D,A, 3,19);
+  function(A,B,C,D, 4, 3);
+  function(D,A,B,C, 5, 7);
+  function(C,D,A,B, 6,11);
+  function(B,C,D,A, 7,19);
+  function(A,B,C,D, 8, 3);
+  function(D,A,B,C, 9, 7);
+  function(C,D,A,B,10,11);
+  function(B,C,D,A,11,19);
+  function(A,B,C,D,12, 3);
+  function(D,A,B,C,13, 7);
+  function(C,D,A,B,14,11);
+  function(B,C,D,A,15,19);
 
 #undef function
 
-    /* Round 2.  */
+  /* Round 2.  */
 #define function(a,b,c,d,k,s) a=rol(a+G(b,c,d)+in[k]+0x5a827999,s);
 
-          function(A,B,C,D, 0, 3);
-          function(D,A,B,C, 4, 5);
-          function(C,D,A,B, 8, 9);
-          function(B,C,D,A,12,13);
-          function(A,B,C,D, 1, 3);
-          function(D,A,B,C, 5, 5);
-          function(C,D,A,B, 9, 9);
-          function(B,C,D,A,13,13);
-          function(A,B,C,D, 2, 3);
-          function(D,A,B,C, 6, 5);
-          function(C,D,A,B,10, 9);
-          function(B,C,D,A,14,13);
-          function(A,B,C,D, 3, 3);
-          function(D,A,B,C, 7, 5);
-          function(C,D,A,B,11, 9);
-          function(B,C,D,A,15,13);
+  function(A,B,C,D, 0, 3);
+  function(D,A,B,C, 4, 5);
+  function(C,D,A,B, 8, 9);
+  function(B,C,D,A,12,13);
+  function(A,B,C,D, 1, 3);
+  function(D,A,B,C, 5, 5);
+  function(C,D,A,B, 9, 9);
+  function(B,C,D,A,13,13);
+  function(A,B,C,D, 2, 3);
+  function(D,A,B,C, 6, 5);
+  function(C,D,A,B,10, 9);
+  function(B,C,D,A,14,13);
+  function(A,B,C,D, 3, 3);
+  function(D,A,B,C, 7, 5);
+  function(C,D,A,B,11, 9);
+  function(B,C,D,A,15,13);
 
 #undef function
 
-    /* Round 3.  */
+  /* Round 3.  */
 #define function(a,b,c,d,k,s) a=rol(a+H(b,c,d)+in[k]+0x6ed9eba1,s);
 
-          function(A,B,C,D, 0, 3);
-          function(D,A,B,C, 8, 9);
-          function(C,D,A,B, 4,11);
-          function(B,C,D,A,12,15);
-          function(A,B,C,D, 2, 3);
-          function(D,A,B,C,10, 9);
-          function(C,D,A,B, 6,11);
-          function(B,C,D,A,14,15);
-          function(A,B,C,D, 1, 3);
-          function(D,A,B,C, 9, 9);
-          function(C,D,A,B, 5,11);
-          function(B,C,D,A,13,15);
-          function(A,B,C,D, 3, 3);
-          function(D,A,B,C,11, 9);
-          function(C,D,A,B, 7,11);
-          function(B,C,D,A,15,15);
-
-
-    /* Put checksum in context given as argument.  */
-    ctx->A += A;
-    ctx->B += B;
-    ctx->C += C;
-    ctx->D += D;
+  function(A,B,C,D, 0, 3);
+  function(D,A,B,C, 8, 9);
+  function(C,D,A,B, 4,11);
+  function(B,C,D,A,12,15);
+  function(A,B,C,D, 2, 3);
+  function(D,A,B,C,10, 9);
+  function(C,D,A,B, 6,11);
+  function(B,C,D,A,14,15);
+  function(A,B,C,D, 1, 3);
+  function(D,A,B,C, 9, 9);
+  function(C,D,A,B, 5,11);
+  function(B,C,D,A,13,15);
+  function(A,B,C,D, 3, 3);
+  function(D,A,B,C,11, 9);
+  function(C,D,A,B, 7,11);
+  function(B,C,D,A,15,15);
+
+
+  /* Put checksum in context given as argument.  */
+  ctx->A += A;
+  ctx->B += B;
+  ctx->C += C;
+  ctx->D += D;
+
+  return /*burn_stack*/ 80+6*sizeof(void*);
 }
 
 
-
-/* The routine updates the message-digest context to
- * account for the presence of each of the characters inBuf[0..inLen-1]
- * in the message whose digest is being computed.
- */
-static void
-md4_write( MD4_CONTEXT *hd, byte *inbuf, size_t inlen)
+static unsigned int
+transform ( void *c, const unsigned char *data, size_t nblks )
 {
-    if( hd->count == 64 ) { /* flush the buffer */
-       transform( hd, hd->buf );
-        burn_stack (80+6*sizeof(void*));
-       hd->count = 0;
-       hd->nblocks++;
-    }
-    if( !inbuf )
-       return;
-    if( hd->count ) {
-       for( ; inlen && hd->count < 64; inlen-- )
-           hd->buf[hd->count++] = *inbuf++;
-       md4_write( hd, NULL, 0 );
-       if( !inlen )
-           return;
-    }
-    burn_stack (80+6*sizeof(void*));
-
-    while( inlen >= 64 ) {
-       transform( hd, inbuf );
-       hd->count = 0;
-       hd->nblocks++;
-       inlen -= 64;
-       inbuf += 64;
+  unsigned int burn;
+
+  do
+    {
+      burn = transform_blk (c, data);
+      data += 64;
     }
-    for( ; inlen && hd->count < 64; inlen-- )
-       hd->buf[hd->count++] = *inbuf++;
+  while (--nblks);
 
+  return burn;
 }
 
 
-
 /* The routine final terminates the message-digest computation and
  * ends with the desired message digest in mdContext->digest[0...15].
  * The handle is prepared for a new MD4 cycle.
@@ -239,157 +206,86 @@ md4_write( MD4_CONTEXT *hd, byte *inbuf, size_t inlen)
  */
 
 static void
-md4_final( MD4_CONTEXT *hd )
+md4_final( void *context )
 {
-    u32 t, msb, lsb;
-    byte *p;
-
-    md4_write(hd, NULL, 0); /* flush */;
-
-    t = hd->nblocks;
-    /* multiply by 64 to make a byte count */
-    lsb = t << 6;
-    msb = t >> 26;
-    /* add the count */
-    t = lsb;
-    if( (lsb += hd->count) < t )
-       msb++;
-    /* multiply by 8 to make a bit count */
-    t = lsb;
-    lsb <<= 3;
-    msb <<= 3;
-    msb |= t >> 29;
-
-    if( hd->count < 56 ) { /* enough room */
-       hd->buf[hd->count++] = 0x80; /* pad */
-       while( hd->count < 56 )
-           hd->buf[hd->count++] = 0;  /* pad */
+  MD4_CONTEXT *hd = context;
+  u32 t, th, msb, lsb;
+  byte *p;
+  unsigned int burn;
+
+  _gcry_md_block_write(hd, NULL, 0); /* flush */;
+
+  t = hd->bctx.nblocks;
+  if (sizeof t == sizeof hd->bctx.nblocks)
+    th = hd->bctx.nblocks_high;
+  else
+    th = hd->bctx.nblocks >> 32;
+
+  /* multiply by 64 to make a byte count */
+  lsb = t << 6;
+  msb = (th << 6) | (t >> 26);
+  /* add the count */
+  t = lsb;
+  if( (lsb += hd->bctx.count) < t )
+    msb++;
+  /* multiply by 8 to make a bit count */
+  t = lsb;
+  lsb <<= 3;
+  msb <<= 3;
+  msb |= t >> 29;
+
+  if( hd->bctx.count < 56 )  /* enough room */
+    {
+      hd->bctx.buf[hd->bctx.count++] = 0x80; /* pad */
+      while( hd->bctx.count < 56 )
+        hd->bctx.buf[hd->bctx.count++] = 0;  /* pad */
     }
-    else { /* need one extra block */
-       hd->buf[hd->count++] = 0x80; /* pad character */
-       while( hd->count < 64 )
-           hd->buf[hd->count++] = 0;
-       md4_write(hd, NULL, 0);  /* flush */;
-       memset(hd->buf, 0, 56 ); /* fill next block with zeroes */
+  else /* need one extra block */
+    {
+      hd->bctx.buf[hd->bctx.count++] = 0x80; /* pad character */
+      while( hd->bctx.count < 64 )
+        hd->bctx.buf[hd->bctx.count++] = 0;
+      _gcry_md_block_write(hd, NULL, 0);  /* flush */;
+      memset(hd->bctx.buf, 0, 56 ); /* fill next block with zeroes */
     }
-    /* append the 64 bit count */
-    hd->buf[56] = lsb     ;
-    hd->buf[57] = lsb >>  8;
-    hd->buf[58] = lsb >> 16;
-    hd->buf[59] = lsb >> 24;
-    hd->buf[60] = msb     ;
-    hd->buf[61] = msb >>  8;
-    hd->buf[62] = msb >> 16;
-    hd->buf[63] = msb >> 24;
-    transform( hd, hd->buf );
-    burn_stack (80+6*sizeof(void*));
-
-    p = hd->buf;
-  #ifdef BIG_ENDIAN_HOST
-    #define X(a) do { *p++ = hd->a      ; *p++ = hd->a >> 8;      \
-                     *p++ = hd->a >> 16; *p++ = hd->a >> 24; } while(0)
-  #else /* little endian */
-    #define X(a) do { *(u32*)p = (*hd).a ; p += 4; } while(0)
-  #endif
-    X(A);
-    X(B);
-    X(C);
-    X(D);
-  #undef X
+  /* append the 64 bit count */
+  buf_put_le32(hd->bctx.buf + 56, lsb);
+  buf_put_le32(hd->bctx.buf + 60, msb);
+  burn = transform ( hd, hd->bctx.buf, 1 );
+  _gcry_burn_stack (burn);
+
+  p = hd->bctx.buf;
+#define X(a) do { buf_put_le32(p, hd->a); p += 4; } while(0)
+  X(A);
+  X(B);
+  X(C);
+  X(D);
+#undef X
 
 }
 
 static byte *
-md4_read( MD4_CONTEXT *hd )
-{
-    return hd->buf;
-}
-
-/****************
- * Return some information about the algorithm.  We need algo here to
- * distinguish different flavors of the algorithm.
- * Returns: A pointer to string describing the algorithm or NULL if
- *         the ALGO is invalid.
- */
-static const char *
-md4_get_info( int algo, size_t *contextsize,
-              byte **r_asnoid, int *r_asnlen, int *r_mdlen,
-              void (**r_init)( void *c ),
-              void (**r_write)( void *c, byte *buf, size_t nbytes ),
-              void (**r_final)( void *c ),
-              byte *(**r_read)( void *c )
-            )
-{
-    static byte asn[18] = /* Object ID is 1.2.840.113549.2.4 */
-                   { 0x30, 0x20, 0x30, 0x0c, 0x06, 0x08, 0x2a, 0x86,0x48,
-                     0x86, 0xf7, 0x0d, 0x02, 0x04, 0x05, 0x00, 0x04, 0x10 };
-
-    if( algo != 301 )
-       return NULL;
-
-    *contextsize = sizeof(MD4_CONTEXT);
-    *r_asnoid = asn;
-    *r_asnlen = DIM(asn);
-    *r_mdlen = 16;
-    *(void  (**)(MD4_CONTEXT *))r_init                = md4_init;
-    *(void  (**)(MD4_CONTEXT *, byte*, size_t))r_write = md4_write;
-    *(void  (**)(MD4_CONTEXT *))r_final               = md4_final;
-    *(byte *(**)(MD4_CONTEXT *))r_read                = md4_read;
-
-    return "MD4";
-}
-
-
-#ifndef IS_MODULE
-static
-#endif
-const char * const gnupgext_version = "MD4 ($Revision$)";
-
-static struct {
-    int class;
-    int version;
-    int  value;
-    void (*func)(void);
-} func_table[] = {
-    { 10, 1, 0, (void(*)(void))md4_get_info },
-    { 11, 1, 301 },
-};
-
-
-#ifndef IS_MODULE
-static
-#endif
-void *
-gnupgext_enum_func( int what, int *sequence, int *class, int *vers )
-{
-    void *ret;
-    int i = *sequence;
-
-    do {
-       if( i >= DIM(func_table) || i < 0 )
-           return NULL;
-       *class = func_table[i].class;
-       *vers  = func_table[i].version;
-       switch( *class ) {
-         case 11: case 21: case 31: ret = &func_table[i].value; break;
-         default:                   ret = func_table[i].func; break;
-       }
-       i++;
-    } while( what && what != *class );
-
-    *sequence = i;
-    return ret;
-}
-
-
-
-
-#ifndef IS_MODULE
-void
-_gcry_md4_constructor(void)
+md4_read (void *context)
 {
-    _gcry_register_internal_cipher_extension( gnupgext_version, gnupgext_enum_func );
+  MD4_CONTEXT *hd = context;
+  return hd->bctx.buf;
 }
-#endif
 
-/* end of file */
+static byte asn[18] = /* Object ID is 1.2.840.113549.2.4 */
+  { 0x30, 0x20, 0x30, 0x0c, 0x06, 0x08, 0x2a, 0x86,0x48,
+    0x86, 0xf7, 0x0d, 0x02, 0x04, 0x05, 0x00, 0x04, 0x10 };
+
+static gcry_md_oid_spec_t oid_spec_md4[] =
+  {
+    /* iso.member-body.us.rsadsi.digestAlgorithm.md4 */
+    { "1.2.840.113549.2.4" },
+    { NULL },
+  };
+
+gcry_md_spec_t _gcry_digest_spec_md4 =
+  {
+    GCRY_MD_MD4, {0, 0},
+    "MD4", asn, DIM (asn), oid_spec_md4,16,
+    md4_init, _gcry_md_block_write, md4_final, md4_read,
+    sizeof (MD4_CONTEXT)
+  };