nettle-benchmark.c 11 KB
Newer Older
1 2 3 4 5 6 7 8
/* nettle-benchmark.c
 *
 * Tries the performance of the various algorithms.
 *
 */
 
/* nettle, low-level cryptographics library
 *
Niels Möller's avatar
Niels Möller committed
9
 * Copyright (C) 2001, 2010 Niels Mller
10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
 *  
 * The nettle library is free software; you can redistribute it and/or modify
 * it under the terms of the GNU Lesser General Public License as published by
 * the Free Software Foundation; either version 2.1 of the License, or (at your
 * option) any later version.
 * 
 * The nettle library is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
 * or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
 * License for more details.
 * 
 * You should have received a copy of the GNU Lesser General Public License
 * along with the nettle library; see the file COPYING.LIB.  If not, write to
 * the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
 * MA 02111-1307, USA.
 */

27 28
#if HAVE_CONFIG_H
# include "config.h"
29 30 31 32
#endif

#include <assert.h>
#include <errno.h>
33
#include <math.h>
34 35 36 37 38
#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#include <time.h>
39

40
#include "aes.h"
Niels Möller's avatar
Niels Möller committed
41 42 43
#include "arcfour.h"
#include "blowfish.h"
#include "cast128.h"
Niels Möller's avatar
Niels Möller committed
44
#include "cbc.h"
Niels Möller's avatar
Niels Möller committed
45
#include "des.h"
46
#include "memxor.h"
Niels Möller's avatar
Niels Möller committed
47
#include "serpent.h"
48
#include "sha.h"
Niels Möller's avatar
Niels Möller committed
49 50
#include "twofish.h"

51 52 53
#include "nettle-meta.h"
#include "nettle-internal.h"

Niels Möller's avatar
Niels Möller committed
54
#include "getopt.h"
55

56
static double frequency = 0.0;
Niels Möller's avatar
Niels Möller committed
57

58
/* Process BENCH_BLOCK bytes at a time, for BENCH_INTERVAL seconds. */
59
#define BENCH_BLOCK 10240
60
#define BENCH_INTERVAL 0.1
Niels Möller's avatar
Niels Möller committed
61

62 63
/* FIXME: Proper configure test for rdtsc? */
#ifndef WITH_CYCLE_COUNTER
64
# if defined(__GNUC__) && (defined(__i386__) || defined(__x86_64__))
65 66 67 68 69 70 71
#  define WITH_CYCLE_COUNTER 1
# else
#  define WITH_CYCLE_COUNTER 0
# endif
#endif

#if WITH_CYCLE_COUNTER
72
# if defined(__i386__)
73 74 75 76 77 78 79 80 81
#define GET_CYCLE_COUNTER(hi, lo)		\
  __asm__("xorl %%eax,%%eax\n"			\
	  "movl %%ebx, %%edi\n"			\
	  "cpuid\n"				\
	  "rdtsc\n"				\
	  "movl %%edi, %%ebx\n"			\
	  : "=a" (lo), "=d" (hi)		\
	  : /* No inputs. */			\
	  : "%edi", "%ecx", "cc")
82 83 84 85 86 87 88 89 90 91 92
# elif defined(__x86_64__)
#define GET_CYCLE_COUNTER(hi, lo)		\
  __asm__("xorl %%eax,%%eax\n"			\
	  "mov %%rbx, %%r10\n"			\
	  "cpuid\n"				\
	  "rdtsc\n"				\
	  "mov %%r10, %%rbx\n"			\
	  : "=a" (lo), "=d" (hi)		\
	  : /* No inputs. */			\
	  : "%r10", "%rcx", "cc")
# endif
93 94 95
#define BENCH_ITERATIONS 10
#endif

