kpqc 0.1.0

data/vendor/NTRUplus/NTRU+1152/poly.c

@@ -0,0 +1,318 @@
+#include <stdint.h>
+#include <string.h>
+#include "params.h"
+#include "poly.h"
+#include "ntt.h"
+#include "symmetric.h"
+
+/*************************************************
+* Name:        crepmod3
+*
+* Description: Compute modulus 3 operation
+*
+* Arguments: - poly *a: pointer to intput integer to be reduced
+*
+* Returns:     integer in {-1,0,1} congruent to a modulo 3.
+**************************************************/
+static inline int16_t crepmod3(int16_t a)
+{
+	int16_t t;
+	const int16_t v = ((1<<15) + 3/2)/3;
+
+	a += (a >> 15) & NTRUPLUS_Q;
+	a -= (NTRUPLUS_Q+1)/2;
+	a += (a >> 15) & NTRUPLUS_Q;
+	a -= (NTRUPLUS_Q-1)/2;
+
+	t  = ((int32_t)v*a + (1<<14)) >> 15;
+	t *= 3;
+	return a - t;
+}
+
+/*************************************************
+* Name:        poly_tobytes
+*
+* Description: Serialization of a polynomial
+*
+* Arguments:   - uint8_t *r: pointer to output byte array
+*                            (needs space for NTRUPLUS_POLYBYTES bytes)
+*              - poly *a:    pointer to input polynomial
+**************************************************/
+void poly_tobytes(uint8_t r[NTRUPLUS_POLYBYTES], const poly *a)
+{
+	int16_t t[2];
+
+	#pragma GCC unroll 2
+	for(size_t i = 0; i < NTRUPLUS_N/2; i++)
+	{
+		t[0] = a->coeffs[2*i];
+		t[0] += (t[0] >> 15) & NTRUPLUS_Q;
+		t[1] = a->coeffs[2*i+1];
+		t[1] += (t[1] >> 15) & NTRUPLUS_Q;
+
+		r[3*i+0] = (t[0] >> 0);
+		r[3*i+1] = (t[0] >> 8) | (t[1] << 4);
+		r[3*i+2] = (t[1] >> 4);
+	}
+}
+
+
+/*************************************************
+* Name:        poly_frombytes
+*
+* Description: De-serialization of a polynomial;
+*              inverse of poly_tobytes
+*
+* Arguments:   - poly *r:          pointer to output polynomial
+*              - const uint8_t *a: pointer to input byte array
+*                                  (of NTRUPLUS_POLYBYTES bytes)
+**************************************************/
+void poly_frombytes(poly *r, const uint8_t a[NTRUPLUS_POLYBYTES])
+{
+	for(size_t i = 0; i < NTRUPLUS_N/2; i++)
+	{
+		r->coeffs[2*i]   = ((a[3*i+0] >> 0) | ((uint16_t)a[3*i+1] << 8)) & 0xFFF;
+		r->coeffs[2*i+1] = ((a[3*i+1] >> 4) | ((uint16_t)a[3*i+2] << 4)) & 0xFFF;
+	}
+}
+
+/*************************************************
+* Name:        poly_cbd1
+*
+* Description: Sample a polynomial deterministically from a random,
+*              with output polynomial close to centered binomial distribution
+*
+* Arguments:   - poly *r: pointer to output polynomial
+*              - const uint8_t *buf: pointer to input random
+*                                     (of length NTRUPLUS_N/4 bytes)
+**************************************************/
+void poly_cbd1(poly *r, const unsigned char buf[NTRUPLUS_N/4])
+{
+	uint8_t t1, t2;
+
+	for(size_t i = 0; i < NTRUPLUS_N / 8; i++)
+	{
+		t1 = buf[i];
+		t2 = buf[i + NTRUPLUS_N / 8];
+
+		for(size_t j = 0; j < 8; j++)
+		{
+			r->coeffs[8*i + j] = (t1 & 0x1) - (t2 & 0x1);
+
+			t1 >>= 1;   
+			t2 >>= 1;
+		}
+	}
+}
+
+/*************************************************
+* Name:        poly_sotp_encode
+*
+* Description: Encode a message deterministically using SOTP and a random,
+			   with output polynomial close to centered binomial distribution
+*
+* Arguments:   - poly *r: pointer to output polynomial
+*              - const uint8_t *msg: pointer to input message
+*              - const uint8_t *buf: pointer to input random
+**************************************************/
+void poly_sotp_encode(poly *r, const uint8_t msg[NTRUPLUS_N/8], const uint8_t buf[NTRUPLUS_N/4])
+{
+    uint8_t tmp[NTRUPLUS_N/4];
+
+    for(int i = 0; i < NTRUPLUS_N/8; i++)
+    {
+         tmp[i] = buf[i]^msg[i];
+    }
+
+    for(int i = NTRUPLUS_N/8; i < NTRUPLUS_N/4; i++)
+    {
+         tmp[i] = buf[i];
+    }
+
+	poly_cbd1(r, tmp);
+}
+
+/*************************************************
+* Name:        poly_sotp_decode
+*
+* Description: Decode a message deterministically using SOTP_INV and a random
+*
+* Arguments:   - uint8_t *msg: pointer to output message
+*              - const poly *a: pointer to iput polynomial
+*              - const uint8_t *buf: pointer to input random
+*
+* Returns 0 (success) or 1 (failure)
+**************************************************/
+int poly_sotp_decode(uint8_t msg[NTRUPLUS_N/8], const poly *a, const uint8_t buf[NTRUPLUS_N/4])
+{
+	uint8_t t1, t2, t3;
+	uint16_t t4;
+	uint32_t r = 0;
+	uint8_t mask;
+
+	for(size_t i = 0; i < NTRUPLUS_N / 8; i++)
+	{
+		t1 = buf[i                 ];
+		t2 = buf[i + NTRUPLUS_N / 8];
+		t3 = 0;
+
+		for(size_t j = 0; j < 8; j++)
+		{
+			t4 = t2 & 0x1;
+			t4 += a->coeffs[8*i + j];
+			r |= t4;
+			t4 = (t4 ^ t1) & 0x1;
+			t3 ^= (uint8_t)(t4 << j);
+
+			t1 >>= 1;
+			t2 >>= 1;
+		}
+
+		msg[i] = t3;
+	}
+
+	r = r >> 1;
+	r = (-(uint32_t)r) >> 31;
+
+	mask = (uint8_t)(r - 1);
+
+	for (size_t i = 0; i < NTRUPLUS_N / 8; i++)
+		msg[i] &= mask;
+
+	return r;
+}
+
+/*************************************************
+* Name:        poly_ntt
+*
+* Description: Computes number-theoretic transform (NTT)
+*
+* Arguments:   - poly *r: pointer to input/output polynomial
+**************************************************/
+void poly_ntt(poly *r)
+{
+	ntt(r->coeffs);
+}
+
+/*************************************************
+* Name:        poly_invntt
+*
+* Description: Computes inverse of number-theoretic transform (NTT)
+*
+* Arguments:   - poly *r: pointer to input/output polynomial
+**************************************************/
+void poly_invntt(poly *r)
+{
+	invntt(r->coeffs);
+}
+
+/*************************************************
+* Name:        poly_baseinv
+*
+* Description: Inverse of polynomial in NTT domain
+*
+* Arguments:   - poly *r:       pointer to output polynomial
+*              - const poly *a: pointer to input polynomial
+**************************************************/
+int poly_baseinv(poly *r, const poly *a)
+{
+	for(size_t i = 0; i < NTRUPLUS_N/8; ++i)
+	{
+		if(baseinv(r->coeffs + 8*i, a->coeffs + 8*i, zetas[144 + i]))
+		{
+			for (size_t j = 0; j < NTRUPLUS_N; ++j)
+				r->coeffs[j] = 0;
+
+			return 1;
+		} 
+
+		if(baseinv(r->coeffs + 8*i + 4, a->coeffs + 8*i + 4, -zetas[144 + i]))
+		{
+			for (size_t j = 0; j < NTRUPLUS_N; ++j)
+				r->coeffs[j] = 0;
+
+			return 1;
+		} 
+	 }
+
+	return 0;
+}
+
+/*************************************************
+* Name:        poly_basemul
+*
+* Description: Multiplication of two polynomials in NTT domain
+*
+* Arguments:   - poly *r:       pointer to output polynomial
+*              - const poly *a: pointer to first input polynomial
+*              - const poly *b: pointer to second input polynomial
+**************************************************/
+void poly_basemul(poly *r, const poly *a, const poly *b)
+{
+	for(int i = 0; i < NTRUPLUS_N/8; ++i)
+	{
+		basemul(r->coeffs + 8*i, a->coeffs + 8*i, b->coeffs + 8*i, zetas[144 + i]);
+		basemul(r->coeffs + 8*i + 4, a->coeffs + 8*i + 4, b->coeffs + 8*i + 4, -zetas[144 + i]);
+	}
+}
+
+/*************************************************
+* Name:        poly_basemul_add
+*
+* Description: Multiplication of two polynomials in NTT domain
+*
+* Arguments:   - poly *r:       pointer to output polynomial
+*              - const poly *a: pointer to first input polynomial
+*              - const poly *b: pointer to second input polynomial
+**************************************************/
+void poly_basemul_add(poly *r, const poly *a, const poly *b, const poly *c)
+{
+	for(int i = 0; i < NTRUPLUS_N/8; ++i)
+	{
+		basemul_add(r->coeffs + 8*i, a->coeffs + 8*i, b->coeffs + 8*i, c->coeffs + 8*i, zetas[144 + i]);
+		basemul_add(r->coeffs + 8*i + 4, a->coeffs + 8*i + 4, b->coeffs + 8*i + 4, c->coeffs + 8*i + 4, -zetas[144 + i]);
+	}
+}
+
+/*************************************************
+* Name:        poly_sub
+*
+* Description: Subtract two polynomials; no modular reduction is performed
+*
+* Arguments: - poly *r: pointer to output polynomial
+*            - const poly *a: pointer to first input polynomial
+*            - const poly *b: pointer to second input polynomial
+**************************************************/
+void poly_sub(poly *r, const poly *a, const poly *b)
+{
+	for(int i = 0; i < NTRUPLUS_N; ++i)
+		r->coeffs[i] = a->coeffs[i] - b->coeffs[i];
+}
+
+/*************************************************
+* Name:        poly_triple
+*
+* Description: Multiply polynomial by 3; no modular reduction is performed
+*
+* Arguments: - poly *r: pointer to output polynomial
+*            - const poly *a: pointer to input polynomial
+**************************************************/
+void poly_triple(poly *r, const poly *a) 
+{
+	for(int i = 0; i < NTRUPLUS_N; ++i)
+		r->coeffs[i] = 3*a->coeffs[i];
+}
+
+/*************************************************
+* Name:        poly_crepmod3
+*
+* Description: Compute modulus 3 operation to polynomial
+*
+* Arguments: - poly *r: pointer to output polynomial
+*            - const poly *a: pointer to input polynomial
+**************************************************/
+void poly_crepmod3(poly *r, const poly *a)
+{
+	for(int i = 0; i < NTRUPLUS_N; i++)
+		r->coeffs[i] = crepmod3(a->coeffs[i]);
+}

