progs/perftest.c: Use from Glibc syscall numbers.

[catacomb] / pub / x448.c

/* -*-c-*-
 *
 * The X448 key-agreement algorithm
 *
 * (c) 2017 Straylight/Edgeware
 */

/*----- Licensing notice --------------------------------------------------*
 *
 * This file is part of Catacomb.
 *
 * Catacomb is free software; you can redistribute it and/or modify
 * it under the terms of the GNU Library General Public License as
 * published by the Free Software Foundation; either version 2 of the
 * License, or (at your option) any later version.
 *
 * Catacomb 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 Library General Public License for more details.
 *
 * You should have received a copy of the GNU Library General Public
 * License along with Catacomb; if not, write to the Free
 * Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
 * MA 02111-1307, USA.
 */

/*----- Header files ------------------------------------------------------*/

#include <mLib/bits.h>

#include "montladder.h"
#include "fgoldi.h"
#include "x448.h"

/*----- Important constants -----------------------------------------------*/

const octet x448_base[56] = { 5, 0, /* ... */ };

#define A0 39081

/*----- Key fetching ------------------------------------------------------*/

const key_fetchdef x448_pubfetch[] = {
  { "pub",	offsetof(x448_pub, pub),	KENC_BINARY,	0 },
  { 0,		0,				0,		0 }
};

static const key_fetchdef priv[] = {
  { "priv",	offsetof(x448_priv, priv),	KENC_BINARY,	0 },
  { 0,		0,				0,		0 }
};

const key_fetchdef x448_privfetch[] = {
  { "pub",	offsetof(x448_priv, pub),	KENC_BINARY,	0 },
  { "private",	0,				KENC_STRUCT,	priv },
  { 0,		0,				0,		0 }
};

/*----- Main code ---------------------------------------------------------*/

/* --- @x448@ --- *
 *
 * Arguments:	@octet zz[X448_OUTSZ]@ = where to put the result
 *		@const octet k[X448_KEYSZ]@ = pointer to private key
 *		@const octet qx[X448_PUBSZ]@ = pointer to public value
 *
 * Returns:	---
 *
 * Use:		Calculates X448 of @k@ and @qx@.
 */

void x448(octet zz[X448_OUTSZ],
	  const octet k[X448_KEYSZ],
	  const octet qx[X448_PUBSZ])
{
  uint32 kw[14];
  fgoldi x1;
  unsigned i;

  /* Load and clamp the key.  The low bits are cleared to kill the small
   * subgroups on the curve and its twist, and a high bit is set to guard
   * against careless implementations, though this isn't one of those.
   */
  for (i = 0; i < 14; i++) kw[i] = LOAD32_L(k + 4*i);
  kw[0] &= 0xfffffffc; kw[13] |= 0x80000000;

  /* And run the ladder. */
  fgoldi_load(&x1, qx);
#define MULA0(z, x) do { fgoldi_mulconst((z), (x), A0); } while (0)
  MONT_LADDER(fgoldi, MULA0, kw, 14, 32, &x1, &x1);
#undef MULA0
  fgoldi_store(zz, &x1);
}

/*----- Test rig ----------------------------------------------------------*/

#ifdef TEST_RIG

#include <mLib/macros.h>
#include <mLib/report.h>
#include <mLib/testrig.h>

#include "ct.h"

static int vrf_x448(dstr dv[])
{
  dstr dz = DSTR_INIT;
  int ok = 1;

  if (dv[0].len != X448_KEYSZ) die(1, "bad key length");
  if (dv[1].len != X448_PUBSZ) die(1, "bad public length");
  if (dv[2].len != X448_OUTSZ) die(1, "bad result length");

  ct_poison(dv[0].buf, dv[0].len);
  dstr_ensure(&dz, X448_OUTSZ); dz.len = X448_OUTSZ;
  x448((octet *)dz.buf,
       (const octet *)dv[0].buf,
       (const octet *)dv[1].buf);
  ct_remedy(dz.buf, dz.len);
  if (MEMCMP(dz.buf, !=, dv[2].buf, X448_OUTSZ)) {
    ok = 0;
    fprintf(stderr, "failed!");
    fprintf(stderr, "\n\t   k = "); type_hex.dump(&dv[0], stderr);
    fprintf(stderr, "\n\t   p = "); type_hex.dump(&dv[1], stderr);
    fprintf(stderr, "\n\twant = "); type_hex.dump(&dv[2], stderr);
    fprintf(stderr, "\n\tcalc = "); type_hex.dump(&dz, stderr);
    fprintf(stderr, "\n");
  }

  dstr_destroy(&dz);
  return (ok);
}

static int vrf_mct(dstr dv[])
{
  octet b0[X448_OUTSZ], b1[X448_OUTSZ], *k = b0, *x = b1, *t;
  unsigned long i, niter;
  dstr d = DSTR_INIT;
  int ok = 1;

  if (dv[0].len != sizeof(b0)) { fprintf(stderr, "k len\n"); exit(2); }
  if (dv[1].len != sizeof(b1)) { fprintf(stderr, "x len\n"); exit(2); }
  if (dv[3].len != sizeof(b0)) { fprintf(stderr, "result len\n"); exit(2); }
  memcpy(b0, dv[0].buf, sizeof(b0));
  memcpy(b1, dv[1].buf, sizeof(b1));
  niter = *(unsigned long *)dv[2].buf;
  dstr_ensure(&d, X448_OUTSZ); d.len = X448_OUTSZ; t = (octet *)d.buf;

  for (i = 0; i < niter; i++) {
    x448(x, k, x);
    t = x; x = k; k = t;
  }
  memcpy(d.buf, k, d.len);

  if (MEMCMP(d.buf, !=, dv[3].buf, d.len)) {
    ok = 0;
    fprintf(stderr, "failed...");
    fprintf(stderr, "\n\tinitial k = "); type_hex.dump(&dv[0], stderr);
    fprintf(stderr, "\n\tinitial x = "); type_hex.dump(&dv[1], stderr);
    fprintf(stderr, "\n\titerations = %lu", niter);
    fprintf(stderr, "\n\texpected = "); type_hex.dump(&dv[3], stderr);
    fprintf(stderr, "\n\tcalculated = "); type_hex.dump(&d, stderr);
    fputc('\n', stderr);
  }

  dstr_destroy(&d);
  return (ok);
}

static test_chunk tests[] = {
  { "x448", vrf_x448, { &type_hex, &type_hex, &type_hex } },
  { "x448-mct", vrf_mct,
    { &type_hex, &type_hex, &type_ulong, &type_hex } },
  { 0, 0, { 0 } }
};

int main(int argc, char *argv[])
{
  test_run(argc, argv, tests, SRCDIR "/t/x448");
  return (0);
}

#endif

/*----- That's all, folks -------------------------------------------------*/