src/TortoiseMerge/libsvn_diff/adler32.c

   1 /* adler32.c -- compute the Adler-32 checksum of a data stream\r
   2  * Copyright (C) 1995-2002 Mark Adler\r
   3  * For conditions of distribution and use, see copyright notice in zlib.h \r
   4  */\r
   5 \r
   6 /* @(#) $Id: adler32.c,v 1.1 2003/05/23 21:05:30 steveking Exp $ */\r
   7 \r
   8 #include "zlib.h"\r
   9 \r
  10 #define BASE 65521L /* largest prime smaller than 65536 */\r
  11 #define NMAX 5552\r
  12 /* NMAX is the largest n such that 255n(n+1)/2 + (n+1)(BASE-1) <= 2^32-1 */\r
  13 \r
  14 #define DO1(buf,i)  {s1 += buf[i]; s2 += s1;}\r
  15 #define DO2(buf,i)  DO1(buf,i); DO1(buf,i+1);\r
  16 #define DO4(buf,i)  DO2(buf,i); DO2(buf,i+2);\r
  17 #define DO8(buf,i)  DO4(buf,i); DO4(buf,i+4);\r
  18 #define DO16(buf)   DO8(buf,0); DO8(buf,8);\r
  19 \r
  20 /* ========================================================================= */\r
  21 uLong ZEXPORT adler32(adler, buf, len)\r
  22     uLong adler;\r
  23     const Bytef *buf;\r
  24     uInt len;\r
  25 {\r
  26     unsigned long s1 = adler & 0xffff;\r
  27     unsigned long s2 = (adler >> 16) & 0xffff;\r
  28     int k;\r
  29 \r
  30     if (buf == Z_NULL) return 1L;\r
  31 \r
  32     while (len > 0) {\r
  33         k = len < NMAX ? len : NMAX;\r
  34         len -= k;\r
  35         while (k >= 16) {\r
  36             DO16(buf);\r
  37             buf += 16;\r
  38             k -= 16;\r
  39         }\r
  40         if (k != 0) do {\r
  41             s1 += *buf++;\r
  42             s2 += s1;\r
  43         } while (--k);\r
  44         s1 %= BASE;\r
  45         s2 %= BASE;\r
  46     }\r
  47     return (s2 << 16) | s1;\r
  48 }\r