data/vendor/NTRUplus/NTRU+1152/poly.h

@@ -0,0 +1,31 @@
+#ifndef POLY_H
+#define POLY_H
+
+#include <stdint.h>
+#include "params.h"
+
+/*
+ * Elements of R_q = Z_q[X]/(X^n - X^n/2 + 1). Represents polynomial
+ * coeffs[0] + X*coeffs[1] + X^2*xoeffs[2] + ... + X^{n-1}*coeffs[n-1]
+ */
+typedef struct{
+	int16_t coeffs[NTRUPLUS_N];
+} poly;
+
+void poly_tobytes(uint8_t r[NTRUPLUS_POLYBYTES], const poly *a);
+void poly_frombytes(poly *r, const uint8_t a[NTRUPLUS_POLYBYTES]);
+
+void poly_cbd1(poly *r, const uint8_t buf[NTRUPLUS_N/4]);
+void poly_sotp_encode(poly *r, const uint8_t msg[NTRUPLUS_N/8], const uint8_t buf[NTRUPLUS_N/4]);
+int  poly_sotp_decode(uint8_t msg[NTRUPLUS_N/8], const poly *a, const uint8_t buf[NTRUPLUS_N/4]);
+
+void poly_ntt(poly *r);
+void poly_invntt(poly *r);
+int  poly_baseinv(poly *r, const poly *a);
+void poly_basemul(poly *r, const poly *a, const poly *b);
+void poly_basemul_add(poly *r, const poly *a, const poly *b, const poly *c);
+void poly_sub(poly *r, const poly *a, const poly *b);
+void poly_triple(poly *r, const poly *a);
+void poly_crepmod3(poly *r, const poly *a);
+
+#endif

