-void keccak_theta(uint64_t *a, uint64_t *b);
-
-static inline
-void keccak_round(uint64_t a[5][5], uint8_t rci){
- uint64_t b[5][5];
- uint8_t i, j;
- union {
- uint64_t v64;
- uint8_t v8[8];
- } t;
- /* theta */
-#if DEBUG & 0
- printf_P(PSTR("\npre-theta(1):"));
- keccak_dump_state(b);
-#endif
- keccak_theta((uint64_t*)a, (uint64_t*)b);
-#if DEBUG
- printf_P(PSTR("\npost-theta(1):"));
- keccak_dump_state(b);
-#endif
- for(i = 0; i < 5; ++i){
- for(j = 0; j < 5; ++j){
- a[j][i] ^= rotate64_1bit_left(b[(i + 1) % 5][0]);
- }
- }
-#if DEBUG
- cli_putstr_P(PSTR("\r\nAfter theta:"));
- keccak_dump_state(a);
-#endif
- /* rho & pi */
- for(i = 0; i < 5; ++i){
- for(j = 0; j < 5; ++j){
- b[(2 * i + 3 * j) % 5][j] = rotate64left_code(a[j][i], pgm_read_byte(&(r[i][j])));
- }
- }
-#if DEBUG & 0
- cli_putstr_P(PSTR("\r\n--- after rho & pi ---"));
- keccak_dump_state(a);
-#endif
- /* chi */
- for(i = 0; i < 5; ++i){
- for(j = 0; j < 5; ++j){
- a[j][i] = b[j][i] ^ ((~(b[j][(i + 1) % 5])) & (b[j][(i + 2) % 5]));
- }
- }
-#if DEBUG & 0
- cli_putstr_P(PSTR("\r\nAfter chi:"));
- keccak_dump_state(a);
-#endif
- /* iota */
-
-// memcpy_P(&t, &(rc_comp[rci]), 8);
- t.v64 = 0;
- t.v8[0] = pgm_read_byte(&(rc_comp[rci]));
- if(t.v8[0] & 0x40){
- t.v8[7] = 0x80;
- }
- if(t.v8[0] & 0x20){
- t.v8[3] = 0x80;
- }
- if(t.v8[0] & 0x10){
- t.v8[1] = 0x80;
- }
- t.v8[0] &= 0x8F;
-
- a[0][0] ^= t.v64;
-#if DEBUG & 0
- cli_putstr_P(PSTR("\r\nAfter iota:"));
- keccak_dump_state(a);
-#endif
-}
-
-void keccak_f1600(uint64_t a[5][5]){
- uint8_t i = 0;
- do {
-#if DEBUG
- cli_putstr_P(PSTR("\r\n\r\n--- Round "));
- cli_hexdump(&i, 1);
- cli_putstr_P(PSTR(" ---"));
-#endif
- keccak_round(a, i);
- } while (++i < 24);
-}