crypto/poly1305/poly1305.c

   1 /*
   2  * Copyright 2015-2016 The OpenSSL Project Authors. All Rights Reserved.
   3  *
   4  * Licensed under the OpenSSL license (the "License").  You may not use
   5  * this file except in compliance with the License.  You can obtain a copy
   6  * in the file LICENSE in the source distribution or at
   7  * https://www.openssl.org/source/license.html
   8  */
   9
  10 #include <stdlib.h>
  11 #include <string.h>
  12 #include <openssl/crypto.h>
  13
  14 #include "internal/poly1305.h"
  15
  16 typedef void (*poly1305_blocks_f) (void *ctx, const unsigned char *inp,
  17                                    size_t len, unsigned int padbit);
  18 typedef void (*poly1305_emit_f) (void *ctx, unsigned char mac[16],
  19                                  const unsigned int nonce[4]);
  20
  21 struct poly1305_context {
  22     double opaque[24];  /* large enough to hold internal state, declared
  23                          * 'double' to ensure at least 64-bit invariant
  24                          * alignment across all platforms and
  25                          * configurations */
  26     unsigned int nonce[4];
  27     unsigned char data[POLY1305_BLOCK_SIZE];
  28     size_t num;
  29     struct {
  30         poly1305_blocks_f blocks;
  31         poly1305_emit_f emit;
  32     } func;
  33 };
  34
  35 size_t Poly1305_ctx_size ()
  36 {
  37     return sizeof(struct poly1305_context);
  38 }
  39
  40 /* pick 32-bit unsigned integer in little endian order */
  41 static unsigned int U8TOU32(const unsigned char *p)
  42 {
  43     return (((unsigned int)(p[0] & 0xff)) |
  44             ((unsigned int)(p[1] & 0xff) << 8) |
  45             ((unsigned int)(p[2] & 0xff) << 16) |
  46             ((unsigned int)(p[3] & 0xff) << 24));
  47 }
  48
  49 /*
  50  * Implementations can be classified by amount of significant bits in
  51  * words making up the multi-precision value, or in other words radix
  52  * or base of numerical representation, e.g. base 2^64, base 2^32,
  53  * base 2^26. Complementary characteristic is how wide is the result of
  54  * multiplication of pair of digits, e.g. it would take 128 bits to
  55  * accommodate multiplication result in base 2^64 case. These are used
  56  * interchangeably. To describe implementation that is. But interface
  57  * is designed to isolate this so that low-level primitives implemented
  58  * in assembly can be self-contained/self-coherent.
  59  */
  60 #ifndef POLY1305_ASM
  61 /*
  62  * Even though there is __int128 reference implementation targeting
  63  * 64-bit platforms provided below, it's not obvious that it's optimal
  64  * choice for every one of them. Depending on instruction set overall
  65  * amount of instructions can be comparable to one in __int64
  66  * implementation. Amount of multiplication instructions would be lower,
  67  * but not necessarily overall. And in out-of-order execution context,
  68  * it is the latter that can be crucial...
  69  *
  70  * On related note. Poly1305 author, D. J. Bernstein, discusses and
  71  * provides floating-point implementations of the algorithm in question.
  72  * It made a lot of sense by the time of introduction, because most
  73  * then-modern processors didn't have pipelined integer multiplier.
  74  * [Not to mention that some had non-constant timing for integer
  75  * multiplications.] Floating-point instructions on the other hand could
  76  * be issued every cycle, which allowed to achieve better performance.
  77  * Nowadays, with SIMD and/or out-or-order execution, shared or
  78  * even emulated FPU, it's more complicated, and floating-point
  79  * implementation is not necessarily optimal choice in every situation,
  80  * rather contrary...
  81  *
  82  *                                              <appro@openssl.org>
  83  */
  84
  85 typedef unsigned int u32;
  86
  87 /*
  88  * poly1305_blocks processes a multiple of POLY1305_BLOCK_SIZE blocks
  89  * of |inp| no longer than |len|. Behaviour for |len| not divisible by
  90  * block size is unspecified in general case, even though in reference
  91  * implementation the trailing chunk is simply ignored. Per algorithm
  92  * specification, every input block, complete or last partial, is to be
  93  * padded with a bit past most significant byte. The latter kind is then
  94  * padded with zeros till block size. This last partial block padding
  95  * is caller(*)'s responsibility, and because of this the last partial
  96  * block is always processed with separate call with |len| set to
  97  * POLY1305_BLOCK_SIZE and |padbit| to 0. In all other cases |padbit|
  98  * should be set to 1 to perform implicit padding with 128th bit.
  99  * poly1305_blocks does not actually check for this constraint though,
 100  * it's caller(*)'s responsibility to comply.
 101  *
 102  * (*)  In the context "caller" is not application code, but higher
 103  *      level Poly1305_* from this very module, so that quirks are
 104  *      handled locally.
 105  */
 106 static void
 107 poly1305_blocks(void *ctx, const unsigned char *inp, size_t len, u32 padbit);
 108
 109 /*
 110  * Type-agnostic "rip-off" from constant_time_locl.h
 111  */
 112 # define CONSTANT_TIME_CARRY(a,b) ( \
 113          (a ^ ((a ^ b) | ((a - b) ^ b))) >> (sizeof(a) * 8 - 1) \
 114          )
 115
 116 # if !defined(PEDANTIC) && \
 117      (defined(__SIZEOF_INT128__) && __SIZEOF_INT128__==16) && \
 118      (defined(__SIZEOF_LONG__) && __SIZEOF_LONG__==8)
 119
 120 typedef unsigned long u64;
 121 typedef unsigned __int128 u128;
 122
 123 typedef struct {
 124     u64 h[3];
 125     u64 r[2];
 126 } poly1305_internal;
 127
 128 /* pick 32-bit unsigned integer in little endian order */
 129 static u64 U8TOU64(const unsigned char *p)
 130 {
 131     return (((u64)(p[0] & 0xff)) |
 132             ((u64)(p[1] & 0xff) << 8) |
 133             ((u64)(p[2] & 0xff) << 16) |
 134             ((u64)(p[3] & 0xff) << 24) |
 135             ((u64)(p[4] & 0xff) << 32) |
 136             ((u64)(p[5] & 0xff) << 40) |
 137             ((u64)(p[6] & 0xff) << 48) |
 138             ((u64)(p[7] & 0xff) << 56));
 139 }
 140
 141 /* store a 32-bit unsigned integer in little endian */
 142 static void U64TO8(unsigned char *p, u64 v)
 143 {
 144     p[0] = (unsigned char)((v) & 0xff);
 145     p[1] = (unsigned char)((v >> 8) & 0xff);
 146     p[2] = (unsigned char)((v >> 16) & 0xff);
 147     p[3] = (unsigned char)((v >> 24) & 0xff);
 148     p[4] = (unsigned char)((v >> 32) & 0xff);
 149     p[5] = (unsigned char)((v >> 40) & 0xff);
 150     p[6] = (unsigned char)((v >> 48) & 0xff);
 151     p[7] = (unsigned char)((v >> 56) & 0xff);
 152 }
 153
 154 static void poly1305_init(void *ctx, const unsigned char key[16])
 155 {
 156     poly1305_internal *st = (poly1305_internal *) ctx;
 157
 158     /* h = 0 */
 159     st->h[0] = 0;
 160     st->h[1] = 0;
 161     st->h[2] = 0;
 162
 163     /* r &= 0xffffffc0ffffffc0ffffffc0fffffff */
 164     st->r[0] = U8TOU64(&key[0]) & 0x0ffffffc0fffffff;
 165     st->r[1] = U8TOU64(&key[8]) & 0x0ffffffc0ffffffc;
 166 }
 167
 168 static void
 169 poly1305_blocks(void *ctx, const unsigned char *inp, size_t len, u32 padbit)
 170 {
 171     poly1305_internal *st = (poly1305_internal *)ctx;
 172     u64 r0, r1;
 173     u64 s1;
 174     u64 h0, h1, h2, c;
 175     u128 d0, d1;
 176
 177     r0 = st->r[0];
 178     r1 = st->r[1];
 179
 180     s1 = r1 + (r1 >> 2);
 181
 182     h0 = st->h[0];
 183     h1 = st->h[1];
 184     h2 = st->h[2];
 185
 186     while (len >= POLY1305_BLOCK_SIZE) {
 187         /* h += m[i] */
 188         h0 = (u64)(d0 = (u128)h0 + U8TOU64(inp + 0));
 189         h1 = (u64)(d1 = (u128)h1 + (d0 >> 64) + U8TOU64(inp + 8));
 190         /*
 191          * padbit can be zero only when original len was
 192          * POLY1306_BLOCK_SIZE, but we don't check
 193          */
 194         h2 += (u64)(d1 >> 64) + padbit;
 195
 196         /* h *= r "%" p, where "%" stands for "partial remainder" */
 197         d0 = ((u128)h0 * r0) +
 198              ((u128)h1 * s1);
 199         d1 = ((u128)h0 * r1) +
 200              ((u128)h1 * r0) +
 201              (h2 * s1);
 202         h2 = (h2 * r0);
 203
 204         /* last reduction step: */
 205         /* a) h2:h0 = h2<<128 + d1<<64 + d0 */
 206         h0 = (u64)d0;
 207         h1 = (u64)(d1 += d0 >> 64);
 208         h2 += (u64)(d1 >> 64);
 209         /* b) (h2:h0 += (h2:h0>>130) * 5) %= 2^130 */
 210         c = (h2 >> 2) + (h2 & ~3UL);
 211         h2 &= 3;
 212         h0 += c;
 213         h1 += (c = CONSTANT_TIME_CARRY(h0,c));
 214         h2 += CONSTANT_TIME_CARRY(h1,c);
 215         /*
 216          * Occasional overflows to 3rd bit of h2 are taken care of
 217          * "naturally". If after this point we end up at the top of
 218          * this loop, then the overflow bit will be accounted for
 219          * in next iteration. If we end up in poly1305_emit, then
 220          * comparison to modulus below will still count as "carry
 221          * into 131st bit", so that properly reduced value will be
 222          * picked in conditional move.
 223          */
 224
 225         inp += POLY1305_BLOCK_SIZE;
 226         len -= POLY1305_BLOCK_SIZE;
 227     }
 228
 229     st->h[0] = h0;
 230     st->h[1] = h1;
 231     st->h[2] = h2;
 232 }
 233
 234 static void poly1305_emit(void *ctx, unsigned char mac[16],
 235                           const u32 nonce[4])
 236 {
 237     poly1305_internal *st = (poly1305_internal *) ctx;
 238     u64 h0, h1, h2;
 239     u64 g0, g1, g2;
 240     u128 t;
 241     u64 mask;
 242
 243     h0 = st->h[0];
 244     h1 = st->h[1];
 245     h2 = st->h[2];
 246
 247     /* compare to modulus by computing h + -p */
 248     g0 = (u64)(t = (u128)h0 + 5);
 249     g1 = (u64)(t = (u128)h1 + (t >> 64));
 250     g2 = h2 + (u64)(t >> 64);
 251
 252     /* if there was carry into 131st bit, h1:h0 = g1:g0 */
 253     mask = 0 - (g2 >> 2);
 254     g0 &= mask;
 255     g1 &= mask;
 256     mask = ~mask;
 257     h0 = (h0 & mask) | g0;
 258     h1 = (h1 & mask) | g1;
 259
 260     /* mac = (h + nonce) % (2^128) */
 261     h0 = (u64)(t = (u128)h0 + nonce[0] + ((u64)nonce[1]<<32));
 262     h1 = (u64)(t = (u128)h1 + nonce[2] + ((u64)nonce[3]<<32) + (t >> 64));
 263
 264     U64TO8(mac + 0, h0);
 265     U64TO8(mac + 8, h1);
 266 }
 267
 268 # else
 269
 270 #  if defined(_WIN32) && !defined(__MINGW32__)
 271 typedef unsigned __int64 u64;
 272 #  elif defined(__arch64__)
 273 typedef unsigned long u64;
 274 #  else
 275 typedef unsigned long long u64;
 276 #  endif
 277
 278 typedef struct {
 279     u32 h[5];
 280     u32 r[4];
 281 } poly1305_internal;
 282
 283 /* store a 32-bit unsigned integer in little endian */
 284 static void U32TO8(unsigned char *p, unsigned int v)
 285 {
 286     p[0] = (unsigned char)((v) & 0xff);
 287     p[1] = (unsigned char)((v >> 8) & 0xff);
 288     p[2] = (unsigned char)((v >> 16) & 0xff);
 289     p[3] = (unsigned char)((v >> 24) & 0xff);
 290 }
 291
 292 static void poly1305_init(void *ctx, const unsigned char key[16])
 293 {
 294     poly1305_internal *st = (poly1305_internal *) ctx;
 295
 296     /* h = 0 */
 297     st->h[0] = 0;
 298     st->h[1] = 0;
 299     st->h[2] = 0;
 300     st->h[3] = 0;
 301     st->h[4] = 0;
 302
 303     /* r &= 0xffffffc0ffffffc0ffffffc0fffffff */
 304     st->r[0] = U8TOU32(&key[0]) & 0x0fffffff;
 305     st->r[1] = U8TOU32(&key[4]) & 0x0ffffffc;
 306     st->r[2] = U8TOU32(&key[8]) & 0x0ffffffc;
 307     st->r[3] = U8TOU32(&key[12]) & 0x0ffffffc;
 308 }
 309
 310 static void
 311 poly1305_blocks(void *ctx, const unsigned char *inp, size_t len, u32 padbit)
 312 {
 313     poly1305_internal *st = (poly1305_internal *)ctx;
 314     u32 r0, r1, r2, r3;
 315     u32 s1, s2, s3;
 316     u32 h0, h1, h2, h3, h4, c;
 317     u64 d0, d1, d2, d3;
 318
 319     r0 = st->r[0];
 320     r1 = st->r[1];
 321     r2 = st->r[2];
 322     r3 = st->r[3];
 323
 324     s1 = r1 + (r1 >> 2);
 325     s2 = r2 + (r2 >> 2);
 326     s3 = r3 + (r3 >> 2);
 327
 328     h0 = st->h[0];
 329     h1 = st->h[1];
 330     h2 = st->h[2];
 331     h3 = st->h[3];
 332     h4 = st->h[4];
 333
 334     while (len >= POLY1305_BLOCK_SIZE) {
 335         /* h += m[i] */
 336         h0 = (u32)(d0 = (u64)h0 + U8TOU32(inp + 0));
 337         h1 = (u32)(d1 = (u64)h1 + (d0 >> 32) + U8TOU32(inp + 4));
 338         h2 = (u32)(d2 = (u64)h2 + (d1 >> 32) + U8TOU32(inp + 8));
 339         h3 = (u32)(d3 = (u64)h3 + (d2 >> 32) + U8TOU32(inp + 12));
 340         h4 += (u32)(d3 >> 32) + padbit;
 341
 342         /* h *= r "%" p, where "%" stands for "partial remainder" */
 343         d0 = ((u64)h0 * r0) +
 344              ((u64)h1 * s3) +
 345              ((u64)h2 * s2) +
 346              ((u64)h3 * s1);
 347         d1 = ((u64)h0 * r1) +
 348              ((u64)h1 * r0) +
 349              ((u64)h2 * s3) +
 350              ((u64)h3 * s2) +
 351              (h4 * s1);
 352         d2 = ((u64)h0 * r2) +
 353              ((u64)h1 * r1) +
 354              ((u64)h2 * r0) +
 355              ((u64)h3 * s3) +
 356              (h4 * s2);
 357         d3 = ((u64)h0 * r3) +
 358              ((u64)h1 * r2) +
 359              ((u64)h2 * r1) +
 360              ((u64)h3 * r0) +
 361              (h4 * s3);
 362         h4 = (h4 * r0);
 363
 364         /* last reduction step: */
 365         /* a) h4:h0 = h4<<128 + d3<<96 + d2<<64 + d1<<32 + d0 */
 366         h0 = (u32)d0;
 367         h1 = (u32)(d1 += d0 >> 32);
 368         h2 = (u32)(d2 += d1 >> 32);
 369         h3 = (u32)(d3 += d2 >> 32);
 370         h4 += (u32)(d3 >> 32);
 371         /* b) (h4:h0 += (h4:h0>>130) * 5) %= 2^130 */
 372         c = (h4 >> 2) + (h4 & ~3U);
 373         h4 &= 3;
 374         h0 += c;
 375         h1 += (c = CONSTANT_TIME_CARRY(h0,c));
 376         h2 += (c = CONSTANT_TIME_CARRY(h1,c));
 377         h3 += (c = CONSTANT_TIME_CARRY(h2,c));
 378         h4 += CONSTANT_TIME_CARRY(h3,c);
 379         /*
 380          * Occasional overflows to 3rd bit of h4 are taken care of
 381          * "naturally". If after this point we end up at the top of
 382          * this loop, then the overflow bit will be accounted for
 383          * in next iteration. If we end up in poly1305_emit, then
 384          * comparison to modulus below will still count as "carry
 385          * into 131st bit", so that properly reduced value will be
 386          * picked in conditional move.
 387          */
 388
 389         inp += POLY1305_BLOCK_SIZE;
 390         len -= POLY1305_BLOCK_SIZE;
 391     }
 392
 393     st->h[0] = h0;
 394     st->h[1] = h1;
 395     st->h[2] = h2;
 396     st->h[3] = h3;
 397     st->h[4] = h4;
 398 }
 399
 400 static void poly1305_emit(void *ctx, unsigned char mac[16],
 401                           const u32 nonce[4])
 402 {
 403     poly1305_internal *st = (poly1305_internal *) ctx;
 404     u32 h0, h1, h2, h3, h4;
 405     u32 g0, g1, g2, g3, g4;
 406     u64 t;
 407     u32 mask;
 408
 409     h0 = st->h[0];
 410     h1 = st->h[1];
 411     h2 = st->h[2];
 412     h3 = st->h[3];
 413     h4 = st->h[4];
 414
 415     /* compare to modulus by computing h + -p */
 416     g0 = (u32)(t = (u64)h0 + 5);
 417     g1 = (u32)(t = (u64)h1 + (t >> 32));
 418     g2 = (u32)(t = (u64)h2 + (t >> 32));
 419     g3 = (u32)(t = (u64)h3 + (t >> 32));
 420     g4 = h4 + (u32)(t >> 32);
 421
 422     /* if there was carry into 131st bit, h3:h0 = g3:g0 */
 423     mask = 0 - (g4 >> 2);
 424     g0 &= mask;
 425     g1 &= mask;
 426     g2 &= mask;
 427     g3 &= mask;
 428     mask = ~mask;
 429     h0 = (h0 & mask) | g0;
 430     h1 = (h1 & mask) | g1;
 431     h2 = (h2 & mask) | g2;
 432     h3 = (h3 & mask) | g3;
 433
 434     /* mac = (h + nonce) % (2^128) */
 435     h0 = (u32)(t = (u64)h0 + nonce[0]);
 436     h1 = (u32)(t = (u64)h1 + (t >> 32) + nonce[1]);
 437     h2 = (u32)(t = (u64)h2 + (t >> 32) + nonce[2]);
 438     h3 = (u32)(t = (u64)h3 + (t >> 32) + nonce[3]);
 439
 440     U32TO8(mac + 0, h0);
 441     U32TO8(mac + 4, h1);
 442     U32TO8(mac + 8, h2);
 443     U32TO8(mac + 12, h3);
 444 }
 445 # endif
 446 #else
 447 int poly1305_init(void *ctx, const unsigned char key[16], void *func);
 448 void poly1305_blocks(void *ctx, const unsigned char *inp, size_t len,
 449                      unsigned int padbit);
 450 void poly1305_emit(void *ctx, unsigned char mac[16],
 451                    const unsigned int nonce[4]);
 452 #endif
 453
 454 void Poly1305_Init(POLY1305 *ctx, const unsigned char key[32])
 455 {
 456     ctx->nonce[0] = U8TOU32(&key[16]);
 457     ctx->nonce[1] = U8TOU32(&key[20]);
 458     ctx->nonce[2] = U8TOU32(&key[24]);
 459     ctx->nonce[3] = U8TOU32(&key[28]);
 460
 461 #ifndef POLY1305_ASM
 462     poly1305_init(ctx->opaque, key);
 463 #else
 464     /*
 465      * Unlike reference poly1305_init assembly counterpart is expected
 466      * to return a value: non-zero if it initializes ctx->func, and zero
 467      * otherwise. Latter is to simplify assembly in cases when there no
 468      * multiple code paths to switch between.
 469      */
 470     if (!poly1305_init(ctx->opaque, key, &ctx->func)) {
 471         ctx->func.blocks = poly1305_blocks;
 472         ctx->func.emit = poly1305_emit;
 473     }
 474 #endif
 475
 476     ctx->num = 0;
 477
 478 }
 479
 480 #ifdef POLY1305_ASM
 481 /*
 482  * This "eclipses" poly1305_blocks and poly1305_emit, but it's
 483  * conscious choice imposed by -Wshadow compiler warnings.
 484  */
 485 # define poly1305_blocks (*poly1305_blocks_p)
 486 # define poly1305_emit   (*poly1305_emit_p)
 487 #endif
 488
 489 void Poly1305_Update(POLY1305 *ctx, const unsigned char *inp, size_t len)
 490 {
 491 #ifdef POLY1305_ASM
 492     /*
 493      * As documented, poly1305_blocks is never called with input
 494      * longer than single block and padbit argument set to 0. This
 495      * property is fluently used in assembly modules to optimize
 496      * padbit handling on loop boundary.
 497      */
 498     poly1305_blocks_f poly1305_blocks_p = ctx->func.blocks;
 499 #endif
 500     size_t rem, num;
 501
 502     if ((num = ctx->num)) {
 503         rem = POLY1305_BLOCK_SIZE - num;
 504         if (len >= rem) {
 505             memcpy(ctx->data + num, inp, rem);
 506             poly1305_blocks(ctx->opaque, ctx->data, POLY1305_BLOCK_SIZE, 1);
 507             inp += rem;
 508             len -= rem;
 509         } else {
 510             /* Still not enough data to process a block. */
 511             memcpy(ctx->data + num, inp, len);
 512             ctx->num = num + len;
 513             return;
 514         }
 515     }
 516
 517     rem = len % POLY1305_BLOCK_SIZE;
 518     len -= rem;
 519
 520     if (len >= POLY1305_BLOCK_SIZE) {
 521         poly1305_blocks(ctx->opaque, inp, len, 1);
 522         inp += len;
 523     }
 524
 525     if (rem)
 526         memcpy(ctx->data, inp, rem);
 527
 528     ctx->num = rem;
 529 }
 530
 531 void Poly1305_Final(POLY1305 *ctx, unsigned char mac[16])
 532 {
 533 #ifdef POLY1305_ASM
 534     poly1305_blocks_f poly1305_blocks_p = ctx->func.blocks;
 535     poly1305_emit_f poly1305_emit_p = ctx->func.emit;
 536 #endif
 537     size_t num;
 538
 539     if ((num = ctx->num)) {
 540         ctx->data[num++] = 1;   /* pad bit */
 541         while (num < POLY1305_BLOCK_SIZE)
 542             ctx->data[num++] = 0;
 543         poly1305_blocks(ctx->opaque, ctx->data, POLY1305_BLOCK_SIZE, 0);
 544     }
 545
 546     poly1305_emit(ctx->opaque, mac, ctx->nonce);
 547
 548     /* zero out the state */
 549     OPENSSL_cleanse(ctx, sizeof(*ctx));
 550 }
 551
 552 #ifdef SELFTEST
 553 #include <stdio.h>
 554
 555 struct poly1305_test {
 556     const char *inputhex;
 557     const char *keyhex;
 558     const char *outhex;
 559 };
 560
 561 static const struct poly1305_test poly1305_tests[] = {
 562     /*
 563      * RFC7539
 564      */
 565     {
 566      "43727970746f6772617068696320466f72756d2052657365617263682047726f"
 567      "7570",
 568      "85d6be7857556d337f4452fe42d506a8""0103808afb0db2fd4abff6af4149f51b",
 569      "a8061dc1305136c6c22b8baf0c0127a9"
 570     },
 571     /*
 572      * test vectors from "The Poly1305-AES message-authentication code"
 573      */
 574     {
 575      "f3f6",
 576      "851fc40c3467ac0be05cc20404f3f700""580b3b0f9447bb1e69d095b5928b6dbc",
 577      "f4c633c3044fc145f84f335cb81953de"
 578     },
 579     {
 580      "",
 581      "a0f3080000f46400d0c7e9076c834403""dd3fab2251f11ac759f0887129cc2ee7",
 582      "dd3fab2251f11ac759f0887129cc2ee7"
 583     },
 584     {
 585      "663cea190ffb83d89593f3f476b6bc24d7e679107ea26adb8caf6652d0656136",
 586      "48443d0bb0d21109c89a100b5ce2c208""83149c69b561dd88298a1798b10716ef",
 587      "0ee1c16bb73f0f4fd19881753c01cdbe"
 588     },
 589     {
 590      "ab0812724a7f1e342742cbed374d94d136c6b8795d45b3819830f2c04491faf0"
 591      "990c62e48b8018b2c3e4a0fa3134cb67fa83e158c994d961c4cb21095c1bf9",
 592      "12976a08c4426d0ce8a82407c4f48207""80f8c20aa71202d1e29179cbcb555a57",
 593      "5154ad0d2cb26e01274fc51148491f1b"
 594     },
 595     /*
 596      * self-generated vectors exercise "significant" lengths, such that
 597      * are handled by different code paths
 598      */
 599     {
 600      "ab0812724a7f1e342742cbed374d94d136c6b8795d45b3819830f2c04491faf0"
 601      "990c62e48b8018b2c3e4a0fa3134cb67fa83e158c994d961c4cb21095c1bf9af",
 602      "12976a08c4426d0ce8a82407c4f48207""80f8c20aa71202d1e29179cbcb555a57",
 603      "812059a5da198637cac7c4a631bee466"
 604     },
 605     {
 606      "ab0812724a7f1e342742cbed374d94d136c6b8795d45b3819830f2c04491faf0"
 607      "990c62e48b8018b2c3e4a0fa3134cb67",
 608      "12976a08c4426d0ce8a82407c4f48207""80f8c20aa71202d1e29179cbcb555a57",
 609      "5b88d7f6228b11e2e28579a5c0c1f761"
 610     },
 611     {
 612      "ab0812724a7f1e342742cbed374d94d136c6b8795d45b3819830f2c04491faf0"
 613      "990c62e48b8018b2c3e4a0fa3134cb67fa83e158c994d961c4cb21095c1bf9af"
 614      "663cea190ffb83d89593f3f476b6bc24d7e679107ea26adb8caf6652d0656136",
 615      "12976a08c4426d0ce8a82407c4f48207""80f8c20aa71202d1e29179cbcb555a57",
 616      "bbb613b2b6d753ba07395b916aaece15"
 617     },
 618     {
 619      "ab0812724a7f1e342742cbed374d94d136c6b8795d45b3819830f2c04491faf0"
 620      "990c62e48b8018b2c3e4a0fa3134cb67fa83e158c994d961c4cb21095c1bf9af"
 621      "48443d0bb0d21109c89a100b5ce2c20883149c69b561dd88298a1798b10716ef"
 622      "663cea190ffb83d89593f3f476b6bc24",
 623      "12976a08c4426d0ce8a82407c4f48207""80f8c20aa71202d1e29179cbcb555a57",
 624      "c794d7057d1778c4bbee0a39b3d97342"
 625     },
 626     {
 627      "ab0812724a7f1e342742cbed374d94d136c6b8795d45b3819830f2c04491faf0"
 628      "990c62e48b8018b2c3e4a0fa3134cb67fa83e158c994d961c4cb21095c1bf9af"
 629      "48443d0bb0d21109c89a100b5ce2c20883149c69b561dd88298a1798b10716ef"
 630      "663cea190ffb83d89593f3f476b6bc24d7e679107ea26adb8caf6652d0656136",
 631      "12976a08c4426d0ce8a82407c4f48207""80f8c20aa71202d1e29179cbcb555a57",
 632      "ffbcb9b371423152d7fca5ad042fbaa9"
 633     },
 634     {
 635      "ab0812724a7f1e342742cbed374d94d136c6b8795d45b3819830f2c04491faf0"
 636      "990c62e48b8018b2c3e4a0fa3134cb67fa83e158c994d961c4cb21095c1bf9af"
 637      "48443d0bb0d21109c89a100b5ce2c20883149c69b561dd88298a1798b10716ef"
 638      "663cea190ffb83d89593f3f476b6bc24d7e679107ea26adb8caf6652d0656136"
 639      "812059a5da198637cac7c4a631bee466",
 640      "12976a08c4426d0ce8a82407c4f48207""80f8c20aa71202d1e29179cbcb555a57",
 641      "069ed6b8ef0f207b3e243bb1019fe632"
 642     },
 643     {
 644      "ab0812724a7f1e342742cbed374d94d136c6b8795d45b3819830f2c04491faf0"
 645      "990c62e48b8018b2c3e4a0fa3134cb67fa83e158c994d961c4cb21095c1bf9af"
 646      "48443d0bb0d21109c89a100b5ce2c20883149c69b561dd88298a1798b10716ef"
 647      "663cea190ffb83d89593f3f476b6bc24d7e679107ea26adb8caf6652d0656136"
 648      "812059a5da198637cac7c4a631bee4665b88d7f6228b11e2e28579a5c0c1f761",
 649      "12976a08c4426d0ce8a82407c4f48207""80f8c20aa71202d1e29179cbcb555a57",
 650      "cca339d9a45fa2368c2c68b3a4179133"
 651     },
 652     {
 653      "ab0812724a7f1e342742cbed374d94d136c6b8795d45b3819830f2c04491faf0"
 654      "990c62e48b8018b2c3e4a0fa3134cb67fa83e158c994d961c4cb21095c1bf9af"
 655      "48443d0bb0d21109c89a100b5ce2c20883149c69b561dd88298a1798b10716ef"
 656      "663cea190ffb83d89593f3f476b6bc24d7e679107ea26adb8caf6652d0656136"
 657      "812059a5da198637cac7c4a631bee4665b88d7f6228b11e2e28579a5c0c1f761"
 658      "ab0812724a7f1e342742cbed374d94d136c6b8795d45b3819830f2c04491faf0"
 659      "990c62e48b8018b2c3e4a0fa3134cb67fa83e158c994d961c4cb21095c1bf9af"
 660      "48443d0bb0d21109c89a100b5ce2c20883149c69b561dd88298a1798b10716ef"
 661      "663cea190ffb83d89593f3f476b6bc24d7e679107ea26adb8caf6652d0656136",
 662      "12976a08c4426d0ce8a82407c4f48207""80f8c20aa71202d1e29179cbcb555a57",
 663      "53f6e828a2f0fe0ee815bf0bd5841a34"
 664     },
 665     {
 666      "ab0812724a7f1e342742cbed374d94d136c6b8795d45b3819830f2c04491faf0"
 667      "990c62e48b8018b2c3e4a0fa3134cb67fa83e158c994d961c4cb21095c1bf9af"
 668      "48443d0bb0d21109c89a100b5ce2c20883149c69b561dd88298a1798b10716ef"
 669      "663cea190ffb83d89593f3f476b6bc24d7e679107ea26adb8caf6652d0656136"
 670      "812059a5da198637cac7c4a631bee4665b88d7f6228b11e2e28579a5c0c1f761"
 671      "ab0812724a7f1e342742cbed374d94d136c6b8795d45b3819830f2c04491faf0"
 672      "990c62e48b8018b2c3e4a0fa3134cb67fa83e158c994d961c4cb21095c1bf9af"
 673      "48443d0bb0d21109c89a100b5ce2c20883149c69b561dd88298a1798b10716ef"
 674      "663cea190ffb83d89593f3f476b6bc24d7e679107ea26adb8caf6652d0656136"
 675      "812059a5da198637cac7c4a631bee4665b88d7f6228b11e2e28579a5c0c1f761",
 676      "12976a08c4426d0ce8a82407c4f48207""80f8c20aa71202d1e29179cbcb555a57",
 677      "b846d44e9bbd53cedffbfbb6b7fa4933"
 678     },
 679     /*
 680      * 4th power of the key spills to 131th bit in SIMD key setup
 681      */
 682     {
 683      "ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff"
 684      "ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff"
 685      "ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff"
 686      "ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff"
 687      "ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff"
 688      "ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff"
 689      "ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff"
 690      "ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff",
 691      "ad628107e8351d0f2c231a05dc4a4106""00000000000000000000000000000000",
 692      "07145a4c02fe5fa32036de68fabe9066"
 693     },
 694     {
 695     /*
 696      * poly1305_ieee754.c failed this in final stage
 697      */
 698      "842364e156336c0998b933a6237726180d9e3fdcbde4cd5d17080fc3beb49614"
 699      "d7122c037463ff104d73f19c12704628d417c4c54a3fe30d3c3d7714382d43b0"
 700      "382a50a5dee54be844b076e8df88201a1cd43b90eb21643fa96f39b518aa8340"
 701      "c942ff3c31baf7c9bdbf0f31ae3fa096bf8c63030609829fe72e179824890bc8"
 702      "e08c315c1cce2a83144dbbff09f74e3efc770b54d0984a8f19b14719e6363564"
 703      "1d6b1eedf63efbf080e1783d32445412114c20de0b837a0dfa33d6b82825fff4"
 704      "4c9a70ea54ce47f07df698e6b03323b53079364a5fc3e9dd034392bdde86dccd"
 705      "da94321c5e44060489336cb65bf3989c36f7282c2f5d2b882c171e74",
 706      "95d5c005503e510d8cd0aa072c4a4d06""6eabc52d11653df47fbf63ab198bcc26",
 707      "f248312e578d9d58f8b7bb4d19105431"
 708     },
 709     /*
 710      * AVX2 in poly1305-x86.pl failed this with 176+32 split
 711      */
 712     {
 713     "248ac31085b6c2adaaa38259a0d7192c5c35d1bb4ef39ad94c38d1c82479e2dd"
 714     "2159a077024b0589bc8a20101b506f0a1ad0bbab76e83a83f1b94be6beae74e8"
 715     "74cab692c5963a75436b776121ec9f62399a3e66b2d22707dae81933b6277f3c"
 716     "8516bcbe26dbbd86f373103d7cf4cad1888c952118fbfbd0d7b4bedc4ae4936a"
 717     "ff91157e7aa47c54442ea78d6ac251d324a0fbe49d89cc3521b66d16e9c66a37"
 718     "09894e4eb0a4eedc4ae19468e66b81f2"
 719     "71351b1d921ea551047abcc6b87a901fde7db79fa1818c11336dbc07244a40eb",
 720     "000102030405060708090a0b0c0d0e0f""00000000000000000000000000000000",
 721     "bc939bc5281480fa99c6d68c258ec42f"
 722     },
 723     /*
 724      * test vectors from Google
 725      */
 726     {
 727      "",
 728      "c8afaac331ee372cd6082de134943b17""4710130e9f6fea8d72293850a667d86c",
 729      "4710130e9f6fea8d72293850a667d86c",
 730     },
 731     {
 732      "48656c6c6f20776f726c6421",
 733      "746869732069732033322d6279746520""6b657920666f7220506f6c7931333035",
 734      "a6f745008f81c916a20dcc74eef2b2f0"
 735     },
 736     {
 737      "0000000000000000000000000000000000000000000000000000000000000000",
 738      "746869732069732033322d6279746520""6b657920666f7220506f6c7931333035",
 739      "49ec78090e481ec6c26b33b91ccc0307"
 740     },
 741     {
 742      "89dab80b7717c1db5db437860a3f70218e93e1b8f461fb677f16f35f6f87e2a9"
 743      "1c99bc3a47ace47640cc95c345be5ecca5a3523c35cc01893af0b64a62033427"
 744      "0372ec12482d1b1e363561698a578b359803495bb4e2ef1930b17a5190b580f1"
 745      "41300df30adbeca28f6427a8bc1a999fd51c554a017d095d8c3e3127daf9f595",
 746      "2d773be37adb1e4d683bf0075e79c4ee""037918535a7f99ccb7040fb5f5f43aea",
 747      "c85d15ed44c378d6b00e23064c7bcd51"
 748     },
 749     {
 750      "000000000000000b1703030200000000"
 751      "06db1f1f368d696a810a349c0c714c9a5e7850c2407d721acded95e018d7a852"
 752      "66a6e1289cdb4aeb18da5ac8a2b0026d24a59ad485227f3eaedbb2e7e35e1c66"
 753      "cd60f9abf716dcc9ac42682dd7dab287a7024c4eefc321cc0574e16793e37cec"
 754      "03c5bda42b54c114a80b57af26416c7be742005e20855c73e21dc8e2edc9d435"
 755      "cb6f6059280011c270b71570051c1c9b3052126620bc1e2730fa066c7a509d53"
 756      "c60e5ae1b40aa6e39e49669228c90eecb4a50db32a50bc49e90b4f4b359a1dfd"
 757      "11749cd3867fcf2fb7bb6cd4738f6a4ad6f7ca5058f7618845af9f020f6c3b96"
 758      "7b8f4cd4a91e2813b507ae66f2d35c18284f7292186062e10fd5510d18775351"
 759      "ef334e7634ab4743f5b68f49adcab384d3fd75f7390f4006ef2a295c8c7a076a"
 760      "d54546cd25d2107fbe1436c840924aaebe5b370893cd63d1325b8616fc481088"
 761      "6bc152c53221b6df373119393255ee72bcaa880174f1717f9184fa91646f17a2"
 762      "4ac55d16bfddca9581a92eda479201f0edbf633600d6066d1ab36d5d2415d713"
 763      "51bbcd608a25108d25641992c1f26c531cf9f90203bc4cc19f5927d834b0a471"
 764      "16d3884bbb164b8ec883d1ac832e56b3918a98601a08d171881541d594db399c"
 765      "6ae6151221745aec814c45b0b05b565436fd6f137aa10a0c0b643761dbd6f9a9"
 766      "dcb99b1a6e690854ce0769cde39761d82fcdec15f0d92d7d8e94ade8eb83fbe0",
 767      "99e5822dd4173c995e3dae0ddefb9774""3fde3b080134b39f76e9bf8d0e88d546",
 768      "2637408fe13086ea73f971e3425e2820"
 769     },
 770     /*
 771      * test vectors from Hanno Böck
 772      */
 773     {
 774      "cccccccccccccccccccccccccccccccccccccccccccccccccc80cccccccccccc"
 775      "cccccccccccccccccccccccccccccccccccccccccccccccccccccccccecccccc"
 776      "ccccccccccccccccccccccccccccccc5cccccccccccccccccccccccccccccccc"
 777      "cccccccccce3cccccccccccccccccccccccccccccccccccccccccccccccccccc"
 778      "ccccccccaccccccccccccccccccccce6cccccccccc000000afcccccccccccccc"
 779      "ccccfffffff50000000000000000000000000000000000000000000000000000"
 780      "00ffffffe7000000000000000000000000000000000000000000000000000000"
 781      "0000000000000000000000000000000000000000000000000000719205a8521d"
 782      "fc",
 783      "7f1b0264000000000000000000000000""0000000000000000cccccccccccccccc",
 784      "8559b876eceed66eb37798c0457baff9"
 785     },
 786     {
 787      "aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa0000000000"
 788      "00000000800264",
 789      "e0001600000000000000000000000000""0000aaaaaaaaaaaaaaaaaaaaaaaaaaaa",
 790      "00bd1258978e205444c9aaaa82006fed"
 791     },
 792     {
 793      "02fc",
 794      "0c0c0c0c0c0c0c0c0c0c0c0c0c0c0c0c""0c0c0c0c0c0c0c0c0c0c0c0c0c0c0c0c",
 795      "06120c0c0c0c0c0c0c0c0c0c0c0c0c0c"
 796     },
 797     {
 798      "7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b"
 799      "7b7b7b7b7b7b7a7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b"
 800      "7b7b5c7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b"
 801      "7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b6e7b007b7b7b7b7b7b7b7b7b"
 802      "7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7a7b7b7b7b7b7b7b7b7b7b7b7b"
 803      "7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b5c7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b"
 804      "7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b7b"
 805      "7b6e7b001300000000b300000000000000000000000000000000000000000000"
 806      "f20000000000000000000000000000000000002000efff000900000000000000"
 807      "0000000000100000000009000000640000000000000000000000001300000000"
 808      "b300000000000000000000000000000000000000000000f20000000000000000"
 809      "000000000000000000002000efff00090000000000000000007a000010000000"
 810      "000900000064000000000000000000000000000000000000000000000000fc",
 811      "00ff0000000000000000000000000000""00000000001e00000000000000007b7b",
 812      "33205bbf9e9f8f7212ab9e2ab9b7e4a5"
 813     },
 814     {
 815      "7777777777777777777777777777777777777777777777777777777777777777"
 816      "7777777777777777777777777777777777777777777777777777777777777777"
 817      "777777777777777777777777ffffffe9e9acacacacacacacacacacac0000acac"
 818      "ec0100acacac2caca2acacacacacacacacacacac64f2",
 819      "0000007f0000007f0100002000000000""0000cf77777777777777777777777777",
 820      "02ee7c8c546ddeb1a467e4c3981158b9"
 821     },
 822     /*
 823      * test vectors from Andrew Moon
 824      */
 825     {   /* nacl */
 826      "8e993b9f48681273c29650ba32fc76ce48332ea7164d96a4476fb8c531a1186a"
 827      "c0dfc17c98dce87b4da7f011ec48c97271d2c20f9b928fe2270d6fb863d51738"
 828      "b48eeee314a7cc8ab932164548e526ae90224368517acfeabd6bb3732bc0e9da"
 829      "99832b61ca01b6de56244a9e88d5f9b37973f622a43d14a6599b1f654cb45a74"
 830      "e355a5",
 831      "eea6a7251c1e72916d11c2cb214d3c25""2539121d8e234e652d651fa4c8cff880",
 832      "f3ffc7703f9400e52a7dfb4b3d3305d9"
 833     },
 834     {   /* wrap 2^130-5 */
 835      "ffffffffffffffffffffffffffffffff",
 836      "02000000000000000000000000000000""00000000000000000000000000000000",
 837      "03000000000000000000000000000000"
 838     },
 839     {   /* wrap 2^128 */
 840      "02000000000000000000000000000000",
 841      "02000000000000000000000000000000""ffffffffffffffffffffffffffffffff",
 842      "03000000000000000000000000000000"
 843     },
 844     {   /* limb carry */
 845      "fffffffffffffffffffffffffffffffff0ffffffffffffffffffffffffffffff"
 846      "11000000000000000000000000000000",
 847      "01000000000000000000000000000000""00000000000000000000000000000000",
 848      "05000000000000000000000000000000"
 849     },
 850     {   /* 2^130-5 */
 851      "fffffffffffffffffffffffffffffffffbfefefefefefefefefefefefefefefe"
 852      "01010101010101010101010101010101",
 853      "01000000000000000000000000000000""00000000000000000000000000000000",
 854      "00000000000000000000000000000000"
 855     },
 856     {   /* 2^130-6 */
 857      "fdffffffffffffffffffffffffffffff",
 858      "02000000000000000000000000000000""00000000000000000000000000000000",
 859      "faffffffffffffffffffffffffffffff"
 860     },
 861     {   /* 5*H+L reduction intermediate */
 862      "e33594d7505e43b900000000000000003394d7505e4379cd0100000000000000"
 863      "0000000000000000000000000000000001000000000000000000000000000000",
 864      "01000000000000000400000000000000""00000000000000000000000000000000",
 865      "14000000000000005500000000000000"
 866     },
 867     {   /* 5*H+L reduction final */
 868      "e33594d7505e43b900000000000000003394d7505e4379cd0100000000000000"
 869      "00000000000000000000000000000000",
 870      "01000000000000000400000000000000""00000000000000000000000000000000",
 871      "13000000000000000000000000000000"
 872     }
 873 };
 874
 875 static unsigned char hex_digit(char h)
 876 {
 877     int i = OPENSSL_hexchar2int(h);
 878
 879     if (i < 0)
 880         abort();
 881     return i;
 882 }
 883
 884 static void hex_decode(unsigned char *out, const char *hex)
 885 {
 886     size_t j = 0;
 887
 888     while (*hex != 0) {
 889         unsigned char v = hex_digit(*hex++);
 890         v <<= 4;
 891         v |= hex_digit(*hex++);
 892         out[j++] = v;
 893     }
 894 }
 895
 896 static void hexdump(unsigned char *a, size_t len)
 897 {
 898     size_t i;
 899
 900     for (i = 0; i < len; i++)
 901         printf("%02x", a[i]);
 902 }
 903
 904 int main()
 905 {
 906     static const unsigned num_tests =
 907         sizeof(poly1305_tests) / sizeof(struct poly1305_test);
 908     unsigned i;
 909     unsigned char key[32], out[16], expected[16];
 910     POLY1305 poly1305;
 911
 912     for (i = 0; i < num_tests; i++) {
 913         const struct poly1305_test *test = &poly1305_tests[i];
 914         unsigned char *in;
 915         size_t inlen = strlen(test->inputhex);
 916
 917         if (strlen(test->keyhex) != sizeof(key) * 2 ||
 918             strlen(test->outhex) != sizeof(out) * 2 || (inlen & 1) == 1)
 919             return 1;
 920
 921         inlen /= 2;
 922
 923         hex_decode(key, test->keyhex);
 924         hex_decode(expected, test->outhex);
 925
 926         in = malloc(inlen);
 927
 928         hex_decode(in, test->inputhex);
 929
 930         Poly1305_Init(&poly1305, key);
 931         Poly1305_Update(&poly1305, in, inlen);
 932         Poly1305_Final(&poly1305, out);
 933
 934         if (memcmp(out, expected, sizeof(expected)) != 0) {
 935             printf("Poly1305 test #%d failed.\n", i);
 936             printf("got:      ");
 937             hexdump(out, sizeof(out));
 938             printf("\nexpected: ");
 939             hexdump(expected, sizeof(expected));
 940             printf("\n");
 941             return 1;
 942         }
 943
 944         if (inlen > 16) {
 945             Poly1305_Init(&poly1305, key);
 946             Poly1305_Update(&poly1305, in, 1);
 947             Poly1305_Update(&poly1305, in+1, inlen-1);
 948             Poly1305_Final(&poly1305, out);
 949
 950             if (memcmp(out, expected, sizeof(expected)) != 0) {
 951                 printf("Poly1305 test #%d/1+(N-1) failed.\n", i);
 952                 printf("got:      ");
 953                 hexdump(out, sizeof(out));
 954                 printf("\nexpected: ");
 955                 hexdump(expected, sizeof(expected));
 956                 printf("\n");
 957                 return 1;
 958             }
 959         }
 960
 961         if (inlen > 32) {
 962             size_t half = inlen / 2;
 963
 964             Poly1305_Init(&poly1305, key);
 965             Poly1305_Update(&poly1305, in, half);
 966             Poly1305_Update(&poly1305, in+half, inlen-half);
 967             Poly1305_Final(&poly1305, out);
 968
 969             if (memcmp(out, expected, sizeof(expected)) != 0) {
 970                 printf("Poly1305 test #%d/2 failed.\n", i);
 971                 printf("got:      ");
 972                 hexdump(out, sizeof(out));
 973                 printf("\nexpected: ");
 974                 hexdump(expected, sizeof(expected));
 975                 printf("\n");
 976                 return 1;
 977             }
 978
 979             for (half = 16; half < inlen; half += 16) {
 980                 Poly1305_Init(&poly1305, key);
 981                 Poly1305_Update(&poly1305, in, half);
 982                 Poly1305_Update(&poly1305, in+half, inlen-half);
 983                 Poly1305_Final(&poly1305, out);
 984
 985                 if (memcmp(out, expected, sizeof(expected)) != 0) {
 986                     printf("Poly1305 test #%d/%d+%d failed.\n",
 987                                            i, half, inlen-half);
 988                     printf("got:      ");
 989                     hexdump(out, sizeof(out));
 990                     printf("\nexpected: ");
 991                     hexdump(expected, sizeof(expected));
 992                     printf("\n");
 993                     return 1;
 994                 }
 995             }
 996         }
 997
 998         free(in);
 999     }
1000
1001     printf("PASS\n");
1002
1003 # ifdef OPENSSL_CPUID_OBJ
1004     {
1005         unsigned char buf[8192];
1006         unsigned long long stopwatch;
1007         unsigned long long OPENSSL_rdtsc();
1008
1009         memset (buf,0x55,sizeof(buf));
1010         memset (key,0xAA,sizeof(key));
1011
1012         Poly1305_Init(&poly1305, key);
1013
1014         for (i=0;i<100000;i++)
1015             Poly1305_Update(&poly1305,buf,sizeof(buf));
1016
1017         stopwatch = OPENSSL_rdtsc();
1018         for (i=0;i<10000;i++)
1019             Poly1305_Update(&poly1305,buf,sizeof(buf));
1020         stopwatch = OPENSSL_rdtsc() - stopwatch;
1021
1022         printf("%g\n",stopwatch/(double)(i*sizeof(buf)));
1023
1024         stopwatch = OPENSSL_rdtsc();
1025         for (i=0;i<10000;i++) {
1026             Poly1305_Init(&poly1305, key);
1027             Poly1305_Update(&poly1305,buf,16);
1028             Poly1305_Final(&poly1305,buf);
1029         }
1030         stopwatch = OPENSSL_rdtsc() - stopwatch;
1031
1032         printf("%g\n",stopwatch/(double)(i));
1033     }
1034 # endif
1035     return 0;
1036 }
1037 #endif