data/vendor/NTRUplus/NTRU+1152/randombytes.h

@@ -0,0 +1,9 @@
+#ifndef RANDOMBYTES_H
+#define RANDOMBYTES_H
+
+#include <stddef.h>
+#include <stdint.h>
+
+int randombytes(uint8_t *out, size_t outlen);
+
+#endif

data/vendor/NTRUplus/NTRU+1152/symmetric.c

@@ -0,0 +1,39 @@
+#include <string.h>
+#include "symmetric.h"
+#include "fips202/fips202.h"
+
+#define HASH_F_INBYTES  (NTRUPLUS_POLYBYTES)
+#define HASH_F_OUTBYTES (32)
+
+#define HASH_G_INBYTES  (NTRUPLUS_POLYBYTES)
+#define HASH_G_OUTBYTES (NTRUPLUS_N / 4)
+
+#define HASH_H_INBYTES  (NTRUPLUS_N / 8 + NTRUPLUS_SYMBYTES)
+#define HASH_H_OUTBYTES (NTRUPLUS_SSBYTES + NTRUPLUS_N / 4)
+
+void hash_f(uint8_t *buf, const uint8_t *msg)
+{
+    uint8_t data[1 + HASH_F_INBYTES];
+
+    data[0] = 0x00;
+    memcpy(data + 1, msg, HASH_F_INBYTES);
+    shake256(buf, HASH_F_OUTBYTES, data, HASH_F_INBYTES + 1);
+}
+
+void hash_g(uint8_t *buf, const uint8_t *msg)
+{
+    uint8_t data[1 + HASH_G_INBYTES];
+
+    data[0] = 0x01;
+    memcpy(data + 1, msg, HASH_G_INBYTES);
+    shake256(buf, HASH_G_OUTBYTES, data, HASH_G_INBYTES + 1);
+}
+
+void hash_h(uint8_t *buf, const uint8_t *msg)
+{
+    uint8_t data[1 + HASH_H_INBYTES];
+
+    data[0] = 0x02;
+    memcpy(data + 1, msg, HASH_H_INBYTES);
+    shake256(buf, HASH_H_OUTBYTES, data, HASH_H_INBYTES + 1);
+}

