1 /* 2 * Copyright (c) 2016 Thomas Pornin <pornin@bolet.org> 3 * 4 * Permission is hereby granted, free of charge, to any person obtaining 5 * a copy of this software and associated documentation files (the 6 * "Software"), to deal in the Software without restriction, including 7 * without limitation the rights to use, copy, modify, merge, publish, 8 * distribute, sublicense, and/or sell copies of the Software, and to 9 * permit persons to whom the Software is furnished to do so, subject to 10 * the following conditions: 11 * 12 * The above copyright notice and this permission notice shall be 13 * included in all copies or substantial portions of the Software. 14 * 15 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, 16 * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF 17 * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND 18 * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS 19 * BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN 20 * ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN 21 * CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE 22 * SOFTWARE. 23 */ 24 25 #include "inner.h" 26 27 /* 28 * Inverse S-box (used in key schedule for decryption). 29 */ 30 static const unsigned char iS[] = { 31 0x52, 0x09, 0x6A, 0xD5, 0x30, 0x36, 0xA5, 0x38, 0xBF, 0x40, 0xA3, 0x9E, 32 0x81, 0xF3, 0xD7, 0xFB, 0x7C, 0xE3, 0x39, 0x82, 0x9B, 0x2F, 0xFF, 0x87, 33 0x34, 0x8E, 0x43, 0x44, 0xC4, 0xDE, 0xE9, 0xCB, 0x54, 0x7B, 0x94, 0x32, 34 0xA6, 0xC2, 0x23, 0x3D, 0xEE, 0x4C, 0x95, 0x0B, 0x42, 0xFA, 0xC3, 0x4E, 35 0x08, 0x2E, 0xA1, 0x66, 0x28, 0xD9, 0x24, 0xB2, 0x76, 0x5B, 0xA2, 0x49, 36 0x6D, 0x8B, 0xD1, 0x25, 0x72, 0xF8, 0xF6, 0x64, 0x86, 0x68, 0x98, 0x16, 37 0xD4, 0xA4, 0x5C, 0xCC, 0x5D, 0x65, 0xB6, 0x92, 0x6C, 0x70, 0x48, 0x50, 38 0xFD, 0xED, 0xB9, 0xDA, 0x5E, 0x15, 0x46, 0x57, 0xA7, 0x8D, 0x9D, 0x84, 39 0x90, 0xD8, 0xAB, 0x00, 0x8C, 0xBC, 0xD3, 0x0A, 0xF7, 0xE4, 0x58, 0x05, 40 0xB8, 0xB3, 0x45, 0x06, 0xD0, 0x2C, 0x1E, 0x8F, 0xCA, 0x3F, 0x0F, 0x02, 41 0xC1, 0xAF, 0xBD, 0x03, 0x01, 0x13, 0x8A, 0x6B, 0x3A, 0x91, 0x11, 0x41, 42 0x4F, 0x67, 0xDC, 0xEA, 0x97, 0xF2, 0xCF, 0xCE, 0xF0, 0xB4, 0xE6, 0x73, 43 0x96, 0xAC, 0x74, 0x22, 0xE7, 0xAD, 0x35, 0x85, 0xE2, 0xF9, 0x37, 0xE8, 44 0x1C, 0x75, 0xDF, 0x6E, 0x47, 0xF1, 0x1A, 0x71, 0x1D, 0x29, 0xC5, 0x89, 45 0x6F, 0xB7, 0x62, 0x0E, 0xAA, 0x18, 0xBE, 0x1B, 0xFC, 0x56, 0x3E, 0x4B, 46 0xC6, 0xD2, 0x79, 0x20, 0x9A, 0xDB, 0xC0, 0xFE, 0x78, 0xCD, 0x5A, 0xF4, 47 0x1F, 0xDD, 0xA8, 0x33, 0x88, 0x07, 0xC7, 0x31, 0xB1, 0x12, 0x10, 0x59, 48 0x27, 0x80, 0xEC, 0x5F, 0x60, 0x51, 0x7F, 0xA9, 0x19, 0xB5, 0x4A, 0x0D, 49 0x2D, 0xE5, 0x7A, 0x9F, 0x93, 0xC9, 0x9C, 0xEF, 0xA0, 0xE0, 0x3B, 0x4D, 50 0xAE, 0x2A, 0xF5, 0xB0, 0xC8, 0xEB, 0xBB, 0x3C, 0x83, 0x53, 0x99, 0x61, 51 0x17, 0x2B, 0x04, 0x7E, 0xBA, 0x77, 0xD6, 0x26, 0xE1, 0x69, 0x14, 0x63, 52 0x55, 0x21, 0x0C, 0x7D 53 }; 54 55 static const uint32_t iSsm0[] = { 56 0x51F4A750, 0x7E416553, 0x1A17A4C3, 0x3A275E96, 0x3BAB6BCB, 0x1F9D45F1, 57 0xACFA58AB, 0x4BE30393, 0x2030FA55, 0xAD766DF6, 0x88CC7691, 0xF5024C25, 58 0x4FE5D7FC, 0xC52ACBD7, 0x26354480, 0xB562A38F, 0xDEB15A49, 0x25BA1B67, 59 0x45EA0E98, 0x5DFEC0E1, 0xC32F7502, 0x814CF012, 0x8D4697A3, 0x6BD3F9C6, 60 0x038F5FE7, 0x15929C95, 0xBF6D7AEB, 0x955259DA, 0xD4BE832D, 0x587421D3, 61 0x49E06929, 0x8EC9C844, 0x75C2896A, 0xF48E7978, 0x99583E6B, 0x27B971DD, 62 0xBEE14FB6, 0xF088AD17, 0xC920AC66, 0x7DCE3AB4, 0x63DF4A18, 0xE51A3182, 63 0x97513360, 0x62537F45, 0xB16477E0, 0xBB6BAE84, 0xFE81A01C, 0xF9082B94, 64 0x70486858, 0x8F45FD19, 0x94DE6C87, 0x527BF8B7, 0xAB73D323, 0x724B02E2, 65 0xE31F8F57, 0x6655AB2A, 0xB2EB2807, 0x2FB5C203, 0x86C57B9A, 0xD33708A5, 66 0x302887F2, 0x23BFA5B2, 0x02036ABA, 0xED16825C, 0x8ACF1C2B, 0xA779B492, 67 0xF307F2F0, 0x4E69E2A1, 0x65DAF4CD, 0x0605BED5, 0xD134621F, 0xC4A6FE8A, 68 0x342E539D, 0xA2F355A0, 0x058AE132, 0xA4F6EB75, 0x0B83EC39, 0x4060EFAA, 69 0x5E719F06, 0xBD6E1051, 0x3E218AF9, 0x96DD063D, 0xDD3E05AE, 0x4DE6BD46, 70 0x91548DB5, 0x71C45D05, 0x0406D46F, 0x605015FF, 0x1998FB24, 0xD6BDE997, 71 0x894043CC, 0x67D99E77, 0xB0E842BD, 0x07898B88, 0xE7195B38, 0x79C8EEDB, 72 0xA17C0A47, 0x7C420FE9, 0xF8841EC9, 0x00000000, 0x09808683, 0x322BED48, 73 0x1E1170AC, 0x6C5A724E, 0xFD0EFFFB, 0x0F853856, 0x3DAED51E, 0x362D3927, 74 0x0A0FD964, 0x685CA621, 0x9B5B54D1, 0x24362E3A, 0x0C0A67B1, 0x9357E70F, 75 0xB4EE96D2, 0x1B9B919E, 0x80C0C54F, 0x61DC20A2, 0x5A774B69, 0x1C121A16, 76 0xE293BA0A, 0xC0A02AE5, 0x3C22E043, 0x121B171D, 0x0E090D0B, 0xF28BC7AD, 77 0x2DB6A8B9, 0x141EA9C8, 0x57F11985, 0xAF75074C, 0xEE99DDBB, 0xA37F60FD, 78 0xF701269F, 0x5C72F5BC, 0x44663BC5, 0x5BFB7E34, 0x8B432976, 0xCB23C6DC, 79 0xB6EDFC68, 0xB8E4F163, 0xD731DCCA, 0x42638510, 0x13972240, 0x84C61120, 80 0x854A247D, 0xD2BB3DF8, 0xAEF93211, 0xC729A16D, 0x1D9E2F4B, 0xDCB230F3, 81 0x0D8652EC, 0x77C1E3D0, 0x2BB3166C, 0xA970B999, 0x119448FA, 0x47E96422, 82 0xA8FC8CC4, 0xA0F03F1A, 0x567D2CD8, 0x223390EF, 0x87494EC7, 0xD938D1C1, 83 0x8CCAA2FE, 0x98D40B36, 0xA6F581CF, 0xA57ADE28, 0xDAB78E26, 0x3FADBFA4, 84 0x2C3A9DE4, 0x5078920D, 0x6A5FCC9B, 0x547E4662, 0xF68D13C2, 0x90D8B8E8, 85 0x2E39F75E, 0x82C3AFF5, 0x9F5D80BE, 0x69D0937C, 0x6FD52DA9, 0xCF2512B3, 86 0xC8AC993B, 0x10187DA7, 0xE89C636E, 0xDB3BBB7B, 0xCD267809, 0x6E5918F4, 87 0xEC9AB701, 0x834F9AA8, 0xE6956E65, 0xAAFFE67E, 0x21BCCF08, 0xEF15E8E6, 88 0xBAE79BD9, 0x4A6F36CE, 0xEA9F09D4, 0x29B07CD6, 0x31A4B2AF, 0x2A3F2331, 89 0xC6A59430, 0x35A266C0, 0x744EBC37, 0xFC82CAA6, 0xE090D0B0, 0x33A7D815, 90 0xF104984A, 0x41ECDAF7, 0x7FCD500E, 0x1791F62F, 0x764DD68D, 0x43EFB04D, 91 0xCCAA4D54, 0xE49604DF, 0x9ED1B5E3, 0x4C6A881B, 0xC12C1FB8, 0x4665517F, 92 0x9D5EEA04, 0x018C355D, 0xFA877473, 0xFB0B412E, 0xB3671D5A, 0x92DBD252, 93 0xE9105633, 0x6DD64713, 0x9AD7618C, 0x37A10C7A, 0x59F8148E, 0xEB133C89, 94 0xCEA927EE, 0xB761C935, 0xE11CE5ED, 0x7A47B13C, 0x9CD2DF59, 0x55F2733F, 95 0x1814CE79, 0x73C737BF, 0x53F7CDEA, 0x5FFDAA5B, 0xDF3D6F14, 0x7844DB86, 96 0xCAAFF381, 0xB968C43E, 0x3824342C, 0xC2A3405F, 0x161DC372, 0xBCE2250C, 97 0x283C498B, 0xFF0D9541, 0x39A80171, 0x080CB3DE, 0xD8B4E49C, 0x6456C190, 98 0x7BCB8461, 0xD532B670, 0x486C5C74, 0xD0B85742 99 }; 100 101 static unsigned 102 mul2(unsigned x) 103 { 104 x <<= 1; 105 return x ^ ((unsigned)(-(int)(x >> 8)) & 0x11B); 106 } 107 108 static unsigned 109 mul9(unsigned x) 110 { 111 return x ^ mul2(mul2(mul2(x))); 112 } 113 114 static unsigned 115 mulb(unsigned x) 116 { 117 unsigned x2; 118 119 x2 = mul2(x); 120 return x ^ x2 ^ mul2(mul2(x2)); 121 } 122 123 static unsigned 124 muld(unsigned x) 125 { 126 unsigned x4; 127 128 x4 = mul2(mul2(x)); 129 return x ^ x4 ^ mul2(x4); 130 } 131 132 static unsigned 133 mule(unsigned x) 134 { 135 unsigned x2, x4; 136 137 x2 = mul2(x); 138 x4 = mul2(x2); 139 return x2 ^ x4 ^ mul2(x4); 140 } 141 142 /* see inner.h */ 143 unsigned 144 br_aes_big_keysched_inv(uint32_t *skey, const void *key, size_t key_len) 145 { 146 unsigned num_rounds; 147 int i, m; 148 149 /* 150 * Sub-keys for decryption are distinct from encryption sub-keys 151 * in that InvMixColumns() is already applied for the inner 152 * rounds. 153 */ 154 num_rounds = br_aes_keysched(skey, key, key_len); 155 m = (int)(num_rounds << 2); 156 for (i = 4; i < m; i ++) { 157 uint32_t p; 158 unsigned p0, p1, p2, p3; 159 uint32_t q0, q1, q2, q3; 160 161 p = skey[i]; 162 p0 = p >> 24; 163 p1 = (p >> 16) & 0xFF; 164 p2 = (p >> 8) & 0xFF; 165 p3 = p & 0xFF; 166 q0 = mule(p0) ^ mulb(p1) ^ muld(p2) ^ mul9(p3); 167 q1 = mul9(p0) ^ mule(p1) ^ mulb(p2) ^ muld(p3); 168 q2 = muld(p0) ^ mul9(p1) ^ mule(p2) ^ mulb(p3); 169 q3 = mulb(p0) ^ muld(p1) ^ mul9(p2) ^ mule(p3); 170 skey[i] = (q0 << 24) | (q1 << 16) | (q2 << 8) | q3; 171 } 172 return num_rounds; 173 } 174 175 static inline uint32_t 176 rotr(uint32_t x, int n) 177 { 178 return (x << (32 - n)) | (x >> n); 179 } 180 181 #define iSboxExt0(x) (iSsm0[x]) 182 #define iSboxExt1(x) (rotr(iSsm0[x], 8)) 183 #define iSboxExt2(x) (rotr(iSsm0[x], 16)) 184 #define iSboxExt3(x) (rotr(iSsm0[x], 24)) 185 186 /* see bearssl.h */ 187 void 188 br_aes_big_decrypt(unsigned num_rounds, const uint32_t *skey, void *data) 189 { 190 unsigned char *buf; 191 uint32_t s0, s1, s2, s3; 192 uint32_t t0, t1, t2, t3; 193 unsigned u; 194 195 buf = data; 196 s0 = br_dec32be(buf); 197 s1 = br_dec32be(buf + 4); 198 s2 = br_dec32be(buf + 8); 199 s3 = br_dec32be(buf + 12); 200 s0 ^= skey[(num_rounds << 2) + 0]; 201 s1 ^= skey[(num_rounds << 2) + 1]; 202 s2 ^= skey[(num_rounds << 2) + 2]; 203 s3 ^= skey[(num_rounds << 2) + 3]; 204 for (u = num_rounds - 1; u > 0; u --) { 205 uint32_t v0 = iSboxExt0(s0 >> 24) 206 ^ iSboxExt1((s3 >> 16) & 0xFF) 207 ^ iSboxExt2((s2 >> 8) & 0xFF) 208 ^ iSboxExt3(s1 & 0xFF); 209 uint32_t v1 = iSboxExt0(s1 >> 24) 210 ^ iSboxExt1((s0 >> 16) & 0xFF) 211 ^ iSboxExt2((s3 >> 8) & 0xFF) 212 ^ iSboxExt3(s2 & 0xFF); 213 uint32_t v2 = iSboxExt0(s2 >> 24) 214 ^ iSboxExt1((s1 >> 16) & 0xFF) 215 ^ iSboxExt2((s0 >> 8) & 0xFF) 216 ^ iSboxExt3(s3 & 0xFF); 217 uint32_t v3 = iSboxExt0(s3 >> 24) 218 ^ iSboxExt1((s2 >> 16) & 0xFF) 219 ^ iSboxExt2((s1 >> 8) & 0xFF) 220 ^ iSboxExt3(s0 & 0xFF); 221 s0 = v0; 222 s1 = v1; 223 s2 = v2; 224 s3 = v3; 225 s0 ^= skey[u << 2]; 226 s1 ^= skey[(u << 2) + 1]; 227 s2 ^= skey[(u << 2) + 2]; 228 s3 ^= skey[(u << 2) + 3]; 229 } 230 t0 = ((uint32_t)iS[s0 >> 24] << 24) 231 | ((uint32_t)iS[(s3 >> 16) & 0xFF] << 16) 232 | ((uint32_t)iS[(s2 >> 8) & 0xFF] << 8) 233 | (uint32_t)iS[s1 & 0xFF]; 234 t1 = ((uint32_t)iS[s1 >> 24] << 24) 235 | ((uint32_t)iS[(s0 >> 16) & 0xFF] << 16) 236 | ((uint32_t)iS[(s3 >> 8) & 0xFF] << 8) 237 | (uint32_t)iS[s2 & 0xFF]; 238 t2 = ((uint32_t)iS[s2 >> 24] << 24) 239 | ((uint32_t)iS[(s1 >> 16) & 0xFF] << 16) 240 | ((uint32_t)iS[(s0 >> 8) & 0xFF] << 8) 241 | (uint32_t)iS[s3 & 0xFF]; 242 t3 = ((uint32_t)iS[s3 >> 24] << 24) 243 | ((uint32_t)iS[(s2 >> 16) & 0xFF] << 16) 244 | ((uint32_t)iS[(s1 >> 8) & 0xFF] << 8) 245 | (uint32_t)iS[s0 & 0xFF]; 246 s0 = t0 ^ skey[0]; 247 s1 = t1 ^ skey[1]; 248 s2 = t2 ^ skey[2]; 249 s3 = t3 ^ skey[3]; 250 br_enc32be(buf, s0); 251 br_enc32be(buf + 4, s1); 252 br_enc32be(buf + 8, s2); 253 br_enc32be(buf + 12, s3); 254 } 255