opencv
1/* adler32.c -- compute the Adler-32 checksum of a data stream
2* Copyright (C) 1995-2011, 2016 Mark Adler
3* For conditions of distribution and use, see copyright notice in zlib.h
4*/
5
6#include "zbuild.h"7#include "functable.h"8#include "adler32_p.h"9
10/* ========================================================================= */
11Z_INTERNAL uint32_t adler32_c(uint32_t adler, const uint8_t *buf, size_t len) {12uint32_t sum2;13unsigned n;14
15/* split Adler-32 into component sums */16sum2 = (adler >> 16) & 0xffff;17adler &= 0xffff;18
19/* in case user likes doing a byte at a time, keep it fast */20if (UNLIKELY(len == 1))21return adler32_len_1(adler, buf, sum2);22
23/* initial Adler-32 value (deferred check for len == 1 speed) */24if (UNLIKELY(buf == NULL))25return 1L;26
27/* in case short lengths are provided, keep it somewhat fast */28if (UNLIKELY(len < 16))29return adler32_len_16(adler, buf, len, sum2);30
31/* do length NMAX blocks -- requires just one modulo operation */32while (len >= NMAX) {33len -= NMAX;34#ifdef UNROLL_MORE35n = NMAX / 16; /* NMAX is divisible by 16 */36#else37n = NMAX / 8; /* NMAX is divisible by 8 */38#endif39do {40#ifdef UNROLL_MORE41DO16(adler, sum2, buf); /* 16 sums unrolled */42buf += 16;43#else44DO8(adler, sum2, buf, 0); /* 8 sums unrolled */45buf += 8;46#endif47} while (--n);48adler %= BASE;49sum2 %= BASE;50}51
52/* do remaining bytes (less than NMAX, still just one modulo) */53return adler32_len_64(adler, buf, len, sum2);54}
55
56#ifdef ZLIB_COMPAT57unsigned long Z_EXPORT PREFIX(adler32_z)(unsigned long adler, const unsigned char *buf, size_t len) {58return (unsigned long)functable.adler32((uint32_t)adler, buf, len);59}
60#else61uint32_t Z_EXPORT PREFIX(adler32_z)(uint32_t adler, const unsigned char *buf, size_t len) {62return functable.adler32(adler, buf, len);63}
64#endif65
66/* ========================================================================= */
67#ifdef ZLIB_COMPAT68unsigned long Z_EXPORT PREFIX(adler32)(unsigned long adler, const unsigned char *buf, unsigned int len) {69return (unsigned long)functable.adler32((uint32_t)adler, buf, len);70}
71#else72uint32_t Z_EXPORT PREFIX(adler32)(uint32_t adler, const unsigned char *buf, uint32_t len) {73return functable.adler32(adler, buf, len);74}
75#endif76
77/* ========================================================================= */
78static uint32_t adler32_combine_(uint32_t adler1, uint32_t adler2, z_off64_t len2) {79uint32_t sum1;80uint32_t sum2;81unsigned rem;82
83/* for negative len, return invalid adler32 as a clue for debugging */84if (len2 < 0)85return 0xffffffff;86
87/* the derivation of this formula is left as an exercise for the reader */88len2 %= BASE; /* assumes len2 >= 0 */89rem = (unsigned)len2;90sum1 = adler1 & 0xffff;91sum2 = rem * sum1;92sum2 %= BASE;93sum1 += (adler2 & 0xffff) + BASE - 1;94sum2 += ((adler1 >> 16) & 0xffff) + ((adler2 >> 16) & 0xffff) + BASE - rem;95if (sum1 >= BASE) sum1 -= BASE;96if (sum1 >= BASE) sum1 -= BASE;97if (sum2 >= ((unsigned long)BASE << 1)) sum2 -= ((unsigned long)BASE << 1);98if (sum2 >= BASE) sum2 -= BASE;99return sum1 | (sum2 << 16);100}
101
102/* ========================================================================= */
103#ifdef ZLIB_COMPAT104unsigned long Z_EXPORT PREFIX(adler32_combine)(unsigned long adler1, unsigned long adler2, z_off_t len2) {105return (unsigned long)adler32_combine_((uint32_t)adler1, (uint32_t)adler2, len2);106}
107
108unsigned long Z_EXPORT PREFIX4(adler32_combine)(unsigned long adler1, unsigned long adler2, z_off64_t len2) {109return (unsigned long)adler32_combine_((uint32_t)adler1, (uint32_t)adler2, len2);110}
111#else112uint32_t Z_EXPORT PREFIX4(adler32_combine)(uint32_t adler1, uint32_t adler2, z_off64_t len2) {113return adler32_combine_(adler1, adler2, len2);114}
115#endif116