checksum.py 5.4 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153
  1. import array
  2. import socket
  3. import struct
  4. # avoid references for performance reasons
  5. unpack = struct.unpack
  6. unpack_word_be = struct.Struct(">H").unpack
  7. array_call = array.array
  8. ntohs = socket.ntohs
  9. # TCP (RFC 793) and UDP (RFC 768) checksum
  10. def in_cksum_add(s, buf):
  11. """Add checksum value to the given value s."""
  12. n = len(buf)
  13. # logger.debug("buflen for checksum: %d" % n)
  14. cnt = int(n / 2) * 2
  15. # logger.debug("slicing at: %d, %s" % (cnt, type(cnt)))
  16. a = array_call("H", buf[:cnt])
  17. # logger.debug("2-byte values: %s" % a)
  18. # logger.debug(buf[-1].to_bytes(1, byteorder='big'))
  19. if cnt != n:
  20. #a.append(unpack_word_be( buf[-1].to_bytes(1, byteorder="big") + b"\x00" )[0])
  21. a.append(unpack_word_be( buf[-1:] + b"\x00" )[0])
  22. return s + sum(a)
  23. def in_cksum_done(s):
  24. """Complete checksum building."""
  25. # add carry to sum itself
  26. s = (s >> 16) + (s & 0xffff)
  27. s += (s >> 16)
  28. # return complement of sums
  29. return ntohs(~s & 0xffff)
  30. def in_cksum(buf):
  31. """Return computed Internet Protocol checksum."""
  32. return in_cksum_done(in_cksum_add(0, buf))
  33. # CRC-32C Checksum
  34. # http://tools.ietf.org/html/rfc3309
  35. crc32c_table = (
  36. 0x00000000, 0xF26B8303, 0xE13B70F7, 0x1350F3F4, 0xC79A971F,
  37. 0x35F1141C, 0x26A1E7E8, 0xD4CA64EB, 0x8AD958CF, 0x78B2DBCC,
  38. 0x6BE22838, 0x9989AB3B, 0x4D43CFD0, 0xBF284CD3, 0xAC78BF27,
  39. 0x5E133C24, 0x105EC76F, 0xE235446C, 0xF165B798, 0x030E349B,
  40. 0xD7C45070, 0x25AFD373, 0x36FF2087, 0xC494A384, 0x9A879FA0,
  41. 0x68EC1CA3, 0x7BBCEF57, 0x89D76C54, 0x5D1D08BF, 0xAF768BBC,
  42. 0xBC267848, 0x4E4DFB4B, 0x20BD8EDE, 0xD2D60DDD, 0xC186FE29,
  43. 0x33ED7D2A, 0xE72719C1, 0x154C9AC2, 0x061C6936, 0xF477EA35,
  44. 0xAA64D611, 0x580F5512, 0x4B5FA6E6, 0xB93425E5, 0x6DFE410E,
  45. 0x9F95C20D, 0x8CC531F9, 0x7EAEB2FA, 0x30E349B1, 0xC288CAB2,
  46. 0xD1D83946, 0x23B3BA45, 0xF779DEAE, 0x05125DAD, 0x1642AE59,
  47. 0xE4292D5A, 0xBA3A117E, 0x4851927D, 0x5B016189, 0xA96AE28A,
  48. 0x7DA08661, 0x8FCB0562, 0x9C9BF696, 0x6EF07595, 0x417B1DBC,
  49. 0xB3109EBF, 0xA0406D4B, 0x522BEE48, 0x86E18AA3, 0x748A09A0,
  50. 0x67DAFA54, 0x95B17957, 0xCBA24573, 0x39C9C670, 0x2A993584,
  51. 0xD8F2B687, 0x0C38D26C, 0xFE53516F, 0xED03A29B, 0x1F682198,
  52. 0x5125DAD3, 0xA34E59D0, 0xB01EAA24, 0x42752927, 0x96BF4DCC,
  53. 0x64D4CECF, 0x77843D3B, 0x85EFBE38, 0xDBFC821C, 0x2997011F,
  54. 0x3AC7F2EB, 0xC8AC71E8, 0x1C661503, 0xEE0D9600, 0xFD5D65F4,
  55. 0x0F36E6F7, 0x61C69362, 0x93AD1061, 0x80FDE395, 0x72966096,
  56. 0xA65C047D, 0x5437877E, 0x4767748A, 0xB50CF789, 0xEB1FCBAD,
  57. 0x197448AE, 0x0A24BB5A, 0xF84F3859, 0x2C855CB2, 0xDEEEDFB1,
  58. 0xCDBE2C45, 0x3FD5AF46, 0x7198540D, 0x83F3D70E, 0x90A324FA,
  59. 0x62C8A7F9, 0xB602C312, 0x44694011, 0x5739B3E5, 0xA55230E6,
  60. 0xFB410CC2, 0x092A8FC1, 0x1A7A7C35, 0xE811FF36, 0x3CDB9BDD,
  61. 0xCEB018DE, 0xDDE0EB2A, 0x2F8B6829, 0x82F63B78, 0x709DB87B,
  62. 0x63CD4B8F, 0x91A6C88C, 0x456CAC67, 0xB7072F64, 0xA457DC90,
  63. 0x563C5F93, 0x082F63B7, 0xFA44E0B4, 0xE9141340, 0x1B7F9043,
  64. 0xCFB5F4A8, 0x3DDE77AB, 0x2E8E845F, 0xDCE5075C, 0x92A8FC17,
  65. 0x60C37F14, 0x73938CE0, 0x81F80FE3, 0x55326B08, 0xA759E80B,
  66. 0xB4091BFF, 0x466298FC, 0x1871A4D8, 0xEA1A27DB, 0xF94AD42F,
  67. 0x0B21572C, 0xDFEB33C7, 0x2D80B0C4, 0x3ED04330, 0xCCBBC033,
  68. 0xA24BB5A6, 0x502036A5, 0x4370C551, 0xB11B4652, 0x65D122B9,
  69. 0x97BAA1BA, 0x84EA524E, 0x7681D14D, 0x2892ED69, 0xDAF96E6A,
  70. 0xC9A99D9E, 0x3BC21E9D, 0xEF087A76, 0x1D63F975, 0x0E330A81,
  71. 0xFC588982, 0xB21572C9, 0x407EF1CA, 0x532E023E, 0xA145813D,
  72. 0x758FE5D6, 0x87E466D5, 0x94B49521, 0x66DF1622, 0x38CC2A06,
  73. 0xCAA7A905, 0xD9F75AF1, 0x2B9CD9F2, 0xFF56BD19, 0x0D3D3E1A,
  74. 0x1E6DCDEE, 0xEC064EED, 0xC38D26C4, 0x31E6A5C7, 0x22B65633,
  75. 0xD0DDD530, 0x0417B1DB, 0xF67C32D8, 0xE52CC12C, 0x1747422F,
  76. 0x49547E0B, 0xBB3FFD08, 0xA86F0EFC, 0x5A048DFF, 0x8ECEE914,
  77. 0x7CA56A17, 0x6FF599E3, 0x9D9E1AE0, 0xD3D3E1AB, 0x21B862A8,
  78. 0x32E8915C, 0xC083125F, 0x144976B4, 0xE622F5B7, 0xF5720643,
  79. 0x07198540, 0x590AB964, 0xAB613A67, 0xB831C993, 0x4A5A4A90,
  80. 0x9E902E7B, 0x6CFBAD78, 0x7FAB5E8C, 0x8DC0DD8F, 0xE330A81A,
  81. 0x115B2B19, 0x020BD8ED, 0xF0605BEE, 0x24AA3F05, 0xD6C1BC06,
  82. 0xC5914FF2, 0x37FACCF1, 0x69E9F0D5, 0x9B8273D6, 0x88D28022,
  83. 0x7AB90321, 0xAE7367CA, 0x5C18E4C9, 0x4F48173D, 0xBD23943E,
  84. 0xF36E6F75, 0x0105EC76, 0x12551F82, 0xE03E9C81, 0x34F4F86A,
  85. 0xC69F7B69, 0xD5CF889D, 0x27A40B9E, 0x79B737BA, 0x8BDCB4B9,
  86. 0x988C474D, 0x6AE7C44E, 0xBE2DA0A5, 0x4C4623A6, 0x5F16D052,
  87. 0xAD7D5351
  88. )
  89. def crc32_add(crc, buf):
  90. # buf = array.array("B", buf)
  91. i = 0
  92. while i < len(buf):
  93. # crc = (crc >> 8) ^ crc32c_table[(crc ^ b) & 0xff]
  94. crc = (crc >> 8) ^ crc32c_table[(crc ^ buf[i]) & 0xff]
  95. i += 1
  96. return crc
  97. def crc32_done(crc):
  98. tmp = ~crc & 0xffffffff
  99. b0 = tmp & 0xff
  100. b1 = (tmp >> 8) & 0xff
  101. b2 = (tmp >> 16) & 0xff
  102. b3 = (tmp >> 24) & 0xff
  103. crc = (b0 << 24) | (b1 << 16) | (b2 << 8) | b3
  104. return crc
  105. def crc32_cksum(buf):
  106. """Return computed CRC-32c checksum."""
  107. return crc32_done(crc32_add(0xffffffff, buf))
  108. def fletcher32(data_to_checksum, amount_words):
  109. # 1 word = 2 Bytes
  110. sum1 = 0xffff
  111. sum2 = 0xffff
  112. datapos = 0
  113. while amount_words > 0:
  114. tlen = 359 if amount_words > 359 else amount_words
  115. amount_words -= tlen
  116. while tlen > 0:
  117. # sum1 += unpack_word_be(data_to_checksum[datapos:datapos+2])[0]
  118. # print("%d" % sum1)
  119. sum1 += unpack_word_be(data_to_checksum[datapos: datapos + 2])[0]
  120. datapos += 2
  121. sum2 += sum1
  122. # print("%d" % sum1)
  123. # print("%d" % sum2)
  124. # print("--")
  125. tlen -= 1
  126. sum1 = (sum1 & 0xffff) + (sum1 >> 16)
  127. sum2 = (sum2 & 0xffff) + (sum2 >> 16)
  128. # Second reduction step to reduce sums to 16 bits
  129. sum1 = (sum1 & 0xffff) + (sum1 >> 16)
  130. sum2 = (sum2 & 0xffff) + (sum2 >> 16)
  131. return (sum2 << 16) | sum1