96
/* Returns second per function call */
97 98 99
static double
time_function(void (*f)(void *arg), void *arg)
{
100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127
  unsigned ncalls;
#if HAVE_CLOCK_GETTIME && defined CLOCK_PROCESS_CPUTIME_ID
  struct timespec before;
  struct timespec after;
  struct timespec done;

  clock_gettime(CLOCK_PROCESS_CPUTIME_ID, &before);
  done = before;
  done.tv_nsec += BENCH_INTERVAL * 1e9;
  if (done.tv_nsec >= 1000000000L)
    {
      done.tv_nsec -= 1000000000L;
      done.tv_sec ++;
    }
  ncalls = 0;

  do 
    {
      f(arg);
      clock_gettime(CLOCK_PROCESS_CPUTIME_ID, &after);
      ncalls++;
    }
  while (after.tv_sec < done.tv_sec
	 || (after.tv_nsec < done.tv_nsec && after.tv_sec == done.tv_sec));
  
  return (after.tv_sec - before.tv_sec
	  + 1e-9 * (after.tv_nsec - before.tv_nsec)) / ncalls;
#else /* !HAVE_CLOCK_GETTIME */
128 129
  clock_t before;
  clock_t after;
130 131
  clock_t done;
  
132
  before = clock();
133
  done = before + BENCH_INTERVAL * CLOCKS_PER_SEC;
134
  ncalls = 0;
135
  
136 137 138 139 140 141 142
  do 
    {
      f(arg);
      after = clock();
      ncalls++;
    }
  while (after < done);
143
  
144
  return ((double)(after - before)) / CLOCKS_PER_SEC / ncalls;
145
#endif /* !HAVE_CLOCK_GETTIME */
146 147
}

148 149 150 151 152 153 154 155 156 157 158 159 160
struct bench_memxor_info
{
  uint8_t *dst;
  const uint8_t *src;
};

static void
bench_memxor(void *arg)
{
  struct bench_memxor_info *info = arg;
  memxor (info->dst, info->src, BENCH_BLOCK);
}

161 162 163
struct bench_hash_info
{
  void *ctx;
164
  nettle_hash_update_func *update;
165 166 167 168 169 170 171 172 173 174
  const uint8_t *data;
};

static void
bench_hash(void *arg)
{
  struct bench_hash_info *info = arg;
  info->update(info->ctx, BENCH_BLOCK, info->data);
}

175 176 177
struct bench_cipher_info
{
  void *ctx;
178
  nettle_crypt_func *crypt;
179 180 181 182 183 184 185
  uint8_t *data;
};

static void
bench_cipher(void *arg)
{
  struct bench_cipher_info *info = arg;
186
  info->crypt(info->ctx, BENCH_BLOCK, info->data, info->data);
187 188 189 190 191
}

struct bench_cbc_info
{
  void *ctx;
192
  nettle_crypt_func *crypt;
193
 
194
  uint8_t *data;
195
  
196 197 198 199 200 201 202 203
  unsigned block_size;
  uint8_t *iv;
};

static void
bench_cbc_encrypt(void *arg)
{
  struct bench_cbc_info *info = arg;
204 205 206
  cbc_encrypt(info->ctx, info->crypt,
	      info->block_size, info->iv,
	      BENCH_BLOCK, info->data, info->data);
207 208 209 210 211 212
}

static void
bench_cbc_decrypt(void *arg)
{
  struct bench_cbc_info *info = arg;
213 214 215
  cbc_decrypt(info->ctx, info->crypt,
	      info->block_size, info->iv,
	      BENCH_BLOCK, info->data, info->data);
216 217 218 219 220 221 222
}

/* Set data[i] = floor(sqrt(i)) */
static void
init_data(uint8_t *data)
{
  unsigned i,j;
Niels Möller's avatar
Niels Möller committed
223
  for (i = j = 0; i<BENCH_BLOCK;  i++)
224 225 226 227 228 229 230 231
    {
      if (j*j < i)
	j++;
      data[i] = j;
    }
}

static void
Niels Möller's avatar
Niels Möller committed
232 233
init_key(unsigned length,
         uint8_t *key)
