Update the imported curve448 code to use OpenSSL copyright headers
[openssl.git] / crypto / ec / curve448 / arch_32 / f_impl.h
1 /*
2  * Copyright 2017 The OpenSSL Project Authors. All Rights Reserved.
3  * Copyright 2014-2016 Cryptography Research, Inc.
4  *
5  * Licensed under the OpenSSL license (the "License").  You may not use
6  * this file except in compliance with the License.  You can obtain a copy
7  * in the file LICENSE in the source distribution or at
8  * https://www.openssl.org/source/license.html
9  *
10  * Originally written by Mike Hamburg
11  */
12 #define GF_HEADROOM 2
13 #define LIMB(x) (x)&((1<<28)-1), (x)>>28
14 #define FIELD_LITERAL(a,b,c,d,e,f,g,h) \
15     {{LIMB(a),LIMB(b),LIMB(c),LIMB(d),LIMB(e),LIMB(f),LIMB(g),LIMB(h)}}
16     
17 #define LIMB_PLACE_VALUE(i) 28
18
19 void gf_add_RAW (gf out, const gf a, const gf b) {
20     unsigned int i;
21
22     for (i=0; i<sizeof(*out)/sizeof(out->limb[0]); i++) {
23         out->limb[i] = a->limb[i] + b->limb[i];
24     }
25 }
26
27 void gf_sub_RAW (gf out, const gf a, const gf b) {
28     unsigned int i;
29
30     for (i=0; i<sizeof(*out)/sizeof(out->limb[0]); i++) {
31         out->limb[i] = a->limb[i] - b->limb[i];
32     }
33 }
34
35 void gf_bias (gf a, int amt) {
36     unsigned int i;
37     uint32_t co1 = ((1<<28)-1)*amt, co2 = co1-amt;
38
39     for (i=0; i<sizeof(*a)/sizeof(a->limb[0]); i++) {
40         a->limb[i] += (i==sizeof(*a)/sizeof(a->limb[0])/2) ? co2 : co1;
41     }
42 }
43
44 void gf_weak_reduce (gf a) {
45     uint32_t mask = (1<<28) - 1;
46     uint32_t tmp = a->limb[15] >> 28;
47     unsigned int i;
48
49     a->limb[8] += tmp;
50     for (i=15; i>0; i--) {
51         a->limb[i] = (a->limb[i] & mask) + (a->limb[i-1]>>28);
52     }
53     a->limb[0] = (a->limb[0] & mask) + tmp;
54 }
55