1*05b00f60SXin Li /*
2*05b00f60SXin Li * Copyright (c) 1998-2006 The TCPDUMP project
3*05b00f60SXin Li *
4*05b00f60SXin Li * Redistribution and use in source and binary forms, with or without
5*05b00f60SXin Li * modification, are permitted provided that: (1) source code
6*05b00f60SXin Li * distributions retain the above copyright notice and this paragraph
7*05b00f60SXin Li * in its entirety, and (2) distributions including binary code include
8*05b00f60SXin Li * the above copyright notice and this paragraph in its entirety in
9*05b00f60SXin Li * the documentation or other materials provided with the distribution.
10*05b00f60SXin Li * THIS SOFTWARE IS PROVIDED ``AS IS'' AND
11*05b00f60SXin Li * WITHOUT ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, WITHOUT
12*05b00f60SXin Li * LIMITATION, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
13*05b00f60SXin Li * FOR A PARTICULAR PURPOSE.
14*05b00f60SXin Li *
15*05b00f60SXin Li * miscellaneous checksumming routines
16*05b00f60SXin Li *
17*05b00f60SXin Li * Original code by Hannes Gredler ([email protected])
18*05b00f60SXin Li */
19*05b00f60SXin Li
20*05b00f60SXin Li #ifdef HAVE_CONFIG_H
21*05b00f60SXin Li #include <config.h>
22*05b00f60SXin Li #endif
23*05b00f60SXin Li
24*05b00f60SXin Li #include "netdissect-stdinc.h"
25*05b00f60SXin Li
26*05b00f60SXin Li #include <stdio.h>
27*05b00f60SXin Li #include <stdlib.h>
28*05b00f60SXin Li #include <string.h>
29*05b00f60SXin Li #include <assert.h>
30*05b00f60SXin Li
31*05b00f60SXin Li #include "netdissect.h"
32*05b00f60SXin Li
33*05b00f60SXin Li /*
34*05b00f60SXin Li * CRC-10 table generated using the following Python snippet:
35*05b00f60SXin Li
36*05b00f60SXin Li import sys
37*05b00f60SXin Li
38*05b00f60SXin Li crc_table = []
39*05b00f60SXin Li for i in range(256):
40*05b00f60SXin Li accum = i << 2
41*05b00f60SXin Li for j in range(8):
42*05b00f60SXin Li accum <<= 1
43*05b00f60SXin Li if accum & 0x400:
44*05b00f60SXin Li accum ^= 0x633
45*05b00f60SXin Li crc_table.append(accum)
46*05b00f60SXin Li
47*05b00f60SXin Li for i in range(len(crc_table)/8):
48*05b00f60SXin Li for j in range(8):
49*05b00f60SXin Li sys.stdout.write("0x%04x, " % crc_table[i*8+j])
50*05b00f60SXin Li sys.stdout.write("\n")
51*05b00f60SXin Li
52*05b00f60SXin Li */
53*05b00f60SXin Li static const uint16_t crc10_table[256] =
54*05b00f60SXin Li {
55*05b00f60SXin Li 0x0000, 0x0233, 0x0255, 0x0066, 0x0299, 0x00aa, 0x00cc, 0x02ff,
56*05b00f60SXin Li 0x0301, 0x0132, 0x0154, 0x0367, 0x0198, 0x03ab, 0x03cd, 0x01fe,
57*05b00f60SXin Li 0x0031, 0x0202, 0x0264, 0x0057, 0x02a8, 0x009b, 0x00fd, 0x02ce,
58*05b00f60SXin Li 0x0330, 0x0103, 0x0165, 0x0356, 0x01a9, 0x039a, 0x03fc, 0x01cf,
59*05b00f60SXin Li 0x0062, 0x0251, 0x0237, 0x0004, 0x02fb, 0x00c8, 0x00ae, 0x029d,
60*05b00f60SXin Li 0x0363, 0x0150, 0x0136, 0x0305, 0x01fa, 0x03c9, 0x03af, 0x019c,
61*05b00f60SXin Li 0x0053, 0x0260, 0x0206, 0x0035, 0x02ca, 0x00f9, 0x009f, 0x02ac,
62*05b00f60SXin Li 0x0352, 0x0161, 0x0107, 0x0334, 0x01cb, 0x03f8, 0x039e, 0x01ad,
63*05b00f60SXin Li 0x00c4, 0x02f7, 0x0291, 0x00a2, 0x025d, 0x006e, 0x0008, 0x023b,
64*05b00f60SXin Li 0x03c5, 0x01f6, 0x0190, 0x03a3, 0x015c, 0x036f, 0x0309, 0x013a,
65*05b00f60SXin Li 0x00f5, 0x02c6, 0x02a0, 0x0093, 0x026c, 0x005f, 0x0039, 0x020a,
66*05b00f60SXin Li 0x03f4, 0x01c7, 0x01a1, 0x0392, 0x016d, 0x035e, 0x0338, 0x010b,
67*05b00f60SXin Li 0x00a6, 0x0295, 0x02f3, 0x00c0, 0x023f, 0x000c, 0x006a, 0x0259,
68*05b00f60SXin Li 0x03a7, 0x0194, 0x01f2, 0x03c1, 0x013e, 0x030d, 0x036b, 0x0158,
69*05b00f60SXin Li 0x0097, 0x02a4, 0x02c2, 0x00f1, 0x020e, 0x003d, 0x005b, 0x0268,
70*05b00f60SXin Li 0x0396, 0x01a5, 0x01c3, 0x03f0, 0x010f, 0x033c, 0x035a, 0x0169,
71*05b00f60SXin Li 0x0188, 0x03bb, 0x03dd, 0x01ee, 0x0311, 0x0122, 0x0144, 0x0377,
72*05b00f60SXin Li 0x0289, 0x00ba, 0x00dc, 0x02ef, 0x0010, 0x0223, 0x0245, 0x0076,
73*05b00f60SXin Li 0x01b9, 0x038a, 0x03ec, 0x01df, 0x0320, 0x0113, 0x0175, 0x0346,
74*05b00f60SXin Li 0x02b8, 0x008b, 0x00ed, 0x02de, 0x0021, 0x0212, 0x0274, 0x0047,
75*05b00f60SXin Li 0x01ea, 0x03d9, 0x03bf, 0x018c, 0x0373, 0x0140, 0x0126, 0x0315,
76*05b00f60SXin Li 0x02eb, 0x00d8, 0x00be, 0x028d, 0x0072, 0x0241, 0x0227, 0x0014,
77*05b00f60SXin Li 0x01db, 0x03e8, 0x038e, 0x01bd, 0x0342, 0x0171, 0x0117, 0x0324,
78*05b00f60SXin Li 0x02da, 0x00e9, 0x008f, 0x02bc, 0x0043, 0x0270, 0x0216, 0x0025,
79*05b00f60SXin Li 0x014c, 0x037f, 0x0319, 0x012a, 0x03d5, 0x01e6, 0x0180, 0x03b3,
80*05b00f60SXin Li 0x024d, 0x007e, 0x0018, 0x022b, 0x00d4, 0x02e7, 0x0281, 0x00b2,
81*05b00f60SXin Li 0x017d, 0x034e, 0x0328, 0x011b, 0x03e4, 0x01d7, 0x01b1, 0x0382,
82*05b00f60SXin Li 0x027c, 0x004f, 0x0029, 0x021a, 0x00e5, 0x02d6, 0x02b0, 0x0083,
83*05b00f60SXin Li 0x012e, 0x031d, 0x037b, 0x0148, 0x03b7, 0x0184, 0x01e2, 0x03d1,
84*05b00f60SXin Li 0x022f, 0x001c, 0x007a, 0x0249, 0x00b6, 0x0285, 0x02e3, 0x00d0,
85*05b00f60SXin Li 0x011f, 0x032c, 0x034a, 0x0179, 0x0386, 0x01b5, 0x01d3, 0x03e0,
86*05b00f60SXin Li 0x021e, 0x002d, 0x004b, 0x0278, 0x0087, 0x02b4, 0x02d2, 0x00e1
87*05b00f60SXin Li };
88*05b00f60SXin Li
89*05b00f60SXin Li static void
init_crc10_table(void)90*05b00f60SXin Li init_crc10_table(void)
91*05b00f60SXin Li {
92*05b00f60SXin Li #define CRC10_POLYNOMIAL 0x633
93*05b00f60SXin Li int i, j;
94*05b00f60SXin Li uint16_t accum;
95*05b00f60SXin Li uint16_t verify_crc10_table[256];
96*05b00f60SXin Li
97*05b00f60SXin Li for ( i = 0; i < 256; i++ )
98*05b00f60SXin Li {
99*05b00f60SXin Li accum = ((unsigned short) i << 2);
100*05b00f60SXin Li for ( j = 0; j < 8; j++ )
101*05b00f60SXin Li {
102*05b00f60SXin Li if ((accum <<= 1) & 0x400) accum ^= CRC10_POLYNOMIAL;
103*05b00f60SXin Li }
104*05b00f60SXin Li verify_crc10_table[i] = accum;
105*05b00f60SXin Li }
106*05b00f60SXin Li assert(memcmp(verify_crc10_table,
107*05b00f60SXin Li crc10_table,
108*05b00f60SXin Li sizeof(verify_crc10_table)) == 0);
109*05b00f60SXin Li #undef CRC10_POLYNOMIAL
110*05b00f60SXin Li }
111*05b00f60SXin Li
112*05b00f60SXin Li uint16_t
verify_crc10_cksum(uint16_t accum,const u_char * p,int length)113*05b00f60SXin Li verify_crc10_cksum(uint16_t accum, const u_char *p, int length)
114*05b00f60SXin Li {
115*05b00f60SXin Li int i;
116*05b00f60SXin Li
117*05b00f60SXin Li for ( i = 0; i < length; i++ )
118*05b00f60SXin Li {
119*05b00f60SXin Li accum = ((accum << 8) & 0x3ff)
120*05b00f60SXin Li ^ crc10_table[( accum >> 2) & 0xff]
121*05b00f60SXin Li ^ *p++;
122*05b00f60SXin Li }
123*05b00f60SXin Li return accum;
124*05b00f60SXin Li }
125*05b00f60SXin Li
126*05b00f60SXin Li /* precompute checksum tables */
127*05b00f60SXin Li void
init_checksum(void)128*05b00f60SXin Li init_checksum(void) {
129*05b00f60SXin Li
130*05b00f60SXin Li init_crc10_table();
131*05b00f60SXin Li
132*05b00f60SXin Li }
133*05b00f60SXin Li
134*05b00f60SXin Li /*
135*05b00f60SXin Li * Creates the OSI Fletcher checksum. See 8473-1, Appendix C, section C.3.
136*05b00f60SXin Li * The checksum field of the passed PDU does not need to be reset to zero.
137*05b00f60SXin Li */
138*05b00f60SXin Li uint16_t
create_osi_cksum(const uint8_t * pptr,int checksum_offset,int length)139*05b00f60SXin Li create_osi_cksum (const uint8_t *pptr, int checksum_offset, int length)
140*05b00f60SXin Li {
141*05b00f60SXin Li
142*05b00f60SXin Li int x;
143*05b00f60SXin Li int y;
144*05b00f60SXin Li uint32_t mul;
145*05b00f60SXin Li uint32_t c0;
146*05b00f60SXin Li uint32_t c1;
147*05b00f60SXin Li uint16_t checksum;
148*05b00f60SXin Li int idx;
149*05b00f60SXin Li
150*05b00f60SXin Li c0 = 0;
151*05b00f60SXin Li c1 = 0;
152*05b00f60SXin Li
153*05b00f60SXin Li for (idx = 0; idx < length; idx++) {
154*05b00f60SXin Li /*
155*05b00f60SXin Li * Ignore the contents of the checksum field.
156*05b00f60SXin Li */
157*05b00f60SXin Li if (idx == checksum_offset ||
158*05b00f60SXin Li idx == checksum_offset+1) {
159*05b00f60SXin Li c1 += c0;
160*05b00f60SXin Li pptr++;
161*05b00f60SXin Li } else {
162*05b00f60SXin Li c0 = c0 + *(pptr++);
163*05b00f60SXin Li c1 += c0;
164*05b00f60SXin Li }
165*05b00f60SXin Li }
166*05b00f60SXin Li
167*05b00f60SXin Li c0 = c0 % 255;
168*05b00f60SXin Li c1 = c1 % 255;
169*05b00f60SXin Li
170*05b00f60SXin Li mul = (length - checksum_offset)*(c0);
171*05b00f60SXin Li
172*05b00f60SXin Li x = mul - c0 - c1;
173*05b00f60SXin Li y = c1 - mul - 1;
174*05b00f60SXin Li
175*05b00f60SXin Li if ( y >= 0 ) y++;
176*05b00f60SXin Li if ( x < 0 ) x--;
177*05b00f60SXin Li
178*05b00f60SXin Li x %= 255;
179*05b00f60SXin Li y %= 255;
180*05b00f60SXin Li
181*05b00f60SXin Li
182*05b00f60SXin Li if (x == 0) x = 255;
183*05b00f60SXin Li if (y == 0) y = 255;
184*05b00f60SXin Li
185*05b00f60SXin Li y &= 0x00FF;
186*05b00f60SXin Li checksum = ((x << 8) | y);
187*05b00f60SXin Li
188*05b00f60SXin Li return checksum;
189*05b00f60SXin Li }
190