Botan  2.19.1
Crypto and TLS for C++11
sm3.cpp
Go to the documentation of this file.
1 /*
2 * SM3
3 * (C) 2017 Ribose Inc.
4 *
5 * Botan is released under the Simplified BSD License (see license.txt)
6 */
7 
8 #include <botan/sm3.h>
9 #include <botan/loadstor.h>
10 #include <botan/rotate.h>
11 
12 namespace Botan {
13 
14 std::unique_ptr<HashFunction> SM3::copy_state() const
15  {
16  return std::unique_ptr<HashFunction>(new SM3(*this));
17  }
18 
19 namespace {
20 
21 const uint32_t SM3_IV[] = {
22  0x7380166fUL, 0x4914b2b9UL, 0x172442d7UL, 0xda8a0600UL,
23  0xa96f30bcUL, 0x163138aaUL, 0xe38dee4dUL, 0xb0fb0e4eUL
24 };
25 
26 inline uint32_t P0(uint32_t X)
27  {
28  return X ^ rotl<9>(X) ^ rotl<17>(X);
29  }
30 
31 inline uint32_t FF1(uint32_t X, uint32_t Y, uint32_t Z)
32  {
33  return (X & Y) | ((X | Y) & Z);
34  //return (X & Y) | (X & Z) | (Y & Z);
35  }
36 
37 inline uint32_t GG1(uint32_t X, uint32_t Y, uint32_t Z)
38  {
39  //return (X & Y) | (~X & Z);
40  return ((Z ^ (X & (Y ^ Z))));
41  }
42 
43 inline void R1(uint32_t A, uint32_t& B, uint32_t C, uint32_t& D,
44  uint32_t E, uint32_t& F, uint32_t G, uint32_t& H,
45  uint32_t TJ, uint32_t Wi, uint32_t Wj)
46  {
47  const uint32_t A12 = rotl<12>(A);
48  const uint32_t SS1 = rotl<7>(A12 + E + TJ);
49  const uint32_t TT1 = (A ^ B ^ C) + D + (SS1 ^ A12) + Wj;
50  const uint32_t TT2 = (E ^ F ^ G) + H + SS1 + Wi;
51 
52  B = rotl<9>(B);
53  D = TT1;
54  F = rotl<19>(F);
55  H = P0(TT2);
56  }
57 
58 inline void R2(uint32_t A, uint32_t& B, uint32_t C, uint32_t& D,
59  uint32_t E, uint32_t& F, uint32_t G, uint32_t& H,
60  uint32_t TJ, uint32_t Wi, uint32_t Wj)
61  {
62  const uint32_t A12 = rotl<12>(A);
63  const uint32_t SS1 = rotl<7>(A12 + E + TJ);
64  const uint32_t TT1 = FF1(A, B, C) + D + (SS1 ^ A12) + Wj;
65  const uint32_t TT2 = GG1(E, F, G) + H + SS1 + Wi;
66 
67  B = rotl<9>(B);
68  D = TT1;
69  F = rotl<19>(F);
70  H = P0(TT2);
71  }
72 
73 inline uint32_t P1(uint32_t X)
74  {
75  return X ^ rotl<15>(X) ^ rotl<23>(X);
76  }
77 
78 inline uint32_t SM3_E(uint32_t W0, uint32_t W7, uint32_t W13, uint32_t W3, uint32_t W10)
79  {
80  return P1(W0 ^ W7 ^ rotl<15>(W13)) ^ rotl<7>(W3) ^ W10;
81  }
82 
83 }
84 
85 /*
86 * SM3 Compression Function
87 */
88 void SM3::compress_n(const uint8_t input[], size_t blocks)
89  {
90  uint32_t A = m_digest[0], B = m_digest[1], C = m_digest[2], D = m_digest[3],
91  E = m_digest[4], F = m_digest[5], G = m_digest[6], H = m_digest[7];
92 
93  for(size_t i = 0; i != blocks; ++i)
94  {
95  uint32_t W00 = load_be<uint32_t>(input, 0);
96  uint32_t W01 = load_be<uint32_t>(input, 1);
97  uint32_t W02 = load_be<uint32_t>(input, 2);
98  uint32_t W03 = load_be<uint32_t>(input, 3);
99  uint32_t W04 = load_be<uint32_t>(input, 4);
100  uint32_t W05 = load_be<uint32_t>(input, 5);
101  uint32_t W06 = load_be<uint32_t>(input, 6);
102  uint32_t W07 = load_be<uint32_t>(input, 7);
103  uint32_t W08 = load_be<uint32_t>(input, 8);
104  uint32_t W09 = load_be<uint32_t>(input, 9);
105  uint32_t W10 = load_be<uint32_t>(input, 10);
106  uint32_t W11 = load_be<uint32_t>(input, 11);
107  uint32_t W12 = load_be<uint32_t>(input, 12);
108  uint32_t W13 = load_be<uint32_t>(input, 13);
109  uint32_t W14 = load_be<uint32_t>(input, 14);
110  uint32_t W15 = load_be<uint32_t>(input, 15);
111 
112  R1(A, B, C, D, E, F, G, H, 0x79CC4519, W00, W00 ^ W04);
113  W00 = SM3_E(W00, W07, W13, W03, W10);
114  R1(D, A, B, C, H, E, F, G, 0xF3988A32, W01, W01 ^ W05);
115  W01 = SM3_E(W01, W08, W14, W04, W11);
116  R1(C, D, A, B, G, H, E, F, 0xE7311465, W02, W02 ^ W06);
117  W02 = SM3_E(W02, W09, W15, W05, W12);
118  R1(B, C, D, A, F, G, H, E, 0xCE6228CB, W03, W03 ^ W07);
119  W03 = SM3_E(W03, W10, W00, W06, W13);
120  R1(A, B, C, D, E, F, G, H, 0x9CC45197, W04, W04 ^ W08);
121  W04 = SM3_E(W04, W11, W01, W07, W14);
122  R1(D, A, B, C, H, E, F, G, 0x3988A32F, W05, W05 ^ W09);
123  W05 = SM3_E(W05, W12, W02, W08, W15);
124  R1(C, D, A, B, G, H, E, F, 0x7311465E, W06, W06 ^ W10);
125  W06 = SM3_E(W06, W13, W03, W09, W00);
126  R1(B, C, D, A, F, G, H, E, 0xE6228CBC, W07, W07 ^ W11);
127  W07 = SM3_E(W07, W14, W04, W10, W01);
128  R1(A, B, C, D, E, F, G, H, 0xCC451979, W08, W08 ^ W12);
129  W08 = SM3_E(W08, W15, W05, W11, W02);
130  R1(D, A, B, C, H, E, F, G, 0x988A32F3, W09, W09 ^ W13);
131  W09 = SM3_E(W09, W00, W06, W12, W03);
132  R1(C, D, A, B, G, H, E, F, 0x311465E7, W10, W10 ^ W14);
133  W10 = SM3_E(W10, W01, W07, W13, W04);
134  R1(B, C, D, A, F, G, H, E, 0x6228CBCE, W11, W11 ^ W15);
135  W11 = SM3_E(W11, W02, W08, W14, W05);
136  R1(A, B, C, D, E, F, G, H, 0xC451979C, W12, W12 ^ W00);
137  W12 = SM3_E(W12, W03, W09, W15, W06);
138  R1(D, A, B, C, H, E, F, G, 0x88A32F39, W13, W13 ^ W01);
139  W13 = SM3_E(W13, W04, W10, W00, W07);
140  R1(C, D, A, B, G, H, E, F, 0x11465E73, W14, W14 ^ W02);
141  W14 = SM3_E(W14, W05, W11, W01, W08);
142  R1(B, C, D, A, F, G, H, E, 0x228CBCE6, W15, W15 ^ W03);
143  W15 = SM3_E(W15, W06, W12, W02, W09);
144  R2(A, B, C, D, E, F, G, H, 0x9D8A7A87, W00, W00 ^ W04);
145  W00 = SM3_E(W00, W07, W13, W03, W10);
146  R2(D, A, B, C, H, E, F, G, 0x3B14F50F, W01, W01 ^ W05);
147  W01 = SM3_E(W01, W08, W14, W04, W11);
148  R2(C, D, A, B, G, H, E, F, 0x7629EA1E, W02, W02 ^ W06);
149  W02 = SM3_E(W02, W09, W15, W05, W12);
150  R2(B, C, D, A, F, G, H, E, 0xEC53D43C, W03, W03 ^ W07);
151  W03 = SM3_E(W03, W10, W00, W06, W13);
152  R2(A, B, C, D, E, F, G, H, 0xD8A7A879, W04, W04 ^ W08);
153  W04 = SM3_E(W04, W11, W01, W07, W14);
154  R2(D, A, B, C, H, E, F, G, 0xB14F50F3, W05, W05 ^ W09);
155  W05 = SM3_E(W05, W12, W02, W08, W15);
156  R2(C, D, A, B, G, H, E, F, 0x629EA1E7, W06, W06 ^ W10);
157  W06 = SM3_E(W06, W13, W03, W09, W00);
158  R2(B, C, D, A, F, G, H, E, 0xC53D43CE, W07, W07 ^ W11);
159  W07 = SM3_E(W07, W14, W04, W10, W01);
160  R2(A, B, C, D, E, F, G, H, 0x8A7A879D, W08, W08 ^ W12);
161  W08 = SM3_E(W08, W15, W05, W11, W02);
162  R2(D, A, B, C, H, E, F, G, 0x14F50F3B, W09, W09 ^ W13);
163  W09 = SM3_E(W09, W00, W06, W12, W03);
164  R2(C, D, A, B, G, H, E, F, 0x29EA1E76, W10, W10 ^ W14);
165  W10 = SM3_E(W10, W01, W07, W13, W04);
166  R2(B, C, D, A, F, G, H, E, 0x53D43CEC, W11, W11 ^ W15);
167  W11 = SM3_E(W11, W02, W08, W14, W05);
168  R2(A, B, C, D, E, F, G, H, 0xA7A879D8, W12, W12 ^ W00);
169  W12 = SM3_E(W12, W03, W09, W15, W06);
170  R2(D, A, B, C, H, E, F, G, 0x4F50F3B1, W13, W13 ^ W01);
171  W13 = SM3_E(W13, W04, W10, W00, W07);
172  R2(C, D, A, B, G, H, E, F, 0x9EA1E762, W14, W14 ^ W02);
173  W14 = SM3_E(W14, W05, W11, W01, W08);
174  R2(B, C, D, A, F, G, H, E, 0x3D43CEC5, W15, W15 ^ W03);
175  W15 = SM3_E(W15, W06, W12, W02, W09);
176  R2(A, B, C, D, E, F, G, H, 0x7A879D8A, W00, W00 ^ W04);
177  W00 = SM3_E(W00, W07, W13, W03, W10);
178  R2(D, A, B, C, H, E, F, G, 0xF50F3B14, W01, W01 ^ W05);
179  W01 = SM3_E(W01, W08, W14, W04, W11);
180  R2(C, D, A, B, G, H, E, F, 0xEA1E7629, W02, W02 ^ W06);
181  W02 = SM3_E(W02, W09, W15, W05, W12);
182  R2(B, C, D, A, F, G, H, E, 0xD43CEC53, W03, W03 ^ W07);
183  W03 = SM3_E(W03, W10, W00, W06, W13);
184  R2(A, B, C, D, E, F, G, H, 0xA879D8A7, W04, W04 ^ W08);
185  W04 = SM3_E(W04, W11, W01, W07, W14);
186  R2(D, A, B, C, H, E, F, G, 0x50F3B14F, W05, W05 ^ W09);
187  W05 = SM3_E(W05, W12, W02, W08, W15);
188  R2(C, D, A, B, G, H, E, F, 0xA1E7629E, W06, W06 ^ W10);
189  W06 = SM3_E(W06, W13, W03, W09, W00);
190  R2(B, C, D, A, F, G, H, E, 0x43CEC53D, W07, W07 ^ W11);
191  W07 = SM3_E(W07, W14, W04, W10, W01);
192  R2(A, B, C, D, E, F, G, H, 0x879D8A7A, W08, W08 ^ W12);
193  W08 = SM3_E(W08, W15, W05, W11, W02);
194  R2(D, A, B, C, H, E, F, G, 0x0F3B14F5, W09, W09 ^ W13);
195  W09 = SM3_E(W09, W00, W06, W12, W03);
196  R2(C, D, A, B, G, H, E, F, 0x1E7629EA, W10, W10 ^ W14);
197  W10 = SM3_E(W10, W01, W07, W13, W04);
198  R2(B, C, D, A, F, G, H, E, 0x3CEC53D4, W11, W11 ^ W15);
199  W11 = SM3_E(W11, W02, W08, W14, W05);
200  R2(A, B, C, D, E, F, G, H, 0x79D8A7A8, W12, W12 ^ W00);
201  W12 = SM3_E(W12, W03, W09, W15, W06);
202  R2(D, A, B, C, H, E, F, G, 0xF3B14F50, W13, W13 ^ W01);
203  W13 = SM3_E(W13, W04, W10, W00, W07);
204  R2(C, D, A, B, G, H, E, F, 0xE7629EA1, W14, W14 ^ W02);
205  W14 = SM3_E(W14, W05, W11, W01, W08);
206  R2(B, C, D, A, F, G, H, E, 0xCEC53D43, W15, W15 ^ W03);
207  W15 = SM3_E(W15, W06, W12, W02, W09);
208  R2(A, B, C, D, E, F, G, H, 0x9D8A7A87, W00, W00 ^ W04);
209  W00 = SM3_E(W00, W07, W13, W03, W10);
210  R2(D, A, B, C, H, E, F, G, 0x3B14F50F, W01, W01 ^ W05);
211  W01 = SM3_E(W01, W08, W14, W04, W11);
212  R2(C, D, A, B, G, H, E, F, 0x7629EA1E, W02, W02 ^ W06);
213  W02 = SM3_E(W02, W09, W15, W05, W12);
214  R2(B, C, D, A, F, G, H, E, 0xEC53D43C, W03, W03 ^ W07);
215  W03 = SM3_E(W03, W10, W00, W06, W13);
216  R2(A, B, C, D, E, F, G, H, 0xD8A7A879, W04, W04 ^ W08);
217  R2(D, A, B, C, H, E, F, G, 0xB14F50F3, W05, W05 ^ W09);
218  R2(C, D, A, B, G, H, E, F, 0x629EA1E7, W06, W06 ^ W10);
219  R2(B, C, D, A, F, G, H, E, 0xC53D43CE, W07, W07 ^ W11);
220  R2(A, B, C, D, E, F, G, H, 0x8A7A879D, W08, W08 ^ W12);
221  R2(D, A, B, C, H, E, F, G, 0x14F50F3B, W09, W09 ^ W13);
222  R2(C, D, A, B, G, H, E, F, 0x29EA1E76, W10, W10 ^ W14);
223  R2(B, C, D, A, F, G, H, E, 0x53D43CEC, W11, W11 ^ W15);
224  R2(A, B, C, D, E, F, G, H, 0xA7A879D8, W12, W12 ^ W00);
225  R2(D, A, B, C, H, E, F, G, 0x4F50F3B1, W13, W13 ^ W01);
226  R2(C, D, A, B, G, H, E, F, 0x9EA1E762, W14, W14 ^ W02);
227  R2(B, C, D, A, F, G, H, E, 0x3D43CEC5, W15, W15 ^ W03);
228 
229  A = (m_digest[0] ^= A);
230  B = (m_digest[1] ^= B);
231  C = (m_digest[2] ^= C);
232  D = (m_digest[3] ^= D);
233  E = (m_digest[4] ^= E);
234  F = (m_digest[5] ^= F);
235  G = (m_digest[6] ^= G);
236  H = (m_digest[7] ^= H);
237 
238  input += hash_block_size();
239  }
240  }
241 
242 /*
243 * Copy out the digest
244 */
245 void SM3::copy_out(uint8_t output[])
246  {
247  copy_out_vec_be(output, output_length(), m_digest);
248  }
249 
250 /*
251 * Clear memory of sensitive data
252 */
254  {
256  std::copy(std::begin(SM3_IV), std::end(SM3_IV), m_digest.begin());
257  }
258 
259 }
size_t output_length() const override
Definition: sm3.h:29
fe X
Definition: ge.cpp:27
SIMD_8x32 H
void copy_out_vec_be(uint8_t out[], size_t out_bytes, const std::vector< T, Alloc > &in)
Definition: loadstor.h:673
void clear() override
Definition: mdx_hash.cpp:41
void BOTAN_FUNC_ISA("avx2") SHACAL2 SIMD_8x32 A
uint32_t load_be< uint32_t >(const uint8_t in[], size_t off)
Definition: loadstor.h:179
fe Y
Definition: ge.cpp:28
SIMD_8x32 D
size_t hash_block_size() const overridefinal
Definition: mdx_hash.h:35
SIMD_8x32 E
SIMD_8x32 G
std::unique_ptr< HashFunction > copy_state() const override
Definition: sm3.cpp:14
SIMD_8x32 B
Definition: alg_id.cpp:13
SIMD_8x32 F
void clear() override
Definition: sm3.cpp:253
SIMD_8x32 C
fe Z
Definition: ge.cpp:29
SM3()
Definition: sm3.h:35