data/vendor/NTRUplus/NTRU+1152/symmetric.h

@@ -0,0 +1,12 @@
+#ifndef SYMMETRIC_H
+#define SYMMETRIC_H
+
+#include <stddef.h>
+#include <stdint.h>
+#include "params.h"
+
+void hash_f(uint8_t *buf, const uint8_t *msg);
+void hash_g(uint8_t *buf, const uint8_t *msg);
+void hash_h(uint8_t *buf, const uint8_t *msg);
+
+#endif /* SYMMETRIC_H */

data/vendor/NTRUplus/NTRU+768/api.h

@@ -0,0 +1,24 @@
+#ifndef API_H
+#define API_H
+
+#include "params.h"
+
+#define CRYPTO_SECRETKEYBYTES  NTRUPLUS_SECRETKEYBYTES
+#define CRYPTO_PUBLICKEYBYTES  NTRUPLUS_PUBLICKEYBYTES
+#define CRYPTO_CIPHERTEXTBYTES NTRUPLUS_CIPHERTEXTBYTES
+#define CRYPTO_BYTES           NTRUPLUS_SSBYTES
+
+#define CRYPTO_ALGNAME NTRUPLUS_ALGNAME
+
+int crypto_kem_keypair(unsigned char *pk,
+                       unsigned char *sk);
+
+int crypto_kem_enc(unsigned char *ct,
+                   unsigned char *ss,
+                   const unsigned char *pk);
+
+int crypto_kem_dec(unsigned char *ss,
+                   const unsigned char *ct,
+                   const unsigned char *sk);
+
+#endif