Line data Source code
1 : /* fd_aes_ref.c was imported from the OpenSSL project circa 2023-Aug.
2 : Original source files: crypto/evp/e_aes.c crypto/modes/gcm128.c */
3 :
4 : #include "fd_aes_gcm.h"
5 :
6 : #include <assert.h>
7 :
8 4349806 : #define fd_gcm_init fd_gcm_init_4bit
9 4001692 : #define fd_gcm_gmult fd_gcm_gmult_4bit
10 14229342 : #define fd_gcm_ghash fd_gcm_ghash_4bit
11 :
12 : static void
13 : fd_aes_gcm_setiv( fd_aes_gcm_ref_t * gcm,
14 4349806 : uchar const iv[ 12 ] ) {
15 :
16 4349806 : uint ctr;
17 4349806 : gcm->len.u[ 0 ] = 0; /* AAD length */
18 4349806 : gcm->len.u[ 1 ] = 0; /* Message length */
19 4349806 : gcm->ares = 0;
20 4349806 : gcm->mres = 0;
21 :
22 4349806 : memcpy( gcm->Yi.c, iv, 12 );
23 4349806 : gcm->Yi.c[12] = 0;
24 4349806 : gcm->Yi.c[13] = 0;
25 4349806 : gcm->Yi.c[14] = 0;
26 4349806 : gcm->Yi.c[15] = 1;
27 4349806 : ctr = 1;
28 :
29 4349806 : gcm->Xi.u[0] = 0;
30 4349806 : gcm->Xi.u[1] = 0;
31 :
32 4349806 : fd_aes_encrypt( gcm->Yi.c, gcm->EK0.c, &gcm->key );
33 4349806 : ctr++;
34 :
35 4349806 : gcm->Yi.d[3] = fd_uint_bswap( ctr );
36 4349806 : }
37 :
38 : void
39 : fd_aes_128_gcm_init_ref( fd_aes_gcm_ref_t * gcm,
40 : uchar const key[ 16 ],
41 4349806 : uchar const iv[ 12 ] ) {
42 :
43 : /* TODO: Check key size */
44 :
45 4349806 : memset( gcm, 0, sizeof(fd_aes_gcm_t) );
46 :
47 4349806 : fd_aes_key_ref_t * ks = &gcm->key;
48 4349806 : fd_aes_set_encrypt_key( key, 128, ks );
49 :
50 4349806 : fd_aes_encrypt( gcm->H.c, gcm->H.c, ks );
51 4349806 : gcm->H.u[ 0 ] = fd_ulong_bswap( gcm->H.u[ 0 ] );
52 4349806 : gcm->H.u[ 1 ] = fd_ulong_bswap( gcm->H.u[ 1 ] );
53 :
54 4349806 : fd_gcm_init( gcm->Htable, gcm->H.u );
55 4349806 : fd_aes_gcm_setiv( gcm, iv );
56 4349806 : }
57 :
58 : static int
59 : fd_gcm128_aad( fd_aes_gcm_ref_t * aes_gcm,
60 : uchar const * aad,
61 4349806 : ulong aad_sz ) {
62 :
63 4349806 : ulong alen = aes_gcm->len.u[ 0 ];
64 :
65 4349806 : if( FD_UNLIKELY( aes_gcm->len.u[ 1 ] ) )
66 0 : return -2;
67 :
68 4349806 : alen += aad_sz;
69 4349806 : if (alen > (1UL<<61) || (sizeof(aad_sz) == 8 && alen < aad_sz))
70 0 : return -1;
71 4349806 : aes_gcm->len.u[0] = alen;
72 :
73 4349806 : uint n = aes_gcm->ares;
74 4349806 : if (n) {
75 0 : while (n && aad_sz) {
76 0 : aes_gcm->Xi.c[n] ^= *(aad++);
77 0 : --aad_sz;
78 0 : n = (n + 1) % 16;
79 0 : }
80 0 : if (n == 0)
81 0 : fd_gcm_gmult( aes_gcm->Xi.u, aes_gcm->Htable );
82 0 : else {
83 0 : aes_gcm->ares = n;
84 0 : return 0;
85 0 : }
86 0 : }
87 4349806 : ulong i;
88 4349806 : if ((i = (aad_sz & (ulong)-16))) {
89 4042035 : fd_gcm_ghash( aes_gcm->Xi.u, aes_gcm->Htable, aad, i );
90 4042035 : aad += i;
91 4042035 : aad_sz -= i;
92 4042035 : }
93 4349806 : if (aad_sz) {
94 4349806 : n = (unsigned int)aad_sz;
95 43557366 : for (i = 0; i < aad_sz; ++i)
96 39207560 : aes_gcm->Xi.c[i] ^= aad[i];
97 4349806 : }
98 :
99 4349806 : aes_gcm->ares = n;
100 4349806 : return 0;
101 4349806 : }
102 :
103 : /* TODO separate reference code and GCM128 */
104 :
105 : static int
106 : fd_gcm128_encrypt( fd_aes_gcm_ref_t * ctx,
107 : uchar const * in,
108 : uchar * out,
109 2175155 : ulong len ) {
110 :
111 2175155 : uint n, ctr, mres;
112 2175155 : ulong i;
113 2175155 : ulong mlen = ctx->len.u[1];
114 2175155 : void *key = &ctx->key;
115 :
116 2175155 : mlen += len;
117 2175155 : if (mlen > ((1UL<<36) - 32) || (sizeof(len) == 8 && mlen < len))
118 0 : return -1;
119 2175155 : ctx->len.u[1] = mlen;
120 :
121 2175155 : mres = ctx->mres;
122 :
123 2175155 : if (ctx->ares) {
124 : /* First call to encrypt finalizes GHASH(AAD) */
125 2175155 : if (len == 0) {
126 2000093 : fd_gcm_gmult( ctx->Xi.u, ctx->Htable );
127 2000093 : ctx->ares = 0;
128 2000093 : return 0;
129 2000093 : }
130 175062 : memcpy(ctx->Xn, ctx->Xi.c, sizeof(ctx->Xi));
131 175062 : ctx->Xi.u[0] = 0;
132 175062 : ctx->Xi.u[1] = 0;
133 175062 : mres = sizeof(ctx->Xi);
134 175062 : ctx->ares = 0;
135 175062 : }
136 :
137 175062 : ctr = fd_uint_bswap( ctx->Yi.d[3] );
138 :
139 175062 : n = mres % 16;
140 136684656 : for (i = 0; i < len; ++i) {
141 136509594 : if (n == 0) {
142 8634154 : fd_aes_encrypt( ctx->Yi.c, ctx->EKi.c, key );
143 8634154 : ++ctr;
144 8634154 : ctx->Yi.d[3] = fd_uint_bswap( ctr );
145 8634154 : }
146 136509594 : ctx->Xn[mres++] = out[i] = in[i] ^ ctx->EKi.c[n];
147 136509594 : n = (n + 1) % 16;
148 136509594 : if (mres == sizeof(ctx->Xn)) {
149 2777040 : fd_gcm_ghash( ctx->Xi.u, ctx->Htable, ctx->Xn, sizeof(ctx->Xn) );
150 2777040 : mres = 0;
151 2777040 : }
152 136509594 : }
153 :
154 175062 : ctx->mres = mres;
155 175062 : return 0;
156 2175155 : }
157 :
158 : static int
159 : fd_gcm128_decrypt( fd_aes_gcm_ref_t * ctx,
160 : uchar const * in,
161 : uchar * out,
162 2174651 : ulong len ) {
163 :
164 2174651 : uint n, ctr, mres;
165 2174651 : ulong i;
166 2174651 : ulong mlen = ctx->len.u[1];
167 2174651 : void * key = &ctx->key;
168 :
169 2174651 : mlen += len;
170 2174651 : if (mlen > ((1UL<<36) - 32) || (sizeof(len) == 8 && mlen < len))
171 0 : return -1;
172 2174651 : ctx->len.u[1] = mlen;
173 :
174 2174651 : mres = ctx->mres;
175 :
176 2174651 : if (ctx->ares) {
177 : /* First call to decrypt finalizes GHASH(AAD) */
178 2174651 : if (len == 0) {
179 2001599 : fd_gcm_gmult( ctx->Xi.u, ctx->Htable );
180 2001599 : ctx->ares = 0;
181 2001599 : return 0;
182 2001599 : }
183 173052 : memcpy(ctx->Xn, ctx->Xi.c, sizeof(ctx->Xi));
184 173052 : ctx->Xi.u[0] = 0;
185 173052 : ctx->Xi.u[1] = 0;
186 173052 : mres = sizeof(ctx->Xi);
187 173052 : ctx->ares = 0;
188 173052 : }
189 :
190 173052 : ctr = fd_uint_bswap( ctx->Yi.d[3] );
191 :
192 173052 : n = mres % 16;
193 137049099 : for (i = 0; i < len; ++i) {
194 136876047 : uchar c;
195 136876047 : if (n == 0) {
196 8652149 : fd_aes_encrypt( ctx->Yi.c, ctx->EKi.c, key );
197 8652149 : ++ctr;
198 8652149 : ctx->Yi.d[3] = fd_uint_bswap( ctr );
199 8652149 : }
200 136876047 : out[i] = (ctx->Xn[mres++] = c = in[i]) ^ ctx->EKi.c[n];
201 136876047 : n = (n + 1) % 16;
202 136876047 : if (mres == sizeof(ctx->Xn)) {
203 2787042 : fd_gcm_ghash( ctx->Xi.u, ctx->Htable, ctx->Xn, sizeof(ctx->Xn) );
204 2787042 : mres = 0;
205 2787042 : }
206 136876047 : }
207 :
208 173052 : ctx->mres = mres;
209 173052 : return 0;
210 2174651 : }
211 :
212 : static void
213 4349806 : fd_gcm128_finish( fd_aes_gcm_ref_t * ctx ) {
214 :
215 4349806 : ulong alen = ctx->len.u[0] << 3; // 176
216 4349806 : ulong clen = ctx->len.u[1] << 3; // 9296
217 :
218 4349806 : struct {
219 4349806 : ulong hi;
220 4349806 : ulong lo;
221 4349806 : } bitlen;
222 4349806 : uint mres = ctx->mres;
223 :
224 4349806 : if( mres ) {
225 347904 : uint blocks = (mres + 15u) & 0xfffffff0u; // 16
226 :
227 347904 : memset(ctx->Xn + mres, 0, blocks - mres);
228 347904 : mres = blocks;
229 347904 : if (mres == sizeof(ctx->Xn)) {
230 273419 : fd_gcm_ghash( ctx->Xi.u, ctx->Htable, ctx->Xn, mres );
231 273419 : mres = 0;
232 273419 : }
233 4001902 : } else if( ctx->ares ) {
234 0 : fd_gcm_gmult( ctx->Xi.u, ctx->Htable );
235 0 : }
236 :
237 4349806 : alen = fd_ulong_bswap( alen );
238 4349806 : clen = fd_ulong_bswap( clen );
239 :
240 4349806 : bitlen.hi = alen;
241 4349806 : bitlen.lo = clen;
242 4349806 : memcpy( ctx->Xn + mres, &bitlen, sizeof(bitlen) );
243 4349806 : mres += (uint)sizeof(bitlen);
244 4349806 : fd_gcm_ghash( ctx->Xi.u, ctx->Htable, ctx->Xn, mres );
245 :
246 4349806 : ctx->Xi.u[0] ^= ctx->EK0.u[0];
247 4349806 : ctx->Xi.u[1] ^= ctx->EK0.u[1];
248 4349806 : }
249 :
250 : void
251 : fd_aes_gcm_encrypt_ref( fd_aes_gcm_ref_t * aes_gcm,
252 : uchar * c,
253 : uchar const * p,
254 : ulong sz,
255 : uchar const * aad,
256 : ulong aad_sz,
257 2175155 : uchar tag[ 16 ] ) {
258 :
259 2175155 : fd_gcm128_aad( aes_gcm, aad, aad_sz );
260 :
261 2175155 : ulong bulk = 0UL;
262 2175155 : assert( 0==fd_gcm128_encrypt( aes_gcm, p+bulk, c+bulk, sz-bulk ) );
263 :
264 : /* CRYPTO_gcm128_tag */
265 2175155 : fd_gcm128_finish( aes_gcm );
266 2175155 : fd_memcpy( tag, aes_gcm->Xi.c, 16 );
267 2175155 : }
268 :
269 : int
270 : fd_aes_gcm_decrypt_ref( fd_aes_gcm_ref_t * aes_gcm,
271 : uchar const * c,
272 : uchar * p,
273 : ulong sz,
274 : uchar const * aad,
275 : ulong aad_sz,
276 2174651 : uchar const tag[ 16 ] ) {
277 :
278 2174651 : fd_gcm128_aad( aes_gcm, aad, aad_sz );
279 :
280 2174651 : ulong bulk = 0UL;
281 2174651 : assert( 0==fd_gcm128_decrypt( aes_gcm, c+bulk, p+bulk, sz-bulk ) );
282 :
283 : /* CRYPTO_gcm128_finish */
284 2174651 : fd_gcm128_finish( aes_gcm );
285 2174651 : return 0==memcmp( aes_gcm->Xi.c, tag, 16 ); /* TODO USE CONSTANT TIME COMPARE */
286 2174651 : }
|