234
{
Niels Möller's avatar
Niels Möller committed
235 236 237
  unsigned i;
  for (i = 0; i<length; i++)
    key[i] = i;
238 239
}

240 241 242 243 244
static void
header(void)
{
  printf("%18s %11s Mbyte/s%s\n",
	 "Algorithm", "mode", 
Niels Möller's avatar
Niels Möller committed
245
	 frequency > 0.0 ? " cycles/byte cycles/block" : "");  
246 247
}

Niels Möller's avatar
Niels Möller committed
248
static void
Niels Möller's avatar
Niels Möller committed
249
display(const char *name, const char *mode, unsigned block_size,
250
	double time)
Niels Möller's avatar
Niels Möller committed
251
{
252
  printf("%18s %11s %7.2f",
Niels Möller's avatar
Niels Möller committed
253
	 name, mode,
254
	 BENCH_BLOCK / (time * 1048576.0));
255
  if (frequency > 0.0)
Niels Möller's avatar
Niels Möller committed
256
    {
257
      printf(" %11.2f", time * frequency / BENCH_BLOCK);
Niels Möller's avatar
Niels Möller committed
258
      if (block_size > 0)
259
	printf(" %12.2f", time * frequency * block_size / BENCH_BLOCK);
Niels Möller's avatar
Niels Möller committed
260
    }
261
  printf("\n");
Niels Möller's avatar
Niels Möller committed
262 263
}

264 265 266 267 268 269 270 271 272 273 274 275 276
static void *
xalloc(size_t size)
{
  void *p = malloc(size);
  if (!p)
    {
      fprintf(stderr, "Virtual memory exhausted.\n");
      abort();
    }

  return p;
}

277 278 279 280 281 282 283 284 285 286
static void
time_memxor(void)
{
  struct bench_memxor_info info;
  uint8_t src[BENCH_BLOCK + 1];
  uint8_t dst[BENCH_BLOCK + 1];

  info.src = src;
  info.dst = dst;

287 288
  display ("memxor", "aligned", sizeof(unsigned long),
	   time_function(bench_memxor, &info));
289
  info.src++;
290 291
  display ("memxor", "unaligned", sizeof(unsigned long),
	   time_function(bench_memxor, &info));  
292 293
}

294 295 296 297 298
static void
time_hash(const struct nettle_hash *hash)
{
  static uint8_t data[BENCH_BLOCK];
  struct bench_hash_info info;
299

300
  info.ctx = xalloc(hash->context_size); 
301 302 303 304 305 306
  info.update = hash->update;
  info.data = data;

  init_data(data);
  hash->init(info.ctx);

Niels Möller's avatar
Niels Möller committed
307
  display(hash->name, "update", hash->block_size,
308
	  time_function(bench_hash, &info));
309 310

  free(info.ctx);
311 312
}

