updated Grøstl for round 3

[avr-crypto-lib.git] / groestl / groestl_large.c

diff --git a/groestl/groestl_large.c b/groestl/groestl_large.c
index bcfd8e51ae67c81f807c466b21d7f20d7db38314..0812683357c2f33fb57307b701e4e5b0629a0c7a 100644 (file)
--- a/groestl/groestl_large.c
+++ b/groestl/groestl_large.c
@@ -65,14 +65,49 @@ static uint8_t matrix[] PROGMEM = {
  2, 3, 4, 5, 3, 5, 7, 2
 };
 
+static
+void shift_columns(uint8_t *a, PGM_VOID_P shifts){
+       uint8_t s;
+       uint8_t tmp[16];
+       uint8_t i,j;
+       for(i=0; i<8; ++i){
+               s = pgm_read_byte(shifts);
+               shifts = (uint8_t*)shifts + 1;
+               if(s==0){
+                       continue;
+               }
+               for(j=0;j<16;++j){
+                       tmp[j] = a[i+j*8];
+               }
+               for(j=0; j<16; ++j){
+                       a[i+((j-s+16)%16)*8] = tmp[j];
+               }
+       }
+}
+
+static const uint8_t p_shifts[] PROGMEM = {
+               0, 1, 2, 3, 4, 5, 6, 11
+};
+
+static const uint8_t q_shifts[] PROGMEM = {
+               1, 3, 5, 11, 0, 2, 4, 6
+};
+
 void groestl_large_rounds(uint8_t *m, uint8_t q){
        uint8_t r,i,j;
        uint8_t tmp[16];
        for(r=0; r<ROUNDS; ++r){
                if(q){
-                       m[7] ^= 0xff ^ r;
+                       for(i=0; i<(8*16); ++i){
+                               m[i] ^= 0xff;
+                       }
+                       for(i=0; i<16; ++i){
+                               m[7+i*8] ^= r ^ (i<<4);
+                       }
                }else{
-                       m[0] ^= r;
+                       for(i=0; i<16; ++i){
+                               m[i*8] ^= r ^ (i<<4);
+                       }
                }       
 #if DEBUG              
                if(r<2){
@@ -83,17 +118,10 @@ void groestl_large_rounds(uint8_t *m, uint8_t q){
                for(i=0;i<16*8; ++i){
                        m[i] = pgm_read_byte(aes_sbox+m[i]);
                }
-               for(i=1; i<7; ++i){
-                       for(j=0; j<16; ++j)
-                               tmp[j] = m[i+8*j];
-                       for(j=0; j<16; ++j){
-                               m[i+((j-i+16)%16)*8] = tmp[j];
-                       }
-               }
-               for(j=0; j<16; ++j)
-                       tmp[j] = m[7+8*j];
-               for(j=0; j<16; ++j){
-                       m[7+((j-11+16)%16)*8] = tmp[j];
+               if(!q){
+                       shift_columns(m, p_shifts);
+               }else{
+                       shift_columns(m, q_shifts);
                }
                
 #if DEBUG