/src/botan/src/lib/block/des/des.cpp

Source (jump to first uncovered line)
/*
* DES
* (C) 1999-2008,2018,2020 Jack Lloyd
*
* Based on a public domain implemenation by Phil Karn (who in turn
* credited Richard Outerbridge and Jim Gillogly)
*
* Botan is released under the Simplified BSD License (see license.txt)
*/

#include <botan/internal/des.h>
#include <botan/internal/loadstor.h>
#include <botan/internal/rotate.h>

namespace Botan {

namespace {

alignas(256) const uint8_t SPBOX_CATS[64*8] = {
   0x54, 0x00, 0x10, 0x55, 0x51, 0x15, 0x01, 0x10, 0x04, 0x54, 0x55, 0x04, 0x45, 0x51, 0x40, 0x01,
   0x05, 0x44, 0x44, 0x14, 0x14, 0x50, 0x50, 0x45, 0x11, 0x41, 0x41, 0x11, 0x00, 0x05, 0x15, 0x40,
   0x10, 0x55, 0x01, 0x50, 0x54, 0x40, 0x40, 0x04, 0x51, 0x10, 0x14, 0x41, 0x04, 0x01, 0x45, 0x15,
   0x55, 0x11, 0x50, 0x45, 0x41, 0x05, 0x15, 0x54, 0x05, 0x44, 0x44, 0x00, 0x11, 0x14, 0x00, 0x51,

   0x55, 0x44, 0x04, 0x15, 0x10, 0x01, 0x51, 0x45, 0x41, 0x55, 0x54, 0x40, 0x44, 0x10, 0x01, 0x51,
   0x14, 0x11, 0x45, 0x00, 0x40, 0x04, 0x15, 0x50, 0x11, 0x41, 0x00, 0x14, 0x05, 0x54, 0x50, 0x05,
   0x00, 0x15, 0x51, 0x10, 0x45, 0x50, 0x54, 0x04, 0x50, 0x44, 0x01, 0x55, 0x15, 0x01, 0x04, 0x40,
   0x05, 0x54, 0x10, 0x41, 0x11, 0x45, 0x41, 0x11, 0x14, 0x00, 0x44, 0x05, 0x40, 0x51, 0x55, 0x14,

   0x09, 0xA8, 0x00, 0xA1, 0x88, 0x00, 0x29, 0x88, 0x21, 0x81, 0x81, 0x20, 0xA9, 0x21, 0xA0, 0x09,
   0x80, 0x01, 0xA8, 0x08, 0x28, 0xA0, 0xA1, 0x29, 0x89, 0x28, 0x20, 0x89, 0x01, 0xA9, 0x08, 0x80,
   0xA8, 0x80, 0x21, 0x09, 0x20, 0xA8, 0x88, 0x00, 0x08, 0x21, 0xA9, 0x88, 0x81, 0x08, 0x00, 0xA1,
   0x89, 0x20, 0x80, 0xA9, 0x01, 0x29, 0x28, 0x81, 0xA0, 0x89, 0x09, 0xA0, 0x29, 0x01, 0xA1, 0x28,

   0x51, 0x15, 0x15, 0x04, 0x54, 0x45, 0x41, 0x11, 0x00, 0x50, 0x50, 0x55, 0x05, 0x00, 0x44, 0x41,
   0x01, 0x10, 0x40, 0x51, 0x04, 0x40, 0x11, 0x14, 0x45, 0x01, 0x14, 0x44, 0x10, 0x54, 0x55, 0x05,
   0x44, 0x41, 0x50, 0x55, 0x05, 0x00, 0x00, 0x50, 0x14, 0x44, 0x45, 0x01, 0x51, 0x15, 0x15, 0x04,
   0x55, 0x05, 0x01, 0x10, 0x41, 0x11, 0x54, 0x45, 0x11, 0x14, 0x40, 0x51, 0x04, 0x40, 0x10, 0x54,

   0x01, 0x29, 0x28, 0xA1, 0x08, 0x01, 0x80, 0x28, 0x89, 0x08, 0x21, 0x89, 0xA1, 0xA8, 0x09, 0x80,
   0x20, 0x88, 0x88, 0x00, 0x81, 0xA9, 0xA9, 0x21, 0xA8, 0x81, 0x00, 0xA0, 0x29, 0x20, 0xA0, 0x09,
   0x08, 0xA1, 0x01, 0x20, 0x80, 0x28, 0xA1, 0x89, 0x21, 0x80, 0xA8, 0x29, 0x89, 0x01, 0x20, 0xA8,
   0xA9, 0x09, 0xA0, 0xA9, 0x28, 0x00, 0x88, 0xA0, 0x09, 0x21, 0x81, 0x08, 0x00, 0x88, 0x29, 0x81,

   0x41, 0x50, 0x04, 0x55, 0x50, 0x01, 0x55, 0x10, 0x44, 0x15, 0x10, 0x41, 0x11, 0x44, 0x40, 0x05,
   0x00, 0x11, 0x45, 0x04, 0x14, 0x45, 0x01, 0x51, 0x51, 0x00, 0x15, 0x54, 0x05, 0x14, 0x54, 0x40,
   0x44, 0x01, 0x51, 0x14, 0x55, 0x10, 0x05, 0x41, 0x10, 0x44, 0x40, 0x05, 0x41, 0x55, 0x14, 0x50,
   0x15, 0x54, 0x00, 0x51, 0x01, 0x04, 0x50, 0x15, 0x04, 0x11, 0x45, 0x00, 0x54, 0x40, 0x11, 0x45,

   0x10, 0x51, 0x45, 0x00, 0x04, 0x45, 0x15, 0x54, 0x55, 0x10, 0x00, 0x41, 0x01, 0x40, 0x51, 0x05,
   0x44, 0x15, 0x11, 0x44, 0x41, 0x50, 0x54, 0x11, 0x50, 0x04, 0x05, 0x55, 0x14, 0x01, 0x40, 0x14,
   0x40, 0x14, 0x10, 0x45, 0x45, 0x51, 0x51, 0x01, 0x11, 0x40, 0x44, 0x10, 0x54, 0x05, 0x15, 0x54,
   0x05, 0x41, 0x55, 0x50, 0x14, 0x00, 0x01, 0x55, 0x00, 0x15, 0x50, 0x04, 0x41, 0x44, 0x04, 0x11,

   0x89, 0x08, 0x20, 0xA9, 0x80, 0x89, 0x01, 0x80, 0x21, 0xA0, 0xA9, 0x28, 0xA8, 0x29, 0x08, 0x01,
   0xA0, 0x81, 0x88, 0x09, 0x28, 0x21, 0xA1, 0xA8, 0x09, 0x00, 0x00, 0xA1, 0x81, 0x88, 0x29, 0x20,
   0x29, 0x20, 0xa8, 0x08, 0x01, 0xA1, 0x08, 0x29, 0x88, 0x01, 0x81, 0xA0, 0xA1, 0x80, 0x20, 0x89,
   0x00, 0xA9, 0x21, 0x81, 0xA0, 0x88, 0x89, 0x00, 0xA9, 0x28, 0x28, 0x09, 0x09, 0x21, 0x80, 0xA8,
};

const uint32_t SPBOX_CAT_0_MUL = 0x70041106;
const uint32_t SPBOX_CAT_1_MUL = 0x02012020;
const uint32_t SPBOX_CAT_2_MUL = 0x00901048;
const uint32_t SPBOX_CAT_3_MUL = 0x8e060221;
const uint32_t SPBOX_CAT_4_MUL = 0x00912140;
const uint32_t SPBOX_CAT_5_MUL = 0x80841018;
const uint32_t SPBOX_CAT_6_MUL = 0xe0120202;
const uint32_t SPBOX_CAT_7_MUL = 0x00212240;

const uint32_t SPBOX_CAT_0_MASK = 0x01010404;
const uint32_t SPBOX_CAT_1_MASK = 0x80108020;
const uint32_t SPBOX_CAT_2_MASK = 0x08020208;
const uint32_t SPBOX_CAT_3_MASK = 0x00802081;
const uint32_t SPBOX_CAT_4_MASK = 0x42080100;
const uint32_t SPBOX_CAT_5_MASK = 0x20404010;
const uint32_t SPBOX_CAT_6_MASK = 0x04200802;
const uint32_t SPBOX_CAT_7_MASK = 0x10041040;

/*
* DES Key Schedule
*/
void des_key_schedule(uint32_t round_key[32], const uint8_t key[8])
   {
   static const uint8_t ROT[16] = { 1, 1, 2, 2, 2, 2, 2, 2,
                                 1, 2, 2, 2, 2, 2, 2, 1 };

   uint32_t C = ((key[7] & 0x80) << 20) | ((key[6] & 0x80) << 19) |
                ((key[5] & 0x80) << 18) | ((key[4] & 0x80) << 17) |
                ((key[3] & 0x80) << 16) | ((key[2] & 0x80) << 15) |
                ((key[1] & 0x80) << 14) | ((key[0] & 0x80) << 13) |
                ((key[7] & 0x40) << 13) | ((key[6] & 0x40) << 12) |
                ((key[5] & 0x40) << 11) | ((key[4] & 0x40) << 10) |
                ((key[3] & 0x40) <<  9) | ((key[2] & 0x40) <<  8) |
                ((key[1] & 0x40) <<  7) | ((key[0] & 0x40) <<  6) |
                ((key[7] & 0x20) <<  6) | ((key[6] & 0x20) <<  5) |
                ((key[5] & 0x20) <<  4) | ((key[4] & 0x20) <<  3) |
                ((key[3] & 0x20) <<  2) | ((key[2] & 0x20) <<  1) |
                ((key[1] & 0x20)      ) | ((key[0] & 0x20) >>  1) |
                ((key[7] & 0x10) >>  1) | ((key[6] & 0x10) >>  2) |
                ((key[5] & 0x10) >>  3) | ((key[4] & 0x10) >>  4);
   uint32_t D = ((key[7] & 0x02) << 26) | ((key[6] & 0x02) << 25) |
                ((key[5] & 0x02) << 24) | ((key[4] & 0x02) << 23) |
                ((key[3] & 0x02) << 22) | ((key[2] & 0x02) << 21) |
                ((key[1] & 0x02) << 20) | ((key[0] & 0x02) << 19) |
                ((key[7] & 0x04) << 17) | ((key[6] & 0x04) << 16) |
                ((key[5] & 0x04) << 15) | ((key[4] & 0x04) << 14) |
                ((key[3] & 0x04) << 13) | ((key[2] & 0x04) << 12) |
                ((key[1] & 0x04) << 11) | ((key[0] & 0x04) << 10) |
                ((key[7] & 0x08) <<  8) | ((key[6] & 0x08) <<  7) |
                ((key[5] & 0x08) <<  6) | ((key[4] & 0x08) <<  5) |
                ((key[3] & 0x08) <<  4) | ((key[2] & 0x08) <<  3) |
                ((key[1] & 0x08) <<  2) | ((key[0] & 0x08) <<  1) |
                ((key[3] & 0x10) >>  1) | ((key[2] & 0x10) >>  2) |
                ((key[1] & 0x10) >>  3) | ((key[0] & 0x10) >>  4);

   for(size_t i = 0; i != 16; ++i)
      {
      C = ((C << ROT[i]) | (C >> (28-ROT[i]))) & 0x0FFFFFFF;
      D = ((D << ROT[i]) | (D >> (28-ROT[i]))) & 0x0FFFFFFF;
      round_key[2*i  ] = ((C & 0x00000010) << 22) | ((C & 0x00000800) << 17) |
                         ((C & 0x00000020) << 16) | ((C & 0x00004004) << 15) |
                         ((C & 0x00000200) << 11) | ((C & 0x00020000) << 10) |
                         ((C & 0x01000000) >>  6) | ((C & 0x00100000) >>  4) |
                         ((C & 0x00010000) <<  3) | ((C & 0x08000000) >>  2) |
                         ((C & 0x00800000) <<  1) | ((D & 0x00000010) <<  8) |
                         ((D & 0x00000002) <<  7) | ((D & 0x00000001) <<  2) |
                         ((D & 0x00000200)      ) | ((D & 0x00008000) >>  2) |
                         ((D & 0x00000088) >>  3) | ((D & 0x00001000) >>  7) |
                         ((D & 0x00080000) >>  9) | ((D & 0x02020000) >> 14) |
                         ((D & 0x00400000) >> 21);
      round_key[2*i+1] = ((C & 0x00000001) << 28) | ((C & 0x00000082) << 18) |
                         ((C & 0x00002000) << 14) | ((C & 0x00000100) << 10) |
                         ((C & 0x00001000) <<  9) | ((C & 0x00040000) <<  6) |
                         ((C & 0x02400000) <<  4) | ((C & 0x00008000) <<  2) |
                         ((C & 0x00200000) >>  1) | ((C & 0x04000000) >> 10) |
                         ((D & 0x00000020) <<  6) | ((D & 0x00000100)      ) |
                         ((D & 0x00000800) >>  1) | ((D & 0x00000040) >>  3) |
                         ((D & 0x00010000) >>  4) | ((D & 0x00000400) >>  5) |
                         ((D & 0x00004000) >> 10) | ((D & 0x04000000) >> 13) |
                         ((D & 0x00800000) >> 14) | ((D & 0x00100000) >> 18) |
                         ((D & 0x01000000) >> 24) | ((D & 0x08000000) >> 26);
      }
   }

inline uint32_t spbox(uint32_t T0, uint32_t T1)
   {
   return
      ((SPBOX_CATS[0*64 + ((T0 >> 24) & 0x3F)] * SPBOX_CAT_0_MUL) & SPBOX_CAT_0_MASK)  ^
      ((SPBOX_CATS[1*64 + ((T1 >> 24) & 0x3F)] * SPBOX_CAT_1_MUL) & SPBOX_CAT_1_MASK)  ^
      ((SPBOX_CATS[2*64 + ((T0 >> 16) & 0x3F)] * SPBOX_CAT_2_MUL) & SPBOX_CAT_2_MASK)  ^
      ((SPBOX_CATS[3*64 + ((T1 >> 16) & 0x3F)] * SPBOX_CAT_3_MUL) & SPBOX_CAT_3_MASK)  ^
      ((SPBOX_CATS[4*64 + ((T0 >>  8) & 0x3F)] * SPBOX_CAT_4_MUL) & SPBOX_CAT_4_MASK)  ^
      ((SPBOX_CATS[5*64 + ((T1 >>  8) & 0x3F)] * SPBOX_CAT_5_MUL) & SPBOX_CAT_5_MASK)  ^
      ((SPBOX_CATS[6*64 + ((T0 >>  0) & 0x3F)] * SPBOX_CAT_6_MUL) & SPBOX_CAT_6_MASK)  ^
      ((SPBOX_CATS[7*64 + ((T1 >>  0) & 0x3F)] * SPBOX_CAT_7_MUL) & SPBOX_CAT_7_MASK);
   }

/*
* DES Encryption
*/
inline void des_encrypt(uint32_t& Lr, uint32_t& Rr,
                        const uint32_t round_key[32])
   {
   uint32_t L = Lr;
   uint32_t R = Rr;
   for(size_t i = 0; i != 16; i += 2)
      {
      L ^= spbox(rotr<4>(R) ^ round_key[2*i  ], R ^ round_key[2*i+1]);
      R ^= spbox(rotr<4>(L) ^ round_key[2*i+2], L ^ round_key[2*i+3]);
      }

   Lr = L;
   Rr = R;
   }

inline void des_encrypt_x2(uint32_t& L0r, uint32_t& R0r,
                           uint32_t& L1r, uint32_t& R1r,
                           const uint32_t round_key[32])
   {
   uint32_t L0 = L0r;
   uint32_t R0 = R0r;
   uint32_t L1 = L1r;
   uint32_t R1 = R1r;

   for(size_t i = 0; i != 16; i += 2)
      {
      L0 ^= spbox(rotr<4>(R0) ^ round_key[2*i  ], R0 ^ round_key[2*i+1]);
      L1 ^= spbox(rotr<4>(R1) ^ round_key[2*i  ], R1 ^ round_key[2*i+1]);

      R0 ^= spbox(rotr<4>(L0) ^ round_key[2*i+2], L0 ^ round_key[2*i+3]);
      R1 ^= spbox(rotr<4>(L1) ^ round_key[2*i+2], L1 ^ round_key[2*i+3]);
      }

   L0r = L0;
   R0r = R0;
   L1r = L1;
   R1r = R1;
   }

/*
* DES Decryption
*/
inline void des_decrypt(uint32_t& Lr, uint32_t& Rr,
                        const uint32_t round_key[32])
   {
   uint32_t L = Lr;
   uint32_t R = Rr;
   for(size_t i = 16; i != 0; i -= 2)
      {
      L ^= spbox(rotr<4>(R) ^ round_key[2*i - 2], R  ^ round_key[2*i - 1]);
      R ^= spbox(rotr<4>(L) ^ round_key[2*i - 4], L  ^ round_key[2*i - 3]);
      }
   Lr = L;
   Rr = R;
   }

inline void des_decrypt_x2(uint32_t& L0r, uint32_t& R0r,
                           uint32_t& L1r, uint32_t& R1r,
                           const uint32_t round_key[32])
   {
   uint32_t L0 = L0r;
   uint32_t R0 = R0r;
   uint32_t L1 = L1r;
   uint32_t R1 = R1r;

   for(size_t i = 16; i != 0; i -= 2)
      {
      L0 ^= spbox(rotr<4>(R0) ^ round_key[2*i - 2], R0  ^ round_key[2*i - 1]);
      L1 ^= spbox(rotr<4>(R1) ^ round_key[2*i - 2], R1  ^ round_key[2*i - 1]);

      R0 ^= spbox(rotr<4>(L0) ^ round_key[2*i - 4], L0  ^ round_key[2*i - 3]);
      R1 ^= spbox(rotr<4>(L1) ^ round_key[2*i - 4], L1  ^ round_key[2*i - 3]);
      }

   L0r = L0;
   R0r = R0;
   L1r = L1;
   R1r = R1;
   }

inline void des_IP(uint32_t& L, uint32_t& R)
   {
   // IP sequence by Wei Dai, taken from public domain Crypto++
   uint32_t T;
   R = rotl<4>(R);
   T = (L ^ R) & 0xF0F0F0F0;
   L ^= T;
   R = rotr<20>(R ^ T);
   T = (L ^ R) & 0xFFFF0000;
   L ^= T;
   R = rotr<18>(R ^ T);
   T = (L ^ R) & 0x33333333;
   L ^= T;
   R = rotr<6>(R ^ T);
   T = (L ^ R) & 0x00FF00FF;
   L ^= T;
   R = rotl<9>(R ^ T);
   T = (L ^ R) & 0xAAAAAAAA;
   L = rotl<1>(L ^ T);
   R ^= T;
   }

inline void des_FP(uint32_t& L, uint32_t& R)
   {
   // FP sequence by Wei Dai, taken from public domain Crypto++
   uint32_t T;

   R = rotr<1>(R);
   T = (L ^ R) & 0xAAAAAAAA;
   R ^= T;
   L = rotr<9>(L ^ T);
   T = (L ^ R) & 0x00FF00FF;
   R ^= T;
   L = rotl<6>(L ^ T);
   T = (L ^ R) & 0x33333333;
   R ^= T;
   L = rotl<18>(L ^ T);
   T = (L ^ R) & 0xFFFF0000;
   R ^= T;
   L = rotl<20>(L ^ T);
   T = (L ^ R) & 0xF0F0F0F0;
   R ^= T;
   L = rotr<4>(L ^ T);
   }

}

/*
* DES Encryption
*/
void DES::encrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
   {
   assert_key_material_set();

   while(blocks >= 2)
      {
      uint32_t L0 = load_be<uint32_t>(in, 0);
      uint32_t R0 = load_be<uint32_t>(in, 1);
      uint32_t L1 = load_be<uint32_t>(in, 2);
      uint32_t R1 = load_be<uint32_t>(in, 3);

      des_IP(L0, R0);
      des_IP(L1, R1);

      des_encrypt_x2(L0, R0, L1, R1, m_round_key.data());

      des_FP(L0, R0);
      des_FP(L1, R1);

      store_be(out, R0, L0, R1, L1);

      in += 2*BLOCK_SIZE;
      out += 2*BLOCK_SIZE;
      blocks -= 2;
      }

   while(blocks > 0)
      {
      uint32_t L0 = load_be<uint32_t>(in, 0);
      uint32_t R0 = load_be<uint32_t>(in, 1);
      des_IP(L0, R0);
      des_encrypt(L0, R0, m_round_key.data());
      des_FP(L0, R0);
      store_be(out, R0, L0);

      in += BLOCK_SIZE;
      out += BLOCK_SIZE;
      blocks -= 1;
      }
   }

/*
* DES Decryption
*/
void DES::decrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
   {
   assert_key_material_set();

   while(blocks >= 2)
      {
      uint32_t L0 = load_be<uint32_t>(in, 0);
      uint32_t R0 = load_be<uint32_t>(in, 1);
      uint32_t L1 = load_be<uint32_t>(in, 2);
      uint32_t R1 = load_be<uint32_t>(in, 3);

      des_IP(L0, R0);
      des_IP(L1, R1);

      des_decrypt_x2(L0, R0, L1, R1, m_round_key.data());

      des_FP(L0, R0);
      des_FP(L1, R1);

      store_be(out, R0, L0, R1, L1);

      in += 2*BLOCK_SIZE;
      out += 2*BLOCK_SIZE;
      blocks -= 2;
      }

   while(blocks > 0)
      {
      uint32_t L0 = load_be<uint32_t>(in, 0);
      uint32_t R0 = load_be<uint32_t>(in, 1);
      des_IP(L0, R0);
      des_decrypt(L0, R0, m_round_key.data());
      des_FP(L0, R0);
      store_be(out, R0, L0);

      in += BLOCK_SIZE;
      out += BLOCK_SIZE;
      blocks -= 1;
      }
   }

bool DES::has_keying_material() const
   {
   return !m_round_key.empty();
   }

/*
* DES Key Schedule
*/
void DES::key_schedule(const uint8_t key[], size_t /*length*/)
   {
   m_round_key.resize(32);
   des_key_schedule(m_round_key.data(), key);
   }

void DES::clear()
   {
   zap(m_round_key);
   }

/*
* TripleDES Encryption
*/
void TripleDES::encrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
   {
   assert_key_material_set();

   while(blocks >= 2)
      {
      uint32_t L0 = load_be<uint32_t>(in, 0);
      uint32_t R0 = load_be<uint32_t>(in, 1);
      uint32_t L1 = load_be<uint32_t>(in, 2);
      uint32_t R1 = load_be<uint32_t>(in, 3);

      des_IP(L0, R0);
      des_IP(L1, R1);

      des_encrypt_x2(L0, R0, L1, R1, &m_round_key[0]);
      des_decrypt_x2(R0, L0, R1, L1, &m_round_key[32]);
      des_encrypt_x2(L0, R0, L1, R1, &m_round_key[64]);

      des_FP(L0, R0);
      des_FP(L1, R1);

      store_be(out, R0, L0, R1, L1);

      in += 2*BLOCK_SIZE;
      out += 2*BLOCK_SIZE;
      blocks -= 2;
      }

   while(blocks > 0)
      {
      uint32_t L0 = load_be<uint32_t>(in, 0);
      uint32_t R0 = load_be<uint32_t>(in, 1);

      des_IP(L0, R0);
      des_encrypt(L0, R0, &m_round_key[0]);
      des_decrypt(R0, L0, &m_round_key[32]);
      des_encrypt(L0, R0, &m_round_key[64]);
      des_FP(L0, R0);

      store_be(out, R0, L0);

      in += BLOCK_SIZE;
      out += BLOCK_SIZE;
      blocks -= 1;
      }
   }

/*
* TripleDES Decryption
*/
void TripleDES::decrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
   {
   assert_key_material_set();

   while(blocks >= 2)
      {
      uint32_t L0 = load_be<uint32_t>(in, 0);
      uint32_t R0 = load_be<uint32_t>(in, 1);
      uint32_t L1 = load_be<uint32_t>(in, 2);
      uint32_t R1 = load_be<uint32_t>(in, 3);

      des_IP(L0, R0);
      des_IP(L1, R1);

      des_decrypt_x2(L0, R0, L1, R1, &m_round_key[64]);
      des_encrypt_x2(R0, L0, R1, L1, &m_round_key[32]);
      des_decrypt_x2(L0, R0, L1, R1, &m_round_key[0]);

      des_FP(L0, R0);
      des_FP(L1, R1);

      store_be(out, R0, L0, R1, L1);

      in += 2*BLOCK_SIZE;
      out += 2*BLOCK_SIZE;
      blocks -= 2;
      }

   while(blocks > 0)
      {
      uint32_t L0 = load_be<uint32_t>(in, 0);
      uint32_t R0 = load_be<uint32_t>(in, 1);

      des_IP(L0, R0);
      des_decrypt(L0, R0, &m_round_key[64]);
      des_encrypt(R0, L0, &m_round_key[32]);
      des_decrypt(L0, R0, &m_round_key[0]);
      des_FP(L0, R0);

      store_be(out, R0, L0);

      in += BLOCK_SIZE;
      out += BLOCK_SIZE;
      blocks -= 1;
      }
   }

bool TripleDES::has_keying_material() const
   {
   return !m_round_key.empty();
   }

/*
* TripleDES Key Schedule
*/
void TripleDES::key_schedule(const uint8_t key[], size_t length)
   {
   m_round_key.resize(3*32);
   des_key_schedule(&m_round_key[0], key);
   des_key_schedule(&m_round_key[32], key + 8);

   if(length == 24)
      des_key_schedule(&m_round_key[64], key + 16);
   else
      copy_mem(&m_round_key[64], &m_round_key[0], 32);
   }

void TripleDES::clear()
   {
   zap(m_round_key);
   }

}