Niels Möller's avatar
Niels Möller committed
313
static void
314
time_cipher(const struct nettle_cipher *cipher)
Niels Möller's avatar
Niels Möller committed
315
{
316 317
  void *ctx = xalloc(cipher->context_size);
  uint8_t *key = xalloc(cipher->key_size);
Niels Möller's avatar
Niels Möller committed
318

319
  static uint8_t data[BENCH_BLOCK];
Niels Möller's avatar
Niels Möller committed
320 321 322 323

  printf("\n");
  
  init_data(data);
324 325

  {
Niels Möller's avatar
Niels Möller committed
326 327 328 329 330
    /* Decent initializers are a GNU extension, so don't use it here. */
    struct bench_cipher_info info;
    info.ctx = ctx;
    info.crypt = cipher->encrypt;
    info.data = data;
331
    
Niels Möller's avatar
Niels Möller committed
332
    init_key(cipher->key_size, key);
333
    cipher->set_encrypt_key(ctx, cipher->key_size, key);
Niels Möller's avatar
Niels Möller committed
334

Niels Möller's avatar
Niels Möller committed
335
    display(cipher->name, "ECB encrypt", cipher->block_size,
Niels Möller's avatar
Niels Möller committed
336
	    time_function(bench_cipher, &info));
337
  }
Niels Möller's avatar
Niels Möller committed
338
  
339
  {
Niels Möller's avatar
Niels Möller committed
340 341 342 343
    struct bench_cipher_info info;
    info.ctx = ctx;
    info.crypt = cipher->decrypt;
    info.data = data;
344
    
Niels Möller's avatar
Niels Möller committed
345
    init_key(cipher->key_size, key);
346
    cipher->set_decrypt_key(ctx, cipher->key_size, key);
Niels Möller's avatar
Niels Möller committed
347

Niels Möller's avatar
Niels Möller committed
348
    display(cipher->name, "ECB decrypt", cipher->block_size,
Niels Möller's avatar
Niels Möller committed
349
	    time_function(bench_cipher, &info));
350 351
  }

Niels Möller's avatar
Niels Möller committed
352 353
  /* Don't use nettle cbc to benchmark openssl ciphers */
  if (cipher->block_size && cipher->name[0] != 'o')
Niels Möller's avatar
Niels Möller committed
354
    {
355
      uint8_t *iv = xalloc(cipher->block_size);
Niels Möller's avatar
Niels Möller committed
356 357 358
      
      /* Do CBC mode */
      {
Niels Möller's avatar
Niels Möller committed
359 360 361 362 363 364
        struct bench_cbc_info info;
	info.ctx = ctx;
	info.crypt = cipher->encrypt;
	info.data = data;
	info.block_size = cipher->block_size;
	info.iv = iv;
365
    
Niels Möller's avatar
Niels Möller committed
366
        memset(iv, 0, sizeof(iv));
367
    
368
        cipher->set_encrypt_key(ctx, cipher->key_size, key);
369

Niels Möller's avatar
Niels Möller committed
370
	display(cipher->name, "CBC encrypt", cipher->block_size,
Niels Möller's avatar
Niels Möller committed
371
		time_function(bench_cbc_encrypt, &info));
Niels Möller's avatar
Niels Möller committed
372
      }
373

Niels Möller's avatar
Niels Möller committed
374
      {
Niels Möller's avatar
Niels Möller committed
375 376 377 378 379 380
        struct bench_cbc_info info;
	info.ctx = ctx;
	info.crypt = cipher->decrypt;
	info.data = data;
	info.block_size = cipher->block_size;
	info.iv = iv;
381
    
Niels Möller's avatar
Niels Möller committed
382
        memset(iv, 0, sizeof(iv));
383

384
        cipher->set_decrypt_key(ctx, cipher->key_size, key);
385

Niels Möller's avatar
Niels Möller committed
386
	display(cipher->name, "CBC decrypt", cipher->block_size,
Niels Möller's avatar
Niels Möller committed
387
		time_function(bench_cbc_decrypt, &info));
Niels Möller's avatar
Niels Möller committed
388
      }
389
      free(iv);
Niels Möller's avatar
Niels Möller committed
390
    }
391 392
  free(ctx);
  free(key);
Niels Möller's avatar
Niels Möller committed
393 394
}

395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442
static int
compare_double(const void *ap, const void *bp)
{
  double a = *(const double *) ap;
  double b = *(const double *) bp;
  if (a < b)
    return -1;
  else if (a > b)
    return 1;
  else
    return 0;
}

