|
Packit |
141393 |
/**
|
|
Packit |
141393 |
* @file p448/f_generic.c
|
|
Packit |
141393 |
* @author Mike Hamburg
|
|
Packit |
141393 |
*
|
|
Packit |
141393 |
* @copyright
|
|
Packit |
141393 |
* Copyright (c) 2015-2016 Cryptography Research, Inc. \n
|
|
Packit |
141393 |
* Released under the MIT License. See LICENSE.txt for license information.
|
|
Packit |
141393 |
*
|
|
Packit |
141393 |
* @brief Generic arithmetic which has to be compiled per field.
|
|
Packit |
141393 |
*
|
|
Packit |
141393 |
* @warning This file was automatically generated in Python.
|
|
Packit |
141393 |
* Please do not edit it.
|
|
Packit |
141393 |
*/
|
|
Packit |
141393 |
#include "field.h"
|
|
Packit |
141393 |
|
|
Packit |
141393 |
static const gf MODULUS = {FIELD_LITERAL(
|
|
Packit |
141393 |
0xffffffffffffff, 0xffffffffffffff, 0xffffffffffffff, 0xffffffffffffff, 0xfffffffffffffe, 0xffffffffffffff, 0xffffffffffffff, 0xffffffffffffff
|
|
Packit |
141393 |
)};
|
|
Packit |
141393 |
|
|
Packit |
141393 |
#if P_MOD_8 == 5
|
|
Packit |
141393 |
const gf SQRT_MINUS_ONE = {FIELD_LITERAL(
|
|
Packit |
141393 |
/* NOPE */
|
|
Packit |
141393 |
)};
|
|
Packit |
141393 |
#endif
|
|
Packit |
141393 |
|
|
Packit |
141393 |
/** Serialize to wire format. */
|
|
Packit |
141393 |
void cryptonite_gf_serialize (uint8_t serial[SER_BYTES], const gf x, int with_hibit) {
|
|
Packit |
141393 |
gf red;
|
|
Packit |
141393 |
cryptonite_gf_copy(red, x);
|
|
Packit |
141393 |
cryptonite_gf_strong_reduce(red);
|
|
Packit |
141393 |
if (!with_hibit) { assert(cryptonite_gf_hibit(red) == 0); }
|
|
Packit |
141393 |
|
|
Packit |
141393 |
unsigned int j=0, fill=0;
|
|
Packit |
141393 |
dword_t buffer = 0;
|
|
Packit |
141393 |
UNROLL for (unsigned int i=0; i<(with_hibit ? X_SER_BYTES : SER_BYTES); i++) {
|
|
Packit |
141393 |
if (fill < 8 && j < NLIMBS) {
|
|
Packit |
141393 |
buffer |= ((dword_t)red->limb[LIMBPERM(j)]) << fill;
|
|
Packit |
141393 |
fill += LIMB_PLACE_VALUE(LIMBPERM(j));
|
|
Packit |
141393 |
j++;
|
|
Packit |
141393 |
}
|
|
Packit |
141393 |
serial[i] = buffer;
|
|
Packit |
141393 |
fill -= 8;
|
|
Packit |
141393 |
buffer >>= 8;
|
|
Packit |
141393 |
}
|
|
Packit |
141393 |
}
|
|
Packit |
141393 |
|
|
Packit |
141393 |
/** Return high bit of x = low bit of 2x mod p */
|
|
Packit |
141393 |
mask_t cryptonite_gf_hibit(const gf x) {
|
|
Packit |
141393 |
gf y;
|
|
Packit |
141393 |
cryptonite_gf_add(y,x,x);
|
|
Packit |
141393 |
cryptonite_gf_strong_reduce(y);
|
|
Packit |
141393 |
return -(y->limb[0]&1;;
|
|
Packit |
141393 |
}
|
|
Packit |
141393 |
|
|
Packit |
141393 |
/** Deserialize from wire format; return -1 on success and 0 on failure. */
|
|
Packit |
141393 |
mask_t cryptonite_gf_deserialize (gf x, const uint8_t serial[SER_BYTES], int with_hibit) {
|
|
Packit |
141393 |
unsigned int j=0, fill=0;
|
|
Packit |
141393 |
dword_t buffer = 0;
|
|
Packit |
141393 |
dsword_t scarry = 0;
|
|
Packit |
141393 |
UNROLL for (unsigned int i=0; i
|
|
Packit |
141393 |
UNROLL while (fill < LIMB_PLACE_VALUE(LIMBPERM(i)) && j < (with_hibit ? X_SER_BYTES : SER_BYTES)) {
|
|
Packit |
141393 |
buffer |= ((dword_t)serial[j]) << fill;
|
|
Packit |
141393 |
fill += 8;
|
|
Packit |
141393 |
j++;
|
|
Packit |
141393 |
}
|
|
Packit |
141393 |
x->limb[LIMBPERM(i)] = (i
|
|
Packit |
141393 |
fill -= LIMB_PLACE_VALUE(LIMBPERM(i));
|
|
Packit |
141393 |
buffer >>= LIMB_PLACE_VALUE(LIMBPERM(i));
|
|
Packit |
141393 |
scarry = (scarry + x->limb[LIMBPERM(i)] - MODULUS->limb[LIMBPERM(i)]) >> (8*sizeof(word_t));
|
|
Packit |
141393 |
}
|
|
Packit |
141393 |
mask_t succ = with_hibit ? -(mask_t)1 : ~cryptonite_gf_hibit(x);
|
|
Packit |
141393 |
return succ & word_is_zero(buffer) & ~word_is_zero(scarry);
|
|
Packit |
141393 |
}
|
|
Packit |
141393 |
|
|
Packit |
141393 |
/** Reduce to canonical form. */
|
|
Packit |
141393 |
void cryptonite_gf_strong_reduce (gf a) {
|
|
Packit |
141393 |
/* first, clear high */
|
|
Packit |
141393 |
cryptonite_gf_weak_reduce(a); /* Determined to have negligible perf impact. */
|
|
Packit |
141393 |
|
|
Packit |
141393 |
/* now the total is less than 2p */
|
|
Packit |
141393 |
|
|
Packit |
141393 |
/* compute total_value - p. No need to reduce mod p. */
|
|
Packit |
141393 |
dsword_t scarry = 0;
|
|
Packit |
141393 |
for (unsigned int i=0; i
|
|
Packit |
141393 |
scarry = scarry + a->limb[LIMBPERM(i)] - MODULUS->limb[LIMBPERM(i)];
|
|
Packit |
141393 |
a->limb[LIMBPERM(i)] = scarry & LIMB_MASK(LIMBPERM(i));
|
|
Packit |
141393 |
scarry >>= LIMB_PLACE_VALUE(LIMBPERM(i));
|
|
Packit |
141393 |
}
|
|
Packit |
141393 |
|
|
Packit |
141393 |
/* uncommon case: it was >= p, so now scarry = 0 and this = x
|
|
Packit |
141393 |
* common case: it was < p, so now scarry = -1 and this = x - p + 2^255
|
|
Packit |
141393 |
* so let's add back in p. will carry back off the top for 2^255.
|
|
Packit |
141393 |
*/
|
|
Packit |
141393 |
assert(word_is_zero(scarry) | word_is_zero(scarry+1));
|
|
Packit |
141393 |
|
|
Packit |
141393 |
word_t scarry_0 = scarry;
|
|
Packit |
141393 |
dword_t carry = 0;
|
|
Packit |
141393 |
|
|
Packit |
141393 |
/* add it back */
|
|
Packit |
141393 |
for (unsigned int i=0; i
|
|
Packit |
141393 |
carry = carry + a->limb[LIMBPERM(i)] + (scarry_0 & MODULUS->limb[LIMBPERM(i)]);
|
|
Packit |
141393 |
a->limb[LIMBPERM(i)] = carry & LIMB_MASK(LIMBPERM(i));
|
|
Packit |
141393 |
carry >>= LIMB_PLACE_VALUE(LIMBPERM(i));
|
|
Packit |
141393 |
}
|
|
Packit |
141393 |
|
|
Packit |
141393 |
assert(word_is_zero(carry + scarry_0));
|
|
Packit |
141393 |
}
|
|
Packit |
141393 |
|
|
Packit |
141393 |
/** Add two gf elements */
|
|
Packit |
141393 |
void cryptonite_gf_sub (gf d, const gf a, const gf b) {
|
|
Packit |
141393 |
cryptonite_gf_sub_RAW ( d, a, b );
|
|
Packit |
141393 |
cryptonite_gf_bias( d, 2 );
|
|
Packit |
141393 |
cryptonite_gf_weak_reduce ( d );
|
|
Packit |
141393 |
}
|
|
Packit |
141393 |
|
|
Packit |
141393 |
/** Subtract d = a-b */
|
|
Packit |
141393 |
void cryptonite_gf_add (gf d, const gf a, const gf b) {
|
|
Packit |
141393 |
cryptonite_gf_add_RAW ( d, a, b );
|
|
Packit |
141393 |
cryptonite_gf_weak_reduce ( d );
|
|
Packit |
141393 |
}
|
|
Packit |
141393 |
|
|
Packit |
141393 |
/** Compare a==b */
|
|
Packit |
141393 |
mask_t cryptonite_gf_eq(const gf a, const gf b) {
|
|
Packit |
141393 |
gf c;
|
|
Packit |
141393 |
cryptonite_gf_sub(c,a,b);
|
|
Packit |
141393 |
cryptonite_gf_strong_reduce(c);
|
|
Packit |
141393 |
mask_t ret=0;
|
|
Packit |
141393 |
for (unsigned int i=0; i
|
|
Packit |
141393 |
ret |= c->limb[LIMBPERM(i)];
|
|
Packit |
141393 |
}
|
|
Packit |
141393 |
|
|
Packit |
141393 |
return word_is_zero(ret);
|
|
Packit |
141393 |
}
|