Coverage Report

Created: 2023-02-22 06:39

Line	Count	Source (jump to first uncovered line)
1		/*
2		* DES
3		* (C) 1999-2008,2018,2020 Jack Lloyd
4		*
5		* Based on a public domain implemenation by Phil Karn (who in turn
6		* credited Richard Outerbridge and Jim Gillogly)
7		*
8		* Botan is released under the Simplified BSD License (see license.txt)
9		*/
10
11		#include <botan/internal/des.h>
12		#include <botan/internal/loadstor.h>
13		#include <botan/internal/rotate.h>
14
15		namespace Botan {
16
17		namespace {
18
19		alignas(256) const uint8_t SPBOX_CATS[64*8] = {
20		0x54, 0x00, 0x10, 0x55, 0x51, 0x15, 0x01, 0x10, 0x04, 0x54, 0x55, 0x04, 0x45, 0x51, 0x40, 0x01,
21		0x05, 0x44, 0x44, 0x14, 0x14, 0x50, 0x50, 0x45, 0x11, 0x41, 0x41, 0x11, 0x00, 0x05, 0x15, 0x40,
22		0x10, 0x55, 0x01, 0x50, 0x54, 0x40, 0x40, 0x04, 0x51, 0x10, 0x14, 0x41, 0x04, 0x01, 0x45, 0x15,
23		0x55, 0x11, 0x50, 0x45, 0x41, 0x05, 0x15, 0x54, 0x05, 0x44, 0x44, 0x00, 0x11, 0x14, 0x00, 0x51,
24
25		0x55, 0x44, 0x04, 0x15, 0x10, 0x01, 0x51, 0x45, 0x41, 0x55, 0x54, 0x40, 0x44, 0x10, 0x01, 0x51,
26		0x14, 0x11, 0x45, 0x00, 0x40, 0x04, 0x15, 0x50, 0x11, 0x41, 0x00, 0x14, 0x05, 0x54, 0x50, 0x05,
27		0x00, 0x15, 0x51, 0x10, 0x45, 0x50, 0x54, 0x04, 0x50, 0x44, 0x01, 0x55, 0x15, 0x01, 0x04, 0x40,
28		0x05, 0x54, 0x10, 0x41, 0x11, 0x45, 0x41, 0x11, 0x14, 0x00, 0x44, 0x05, 0x40, 0x51, 0x55, 0x14,
29
30		0x09, 0xA8, 0x00, 0xA1, 0x88, 0x00, 0x29, 0x88, 0x21, 0x81, 0x81, 0x20, 0xA9, 0x21, 0xA0, 0x09,
31		0x80, 0x01, 0xA8, 0x08, 0x28, 0xA0, 0xA1, 0x29, 0x89, 0x28, 0x20, 0x89, 0x01, 0xA9, 0x08, 0x80,
32		0xA8, 0x80, 0x21, 0x09, 0x20, 0xA8, 0x88, 0x00, 0x08, 0x21, 0xA9, 0x88, 0x81, 0x08, 0x00, 0xA1,
33		0x89, 0x20, 0x80, 0xA9, 0x01, 0x29, 0x28, 0x81, 0xA0, 0x89, 0x09, 0xA0, 0x29, 0x01, 0xA1, 0x28,
34
35		0x51, 0x15, 0x15, 0x04, 0x54, 0x45, 0x41, 0x11, 0x00, 0x50, 0x50, 0x55, 0x05, 0x00, 0x44, 0x41,
36		0x01, 0x10, 0x40, 0x51, 0x04, 0x40, 0x11, 0x14, 0x45, 0x01, 0x14, 0x44, 0x10, 0x54, 0x55, 0x05,
37		0x44, 0x41, 0x50, 0x55, 0x05, 0x00, 0x00, 0x50, 0x14, 0x44, 0x45, 0x01, 0x51, 0x15, 0x15, 0x04,
38		0x55, 0x05, 0x01, 0x10, 0x41, 0x11, 0x54, 0x45, 0x11, 0x14, 0x40, 0x51, 0x04, 0x40, 0x10, 0x54,
39
40		0x01, 0x29, 0x28, 0xA1, 0x08, 0x01, 0x80, 0x28, 0x89, 0x08, 0x21, 0x89, 0xA1, 0xA8, 0x09, 0x80,
41		0x20, 0x88, 0x88, 0x00, 0x81, 0xA9, 0xA9, 0x21, 0xA8, 0x81, 0x00, 0xA0, 0x29, 0x20, 0xA0, 0x09,
42		0x08, 0xA1, 0x01, 0x20, 0x80, 0x28, 0xA1, 0x89, 0x21, 0x80, 0xA8, 0x29, 0x89, 0x01, 0x20, 0xA8,
43		0xA9, 0x09, 0xA0, 0xA9, 0x28, 0x00, 0x88, 0xA0, 0x09, 0x21, 0x81, 0x08, 0x00, 0x88, 0x29, 0x81,
44
45		0x41, 0x50, 0x04, 0x55, 0x50, 0x01, 0x55, 0x10, 0x44, 0x15, 0x10, 0x41, 0x11, 0x44, 0x40, 0x05,
46		0x00, 0x11, 0x45, 0x04, 0x14, 0x45, 0x01, 0x51, 0x51, 0x00, 0x15, 0x54, 0x05, 0x14, 0x54, 0x40,
47		0x44, 0x01, 0x51, 0x14, 0x55, 0x10, 0x05, 0x41, 0x10, 0x44, 0x40, 0x05, 0x41, 0x55, 0x14, 0x50,
48		0x15, 0x54, 0x00, 0x51, 0x01, 0x04, 0x50, 0x15, 0x04, 0x11, 0x45, 0x00, 0x54, 0x40, 0x11, 0x45,
49
50		0x10, 0x51, 0x45, 0x00, 0x04, 0x45, 0x15, 0x54, 0x55, 0x10, 0x00, 0x41, 0x01, 0x40, 0x51, 0x05,
51		0x44, 0x15, 0x11, 0x44, 0x41, 0x50, 0x54, 0x11, 0x50, 0x04, 0x05, 0x55, 0x14, 0x01, 0x40, 0x14,
52		0x40, 0x14, 0x10, 0x45, 0x45, 0x51, 0x51, 0x01, 0x11, 0x40, 0x44, 0x10, 0x54, 0x05, 0x15, 0x54,
53		0x05, 0x41, 0x55, 0x50, 0x14, 0x00, 0x01, 0x55, 0x00, 0x15, 0x50, 0x04, 0x41, 0x44, 0x04, 0x11,
54
55		0x89, 0x08, 0x20, 0xA9, 0x80, 0x89, 0x01, 0x80, 0x21, 0xA0, 0xA9, 0x28, 0xA8, 0x29, 0x08, 0x01,
56		0xA0, 0x81, 0x88, 0x09, 0x28, 0x21, 0xA1, 0xA8, 0x09, 0x00, 0x00, 0xA1, 0x81, 0x88, 0x29, 0x20,
57		0x29, 0x20, 0xa8, 0x08, 0x01, 0xA1, 0x08, 0x29, 0x88, 0x01, 0x81, 0xA0, 0xA1, 0x80, 0x20, 0x89,
58		0x00, 0xA9, 0x21, 0x81, 0xA0, 0x88, 0x89, 0x00, 0xA9, 0x28, 0x28, 0x09, 0x09, 0x21, 0x80, 0xA8,
59		};
60
61		const uint32_t SPBOX_CAT_0_MUL = 0x70041106;
62		const uint32_t SPBOX_CAT_1_MUL = 0x02012020;
63		const uint32_t SPBOX_CAT_2_MUL = 0x00901048;
64		const uint32_t SPBOX_CAT_3_MUL = 0x8e060221;
65		const uint32_t SPBOX_CAT_4_MUL = 0x00912140;
66		const uint32_t SPBOX_CAT_5_MUL = 0x80841018;
67		const uint32_t SPBOX_CAT_6_MUL = 0xe0120202;
68		const uint32_t SPBOX_CAT_7_MUL = 0x00212240;
69
70		const uint32_t SPBOX_CAT_0_MASK = 0x01010404;
71		const uint32_t SPBOX_CAT_1_MASK = 0x80108020;
72		const uint32_t SPBOX_CAT_2_MASK = 0x08020208;
73		const uint32_t SPBOX_CAT_3_MASK = 0x00802081;
74		const uint32_t SPBOX_CAT_4_MASK = 0x42080100;
75		const uint32_t SPBOX_CAT_5_MASK = 0x20404010;
76		const uint32_t SPBOX_CAT_6_MASK = 0x04200802;
77		const uint32_t SPBOX_CAT_7_MASK = 0x10041040;
78
79		/*
80		* DES Key Schedule
81		*/
82		void des_key_schedule(uint32_t round_key[32], const uint8_t key[8])
83	0	{
84	0	static const uint8_t ROT[16] = { 1, 1, 2, 2, 2, 2, 2, 2,
85	0	1, 2, 2, 2, 2, 2, 2, 1 };
86
87	0	uint32_t C = ((key[7] & 0x80) << 20) \| ((key[6] & 0x80) << 19) \|
88	0	((key[5] & 0x80) << 18) \| ((key[4] & 0x80) << 17) \|
89	0	((key[3] & 0x80) << 16) \| ((key[2] & 0x80) << 15) \|
90	0	((key[1] & 0x80) << 14) \| ((key[0] & 0x80) << 13) \|
91	0	((key[7] & 0x40) << 13) \| ((key[6] & 0x40) << 12) \|
92	0	((key[5] & 0x40) << 11) \| ((key[4] & 0x40) << 10) \|
93	0	((key[3] & 0x40) << 9) \| ((key[2] & 0x40) << 8) \|
94	0	((key[1] & 0x40) << 7) \| ((key[0] & 0x40) << 6) \|
95	0	((key[7] & 0x20) << 6) \| ((key[6] & 0x20) << 5) \|
96	0	((key[5] & 0x20) << 4) \| ((key[4] & 0x20) << 3) \|
97	0	((key[3] & 0x20) << 2) \| ((key[2] & 0x20) << 1) \|
98	0	((key[1] & 0x20) ) \| ((key[0] & 0x20) >> 1) \|
99	0	((key[7] & 0x10) >> 1) \| ((key[6] & 0x10) >> 2) \|
100	0	((key[5] & 0x10) >> 3) \| ((key[4] & 0x10) >> 4);
101	0	uint32_t D = ((key[7] & 0x02) << 26) \| ((key[6] & 0x02) << 25) \|
102	0	((key[5] & 0x02) << 24) \| ((key[4] & 0x02) << 23) \|
103	0	((key[3] & 0x02) << 22) \| ((key[2] & 0x02) << 21) \|
104	0	((key[1] & 0x02) << 20) \| ((key[0] & 0x02) << 19) \|
105	0	((key[7] & 0x04) << 17) \| ((key[6] & 0x04) << 16) \|
106	0	((key[5] & 0x04) << 15) \| ((key[4] & 0x04) << 14) \|
107	0	((key[3] & 0x04) << 13) \| ((key[2] & 0x04) << 12) \|
108	0	((key[1] & 0x04) << 11) \| ((key[0] & 0x04) << 10) \|
109	0	((key[7] & 0x08) << 8) \| ((key[6] & 0x08) << 7) \|
110	0	((key[5] & 0x08) << 6) \| ((key[4] & 0x08) << 5) \|
111	0	((key[3] & 0x08) << 4) \| ((key[2] & 0x08) << 3) \|
112	0	((key[1] & 0x08) << 2) \| ((key[0] & 0x08) << 1) \|
113	0	((key[3] & 0x10) >> 1) \| ((key[2] & 0x10) >> 2) \|
114	0	((key[1] & 0x10) >> 3) \| ((key[0] & 0x10) >> 4);
115
116	0	for(size_t i = 0; i != 16; ++i)
117	0	{
118	0	C = ((C << ROT[i]) \| (C >> (28-ROT[i]))) & 0x0FFFFFFF;
119	0	D = ((D << ROT[i]) \| (D >> (28-ROT[i]))) & 0x0FFFFFFF;
120	0	round_key[2*i ] = ((C & 0x00000010) << 22) \| ((C & 0x00000800) << 17) \|
121	0	((C & 0x00000020) << 16) \| ((C & 0x00004004) << 15) \|
122	0	((C & 0x00000200) << 11) \| ((C & 0x00020000) << 10) \|
123	0	((C & 0x01000000) >> 6) \| ((C & 0x00100000) >> 4) \|
124	0	((C & 0x00010000) << 3) \| ((C & 0x08000000) >> 2) \|
125	0	((C & 0x00800000) << 1) \| ((D & 0x00000010) << 8) \|
126	0	((D & 0x00000002) << 7) \| ((D & 0x00000001) << 2) \|
127	0	((D & 0x00000200) ) \| ((D & 0x00008000) >> 2) \|
128	0	((D & 0x00000088) >> 3) \| ((D & 0x00001000) >> 7) \|
129	0	((D & 0x00080000) >> 9) \| ((D & 0x02020000) >> 14) \|
130	0	((D & 0x00400000) >> 21);
131	0	round_key[2*i+1] = ((C & 0x00000001) << 28) \| ((C & 0x00000082) << 18) \|
132	0	((C & 0x00002000) << 14) \| ((C & 0x00000100) << 10) \|
133	0	((C & 0x00001000) << 9) \| ((C & 0x00040000) << 6) \|
134	0	((C & 0x02400000) << 4) \| ((C & 0x00008000) << 2) \|
135	0	((C & 0x00200000) >> 1) \| ((C & 0x04000000) >> 10) \|
136	0	((D & 0x00000020) << 6) \| ((D & 0x00000100) ) \|
137	0	((D & 0x00000800) >> 1) \| ((D & 0x00000040) >> 3) \|
138	0	((D & 0x00010000) >> 4) \| ((D & 0x00000400) >> 5) \|
139	0	((D & 0x00004000) >> 10) \| ((D & 0x04000000) >> 13) \|
140	0	((D & 0x00800000) >> 14) \| ((D & 0x00100000) >> 18) \|
141	0	((D & 0x01000000) >> 24) \| ((D & 0x08000000) >> 26);
142	0	}
143	0	}
144
145		inline uint32_t spbox(uint32_t T0, uint32_t T1)
146	0	{
147	0	return
148	0	((SPBOX_CATS[064 + ((T0 >> 24) & 0x3F)] SPBOX_CAT_0_MUL) & SPBOX_CAT_0_MASK) ^
149	0	((SPBOX_CATS[164 + ((T1 >> 24) & 0x3F)] SPBOX_CAT_1_MUL) & SPBOX_CAT_1_MASK) ^
150	0	((SPBOX_CATS[264 + ((T0 >> 16) & 0x3F)] SPBOX_CAT_2_MUL) & SPBOX_CAT_2_MASK) ^
151	0	((SPBOX_CATS[364 + ((T1 >> 16) & 0x3F)] SPBOX_CAT_3_MUL) & SPBOX_CAT_3_MASK) ^
152	0	((SPBOX_CATS[464 + ((T0 >> 8) & 0x3F)] SPBOX_CAT_4_MUL) & SPBOX_CAT_4_MASK) ^
153	0	((SPBOX_CATS[564 + ((T1 >> 8) & 0x3F)] SPBOX_CAT_5_MUL) & SPBOX_CAT_5_MASK) ^
154	0	((SPBOX_CATS[664 + ((T0 >> 0) & 0x3F)] SPBOX_CAT_6_MUL) & SPBOX_CAT_6_MASK) ^
155	0	((SPBOX_CATS[764 + ((T1 >> 0) & 0x3F)] SPBOX_CAT_7_MUL) & SPBOX_CAT_7_MASK);
156	0	}
157
158		/*
159		* DES Encryption
160		*/
161		inline void des_encrypt(uint32_t& Lr, uint32_t& Rr,
162		const uint32_t round_key[32])
163	0	{
164	0	uint32_t L = Lr;
165	0	uint32_t R = Rr;
166	0	for(size_t i = 0; i != 16; i += 2)
167	0	{
168	0	L ^= spbox(rotr<4>(R) ^ round_key[2i ], R ^ round_key[2i+1]);
169	0	R ^= spbox(rotr<4>(L) ^ round_key[2i+2], L ^ round_key[2i+3]);
170	0	}
171
172	0	Lr = L;
173	0	Rr = R;
174	0	}
175
176		inline void des_encrypt_x2(uint32_t& L0r, uint32_t& R0r,
177		uint32_t& L1r, uint32_t& R1r,
178		const uint32_t round_key[32])
179	0	{
180	0	uint32_t L0 = L0r;
181	0	uint32_t R0 = R0r;
182	0	uint32_t L1 = L1r;
183	0	uint32_t R1 = R1r;
184
185	0	for(size_t i = 0; i != 16; i += 2)
186	0	{
187	0	L0 ^= spbox(rotr<4>(R0) ^ round_key[2i ], R0 ^ round_key[2i+1]);
188	0	L1 ^= spbox(rotr<4>(R1) ^ round_key[2i ], R1 ^ round_key[2i+1]);
189
190	0	R0 ^= spbox(rotr<4>(L0) ^ round_key[2i+2], L0 ^ round_key[2i+3]);
191	0	R1 ^= spbox(rotr<4>(L1) ^ round_key[2i+2], L1 ^ round_key[2i+3]);
192	0	}
193
194	0	L0r = L0;
195	0	R0r = R0;
196	0	L1r = L1;
197	0	R1r = R1;
198	0	}
199
200		/*
201		* DES Decryption
202		*/
203		inline void des_decrypt(uint32_t& Lr, uint32_t& Rr,
204		const uint32_t round_key[32])
205	0	{
206	0	uint32_t L = Lr;
207	0	uint32_t R = Rr;
208	0	for(size_t i = 16; i != 0; i -= 2)
209	0	{
210	0	L ^= spbox(rotr<4>(R) ^ round_key[2i - 2], R ^ round_key[2i - 1]);
211	0	R ^= spbox(rotr<4>(L) ^ round_key[2i - 4], L ^ round_key[2i - 3]);
212	0	}
213	0	Lr = L;
214	0	Rr = R;
215	0	}
216
217		inline void des_decrypt_x2(uint32_t& L0r, uint32_t& R0r,
218		uint32_t& L1r, uint32_t& R1r,
219		const uint32_t round_key[32])
220	0	{
221	0	uint32_t L0 = L0r;
222	0	uint32_t R0 = R0r;
223	0	uint32_t L1 = L1r;
224	0	uint32_t R1 = R1r;
225
226	0	for(size_t i = 16; i != 0; i -= 2)
227	0	{
228	0	L0 ^= spbox(rotr<4>(R0) ^ round_key[2i - 2], R0 ^ round_key[2i - 1]);
229	0	L1 ^= spbox(rotr<4>(R1) ^ round_key[2i - 2], R1 ^ round_key[2i - 1]);
230
231	0	R0 ^= spbox(rotr<4>(L0) ^ round_key[2i - 4], L0 ^ round_key[2i - 3]);
232	0	R1 ^= spbox(rotr<4>(L1) ^ round_key[2i - 4], L1 ^ round_key[2i - 3]);
233	0	}
234
235	0	L0r = L0;
236	0	R0r = R0;
237	0	L1r = L1;
238	0	R1r = R1;
239	0	}
240
241		inline void des_IP(uint32_t& L, uint32_t& R)
242	0	{
243		// IP sequence by Wei Dai, taken from public domain Crypto++
244	0	uint32_t T;
245	0	R = rotl<4>(R);
246	0	T = (L ^ R) & 0xF0F0F0F0;
247	0	L ^= T;
248	0	R = rotr<20>(R ^ T);
249	0	T = (L ^ R) & 0xFFFF0000;
250	0	L ^= T;
251	0	R = rotr<18>(R ^ T);
252	0	T = (L ^ R) & 0x33333333;
253	0	L ^= T;
254	0	R = rotr<6>(R ^ T);
255	0	T = (L ^ R) & 0x00FF00FF;
256	0	L ^= T;
257	0	R = rotl<9>(R ^ T);
258	0	T = (L ^ R) & 0xAAAAAAAA;
259	0	L = rotl<1>(L ^ T);
260	0	R ^= T;
261	0	}
262
263		inline void des_FP(uint32_t& L, uint32_t& R)
264	0	{
265		// FP sequence by Wei Dai, taken from public domain Crypto++
266	0	uint32_t T;
267
268	0	R = rotr<1>(R);
269	0	T = (L ^ R) & 0xAAAAAAAA;
270	0	R ^= T;
271	0	L = rotr<9>(L ^ T);
272	0	T = (L ^ R) & 0x00FF00FF;
273	0	R ^= T;
274	0	L = rotl<6>(L ^ T);
275	0	T = (L ^ R) & 0x33333333;
276	0	R ^= T;
277	0	L = rotl<18>(L ^ T);
278	0	T = (L ^ R) & 0xFFFF0000;
279	0	R ^= T;
280	0	L = rotl<20>(L ^ T);
281	0	T = (L ^ R) & 0xF0F0F0F0;
282	0	R ^= T;
283	0	L = rotr<4>(L ^ T);
284	0	}
285
286		}
287
288		/*
289		* DES Encryption
290		*/
291		void DES::encrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
292	0	{
293	0	assert_key_material_set();
294
295	0	while(blocks >= 2)
296	0	{
297	0	uint32_t L0 = load_be<uint32_t>(in, 0);
298	0	uint32_t R0 = load_be<uint32_t>(in, 1);
299	0	uint32_t L1 = load_be<uint32_t>(in, 2);
300	0	uint32_t R1 = load_be<uint32_t>(in, 3);
301
302	0	des_IP(L0, R0);
303	0	des_IP(L1, R1);
304
305	0	des_encrypt_x2(L0, R0, L1, R1, m_round_key.data());
306
307	0	des_FP(L0, R0);
308	0	des_FP(L1, R1);
309
310	0	store_be(out, R0, L0, R1, L1);
311
312	0	in += 2*BLOCK_SIZE;
313	0	out += 2*BLOCK_SIZE;
314	0	blocks -= 2;
315	0	}
316
317	0	while(blocks > 0)
318	0	{
319	0	uint32_t L0 = load_be<uint32_t>(in, 0);
320	0	uint32_t R0 = load_be<uint32_t>(in, 1);
321	0	des_IP(L0, R0);
322	0	des_encrypt(L0, R0, m_round_key.data());
323	0	des_FP(L0, R0);
324	0	store_be(out, R0, L0);
325
326	0	in += BLOCK_SIZE;
327	0	out += BLOCK_SIZE;
328	0	blocks -= 1;
329	0	}
330	0	}
331
332		/*
333		* DES Decryption
334		*/
335		void DES::decrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
336	0	{
337	0	assert_key_material_set();
338
339	0	while(blocks >= 2)
340	0	{
341	0	uint32_t L0 = load_be<uint32_t>(in, 0);
342	0	uint32_t R0 = load_be<uint32_t>(in, 1);
343	0	uint32_t L1 = load_be<uint32_t>(in, 2);
344	0	uint32_t R1 = load_be<uint32_t>(in, 3);
345
346	0	des_IP(L0, R0);
347	0	des_IP(L1, R1);
348
349	0	des_decrypt_x2(L0, R0, L1, R1, m_round_key.data());
350
351	0	des_FP(L0, R0);
352	0	des_FP(L1, R1);
353
354	0	store_be(out, R0, L0, R1, L1);
355
356	0	in += 2*BLOCK_SIZE;
357	0	out += 2*BLOCK_SIZE;
358	0	blocks -= 2;
359	0	}
360
361	0	while(blocks > 0)
362	0	{
363	0	uint32_t L0 = load_be<uint32_t>(in, 0);
364	0	uint32_t R0 = load_be<uint32_t>(in, 1);
365	0	des_IP(L0, R0);
366	0	des_decrypt(L0, R0, m_round_key.data());
367	0	des_FP(L0, R0);
368	0	store_be(out, R0, L0);
369
370	0	in += BLOCK_SIZE;
371	0	out += BLOCK_SIZE;
372	0	blocks -= 1;
373	0	}
374	0	}
375
376		bool DES::has_keying_material() const
377	0	{
378	0	return !m_round_key.empty();
379	0	}
380
381		/*
382		* DES Key Schedule
383		*/
384		void DES::key_schedule(const uint8_t key[], size_t /length/)
385	0	{
386	0	m_round_key.resize(32);
387	0	des_key_schedule(m_round_key.data(), key);
388	0	}
389
390		void DES::clear()
391	0	{
392	0	zap(m_round_key);
393	0	}
394
395		/*
396		* TripleDES Encryption
397		*/
398		void TripleDES::encrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
399	0	{
400	0	assert_key_material_set();
401
402	0	while(blocks >= 2)
403	0	{
404	0	uint32_t L0 = load_be<uint32_t>(in, 0);
405	0	uint32_t R0 = load_be<uint32_t>(in, 1);
406	0	uint32_t L1 = load_be<uint32_t>(in, 2);
407	0	uint32_t R1 = load_be<uint32_t>(in, 3);
408
409	0	des_IP(L0, R0);
410	0	des_IP(L1, R1);
411
412	0	des_encrypt_x2(L0, R0, L1, R1, &m_round_key[0]);
413	0	des_decrypt_x2(R0, L0, R1, L1, &m_round_key[32]);
414	0	des_encrypt_x2(L0, R0, L1, R1, &m_round_key[64]);
415
416	0	des_FP(L0, R0);
417	0	des_FP(L1, R1);
418
419	0	store_be(out, R0, L0, R1, L1);
420
421	0	in += 2*BLOCK_SIZE;
422	0	out += 2*BLOCK_SIZE;
423	0	blocks -= 2;
424	0	}
425
426	0	while(blocks > 0)
427	0	{
428	0	uint32_t L0 = load_be<uint32_t>(in, 0);
429	0	uint32_t R0 = load_be<uint32_t>(in, 1);
430
431	0	des_IP(L0, R0);
432	0	des_encrypt(L0, R0, &m_round_key[0]);
433	0	des_decrypt(R0, L0, &m_round_key[32]);
434	0	des_encrypt(L0, R0, &m_round_key[64]);
435	0	des_FP(L0, R0);
436
437	0	store_be(out, R0, L0);
438
439	0	in += BLOCK_SIZE;
440	0	out += BLOCK_SIZE;
441	0	blocks -= 1;
442	0	}
443	0	}
444
445		/*
446		* TripleDES Decryption
447		*/
448		void TripleDES::decrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
449	0	{
450	0	assert_key_material_set();
451
452	0	while(blocks >= 2)
453	0	{
454	0	uint32_t L0 = load_be<uint32_t>(in, 0);
455	0	uint32_t R0 = load_be<uint32_t>(in, 1);
456	0	uint32_t L1 = load_be<uint32_t>(in, 2);
457	0	uint32_t R1 = load_be<uint32_t>(in, 3);
458
459	0	des_IP(L0, R0);
460	0	des_IP(L1, R1);
461
462	0	des_decrypt_x2(L0, R0, L1, R1, &m_round_key[64]);
463	0	des_encrypt_x2(R0, L0, R1, L1, &m_round_key[32]);
464	0	des_decrypt_x2(L0, R0, L1, R1, &m_round_key[0]);
465
466	0	des_FP(L0, R0);
467	0	des_FP(L1, R1);
468
469	0	store_be(out, R0, L0, R1, L1);
470
471	0	in += 2*BLOCK_SIZE;
472	0	out += 2*BLOCK_SIZE;
473	0	blocks -= 2;
474	0	}
475
476	0	while(blocks > 0)
477	0	{
478	0	uint32_t L0 = load_be<uint32_t>(in, 0);
479	0	uint32_t R0 = load_be<uint32_t>(in, 1);
480
481	0	des_IP(L0, R0);
482	0	des_decrypt(L0, R0, &m_round_key[64]);
483	0	des_encrypt(R0, L0, &m_round_key[32]);
484	0	des_decrypt(L0, R0, &m_round_key[0]);
485	0	des_FP(L0, R0);
486
487	0	store_be(out, R0, L0);
488
489	0	in += BLOCK_SIZE;
490	0	out += BLOCK_SIZE;
491	0	blocks -= 1;
492	0	}
493	0	}
494
495		bool TripleDES::has_keying_material() const
496	0	{
497	0	return !m_round_key.empty();
498	0	}
499
500		/*
501		* TripleDES Key Schedule
502		*/
503		void TripleDES::key_schedule(const uint8_t key[], size_t length)
504	0	{
505	0	m_round_key.resize(3*32);
506	0	des_key_schedule(&m_round_key[0], key);
507	0	des_key_schedule(&m_round_key[32], key + 8);
508
509	0	if(length == 24)
510	0	des_key_schedule(&m_round_key[64], key + 16);
511	0	else
512	0	copy_mem(&m_round_key[64], &m_round_key[0], 32);
513	0	}
514
515		void TripleDES::clear()
516	0	{
517	0	zap(m_round_key);
518	0	}
519
520		}