oqs-sys 0.11.0+liboqs-0.13.0

Bindings to liboqs
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68

// SPDX-License-Identifier: Apache-2.0

#ifndef ARITHMETIC_DYNAMIC_H
#define ARITHMETIC_DYNAMIC_H

#include <stdint.h>
#include <simple_arithmetic.h>

static
inline void m_vec_copy(int m_vec_limbs, const uint64_t *in, uint64_t *out) {
    for (int i = 0; i < m_vec_limbs; i++) {
        out[i] = in[i];
    }
}

static
inline void m_vec_add(int m_vec_limbs, const uint64_t *in, uint64_t *acc) {
    for (int i = 0; i < m_vec_limbs; i++) {
        acc[i] ^= in[i];
    }
}

static inline
void m_vec_mul_add(int m_vec_limbs, const uint64_t *in, unsigned char a, uint64_t *acc) {
    for(int i=0; i < m_vec_limbs; i++){
        acc[i] ^= gf16v_mul_u64(in[i], a);
    }
}

inline
static void m_vec_mul_add_x (int m_vec_limbs, const uint64_t *in, uint64_t *acc) {
    uint64_t mask_msb = 0x8888888888888888ULL;
    for(int i=0; i < m_vec_limbs; i++){
        uint64_t t = in[i] & mask_msb;
        acc[i] ^= ((in[i] ^ t) << 1) ^ ((t >> 3) * 3);
    }
}

inline
static void m_vec_mul_add_x_inv (int m_vec_limbs, const uint64_t *in, uint64_t *acc) {
    uint64_t mask_lsb = 0x1111111111111111ULL;
    for(int i=0; i < m_vec_limbs; i++){
        uint64_t t = in[i] & mask_lsb;
        acc[i] ^= ((in[i] ^ t) >> 1) ^ (t * 9);
    }
}

static
inline void m_vec_multiply_bins (int m_vec_limbs, uint64_t *bins, uint64_t *out) {

    m_vec_mul_add_x_inv (m_vec_limbs, bins +  5 * m_vec_limbs, bins +  10 * m_vec_limbs);
    m_vec_mul_add_x (m_vec_limbs, bins + 11 * m_vec_limbs, bins + 12 * m_vec_limbs);
    m_vec_mul_add_x_inv (m_vec_limbs, bins +  10 * m_vec_limbs, bins +  7 * m_vec_limbs);
    m_vec_mul_add_x (m_vec_limbs, bins + 12 * m_vec_limbs, bins +  6 * m_vec_limbs);
    m_vec_mul_add_x_inv (m_vec_limbs, bins +  7 * m_vec_limbs, bins +  14 * m_vec_limbs);
    m_vec_mul_add_x (m_vec_limbs, bins +  6 * m_vec_limbs, bins +  3 * m_vec_limbs);
    m_vec_mul_add_x_inv (m_vec_limbs, bins +  14 * m_vec_limbs, bins +  15 * m_vec_limbs);
    m_vec_mul_add_x (m_vec_limbs, bins +  3 * m_vec_limbs, bins +  8 * m_vec_limbs);
    m_vec_mul_add_x_inv (m_vec_limbs, bins +  15 * m_vec_limbs, bins +  13 * m_vec_limbs);
    m_vec_mul_add_x (m_vec_limbs, bins +  8 * m_vec_limbs, bins +  4 * m_vec_limbs);
    m_vec_mul_add_x_inv (m_vec_limbs, bins +  13 * m_vec_limbs, bins +  9 * m_vec_limbs);
    m_vec_mul_add_x (m_vec_limbs, bins +  4 * m_vec_limbs, bins +  2 * m_vec_limbs);
    m_vec_mul_add_x_inv (m_vec_limbs, bins +   9 * m_vec_limbs, bins +  1 * m_vec_limbs);
    m_vec_mul_add_x (m_vec_limbs, bins +  2 * m_vec_limbs, bins +  1 * m_vec_limbs);
    m_vec_copy (m_vec_limbs, bins + m_vec_limbs, out);
}

#endif