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

Source (jump to first uncovered line)
/*
* DES
* (C) 1999-2008,2018 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/des.h>
#include <botan/loadstor.h>
#include <botan/rotate.h>

namespace Botan {

namespace {

/*
* 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 DES_SPBOX1[get_byte(0, T0)] ^ DES_SPBOX2[get_byte(0, T1)] ^
          DES_SPBOX3[get_byte(1, T0)] ^ DES_SPBOX4[get_byte(1, T1)] ^
          DES_SPBOX5[get_byte(2, T0)] ^ DES_SPBOX6[get_byte(2, T1)] ^
          DES_SPBOX7[get_byte(3, T0)] ^ DES_SPBOX8[get_byte(3, T1)];
   }

/*
* 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, const uint8_t block[])
   {
   // IP sequence by Wei Dai, taken from public domain Crypto++
   L = load_be<uint32_t>(block, 0);
   R = load_be<uint32_t>(block, 1);

   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, uint8_t out[])
   {
   // 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);

   store_be(out, R, L);
   }

}

/*
* DES Encryption
*/
void DES::encrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
   {
   verify_key_set(m_round_key.empty() == false);

   while(blocks >= 2)
      {
      uint32_t L0, R0;
      uint32_t L1, R1;

      des_IP(L0, R0, in);
      des_IP(L1, R1, in + BLOCK_SIZE);

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

      des_FP(L0, R0, out);
      des_FP(L1, R1, out + BLOCK_SIZE);

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

   for(size_t i = 0; i < blocks; ++i)
      {
      uint32_t L, R;
      des_IP(L, R, in + BLOCK_SIZE*i);
      des_encrypt(L, R, m_round_key.data());
      des_FP(L, R, out + BLOCK_SIZE*i);
      }
   }

/*
* DES Decryption
*/
void DES::decrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
   {
   verify_key_set(m_round_key.empty() == false);

   while(blocks >= 2)
      {
      uint32_t L0, R0;
      uint32_t L1, R1;

      des_IP(L0, R0, in);
      des_IP(L1, R1, in + BLOCK_SIZE);

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

      des_FP(L0, R0, out);
      des_FP(L1, R1, out + BLOCK_SIZE);

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

   for(size_t i = 0; i < blocks; ++i)
      {
      uint32_t L, R;
      des_IP(L, R, in + BLOCK_SIZE*i);
      des_decrypt(L, R, m_round_key.data());
      des_FP(L, R, out + BLOCK_SIZE*i);
      }
   }

/*
* DES Key Schedule
*/
void DES::key_schedule(const uint8_t key[], size_t)
   {
   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
   {
   verify_key_set(m_round_key.empty() == false);

   while(blocks >= 2)
      {
      uint32_t L0, R0;
      uint32_t L1, R1;

      des_IP(L0, R0, in);
      des_IP(L1, R1, in + BLOCK_SIZE);

      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, out);
      des_FP(L1, R1, out + BLOCK_SIZE);

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

   for(size_t i = 0; i != blocks; ++i)
      {
      uint32_t L, R;
      des_IP(L, R, in + BLOCK_SIZE*i);

      des_encrypt(L, R, &m_round_key[0]);
      des_decrypt(R, L, &m_round_key[32]);
      des_encrypt(L, R, &m_round_key[64]);

      des_FP(L, R, out + BLOCK_SIZE*i);
      }
   }

/*
* TripleDES Decryption
*/
void TripleDES::decrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
   {
   verify_key_set(m_round_key.empty() == false);

   while(blocks >= 2)
      {
      uint32_t L0, R0;
      uint32_t L1, R1;

      des_IP(L0, R0, in);
      des_IP(L1, R1, in + BLOCK_SIZE);

      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, out);
      des_FP(L1, R1, out + BLOCK_SIZE);

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

   for(size_t i = 0; i != blocks; ++i)
      {
      uint32_t L, R;
      des_IP(L, R, in + BLOCK_SIZE*i);

      des_decrypt(L, R, &m_round_key[64]);
      des_encrypt(R, L, &m_round_key[32]);
      des_decrypt(L, R, &m_round_key[0]);

      des_FP(L, R, out + BLOCK_SIZE*i);
      }
   }

/*
* 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: 2020-05-23 13:54

Line	Count	Source (jump to first uncovered line)
1		/*
2		* DES
3		* (C) 1999-2008,2018 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/des.h>
12		#include <botan/loadstor.h>
13		#include <botan/rotate.h>
14
15		namespace Botan {
16
17		namespace {
18
19		/*
20		* DES Key Schedule
21		*/
22		void des_key_schedule(uint32_t round_key[32], const uint8_t key[8])
23	885	{
24	885	static const uint8_t ROT[16] = { 1, 1, 2, 2, 2, 2, 2, 2,
25	885	1, 2, 2, 2, 2, 2, 2, 1 };
26	885
27	885	uint32_t C = ((key[7] & 0x80) << 20) \| ((key[6] & 0x80) << 19) \|
28	885	((key[5] & 0x80) << 18) \| ((key[4] & 0x80) << 17) \|
29	885	((key[3] & 0x80) << 16) \| ((key[2] & 0x80) << 15) \|
30	885	((key[1] & 0x80) << 14) \| ((key[0] & 0x80) << 13) \|
31	885	((key[7] & 0x40) << 13) \| ((key[6] & 0x40) << 12) \|
32	885	((key[5] & 0x40) << 11) \| ((key[4] & 0x40) << 10) \|
33	885	((key[3] & 0x40) << 9) \| ((key[2] & 0x40) << 8) \|
34	885	((key[1] & 0x40) << 7) \| ((key[0] & 0x40) << 6) \|
35	885	((key[7] & 0x20) << 6) \| ((key[6] & 0x20) << 5) \|
36	885	((key[5] & 0x20) << 4) \| ((key[4] & 0x20) << 3) \|
37	885	((key[3] & 0x20) << 2) \| ((key[2] & 0x20) << 1) \|
38	885	((key[1] & 0x20) ) \| ((key[0] & 0x20) >> 1) \|
39	885	((key[7] & 0x10) >> 1) \| ((key[6] & 0x10) >> 2) \|
40	885	((key[5] & 0x10) >> 3) \| ((key[4] & 0x10) >> 4);
41	885	uint32_t D = ((key[7] & 0x02) << 26) \| ((key[6] & 0x02) << 25) \|
42	885	((key[5] & 0x02) << 24) \| ((key[4] & 0x02) << 23) \|
43	885	((key[3] & 0x02) << 22) \| ((key[2] & 0x02) << 21) \|
44	885	((key[1] & 0x02) << 20) \| ((key[0] & 0x02) << 19) \|
45	885	((key[7] & 0x04) << 17) \| ((key[6] & 0x04) << 16) \|
46	885	((key[5] & 0x04) << 15) \| ((key[4] & 0x04) << 14) \|
47	885	((key[3] & 0x04) << 13) \| ((key[2] & 0x04) << 12) \|
48	885	((key[1] & 0x04) << 11) \| ((key[0] & 0x04) << 10) \|
49	885	((key[7] & 0x08) << 8) \| ((key[6] & 0x08) << 7) \|
50	885	((key[5] & 0x08) << 6) \| ((key[4] & 0x08) << 5) \|
51	885	((key[3] & 0x08) << 4) \| ((key[2] & 0x08) << 3) \|
52	885	((key[1] & 0x08) << 2) \| ((key[0] & 0x08) << 1) \|
53	885	((key[3] & 0x10) >> 1) \| ((key[2] & 0x10) >> 2) \|
54	885	((key[1] & 0x10) >> 3) \| ((key[0] & 0x10) >> 4);
55	885
56	15.0k	for(size_t i = 0; i != 16; ++i)
57	14.1k	{
58	14.1k	C = ((C << ROT[i]) \| (C >> (28-ROT[i]))) & 0x0FFFFFFF;
59	14.1k	D = ((D << ROT[i]) \| (D >> (28-ROT[i]))) & 0x0FFFFFFF;
60	14.1k	round_key[2*i ] = ((C & 0x00000010) << 22) \| ((C & 0x00000800) << 17) \|
61	14.1k	((C & 0x00000020) << 16) \| ((C & 0x00004004) << 15) \|
62	14.1k	((C & 0x00000200) << 11) \| ((C & 0x00020000) << 10) \|
63	14.1k	((C & 0x01000000) >> 6) \| ((C & 0x00100000) >> 4) \|
64	14.1k	((C & 0x00010000) << 3) \| ((C & 0x08000000) >> 2) \|
65	14.1k	((C & 0x00800000) << 1) \| ((D & 0x00000010) << 8) \|
66	14.1k	((D & 0x00000002) << 7) \| ((D & 0x00000001) << 2) \|
67	14.1k	((D & 0x00000200) ) \| ((D & 0x00008000) >> 2) \|
68	14.1k	((D & 0x00000088) >> 3) \| ((D & 0x00001000) >> 7) \|
69	14.1k	((D & 0x00080000) >> 9) \| ((D & 0x02020000) >> 14) \|
70	14.1k	((D & 0x00400000) >> 21);
71	14.1k	round_key[2*i+1] = ((C & 0x00000001) << 28) \| ((C & 0x00000082) << 18) \|
72	14.1k	((C & 0x00002000) << 14) \| ((C & 0x00000100) << 10) \|
73	14.1k	((C & 0x00001000) << 9) \| ((C & 0x00040000) << 6) \|
74	14.1k	((C & 0x02400000) << 4) \| ((C & 0x00008000) << 2) \|
75	14.1k	((C & 0x00200000) >> 1) \| ((C & 0x04000000) >> 10) \|
76	14.1k	((D & 0x00000020) << 6) \| ((D & 0x00000100) ) \|
77	14.1k	((D & 0x00000800) >> 1) \| ((D & 0x00000040) >> 3) \|
78	14.1k	((D & 0x00010000) >> 4) \| ((D & 0x00000400) >> 5) \|
79	14.1k	((D & 0x00004000) >> 10) \| ((D & 0x04000000) >> 13) \|
80	14.1k	((D & 0x00800000) >> 14) \| ((D & 0x00100000) >> 18) \|
81	14.1k	((D & 0x01000000) >> 24) \| ((D & 0x08000000) >> 26);
82	14.1k	}
83	885	}
84
85		inline uint32_t spbox(uint32_t T0, uint32_t T1)
86	795k	{
87	795k	return DES_SPBOX1[get_byte(0, T0)] ^ DES_SPBOX2[get_byte(0, T1)] ^
88	795k	DES_SPBOX3[get_byte(1, T0)] ^ DES_SPBOX4[get_byte(1, T1)] ^
89	795k	DES_SPBOX5[get_byte(2, T0)] ^ DES_SPBOX6[get_byte(2, T1)] ^
90	795k	DES_SPBOX7[get_byte(3, T0)] ^ DES_SPBOX8[get_byte(3, T1)];
91	795k	}
92
93		/*
94		* DES Encryption
95		*/
96		inline void des_encrypt(uint32_t& Lr, uint32_t& Rr,
97		const uint32_t round_key[32])
98	1.89k	{
99	1.89k	uint32_t L = Lr;
100	1.89k	uint32_t R = Rr;
101	17.0k	for(size_t i = 0; i != 16; i += 2)
102	15.1k	{
103	15.1k	L ^= spbox(rotr<4>(R) ^ round_key[2i ], R ^ round_key[2i+1]);
104	15.1k	R ^= spbox(rotr<4>(L) ^ round_key[2i+2], L ^ round_key[2i+3]);
105	15.1k	}
106	1.89k
107	1.89k	Lr = L;
108	1.89k	Rr = R;
109	1.89k	}
110
111		inline void des_encrypt_x2(uint32_t& L0r, uint32_t& R0r,
112		uint32_t& L1r, uint32_t& R1r,
113		const uint32_t round_key[32])
114	7.79k	{
115	7.79k	uint32_t L0 = L0r;
116	7.79k	uint32_t R0 = R0r;
117	7.79k	uint32_t L1 = L1r;
118	7.79k	uint32_t R1 = R1r;
119	7.79k
120	70.1k	for(size_t i = 0; i != 16; i += 2)
121	62.3k	{
122	62.3k	L0 ^= spbox(rotr<4>(R0) ^ round_key[2i ], R0 ^ round_key[2i+1]);
123	62.3k	L1 ^= spbox(rotr<4>(R1) ^ round_key[2i ], R1 ^ round_key[2i+1]);
124	62.3k
125	62.3k	R0 ^= spbox(rotr<4>(L0) ^ round_key[2i+2], L0 ^ round_key[2i+3]);
126	62.3k	R1 ^= spbox(rotr<4>(L1) ^ round_key[2i+2], L1 ^ round_key[2i+3]);
127	62.3k	}
128	7.79k
129	7.79k	L0r = L0;
130	7.79k	R0r = R0;
131	7.79k	L1r = L1;
132	7.79k	R1r = R1;
133	7.79k	}
134
135		/*
136		* DES Decryption
137		*/
138		inline void des_decrypt(uint32_t& Lr, uint32_t& Rr,
139		const uint32_t round_key[32])
140	1.05k	{
141	1.05k	uint32_t L = Lr;
142	1.05k	uint32_t R = Rr;
143	9.51k	for(size_t i = 16; i != 0; i -= 2)
144	8.45k	{
145	8.45k	L ^= spbox(rotr<4>(R) ^ round_key[2i - 2], R ^ round_key[2i - 1]);
146	8.45k	R ^= spbox(rotr<4>(L) ^ round_key[2i - 4], L ^ round_key[2i - 3]);
147	8.45k	}
148	1.05k	Lr = L;
149	1.05k	Rr = R;
150	1.05k	}
151
152		inline void des_decrypt_x2(uint32_t& L0r, uint32_t& R0r,
153		uint32_t& L1r, uint32_t& R1r,
154		const uint32_t round_key[32])
155	15.5k	{
156	15.5k	uint32_t L0 = L0r;
157	15.5k	uint32_t R0 = R0r;
158	15.5k	uint32_t L1 = L1r;
159	15.5k	uint32_t R1 = R1r;
160	15.5k
161	140k	for(size_t i = 16; i != 0; i -= 2)
162	124k	{
163	124k	L0 ^= spbox(rotr<4>(R0) ^ round_key[2i - 2], R0 ^ round_key[2i - 1]);
164	124k	L1 ^= spbox(rotr<4>(R1) ^ round_key[2i - 2], R1 ^ round_key[2i - 1]);
165	124k
166	124k	R0 ^= spbox(rotr<4>(L0) ^ round_key[2i - 4], L0 ^ round_key[2i - 3]);
167	124k	R1 ^= spbox(rotr<4>(L1) ^ round_key[2i - 4], L1 ^ round_key[2i - 3]);
168	124k	}
169	15.5k
170	15.5k	L0r = L0;
171	15.5k	R0r = R0;
172	15.5k	L1r = L1;
173	15.5k	R1r = R1;
174	15.5k	}
175
176		inline void des_IP(uint32_t& L, uint32_t& R, const uint8_t block[])
177	16.5k	{
178	16.5k	// IP sequence by Wei Dai, taken from public domain Crypto++
179	16.5k	L = load_be<uint32_t>(block, 0);
180	16.5k	R = load_be<uint32_t>(block, 1);
181	16.5k
182	16.5k	uint32_t T;
183	16.5k	R = rotl<4>(R);
184	16.5k	T = (L ^ R) & 0xF0F0F0F0;
185	16.5k	L ^= T;
186	16.5k	R = rotr<20>(R ^ T);
187	16.5k	T = (L ^ R) & 0xFFFF0000;
188	16.5k	L ^= T;
189	16.5k	R = rotr<18>(R ^ T);
190	16.5k	T = (L ^ R) & 0x33333333;
191	16.5k	L ^= T;
192	16.5k	R = rotr<6>(R ^ T);
193	16.5k	T = (L ^ R) & 0x00FF00FF;
194	16.5k	L ^= T;
195	16.5k	R = rotl<9>(R ^ T);
196	16.5k	T = (L ^ R) & 0xAAAAAAAA;
197	16.5k	L = rotl<1>(L ^ T);
198	16.5k	R ^= T;
199	16.5k	}
200
201		inline void des_FP(uint32_t L, uint32_t R, uint8_t out[])
202	16.5k	{
203	16.5k	// FP sequence by Wei Dai, taken from public domain Crypto++
204	16.5k	uint32_t T;
205	16.5k
206	16.5k	R = rotr<1>(R);
207	16.5k	T = (L ^ R) & 0xAAAAAAAA;
208	16.5k	R ^= T;
209	16.5k	L = rotr<9>(L ^ T);
210	16.5k	T = (L ^ R) & 0x00FF00FF;
211	16.5k	R ^= T;
212	16.5k	L = rotl<6>(L ^ T);
213	16.5k	T = (L ^ R) & 0x33333333;
214	16.5k	R ^= T;
215	16.5k	L = rotl<18>(L ^ T);
216	16.5k	T = (L ^ R) & 0xFFFF0000;
217	16.5k	R ^= T;
218	16.5k	L = rotl<20>(L ^ T);
219	16.5k	T = (L ^ R) & 0xF0F0F0F0;
220	16.5k	R ^= T;
221	16.5k	L = rotr<4>(L ^ T);
222	16.5k
223	16.5k	store_be(out, R, L);
224	16.5k	}
225
226		}
227
228		/*
229		* DES Encryption
230		*/
231		void DES::encrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
232	0	{
233	0	verify_key_set(m_round_key.empty() == false);
234	0
235	0	while(blocks >= 2)
236	0	{
237	0	uint32_t L0, R0;
238	0	uint32_t L1, R1;
239	0
240	0	des_IP(L0, R0, in);
241	0	des_IP(L1, R1, in + BLOCK_SIZE);
242	0
243	0	des_encrypt_x2(L0, R0, L1, R1, m_round_key.data());
244	0
245	0	des_FP(L0, R0, out);
246	0	des_FP(L1, R1, out + BLOCK_SIZE);
247	0
248	0	in += 2*BLOCK_SIZE;
249	0	out += 2*BLOCK_SIZE;
250	0	blocks -= 2;
251	0	}
252	0
253	0	for(size_t i = 0; i < blocks; ++i)
254	0	{
255	0	uint32_t L, R;
256	0	des_IP(L, R, in + BLOCK_SIZE*i);
257	0	des_encrypt(L, R, m_round_key.data());
258	0	des_FP(L, R, out + BLOCK_SIZE*i);
259	0	}
260	0	}
261
262		/*
263		* DES Decryption
264		*/
265		void DES::decrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
266	0	{
267	0	verify_key_set(m_round_key.empty() == false);
268	0
269	0	while(blocks >= 2)
270	0	{
271	0	uint32_t L0, R0;
272	0	uint32_t L1, R1;
273	0
274	0	des_IP(L0, R0, in);
275	0	des_IP(L1, R1, in + BLOCK_SIZE);
276	0
277	0	des_decrypt_x2(L0, R0, L1, R1, m_round_key.data());
278	0
279	0	des_FP(L0, R0, out);
280	0	des_FP(L1, R1, out + BLOCK_SIZE);
281	0
282	0	in += 2*BLOCK_SIZE;
283	0	out += 2*BLOCK_SIZE;
284	0	blocks -= 2;
285	0	}
286	0
287	0	for(size_t i = 0; i < blocks; ++i)
288	0	{
289	0	uint32_t L, R;
290	0	des_IP(L, R, in + BLOCK_SIZE*i);
291	0	des_decrypt(L, R, m_round_key.data());
292	0	des_FP(L, R, out + BLOCK_SIZE*i);
293	0	}
294	0	}
295
296		/*
297		* DES Key Schedule
298		*/
299		void DES::key_schedule(const uint8_t key[], size_t)
300	0	{
301	0	m_round_key.resize(32);
302	0	des_key_schedule(m_round_key.data(), key);
303	0	}
304
305		void DES::clear()
306	0	{
307	0	zap(m_round_key);
308	0	}
309
310		/*
311		* TripleDES Encryption
312		*/
313		void TripleDES::encrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
314	909	{
315	909	verify_key_set(m_round_key.empty() == false);
316	909
317	909	while(blocks >= 2)
318	0	{
319	0	uint32_t L0, R0;
320	0	uint32_t L1, R1;
321	0
322	0	des_IP(L0, R0, in);
323	0	des_IP(L1, R1, in + BLOCK_SIZE);
324	0
325	0	des_encrypt_x2(L0, R0, L1, R1, &m_round_key[0]);
326	0	des_decrypt_x2(R0, L0, R1, L1, &m_round_key[32]);
327	0	des_encrypt_x2(L0, R0, L1, R1, &m_round_key[64]);
328	0
329	0	des_FP(L0, R0, out);
330	0	des_FP(L1, R1, out + BLOCK_SIZE);
331	0
332	0	in += 2*BLOCK_SIZE;
333	0	out += 2*BLOCK_SIZE;
334	0	blocks -= 2;
335	0	}
336	909
337	1.81k	for(size_t i = 0; i != blocks; ++i)
338	909	{
339	909	uint32_t L, R;
340	909	des_IP(L, R, in + BLOCK_SIZE*i);
341	909
342	909	des_encrypt(L, R, &m_round_key[0]);
343	909	des_decrypt(R, L, &m_round_key[32]);
344	909	des_encrypt(L, R, &m_round_key[64]);
345	909
346	909	des_FP(L, R, out + BLOCK_SIZE*i);
347	909	}
348	909	}
349
350		/*
351		* TripleDES Decryption
352		*/
353		void TripleDES::decrypt_n(const uint8_t in[], uint8_t out[], size_t blocks) const
354	3.96k	{
355	3.96k	verify_key_set(m_round_key.empty() == false);
356	3.96k
357	11.7k	while(blocks >= 2)
358	7.79k	{
359	7.79k	uint32_t L0, R0;
360	7.79k	uint32_t L1, R1;
361	7.79k
362	7.79k	des_IP(L0, R0, in);
363	7.79k	des_IP(L1, R1, in + BLOCK_SIZE);
364	7.79k
365	7.79k	des_decrypt_x2(L0, R0, L1, R1, &m_round_key[64]);
366	7.79k	des_encrypt_x2(R0, L0, R1, L1, &m_round_key[32]);
367	7.79k	des_decrypt_x2(L0, R0, L1, R1, &m_round_key[0]);
368	7.79k
369	7.79k	des_FP(L0, R0, out);
370	7.79k	des_FP(L1, R1, out + BLOCK_SIZE);
371	7.79k
372	7.79k	in += 2*BLOCK_SIZE;
373	7.79k	out += 2*BLOCK_SIZE;
374	7.79k	blocks -= 2;
375	7.79k	}
376	3.96k
377	4.03k	for(size_t i = 0; i != blocks; ++i)
378	74	{
379	74	uint32_t L, R;
380	74	des_IP(L, R, in + BLOCK_SIZE*i);
381	74
382	74	des_decrypt(L, R, &m_round_key[64]);
383	74	des_encrypt(R, L, &m_round_key[32]);
384	74	des_decrypt(L, R, &m_round_key[0]);
385	74
386	74	des_FP(L, R, out + BLOCK_SIZE*i);
387	74	}
388	3.96k	}
389
390		/*
391		* TripleDES Key Schedule
392		*/
393		void TripleDES::key_schedule(const uint8_t key[], size_t length)
394	295	{
395	295	m_round_key.resize(3*32);
396	295	des_key_schedule(&m_round_key[0], key);
397	295	des_key_schedule(&m_round_key[32], key + 8);
398	295
399	295	if(length == 24)
400	295	des_key_schedule(&m_round_key[64], key + 16);
401	0	else
402	0	copy_mem(&m_round_key[64], &m_round_key[0], 32);
403	295	}
404
405		void TripleDES::clear()
406	0	{
407	0	zap(m_round_key);
408	0	}
409
410		}