michael@0: /*
michael@0:  *  mksp.c
michael@0:  *
michael@0:  *  Generate SP tables for DES-150 library
michael@0:  *
michael@0:  * This Source Code Form is subject to the terms of the Mozilla Public
michael@0:  * License, v. 2.0. If a copy of the MPL was not distributed with this
michael@0:  * file, You can obtain one at http://mozilla.org/MPL/2.0/. */
michael@0: 
michael@0: #include <stdio.h>
michael@0: 
michael@0: /*
michael@0:  * sboxes - the tables for the s-box functions
michael@0:  *        from FIPS 46, pages 15-16.
michael@0:  */
michael@0: unsigned char S[8][64] = {
michael@0: /* Func S1 = */ {
michael@0: 	14,  0,  4, 15, 13,  7,  1,  4,  2, 14, 15,  2, 11, 13,  8,  1,
michael@0: 	 3, 10, 10,  6,  6, 12, 12, 11,  5,  9,  9,  5,  0,  3,  7,  8,
michael@0: 	 4, 15,  1, 12, 14,  8,  8,  2, 13,  4,  6,  9,  2,  1, 11,  7,
michael@0: 	15,  5, 12, 11,  9,  3,  7, 14,  3, 10, 10,  0,  5,  6,  0, 13
michael@0:     },
michael@0: /* Func S2 = */ {
michael@0: 	15,  3,  1, 13,  8,  4, 14,  7,  6, 15, 11,  2,  3,  8,  4, 14,
michael@0: 	 9, 12,  7,  0,  2,  1, 13, 10, 12,  6,  0,  9,  5, 11, 10,  5,
michael@0: 	 0, 13, 14,  8,  7, 10, 11,  1, 10,  3,  4, 15, 13,  4,  1,  2,
michael@0: 	 5, 11,  8,  6, 12,  7,  6, 12,  9,  0,  3,  5,  2, 14, 15,  9
michael@0:     },
michael@0: /* Func S3 = */ {
michael@0: 	10, 13,  0,  7,  9,  0, 14,  9,  6,  3,  3,  4, 15,  6,  5, 10,
michael@0: 	 1,  2, 13,  8, 12,  5,  7, 14, 11, 12,  4, 11,  2, 15,  8,  1,
michael@0: 	13,  1,  6, 10,  4, 13,  9,  0,  8,  6, 15,  9,  3,  8,  0,  7,
michael@0: 	11,  4,  1, 15,  2, 14, 12,  3,  5, 11, 10,  5, 14,  2,  7, 12
michael@0:     },
michael@0: /* Func S4 = */ {
michael@0: 	 7, 13, 13,  8, 14, 11,  3,  5,  0,  6,  6, 15,  9,  0, 10,  3,
michael@0: 	 1,  4,  2,  7,  8,  2,  5, 12, 11,  1, 12, 10,  4, 14, 15,  9,
michael@0: 	10,  3,  6, 15,  9,  0,  0,  6, 12, 10, 11,  1,  7, 13, 13,  8,
michael@0: 	15,  9,  1,  4,  3,  5, 14, 11,  5, 12,  2,  7,  8,  2,  4, 14
michael@0:     },
michael@0: /* Func S5 = */ {
michael@0: 	 2, 14, 12, 11,  4,  2,  1, 12,  7,  4, 10,  7, 11, 13,  6,  1,
michael@0: 	 8,  5,  5,  0,  3, 15, 15, 10, 13,  3,  0,  9, 14,  8,  9,  6,
michael@0: 	 4, 11,  2,  8,  1, 12, 11,  7, 10,  1, 13, 14,  7,  2,  8, 13,
michael@0: 	15,  6,  9, 15, 12,  0,  5,  9,  6, 10,  3,  4,  0,  5, 14,  3
michael@0:     },
michael@0: /* Func S6 = */ {
michael@0: 	12, 10,  1, 15, 10,  4, 15,  2,  9,  7,  2, 12,  6,  9,  8,  5,
michael@0: 	 0,  6, 13,  1,  3, 13,  4, 14, 14,  0,  7, 11,  5,  3, 11,  8,
michael@0: 	 9,  4, 14,  3, 15,  2,  5, 12,  2,  9,  8,  5, 12, 15,  3, 10,
michael@0: 	 7, 11,  0, 14,  4,  1, 10,  7,  1,  6, 13,  0, 11,  8,  6, 13
michael@0:     },
michael@0: /* Func S7 = */ {
michael@0: 	 4, 13, 11,  0,  2, 11, 14,  7, 15,  4,  0,  9,  8,  1, 13, 10,
michael@0: 	 3, 14, 12,  3,  9,  5,  7, 12,  5,  2, 10, 15,  6,  8,  1,  6,
michael@0: 	 1,  6,  4, 11, 11, 13, 13,  8, 12,  1,  3,  4,  7, 10, 14,  7,
michael@0: 	10,  9, 15,  5,  6,  0,  8, 15,  0, 14,  5,  2,  9,  3,  2, 12
michael@0:     },
michael@0: /* Func S8 = */ {
michael@0: 	13,  1,  2, 15,  8, 13,  4,  8,  6, 10, 15,  3, 11,  7,  1,  4,
michael@0: 	10, 12,  9,  5,  3,  6, 14, 11,  5,  0,  0, 14, 12,  9,  7,  2,
michael@0: 	 7,  2, 11,  1,  4, 14,  1,  7,  9,  4, 12, 10, 14,  8,  2, 13,
michael@0: 	 0, 15,  6, 12, 10,  9, 13,  0, 15,  3,  3,  5,  5,  6,  8, 11
michael@0:     }
michael@0: };
michael@0: 
michael@0: /*
michael@0:  * Permutation function for results from s-boxes
michael@0:  *   from FIPS 46 pages 12 and 16.
michael@0:  * P = 
michael@0:  */
michael@0: unsigned char P[32] = {
michael@0: 	16,   7,  20,  21,  29,  12,  28,  17,
michael@0: 	 1,  15,  23,  26,   5,  18,  31,  10,
michael@0: 	 2,   8,  24,  14,  32,  27,   3,   9,
michael@0: 	19,  13,  30,   6,  22,  11,   4,  25
michael@0: };
michael@0: 
michael@0: unsigned int Pinv[32];
michael@0: unsigned int SP[8][64];
michael@0: 
michael@0: void
michael@0: makePinv(void)
michael@0: {
michael@0:     int i;
michael@0:     unsigned int Pi = 0x80000000;
michael@0:     for (i = 0; i < 32; ++i) {
michael@0:     	int j = 32 - P[i];
michael@0: 	Pinv[j] = Pi;
michael@0: 	Pi >>= 1;
michael@0:     }
michael@0: }
michael@0: 
michael@0: void
michael@0: makeSP(void)
michael@0: {
michael@0:     int box;
michael@0:     for (box = 0; box < 8; ++box) {
michael@0: 	int item;
michael@0: 	printf("/* box S%d */ {\n", box + 1);
michael@0:     	for (item = 0; item < 64; ++item ) {
michael@0: 	    unsigned int s = S[box][item];
michael@0: 	    unsigned int val = 0;
michael@0: 	    unsigned int bitnum = (7-box) * 4;
michael@0: 	    for (; s; s >>= 1, ++bitnum) {
michael@0: 		if (s & 1) {
michael@0: 		    val |= Pinv[bitnum];
michael@0: 		}
michael@0: 	    }
michael@0: 	    val = (val << 3) | (val >> 29);
michael@0: 	    SP[box][item] = val;
michael@0: 	}
michael@0: 	for (item = 0; item < 64; item += 4) {
michael@0: 	    printf("\t0x%08x, 0x%08x, 0x%08x, 0x%08x,\n",
michael@0: 	    SP[box][item], SP[box][item+1], SP[box][item+2], SP[box][item+3]);
michael@0: 	}
michael@0: 	printf("    },\n");
michael@0:     }
michael@0: }
michael@0: 
michael@0: int
michael@0: main()
michael@0: {
michael@0:     makePinv();
michael@0:     makeSP();
michael@0:     return 0;
michael@0: }