/* Try to get accurate cycle times for assembler functions. */
static void
bench_sha1_compress(void)
{
#if WITH_CYCLE_COUNTER
  uint32_t state[_SHA1_DIGEST_LENGTH];
  uint8_t data[BENCH_ITERATIONS * SHA1_DATA_SIZE];
  uint32_t start_lo, start_hi, end_lo, end_hi;

  double count[5];
  
  uint8_t *p;
  unsigned i, j;

  for (j = 0; j < 5; j++)
    {
      i = 0;
      p = data;
      GET_CYCLE_COUNTER(start_hi, start_lo);
      for (; i < BENCH_ITERATIONS; i++, p += SHA1_DATA_SIZE)
	_nettle_sha1_compress(state, p);

      GET_CYCLE_COUNTER(end_hi, end_lo);

      end_hi -= (start_hi + (start_lo > end_lo));
      end_lo -= start_lo;

      count[j] = ldexp(end_hi, 32) + end_lo;
    }

  qsort(count, 5, sizeof(double), compare_double);
  printf("sha1_compress: %.2f cycles\n\n", count[2] / BENCH_ITERATIONS);  
#endif
}

Niels Möller's avatar
Niels Möller committed
443
#if WITH_OPENSSL
444 445 446 447
# define OPENSSL(x) x,
#else
# define OPENSSL(x)
#endif
Niels Möller's avatar
Niels Möller committed
448 449

int
450
main(int argc, char **argv)
Niels Möller's avatar
Niels Möller committed
451 452
{
  unsigned i;
453
  int c;
454
  const char *alg;
455 456 457 458

  const struct nettle_hash *hashes[] =
    {
      &nettle_md2, &nettle_md4, &nettle_md5,
459
      OPENSSL(&nettle_openssl_md5)
460
      &nettle_sha1, OPENSSL(&nettle_openssl_sha1)
Niels Möller's avatar
Niels Möller committed
461 462
      &nettle_sha224, &nettle_sha256,
      &nettle_sha384, &nettle_sha512,
463 464 465
      NULL
    };

466
  const struct nettle_cipher *ciphers[] =
Niels Möller's avatar
Niels Möller committed
467
    {
468
      &nettle_aes128, &nettle_aes192, &nettle_aes256,
Niels Möller's avatar
Niels Möller committed
469 470 471 472 473
      OPENSSL(&nettle_openssl_aes128)
      OPENSSL(&nettle_openssl_aes192)
      OPENSSL(&nettle_openssl_aes256)
      &nettle_arcfour128, OPENSSL(&nettle_openssl_arcfour128)
      &nettle_blowfish128, OPENSSL(&nettle_openssl_blowfish128)
Niels Möller's avatar
Niels Möller committed
474
      &nettle_camellia128, &nettle_camellia192, &nettle_camellia256,
475 476 477
      &nettle_cast128, OPENSSL(&nettle_openssl_cast128)
      &nettle_des, OPENSSL(&nettle_openssl_des)
      &nettle_des3,
478 479
      &nettle_serpent256,
      &nettle_twofish128, &nettle_twofish192, &nettle_twofish256,
480
      NULL
Niels Möller's avatar
Niels Möller committed
481
    };
482

483 484 485 486 487 488 489 490 491
  while ( (c = getopt(argc, argv, "f:")) != -1)
    switch (c)
      {
      case 'f':
	frequency = atof(optarg);
	if (frequency > 0.0)
	  break;

      case ':': case '?':
492
	fprintf(stderr, "Usage: nettle-benchmark [-f clock frequency] [alg]\n");
493 494 495 496 497 498
	return EXIT_FAILURE;

      default:
	abort();
    }

499 500
  alg = argv[optind];
  
501 502
  bench_sha1_compress();

503 504
  header();

505 506 507 508 509
  if (!alg || strstr ("memxor", alg))
    {
      time_memxor();
      printf("\n");
    }
510
  
511 512 513 514 515
  for (i = 0; hashes[i]; i++)
    {
      if (!alg || strstr(hashes[i]->name, alg))
	time_hash(hashes[i]);
    }
516
  for (i = 0; ciphers[i]; i++)
517 518 519
    if (!alg || strstr(ciphers[i]->name, alg))
      time_cipher(ciphers[i]);

520 521
  return 0;
}