WebSVN – Kolibri OS – Diff – /programs/fs/kfar/trunk/zlib/trees.asm


; trees.asm -- output deflated data using Huffman coding
; Copyright (C) 1995-2012 Jean-loup Gailly
; detect_data_type() function provided freely by Cosmin Truta, 2006
; For conditions of distribution and use, see copyright notice in zlib.h
 
;  ALGORITHM
 
;      The "deflation" process uses several Huffman trees. The more
;      common source values are represented by shorter bit sequences.
 
;      Each code tree is stored in a compressed form which is itself
; a Huffman encoding of the lengths of all the code strings (in
; ascending order by source values).  The actual code strings are
; reconstructed from the lengths in the inflate process, as described
; in the deflate specification.
 
;  REFERENCES
 
;      Deutsch, L.P.,"'Deflate' Compressed Data Format Specification".
;      Available in ftp.uu.net:/pub/archiving/zip/doc/deflate-1.1.doc
 
;      Storer, James A.
;          Data Compression:  Methods and Theory, pp. 49-50.
;          Computer Science Press, 1988.  ISBN 0-7167-8156-5.
 
;      Sedgewick, R.
;          Algorithms, p290.
;          Addison-Wesley, 1983. ISBN 0-201-06672-6.
 
; ===========================================================================
; Constants
 
 
MAX_BL_BITS equ 7
; Bit length codes must not exceed MAX_BL_BITS bits
 
END_BLOCK equ 256
; end of block literal code
 
REP_3_6     equ 16
; repeat previous bit length 3-6 times (2 bits of repeat count)
 
REPZ_3_10   equ 17
; repeat a zero length 3-10 times  (3 bits of repeat count)
 
REPZ_11_138 equ 18
; repeat a zero length 11-138 times  (7 bits of repeat count)
 
align 4
extra_lbits dd \ ;int [LENGTH_CODES] ;extra bits for each length code
	0,0,0,0,0,0,0,0,1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4,5,5,5,5,0
 
align 4
extra_dbits dd \ ;int [D_CODES] ;extra bits for each distance code
	0,0,0,0,1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,11,11,12,12,13,13
 
align 4
extra_blbits dd \ ;int [BL_CODES] ;extra bits for each bit length code
	0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,3,7
 
align 4
bl_order db 16,17,18,0,8,7,9,6,10,5,11,4,12,3,13,2,14,1,15
; The lengths of the bit length codes are sent in order of decreasing
; probability, to avoid transmitting the lengths for unused bit length codes.
 
 
; ===========================================================================
; Local data. These are initialized only once.
 
 
DIST_CODE_LEN equ 512 ;see definition of array dist_code below
 
if GEN_TREES_H eq 1 ;| !(STDC)
; non ANSI compilers may not accept trees.inc
 
align 4
static_ltree rb sizeof.ct_data * (L_CODES+2)
; The static literal tree. Since the bit lengths are imposed, there is no
; need for the L_CODES extra codes used during heap construction. However
; The codes 286 and 287 are needed to build a canonical tree (see _tr_init
; below).
 
align 4
static_dtree rb sizeof.ct_data * D_CODES
; The static distance tree. (Actually a trivial tree since all codes use
; 5 bits.)
 
align 4
_dist_code rb DIST_CODE_LEN ;uch[]
; Distance codes. The first 256 values correspond to the distances
; 3 .. 258, the last 256 values correspond to the top 8 bits of
; the 15 bit distances.
 
align 4
_length_code rb MAX_MATCH-MIN_MATCH+1 ;uch[]
; length code for each normalized match length (0 == MIN_MATCH)
 
align 4
base_length rd LENGTH_CODES ;int[]
; First normalized length for each code (0 = MIN_MATCH)
 
align 4
base_dist rd D_CODES ;int[]
; First normalized distance for each code (0 = distance of 1)
 
else
include 'trees.inc'
end if ;GEN_TREES_H
 
struct static_tree_desc ;_s
	static_tree dd ? ;const ct_data * ;static tree or NULL
	extra_bits  dd ? ;const intf * ;extra bits for each code or NULL
	extra_base  dd ? ;int ;base index for extra_bits
	elems       dd ? ;int ;max number of elements in the tree
	max_length  dd ? ;int ;max bit length for the codes
ends
 
align 4
static_l_desc static_tree_desc static_ltree, extra_lbits, LITERALS+1, L_CODES, MAX_BITS
 
align 4
static_d_desc static_tree_desc static_dtree, extra_dbits, 0, D_CODES, MAX_BITS
 
align 4
static_bl_desc static_tree_desc 0, extra_blbits, 0, BL_CODES, MAX_BL_BITS
 
; ===========================================================================
; Local (static) routines in this file.
 
 
macro send_code s, c, tree
{
if DEBUG eq 1
;	if (z_verbose>2) fprintf(stderr,"\ncd %3d ",(c))
end if
push eax ebx
if c eq eax
else
	mov eax,c
end if
	imul eax,sizeof.ct_data
	add eax,tree
	movzx ebx,word[eax+Len]
	push ebx
	movzx ebx,word[eax+Code]
	push ebx
	stdcall send_bits, s ;tree[c].Code, tree[c].Len
pop ebx eax
}
; Send a code of the given tree[c] and tree must not have side effects
 
; ===========================================================================
; Output a short LSB first on the stream.
; IN assertion: there is enough room in pendingBuf.
 
macro put_short s, w
{
	mov eax,[s+deflate_state.pending]
	add eax,[s+deflate_state.pending_buf]
	mov word[eax],w
	add dword[s+deflate_state.pending],2
}
 
; ===========================================================================
; Send a value on a given number of bits.
; IN assertion: length <= 16 and value fits in length bits.
 
;void (s, value, length)
;    deflate_state* s
;    int value  ;value to send
;    int length ;number of bits
align 4
proc send_bits uses eax ecx edi, s:dword, value:dword, length:dword
;    Tracevv((stderr," l %2d v %4x ", length, value));
	zlib_debug 'send_bits value = %d',[value]
;if DEBUG eq 1
	mov eax,[length]
	cmp eax,0
	jle @f
	cmp eax,15
	jle .end1
	@@:
		zlib_assert 'invalid length' ;Assert(..>0 && ..<=15)
	.end1:
	mov edi,[s]
	add [edi+deflate_state.bits_sent],eax
 
	; If not enough room in bi_buf, use (valid) bits from bi_buf and
	; (16 - bi_valid) bits from value, leaving (width - (16-bi_valid))
	; unused bits in value.
 
	mov ecx,Buf_size
	sub ecx,eax
	cmp [edi+deflate_state.bi_valid],ecx
	jle @f ;if (..>..)
		mov eax,[value]
		mov ecx,[edi+deflate_state.bi_valid]
		shl eax,cl
		or [edi+deflate_state.bi_buf],ax
		mov cx,[edi+deflate_state.bi_buf]
		put_short edi, cx
		mov eax,[value]
		mov ecx,Buf_size
		sub ecx,[edi+deflate_state.bi_valid]
		shr eax,cl
		mov [edi+deflate_state.bi_buf],ax
		mov eax,[length]
		sub eax,Buf_size
		jmp .end0
	@@: ;else
		mov eax,[value]
		mov ecx,[edi+deflate_state.bi_valid]
		shl eax,cl
		or [edi+deflate_state.bi_buf],ax
		mov eax,[length]
	.end0:
	add [edi+deflate_state.bi_valid],eax
;else ;!DEBUG
 
;{ int len = length;
;  if (s->bi_valid > (int)Buf_size - len) {
;    int val = value;
;    s->bi_buf |= (uint_16)val << s->bi_valid;
;    put_short(s, s->bi_buf);
;    s->bi_buf = (uint_16)val >> (Buf_size - s->bi_valid);
;    s->bi_valid += len - Buf_size;
;  } else {
;    s->bi_buf |= (uint_16)(value) << s->bi_valid;
;    s->bi_valid += len;
;  }
;}
;end if ;DEBUG
	ret
endp
 
; the arguments must not have side effects
 
; ===========================================================================
; Initialize the various 'constant' tables.
 
;int static_init_done = 0
 
;void ()
align 4
proc tr_static_init
if GEN_TREES_H eq 1
 
;    int n      ;iterates over tree elements
;    int bits   ;bit counter
;    int length ;length value
;    int code   ;code value
;    int dist   ;distance index
;    uint_16 bl_count[MAX_BITS+1];
	; number of codes at each bit length for an optimal tree
 
;    if (static_init_done) return;
 
	; For some embedded targets, global variables are not initialized:
;if NO_INIT_GLOBAL_POINTERS
;    static_l_desc.static_tree = static_ltree;
;    static_l_desc.extra_bits = extra_lbits;
;    static_d_desc.static_tree = static_dtree;
;    static_d_desc.extra_bits = extra_dbits;
;    static_bl_desc.extra_bits = extra_blbits;
;end if
 
	; Initialize the mapping length (0..255) -> length code (0..28)
;    length = 0;
;    for (code = 0; code < LENGTH_CODES-1; code++) {
;        base_length[code] = length;
;        for (n = 0; n < (1<
;            _length_code[length++] = (uch)code;
;        }
;    }
;    Assert (length == 256, "tr_static_init: length != 256");
	; Note that the length 255 (match length 258) can be represented
	; in two different ways: code 284 + 5 bits or code 285, so we
	; overwrite length_code[255] to use the best encoding:
 
;    _length_code[length-1] = (uch)code;
 
	; Initialize the mapping dist (0..32K) -> dist code (0..29)
;    dist = 0;
;    for (code = 0 ; code < 16; code++) {
;        base_dist[code] = dist;
;        for (n = 0; n < (1<
;            _dist_code[dist++] = (uch)code;
;        }
;    }
;    Assert (dist == 256, "tr_static_init: dist != 256");
;    dist >>= 7; /* from now on, all distances are divided by 128 */
;    for ( ; code < D_CODES; code++) {
;        base_dist[code] = dist << 7;
;        for (n = 0; n < (1<<(extra_dbits[code]-7)); n++) {
;            _dist_code[256 + dist++] = (uch)code;
;        }
;    }
;    Assert (dist == 256, "tr_static_init: 256+dist != 512");
 
	; Construct the codes of the static literal tree
;    for (bits = 0; bits <= MAX_BITS; bits++) bl_count[bits] = 0;
;    n = 0;
;    while (n <= 143) static_ltree[n++].Len = 8, bl_count[8]++;
;    while (n <= 255) static_ltree[n++].Len = 9, bl_count[9]++;
;    while (n <= 279) static_ltree[n++].Len = 7, bl_count[7]++;
;    while (n <= 287) static_ltree[n++].Len = 8, bl_count[8]++;
	; Codes 286 and 287 do not exist, but we must include them in the
	; tree construction to get a canonical Huffman tree (longest code
	; all ones)
 
;    gen_codes((ct_data *)static_ltree, L_CODES+1, bl_count);
 
	; The static distance tree is trivial:
;    for (n = 0; n < D_CODES; n++) {
;        static_dtree[n].Len = 5;
;        static_dtree[n].Code = bi_reverse((unsigned)n, 5);
;    }
;    static_init_done = 1;
 
if GEN_TREES_H eq 1
	call gen_trees_header
end if
end if ;(GEN_TREES_H) | !(STDC)
	ret
endp
 
; ===========================================================================
; Genererate the file trees.h describing the static trees.
 
;#  define SEPARATOR(i, last, width) \
;      ((i) == (last)? "\n};\n\n" :    \
;       ((i) % (width) == (width)-1 ? ",\n" : ", "))
 
;void ()
align 4
proc gen_trees_header
;    FILE *header = fopen("trees.inc", "w");
;    int i;
 
;    Assert (header != NULL, "Can't open trees.inc");
;    fprintf(header,
;            "/* header created automatically with -DGEN_TREES_H */\n\n");
 
;    fprintf(header, "local const ct_data static_ltree[L_CODES+2] = {\n");
;    for (i = 0; i < L_CODES+2; i++) {
;        fprintf(header, "{{%3u},{%3u}}%s", static_ltree[i].Code,
;                static_ltree[i].Len, SEPARATOR(i, L_CODES+1, 5));
;    }
 
;    fprintf(header, "local const ct_data static_dtree[D_CODES] = {\n");
;    for (i = 0; i < D_CODES; i++) {
;        fprintf(header, "{{%2u},{%2u}}%s", static_dtree[i].Code,
;                static_dtree[i].Len, SEPARATOR(i, D_CODES-1, 5));
;    }
 
;    fprintf(header, "const uch ZLIB_INTERNAL _dist_code[DIST_CODE_LEN] = {\n");
;    for (i = 0; i < DIST_CODE_LEN; i++) {
;        fprintf(header, "%2u%s", _dist_code[i],
;                SEPARATOR(i, DIST_CODE_LEN-1, 20));
;    }
 
;    fprintf(header,
;        "const uch ZLIB_INTERNAL _length_code[MAX_MATCH-MIN_MATCH+1]= {\n");
;    for (i = 0; i < MAX_MATCH-MIN_MATCH+1; i++) {
;        fprintf(header, "%2u%s", _length_code[i],
;                SEPARATOR(i, MAX_MATCH-MIN_MATCH, 20));
;    }
 
;    fprintf(header, "local const int base_length[LENGTH_CODES] = {\n");
;    for (i = 0; i < LENGTH_CODES; i++) {
;        fprintf(header, "%1u%s", base_length[i],
;                SEPARATOR(i, LENGTH_CODES-1, 20));
;    }
 
;    fprintf(header, "local const int base_dist[D_CODES] = {\n");
;    for (i = 0; i < D_CODES; i++) {
;        fprintf(header, "%5u%s", base_dist[i],
;                SEPARATOR(i, D_CODES-1, 10));
;    }
 
;    fclose(header);
	ret
endp
 
; ===========================================================================
; Initialize the tree data structures for a new zlib stream.
 
;void (s)
;    deflate_state* s
align 4
proc _tr_init uses eax edi, s:dword
	mov edi,[s]
	zlib_debug '_tr_init'
	call tr_static_init
 
	mov eax,edi
	add eax,deflate_state.dyn_ltree
	mov [edi+deflate_state.l_desc.dyn_tree],eax
	mov [edi+deflate_state.l_desc.stat_desc],static_l_desc
 
	add eax,deflate_state.dyn_dtree-deflate_state.dyn_ltree
	mov [edi+deflate_state.d_desc.dyn_tree],eax
	mov [edi+deflate_state.d_desc.stat_desc],static_d_desc
 
	add eax,deflate_state.bl_tree-deflate_state.dyn_dtree
	mov [edi+deflate_state.bl_desc.dyn_tree],eax
	mov [edi+deflate_state.bl_desc.stat_desc],static_bl_desc;
 
	mov word[edi+deflate_state.bi_buf],0
	mov dword[edi+deflate_state.bi_valid],0
if DEBUG eq 1
	mov dword[edi+deflate_state.compressed_len],0
	mov dword[edi+deflate_state.bits_sent],0
end if
 
	; Initialize the first block of the first file:
	stdcall init_block,edi
	ret
endp
 
; ===========================================================================
; Initialize a new block.
 
;void (s)
;    deflate_state* s
align 4
proc init_block uses eax ecx edi, s:dword
	mov edi,[s]
 
	; Initialize the trees.
	mov eax,edi
	add eax,deflate_state.dyn_ltree+Freq
	mov ecx,L_CODES
	@@:
		mov word[eax],0
		add eax,sizeof.ct_data
		loop @b
	mov eax,edi
	add eax,deflate_state.dyn_dtree+Freq
	mov ecx,D_CODES
	@@:
		mov word[eax],0
		add eax,sizeof.ct_data
		loop @b
	mov eax,edi
	add eax,deflate_state.bl_tree+Freq
	mov ecx,BL_CODES
	@@:
		mov word[eax],0
		add eax,sizeof.ct_data
		loop @b
 
	mov word[edi+sizeof.ct_data*END_BLOCK+deflate_state.dyn_ltree+Freq],1
	mov dword[edi+deflate_state.static_len],0
	mov dword[edi+deflate_state.opt_len],0
	mov dword[edi+deflate_state.matches],0
	mov dword[edi+deflate_state.last_lit],0
	ret
endp
 
SMALLEST equ 1
; Index within the heap array of least frequent node in the Huffman tree
 
 
; ===========================================================================
; Remove the smallest element from the heap and recreate the heap with
; one less element. Updates heap and heap_len.
 
macro pqremove s, tree, top
{
	mov eax,s
	add eax,deflate_state.heap+2*SMALLEST
	movzx top,word[eax]
push ebx
	mov ebx,[s+deflate_state.heap_len]
	mov bx,[s+deflate_state.heap+2*ebx]
	mov word[eax],bx
	dec dword[s+deflate_state.heap_len]
pop ebx
	stdcall pqdownheap, s, tree, SMALLEST
}
 
; ===========================================================================
; Compares to subtrees, using the tree depth as tie breaker when
; the subtrees have equal frequency. This minimizes the worst case length.
 
macro smaller tree, n, m, depth, m_end
{
;if (..<.. || (..==.. && depth[n] <= depth[m]))
local .end0
	movzx eax,n
	imul eax,sizeof.ct_data
	add eax,tree
	mov ax,word[eax+Freq]
	movzx ebx,m
	imul ebx,sizeof.ct_data
	add ebx,tree
	mov bx,word[ebx+Freq]
	cmp ax,bx
	jl .end0
	jne m_end
	movzx eax,n
	add eax,depth
	mov al,byte[eax]
	movzx ebx,m
	add ebx,depth
	mov bl,byte[ebx]
	cmp al,bl
	jg m_end
	.end0:
}
 
; ===========================================================================
; Restore the heap property by moving down the tree starting at node k,
; exchanging a node with the smallest of its two sons if necessary, stopping
; when the heap property is re-established (each father smaller than its
; two sons).
 
;void (s, tree, k)
;    deflate_state* s
;    ct_data* tree ;the tree to restore
;    int      k    ;node to move down
align 4
proc pqdownheap, s:dword, tree:dword, k:dword
locals
	v dw ?
endl
pushad
	mov edi,[s]
	mov eax,[k]
	zlib_debug 'pqdownheap k = %d',eax
	mov esi,eax
	shl esi,1
	mov ax,[edi+deflate_state.heap+2*eax]
	mov [v],ax
	;esi = j ;left son of k
	.cycle0: ;while (..<=..)
		cmp esi,[edi+deflate_state.heap_len]
		jg .cycle0end
		; Set j to the smallest of the two sons:
		;;cmp esi,[edi+deflate_state.heap_len]
		jge .end1 ;if (..<.. &&
		mov ecx,edi
		add ecx,deflate_state.depth
		mov edx,esi
		shl edx,1
		add edx,edi
		add edx,deflate_state.heap
		smaller [tree], word[edx+2], word[edx], ecx, .end1
			inc esi
		.end1:
		; Exit if v is smaller than both sons
		mov ecx,edi
		add ecx,deflate_state.depth
		mov dx,[edi+deflate_state.heap+2*esi]
		smaller [tree], [v], dx, ecx, .end2
			jmp .cycle0end ;break
		.end2:
		; Exchange v with the smallest son
		mov dx,[edi+deflate_state.heap+2*esi]
		mov eax,[k]
		mov [edi+deflate_state.heap+2*eax],dx
		mov [k],esi
		; And continue down the tree, setting j to the left son of k
		shl esi,1
		jmp .cycle0
	.cycle0end:
	mov eax,[k]
	mov bx,[v]
	mov [edi+deflate_state.heap+2*eax],bx
popad
	ret
endp
 
; ===========================================================================
; Compute the optimal bit lengths for a tree and update the total bit length
; for the current block.
; IN assertion: the fields freq and dad are set, heap[heap_max] and
;    above are the tree nodes sorted by increasing frequency.
; OUT assertions: the field len is set to the optimal bit length, the
;     array bl_count contains the frequencies for each bit length.
;     The length opt_len is updated; static_len is also updated if stree is
;     not null.
 
;void (s, desc)
;    deflate_state* s
;    tree_desc* desc ;the tree descriptor
align 4
proc gen_bitlen, s:dword, desc:dword
locals
	tree  dd ? ;ct_data* ;= desc.dyn_tree
	max_code dd ? ;int   ;= desc.max_code
	stree dd ? ;ct_data* ;= desc.stat_desc.static_tree
	extra dd ? ;intf*    ;= desc.stat_desc.extra_bits
	base  dd ? ;int      ;= desc.stat_desc.extra_base
	max_length dd ? ;int ;= desc.stat_desc.max_length
	h     dd ? ;int ;heap index
	m     dd ? ;int ;iterate over the tree elements
	bits  dd ? ;int ;bit length
	xbits dd ? ;int ;extra bits
	f     dw ? ;uint_16 ;frequency
	overflow dd 0 ;int ;number of elements with bit length too large
endl
pushad
	zlib_debug 'gen_bitlen'
	mov edi,[s]
	mov edx,[desc]
	mov eax,[edx+tree_desc.dyn_tree]
	mov [tree],eax
	mov eax,[edx+tree_desc.max_code]
	mov [max_code],eax
	mov ebx,[edx+tree_desc.stat_desc]
	mov eax,[ebx+static_tree_desc.static_tree]
	mov [stree],eax
	mov eax,[ebx+static_tree_desc.extra_bits]
	mov [extra],eax
	mov eax,[ebx+static_tree_desc.extra_base]
	mov [base],eax
	mov eax,[ebx+static_tree_desc.max_length]
	mov [max_length],eax
 
	xor ecx,ecx
	.cycle0:
	cmp ecx,MAX_BITS
	jg .cycle0end ;for (..;..<=..;..)
		mov word[edi+deflate_state.bl_count+2*ecx],0
		inc ecx
		jmp .cycle0
align 4
	.cycle0end:
 
	; In a first pass, compute the optimal bit lengths (which may
	; overflow in the case of the bit length tree).
 
	mov eax,[edi+deflate_state.heap_max]
	movzx eax,word[edi+deflate_state.heap+2*eax]
	imul eax,sizeof.ct_data
	add eax,[tree]
	mov word[eax+Len],0 ;root of the heap
 
	mov eax,[edi+deflate_state.heap_max]
	inc eax
	mov [h],eax
	.cycle1:
	cmp dword[h],HEAP_SIZE
	jge .cycle1end ;for (..;..<..;..)
		mov eax,[h]
		movzx ecx,word[edi+deflate_state.heap+2*eax]
		;ecx = n
		mov eax,sizeof.ct_data
		imul eax,ecx
		add eax,[tree]
		movzx eax,word[eax+Dad]
		imul eax,sizeof.ct_data
		add eax,[tree]
		movzx eax,word[eax+Len]
		inc eax
		mov [bits],eax ;bits = tree[tree[n].Dad].Len + 1
		mov eax,[max_length]
		cmp [bits],eax
		jle @f ;if (..>..)
			mov [bits],eax
			inc dword[overflow]
		@@:
		mov esi,[bits]
		mov eax,sizeof.ct_data
		imul eax,ecx
		add eax,[tree]
		mov word[eax+Len],si
		; We overwrite tree[n].Dad which is no longer needed
 
		cmp ecx,[max_code]
		jle @f
			inc dword[h]
			jmp .cycle1 ;if (..>..) continue ;not a leaf node
		@@:
 
		mov eax,[bits]
		shl eax,1 ;*= sizeof.uint_16
		inc word[eax+edi+deflate_state.bl_count]
		mov dword[xbits],0
		cmp ecx,[base]
		jl @f ;if (..>=..)
			mov eax,ecx
			sub eax,[base]
			shl eax,2 ;*= sizeof.dd
			add eax,[extra]
			mov eax,[eax]
			mov [xbits],eax
		@@:
		mov eax,sizeof.ct_data
		imul eax,ecx
		add eax,[tree]
		movzx eax,word[eax+Freq]
		mov [f],ax
		mov esi,[bits]
		add esi,[xbits]
		imul eax,esi
		add [edi+deflate_state.opt_len],eax
		cmp dword[stree],0
		je @f ;if (..)
			movzx eax,word[f]
			mov esi,sizeof.ct_data
			imul esi,ecx
			add esi,[tree]
			movzx esi,word[esi+Len]
			add esi,[xbits]
			imul eax,esi
			add [edi+deflate_state.static_len],eax
		@@:
		inc dword[h]
		jmp .cycle1
align 4
	.cycle1end:
	cmp dword[overflow],0
	je .end_f ;if (..==0) return
 
;    Trace((stderr,"\nbit length overflow\n"));
	; This happens for example on obj2 and pic of the Calgary corpus
 
	; Find the first bit length which could increase:
	.cycle2: ;do
		mov eax,[max_length]
		dec eax
		mov [bits],eax
		shl eax,1 ;*= sizeof.dw
		add eax,edi
		add eax,deflate_state.bl_count
		@@:
		cmp word[eax],0
		jne @f ;while (..==0) bits--
			dec dword[bits]
			sub eax,2
			jmp @b
align 4
		@@:
		dec word[eax]     ;move one leaf down the tree
		add word[eax+2],2 ;move one overflow item as its brother
		mov eax,[max_length]
		dec word[edi+deflate_state.bl_count+2*eax]
		; The brother of the overflow item also moves one step up,
		; but this does not affect bl_count[max_length]
 
		sub dword[overflow],2
		cmp dword[overflow],0
		jg .cycle2 ;while (..>0)
 
	; Now recompute all bit lengths, scanning in increasing frequency.
	; h is still equal to HEAP_SIZE. (It is simpler to reconstruct all
	; lengths instead of fixing only the wrong ones. This idea is taken
	; from 'ar' written by Haruhiko Okumura.)
 
	mov eax,[max_length]
	mov [bits],eax
	.cycle3:
	cmp dword[bits],0
	je .end_f ;for (..;..!=0;..)
		mov eax,[bits]
		shl eax,1 ;*= sizeof.dw
		movzx ecx,word[eax+edi+deflate_state.bl_count]
		.cycle4: ;while (..!=0)
		cmp ecx,0
		je .cycle4end
			dec dword[h]
			mov eax,[h]
			movzx eax,word[edi+deflate_state.heap+2*eax]
			mov [m],eax ;m = s.heap[--h]
			cmp eax,[max_code]
			jg .cycle4 ;if (..>..) continue
 
 
 
			mov esi,[m]
			imul esi,sizeof.ct_data
			add esi,[tree] ;esi = &tree[m]
			mov eax,[bits]
			cmp word[esi+Len],ax
			je @f ;if (..!=..)
;                Trace((stderr,"code %d bits %d->%d\n", m, tree[m].Len, bits));
				movzx ebx,word[esi+Len]
				sub eax,ebx
				movzx ebx,word[esi+Freq]
				imul eax,ebx ;eax = (bits - tree[m].Len) * tree[m].Freq
				add [edi+deflate_state.opt_len],eax
				mov eax,[bits]
				mov word[esi+Len],ax
			@@:
			dec ecx
			jmp .cycle4
align 4
		.cycle4end:
		dec dword[bits]
		jmp .cycle3
align 4
.end_f:
popad
	ret
endp
 
; ===========================================================================
; Generate the codes for a given tree and bit counts (which need not be
; optimal).
; IN assertion: the array bl_count contains the bit length statistics for
; the given tree and the field len is set for all tree elements.
; OUT assertion: the field code is set for all tree elements of non
;     zero code length.
 
;void (tree, max_code, bl_count)
;    ct_data *tree     ;the tree to decorate
;    int max_code      ;largest code with non zero frequency
;    uint_16p bl_count ;number of codes at each bit length
align 4
proc gen_codes uses eax ebx ecx edx edi, tree:dword, max_code:dword, bl_count:dword
locals
	u_code dw 0 ;uint_16 ;running code value
	bits   dd 1 ;int ;bit index
	next_code rw MAX_BITS+1 ;uint_16[] ;next code value for each bit length
endl
	; The distribution counts are first used to generate the code values
	; without bit reversal.
	zlib_debug 'gen_codes'
	mov ebx,ebp
	sub ebx,2*(MAX_BITS+1)
 
	.cycle0: ;for (..;..<=..;..)
	cmp dword[bits],MAX_BITS
	jg .cycle0end
		mov eax,[bits]
		dec eax
		shl eax,1
		add eax,[bl_count]
		mov ax,word[eax]
		add ax,[u_code]
		shl ax,1 ;ax = (u_code + bl_count[bits-1]) << 1
		mov [u_code],ax
		mov ecx,[bits]
		mov word[ebx+2*ecx],ax ;next_code[bits] = u_code
		inc dword[bits]
		jmp .cycle0
	.cycle0end:
	; Check that the bit counts in bl_count are consistent. The last code
	; must be all ones.
 
	mov eax,[bl_count]
	mov ax,word[eax+2*MAX_BITS]
	add ax,[u_code]
	dec ax
	cmp ax,(1 shl MAX_BITS)-1
	je @f
		zlib_assert 'inconsistent bit counts' ;Assert(..==..)
	@@:
;    Tracev((stderr,"\ngen_codes: max_code %d ", max_code));
 
	xor ecx,ecx ;n = 0
	.cycle1: ;for (..;..<=..;..)
	cmp ecx,[max_code]
	jg .cycle1end
		mov edx,sizeof.ct_data
		imul edx,ecx
		add edx,[tree] ;edx = &tree[n]
		movzx edi,word[edx+Len]
		cmp edi,0
		jne @f ;if (..==0) continue
			inc ecx
			jmp .cycle1
		@@:
		; Now reverse the bits
		movzx eax,word[ebx+2*edi]
		stdcall bi_reverse, eax, edi
		mov word[edx+Code],ax
		inc word[ebx+2*edi]
 
;        Tracecv(tree != static_ltree, (stderr,"\nn %3d %c l %2d c %4x (%x) ",
;             n, (isgraph(n) ? n : ' '), len, tree[n].Code, next_code[len]-1));
		inc ecx
		jmp .cycle1
	.cycle1end:
	ret
endp
 
; ===========================================================================
; Construct one Huffman tree and assigns the code bit strings and lengths.
; Update the total bit length for the current block.
; IN assertion: the field freq is set for all tree elements.
; OUT assertions: the fields len and code are set to the optimal bit length
;     and corresponding code. The length opt_len is updated; static_len is
;     also updated if stree is not null. The field max_code is set.
 
;void (s, desc)
;    deflate_state* s
;    tree_desc *desc ;the tree descriptor
align 4
proc build_tree uses eax ebx ecx edx edi, s:dword, desc:dword
locals
	tree     dd  ? ;ct_data* ;= desc.dyn_tree
	stree    dd  ? ;ct_data* ;= desc.stat_desc.static_tree
	elems    dd  ? ;int      ;= desc.stat_desc.elems
	m        dd  ? ;int ;iterate over heap elements
	max_code dd -1 ;int ;largest code with non zero frequency
	node     dd  ? ;int ;new node being created
endl
	; Construct the initial heap, with least frequent element in
	; heap[SMALLEST]. The sons of heap[n] are heap[2*n] and heap[2*n+1].
	; heap[0] is not used.
	mov ebx,[desc]
	mov eax,[ebx+tree_desc.dyn_tree]
	mov [tree],eax
	mov ecx,[ebx+tree_desc.stat_desc]
	mov eax,[ecx+static_tree_desc.static_tree]
	mov [stree],eax
	mov ecx,[ecx+static_tree_desc.elems]
	mov [elems],ecx
	mov edi,[s]
	zlib_debug 'build_tree cycle0 ecx = %d',ecx
 
	mov dword[edi+deflate_state.heap_len],0
	mov dword[edi+deflate_state.heap_max],HEAP_SIZE
 
 
	mov edx,[tree]
	xor ecx,ecx
	.cycle0: ;for (..;..<..;..)
 
	cmp ecx,[elems]
	jge .cycle0end
		cmp word[edx+Freq],0
		je @f ;if (..!=0)
			inc dword[edi+deflate_state.heap_len]
			mov eax,[edi+deflate_state.heap_len]
			mov [max_code],ecx
			mov [edi+deflate_state.heap+2*eax],cx
 
 
			mov byte[edi+deflate_state.depth+ecx],0
 
			jmp .end0
align 4
		@@: ;else
			mov word[edx+Len],0
		.end0:
		add edx,sizeof.ct_data
		inc ecx
		jmp .cycle0
align 4
	.cycle0end:
 
	; The pkzip format requires that at least one distance code exists,
	; and that at least one bit should be sent even if there is only one
	; possible code. So to avoid special checks later on we force at least
	; two codes of non zero frequency.
 
	.cycle1: ;while (..<..)
		cmp dword[edi+deflate_state.heap_len],2
		jge .cycle1end
		inc dword[edi+deflate_state.heap_len]
		xor eax,eax
		cmp dword[max_code],2
		jge @f
			inc dword[max_code]
			mov eax,[max_code]
		@@:
		mov ecx,[edi+deflate_state.heap_len]
		mov [edi+deflate_state.heap+2*ecx],ax
		mov [node],eax
		imul eax,sizeof.ct_data
		add eax,[tree]
		mov word[eax+Freq],1
		mov eax,[node]
		mov byte[edi+deflate_state.depth+eax],0
		dec dword[edi+deflate_state.opt_len]
		cmp dword[stree],0
		je .cycle1 ;if (..)
			mov eax,[node]
			imul eax,sizeof.ct_data
			add eax,[stree]
			movzx eax,word[eax+Len]
			sub [edi+deflate_state.static_len],eax
		; node is 0 or 1 so it does not have extra bits
		jmp .cycle1
align 4
	.cycle1end:
	mov eax,[max_code]
	mov [ebx+tree_desc.max_code],eax
 
	; The elements heap[heap_len/2+1 .. heap_len] are leaves of the tree,
	; establish sub-heaps of increasing lengths:
 
	mov ecx,[edi+deflate_state.heap_len]
	shr ecx,1
	.cycle2: ;for (..;..>=..;..)
		cmp ecx,1
		jl .cycle2end
		stdcall pqdownheap, edi, [tree], ecx
		dec ecx
		jmp .cycle2
align 4
	.cycle2end:
 
	; Construct the Huffman tree by repeatedly combining the least two
	; frequent nodes.
 
	mov eax,[elems]
	mov [node],eax ;next internal node of the tree
	.cycle3: ;do
		pqremove edi, [tree], ecx ;n = node of least frequency
		movzx edx,word[eax]
		mov [m],edx ;m = node of next least frequency
 
		mov eax,[edi+deflate_state.heap_max]
		dec eax
		mov [edi+deflate_state.heap+2*eax],cx ;keep the nodes sorted by frequency
		dec eax
		mov [edi+deflate_state.heap_max],eax
		mov [edi+deflate_state.heap+2*eax],dx
 
		; Create a new node father of n and m
		;;mov edx,[m]
		imul edx,sizeof.ct_data
		add edx,[tree]
		mov ax,word[edx+Freq]
		mov edx,ecx
		imul edx,sizeof.ct_data
		add edx,[tree]
		add ax,word[edx+Freq]
		mov edx,[node]
		imul edx,sizeof.ct_data
		add edx,[tree]
		mov word[edx+Freq],ax
 
		mov eax,ecx
		add eax,edi
		mov al,byte[eax+deflate_state.depth]
 
		mov edx,[m]
		add edx,edi
		mov ah,byte[edx+deflate_state.depth]
 
		cmp al,ah
		jl @f ;if (al>=ah) al=al : al=ah
			mov al,ah
		@@:
		inc al
		mov edx,[node]
		add edx,edi
		mov byte[edx+deflate_state.depth],al
 
 
		mov eax,[node]
		mov edx,[m]
		imul edx,sizeof.ct_data
		add edx,[tree]
		mov [edx+Dad],ax
		mov edx,ecx
		imul edx,sizeof.ct_data
		add edx,[tree]
		mov [edx+Dad],ax
;if DUMP_BL_TREE eq 1
;        if (tree == s->bl_tree) {
;            fprintf(stderr,"\nnode %d(%d), sons %d(%d) %d(%d)",
;                    node, tree[node].Freq, n, tree[n].Freq, m, tree[m].Freq);
;        }
;end if
		; and insert the new node in the heap
		mov ecx,[node]
		mov [edi+deflate_state.heap+2*SMALLEST],cx
		inc dword[node]
		stdcall pqdownheap, edi, [tree], SMALLEST
		cmp dword[edi+deflate_state.heap_len],2
		jge .cycle3 ;while (..>=..)
 
	mov cx,[edi+deflate_state.heap+2*SMALLEST]
	dec dword[edi+deflate_state.heap_max]
	mov eax,[edi+deflate_state.heap_max]
	mov [edi+deflate_state.heap+2*eax],cx
 
	; At this point, the fields freq and dad are set. We can now
	; generate the bit lengths.
 
	stdcall gen_bitlen, edi, [desc]
 
	; The field len is now set, we can generate the bit codes
	mov eax,edi
	add eax,deflate_state.bl_count
	stdcall gen_codes, [tree], [max_code], eax
	ret
endp
 
; ===========================================================================
; Scan a literal or distance tree to determine the frequencies of the codes
; in the bit length tree.
 
;void (s, tree, max_code)
;    deflate_state* s
;    ct_data *tree ;the tree to be scanned
;    int max_code  ;and its largest code of non zero frequency
align 4
proc scan_tree uses eax ebx ecx edi, s:dword, tree:dword, max_code:dword
locals
	n dd ? ;int ;iterates over all tree elements
	prevlen  dd -1 ;int ;last emitted length
	curlen    dd ? ;int ;length of current code
	nextlen   dd ? ;int ;= tree[0].Len ;length of next code
	count     dd 0 ;int ;repeat count of the current code
	max_count dd 7 ;int ;max repeat count
	min_count dd 4 ;int ;min repeat count
endl
	mov edi,[s]
	zlib_debug 'scan_tree'
	mov eax,[tree]
	movzx eax,word[eax+Len]
	mov [nextlen],eax
	cmp eax,0
	jne @f ;if (..==0)
		mov dword[max_count],138
		mov dword[min_count],3
	@@:
	mov eax,[max_code]
	inc eax
	imul eax,sizeof.ct_data
	add eax,[tree]
	mov word[eax+Len],0xffff ;guard
 
	xor ecx,ecx
	.cycle0:
		cmp ecx,[max_code]
		jg .cycle0end ;for (..;..<=..;..)
		mov eax,[nextlen]
		mov [curlen],eax
		mov eax,ecx
		inc eax
		imul eax,sizeof.ct_data
		add eax,[tree]
		movzx eax,word[eax+Len]
		mov [nextlen],eax
		inc dword[count]
		mov ebx,[count]
		cmp ebx,[max_count]
		jge .end0
		mov eax,[nextlen]
		cmp [curlen],eax
		jne .end0 ;if (..<.. && ..==..)
			inc ecx
			jmp .cycle0 ;continue
		.end0:
		cmp ebx,[min_count]
		jge .end1 ;else if (..<..)
			mov eax,[curlen]
			imul eax,sizeof.ct_data
			add eax,edi
			add word[eax+deflate_state.bl_tree+Freq],bx
 
			jmp .end4
		.end1:
		cmp dword[curlen],0
		je .end2 ;else if (..!=0)
			mov eax,[curlen]
			cmp eax,[prevlen]
			je @f ;if (..!=..)
				imul eax,sizeof.ct_data
				add eax,edi
				inc word[eax+deflate_state.bl_tree+Freq]
 
			@@:
			mov eax,REP_3_6
			imul eax,sizeof.ct_data
			add eax,edi
			inc word[eax+deflate_state.bl_tree+Freq]
 
			jmp .end4
		.end2:
		cmp ebx,10
		jg .end3 ;else if (..<=..)
			mov eax,REPZ_3_10
			imul eax,sizeof.ct_data
			add eax,edi
			inc word[eax+deflate_state.bl_tree+Freq]
 
			jmp .end4
		.end3: ;else
			mov eax,REPZ_11_138
			imul eax,sizeof.ct_data
			add eax,edi
			inc word[eax+deflate_state.bl_tree+Freq]
 
		.end4:
		mov dword[curlen],0
		mov eax,[curlen]
		mov [prevlen],eax
		mov [nextlen],eax
		cmp eax,0
		jne .end5 ;if (..==0)
			mov dword[max_count],138
			mov dword[min_count],3
			jmp .end7
		.end5:
		mov eax,[curlen]
		cmp eax,[nextlen]
		jne .end6 ;else if (..==..)
			mov dword[max_count],6
			mov dword[min_count],3
			jmp .end7
		.end6: ;else
			mov dword[max_count],7
			mov dword[min_count],4
		.end7:
		inc ecx
		jmp .cycle0
	.cycle0end:
	ret
endp
 
; ===========================================================================
; Send a literal or distance tree in compressed form, using the codes in
; bl_tree.
 
;void (s, tree, max_code)
;    deflate_state* s
;    ct_data *tree ;the tree to be scanned
;    int max_code  ;and its largest code of non zero frequency
align 4
proc send_tree uses eax ebx ecx edi, s:dword, tree:dword, max_code:dword
locals
	n dd ? ;int ;iterates over all tree elements
	prevlen  dd -1 ;int ;last emitted length
	curlen    dd ? ;int ;length of current code
	nextlen   dd ? ;int ;= tree[0].Len ;length of next code
	count     dd 0 ;int ;repeat count of the current code
	max_count dd 7 ;int ;max repeat count
	min_count dd 4 ;int ;min repeat count
endl
	mov edi,[s]
	zlib_debug 'send_tree'
	; *** tree[max_code+1].Len = -1 ;guard already set
	mov eax,[tree]
	movzx eax,word[eax+Len]
	mov [nextlen],eax
	xor ecx,ecx
	test eax,eax
	jnz .cycle0 ;if (..==0)
		mov dword[max_count],138
		mov dword[min_count],3
 
 
align 4
	.cycle0: ;for (..;..<=..;..)
	cmp ecx,[max_code]
	jg .cycle0end
		mov eax,[nextlen]
		mov [curlen],eax
		mov eax,ecx
		inc eax
		imul eax,sizeof.ct_data
		add eax,[tree]
		movzx eax,word[eax+Len]
		mov [nextlen],eax
		inc dword[count]
		mov ebx,[count]
		cmp ebx,[max_count]
		jge .end0
		mov eax,[nextlen]
		cmp [curlen],eax
		jne .end0 ;if (..<.. && ..==..)
			inc ecx
			jmp .cycle0 ;continue
align 4
		.end0:
		cmp ebx,[min_count]
		jge .end1 ;else if (..<..)
			@@: ;do
				mov ebx,edi
				add ebx,deflate_state.bl_tree
				send_code edi, [curlen], ebx
				dec dword[count]
 
				jnz @b ;while (..!=0)
			jmp .end4
align 4
		.end1:
		cmp dword[curlen],0
		je .end2 ;else if (..!=0)
			mov eax,[curlen]
			cmp eax,[prevlen]
			je @f ;if (..!=..)
				mov ebx,edi
				add ebx,deflate_state.bl_tree
				send_code edi, eax, ebx
				dec dword[count]
			@@:
			cmp dword[count],3
			jl @f
			cmp dword[count],6
			jle .end8
			@@:
				zlib_assert ' 3_6?' ;Assert(..>=.. && ..<=..)
			.end8:
			mov ebx,edi
			add ebx,deflate_state.bl_tree
			send_code edi, REP_3_6, ebx
			mov ebx,[count]
			sub ebx,3
			stdcall send_bits, edi, ebx, 2
			jmp .end4
		.end2:
		cmp ebx,10
		jg .end3 ;else if (..<=..)
			mov ebx,edi
			add ebx,deflate_state.bl_tree
			send_code edi, REPZ_3_10, ebx
			mov ebx,[count]
			sub ebx,3
			stdcall send_bits, edi, ebx, 3
			jmp .end4
		.end3: ;else
			mov ebx,edi
			add ebx,deflate_state.bl_tree
			send_code edi, REPZ_11_138, ebx
			mov ebx,[count]
			sub ebx,11
			stdcall send_bits, edi, ebx, 7
		.end4:
		mov dword[curlen],0
		mov eax,[curlen]
		mov [prevlen],eax
		mov [nextlen],eax
		cmp eax,0
		jne .end5 ;if (..==0)
			mov dword[max_count],138
			mov dword[min_count],3
			jmp .end7
		.end5:
		mov eax,[curlen]
		cmp eax,[nextlen]
		jne .end6 ;else if (..==..)
			mov dword[max_count],6
			mov dword[min_count],3
			jmp .end7
		.end6: ;else
			mov dword[max_count],7
			mov dword[min_count],4
		.end7:
		inc ecx
		jmp .cycle0
align 4
	.cycle0end:
	ret
endp
 
; ===========================================================================
; Construct the Huffman tree for the bit lengths and return the index in
; bl_order of the last bit length code to send.
 
;int (s)
;    deflate_state* s
align 4
proc build_bl_tree uses edi, s:dword
locals
	max_blindex dd ? ;int ;index of last bit length code of non zero freq
endl
	mov edi,[s]
	; Determine the bit length frequencies for literal and distance trees
	mov eax,edi
	add eax,deflate_state.dyn_ltree
	stdcall scan_tree, edi, eax, [edi+deflate_state.l_desc.max_code]
	add eax,deflate_state.dyn_dtree-deflate_state.dyn_ltree
	stdcall scan_tree, edi, eax, [edi+deflate_state.d_desc.max_code]
 
	; Build the bit length tree:
	add eax,deflate_state.bl_desc-deflate_state.dyn_dtree
	stdcall build_tree, edi, eax
	; opt_len now includes the length of the tree representations, except
	; the lengths of the bit lengths codes and the 5+5+4 bits for the counts.
 
	; Determine the number of bit length codes to send. The pkzip format
	; requires that at least 4 bit length codes be sent. (appnote.txt says
	; 3 but the actual value used is 4.)
 
	mov dword[max_blindex],BL_CODES-1
	.cycle0: ;for (..;..>=..;..)
		cmp dword[max_blindex],3
		jl .cycle0end
		dec dword[max_blindex]
		mov eax,[max_blindex]
		add eax,bl_order
		movzx eax,byte[eax]
		imul eax,sizeof.ct_data
		add eax,edi
		cmp word[eax+deflate_state.bl_tree+Len],0
 
		jne .cycle0end ;if (..!=0) break
		jmp .cycle0
align 4
	.cycle0end:
	; Update opt_len to include the bit length tree and counts
	mov eax,[max_blindex]
	inc eax
	imul eax,3
	add eax,5+5+4
	add [edi+deflate_state.opt_len],eax
;    Tracev((stderr, "\ndyn trees: dyn %ld, stat %ld", s->opt_len, s->static_len));
 
	mov eax,[max_blindex]
	ret
endp
 
; ===========================================================================
; Send the header for a block using dynamic Huffman trees: the counts, the
; lengths of the bit length codes, the literal tree and the distance tree.
; IN assertion: lcodes >= 257, dcodes >= 1, blcodes >= 4.
 
;void (s, lcodes, dcodes, blcodes)
;    deflate_state* s
;    int lcodes, dcodes, blcodes ;number of codes for each tree
align 4
proc send_all_trees uses eax ebx ecx edi, s:dword, lcodes:dword, dcodes:dword, blcodes:dword
;ecx = index in bl_order
	zlib_debug 'send_all_trees'
	cmp dword[lcodes],257
	jl @f
	cmp dword[dcodes],1
	jl @f
	cmp dword[blcodes],4
	jge .end0
	@@:
		zlib_assert 'not enough codes' ;Assert(..>=.. && ..>=.. && ..>=..)
	.end0:
	cmp dword[lcodes],L_CODES
	jg @f
	cmp dword[dcodes],D_CODES
	jg @f
	cmp dword[blcodes],BL_CODES
	jle .end1
	@@:
		zlib_assert 'too many codes' ;Assert(..<=.. && ..<=.. && ..<=..)
	.end1:
;    Tracev((stderr, "\nbl counts: "));
	mov edi,[s]
	mov eax,[lcodes]
	sub eax,257
	stdcall send_bits, edi, eax, 5 ;not +255 as stated in appnote.txt
	mov eax,[dcodes]
	dec eax
	stdcall send_bits, edi, eax, 5
	mov eax,[blcodes]
	sub eax,4
	stdcall send_bits, edi, eax, 4 ;not -3 as stated in appnote.txt
	xor ecx,ecx
	.cycle0:
		cmp ecx,[blcodes]
		jge .cycle0end ;for (..;..<..;..)
;        Tracev((stderr, "\nbl code %2d ", bl_order[ecx]));
		mov eax,ecx
		add eax,bl_order
		movzx eax,byte[eax]
		imul eax,sizeof.ct_data
		add eax,edi
		movzx eax,word[eax+deflate_state.bl_tree+Len]
 
		stdcall send_bits, edi, eax, 3
		inc ecx
		jmp .cycle0
align 4
	.cycle0end:
;    Tracev((stderr, "\nbl tree: sent %ld", s->bits_sent));
 
	mov ebx,[lcodes]
	dec ebx
	mov eax,edi
	add eax,deflate_state.dyn_ltree
	stdcall send_tree, edi, eax, ebx ;literal tree
;    Tracev((stderr, "\nlit tree: sent %ld", s->bits_sent));
 
	mov ebx,[dcodes]
	dec ebx
	add eax,deflate_state.dyn_dtree-deflate_state.dyn_ltree
	stdcall send_tree, edi, eax, ebx ;distance tree
;    Tracev((stderr, "\ndist tree: sent %ld", s->bits_sent));
	ret
endp
 
; ===========================================================================
; Send a stored block
 
;void (s, buf, stored_len, last)
;    deflate_state* s
;    charf *buf     ;input block
;    ulg stored_len ;length of input block
;    int last       ;one if this is the last block for a file
align 4
proc _tr_stored_block uses eax edi, s:dword, buf:dword, stored_len:dword, last:dword
	mov edi,[s]
	mov eax,[last]
	add eax,STORED_BLOCK shl 1
	stdcall send_bits, edi, eax, 3 ;send block type
if DEBUG eq 1
	mov eax,[edi+deflate_state.compressed_len]
	add eax,3+7
	and eax,not 7
	mov [edi+deflate_state.compressed_len],eax
	mov eax,[stored_len]
	add eax,4
	shl eax,3
	add [edi+deflate_state.compressed_len],eax
end if
	stdcall copy_block, edi, [buf], [stored_len], 1 ;with header
	ret
endp
 
; ===========================================================================
; Flush the bits in the bit buffer to pending output (leaves at most 7 bits)
 
;void (s)
;    deflate_state* s
align 4
proc _tr_flush_bits, s:dword
	stdcall bi_flush, [s]
	ret
endp
 
; ===========================================================================
; Send one empty static block to give enough lookahead for inflate.
; This takes 10 bits, of which 7 may remain in the bit buffer.
 
;void (s)
;    deflate_state* s
align 4
proc _tr_align uses edi, s:dword
	mov edi,[s]
	stdcall send_bits, edi, STATIC_TREES shl 1, 3
	send_code edi, END_BLOCK, static_ltree
if DEBUG eq 1
	add [edi+deflate_state.compressed_len],10 ;3 for block type, 7 for EOB
end if
	stdcall bi_flush, edi
	ret
endp
 
; ===========================================================================
; Determine the best encoding for the current block: dynamic trees, static
; trees or store, and output the encoded block to the zip file.
 
;void (s, buf, stored_len, last)
;    deflate_state* s
;    charf *buf     ;input block, or NULL if too old
;    ulg stored_len ;length of input block
;    int last       ;one if this is the last block for a file
align 4
proc _tr_flush_block uses eax ebx edi, s:dword, buf:dword, stored_len:dword, last:dword
locals
	opt_lenb dd ? ;ulg
	static_lenb dd ? ;opt_len and static_len in bytes
	max_blindex dd 0 ;int ;index of last bit length code of non zero freq
endl
	; Build the Huffman trees unless a stored block is forced
	mov edi,[s]
	zlib_debug '_tr_flush_block'
	cmp word[edi+deflate_state.level],0
	jle .end0 ;if (..>0)
 
		; Check if the file is binary or text
		mov ebx,[edi+deflate_state.strm]
		cmp dword[ebx+z_stream.data_type],Z_UNKNOWN
		jne @f ;if (..==..)
			stdcall detect_data_type, edi
			mov [ebx+z_stream.data_type],eax
		@@:
 
		; Construct the literal and distance trees
		mov eax,edi
		add eax,deflate_state.l_desc
		stdcall build_tree, edi, eax
;        Tracev((stderr, "\nlit data: dyn %ld, stat %ld", s->opt_len, s->static_len));
 
		mov eax,edi
		add eax,deflate_state.d_desc
		stdcall build_tree, edi, eax
;        Tracev((stderr, "\ndist data: dyn %ld, stat %ld", s->opt_len, s->static_len));
		; At this point, opt_len and static_len are the total bit lengths of
		; the compressed block data, excluding the tree representations.
 
		; Build the bit length tree for the above two trees, and get the index
		; in bl_order of the last bit length code to send.
 
		stdcall build_bl_tree, edi
		mov [max_blindex],eax
 
		; Determine the best encoding. Compute the block lengths in bytes.
		mov eax,[edi+deflate_state.opt_len]
		add eax,3+7
		shr eax,3
		mov [opt_lenb],eax
		mov eax,[edi+deflate_state.static_len]
		add eax,3+7
		shr eax,3
		mov [static_lenb],eax
 
;        Tracev((stderr, "\nopt %lu(%lu) stat %lu(%lu) stored %lu lit %u ",
;                opt_lenb, s->opt_len, static_lenb, s->static_len, stored_len,
;                s->last_lit));
 
		cmp eax,[opt_lenb]
		jg .end1 ;if (..<=..)
			mov [opt_lenb],eax
		jmp .end1
	.end0: ;else
		cmp dword[buf],0
		jne @f
			zlib_assert 'lost buf' ;Assert(..!=0)
		@@:
		mov eax,[stored_len]
		add eax,5
		mov [static_lenb],eax
		mov [opt_lenb],eax ;force a stored block
	.end1:
 
if FORCE_STORED eq 1
	cmp dword[buf],0
	je .end2 ;if (..!=0) ;force stored block
else
	mov eax,[stored_len]
	add eax,4
	cmp eax,[opt_lenb]
	jg .end2
	cmp dword[buf],0
	je .end2 ;if (..<=.. && ..!=0)
		;4: two words for the lengths
end if
		; The test buf != NULL is only necessary if LIT_BUFSIZE > WSIZE.
		; Otherwise we can't have processed more than WSIZE input bytes since
		; the last block flush, because compression would have been
		; successful. If LIT_BUFSIZE <= WSIZE, it is never too late to
		; transform a block into a stored block.
 
		stdcall _tr_stored_block, edi, [buf], [stored_len], [last]
		jmp .end4
	.end2:
if FORCE_STATIC eq 1
	cmp dword[static_lenb],0
	jl .end3 ;else if (..>=0) ;force static trees
else
	cmp word[edi+deflate_state.strategy],Z_FIXED
	je @f
	mov eax,[opt_lenb]
	cmp [static_lenb],eax
	je @f ;else if (..==.. || ..==..)
		jmp .end3
	@@:
end if
		mov eax,STATIC_TREES shl 1
		add eax,[last]
		stdcall send_bits, edi, eax, 3
		stdcall compress_block, edi, static_ltree, static_dtree
if DEBUG eq 1
		mov eax,[edi+deflate_state.static_len]
		add eax,3
		add [edi+deflate_state.compressed_len],eax
end if
		jmp .end4
	.end3: ;else
		mov eax,DYN_TREES shl 1
		add eax,[last]
		stdcall send_bits, edi, eax, 3
		mov eax,[max_blindex]
		inc eax
		push eax
		mov eax,[edi+deflate_state.d_desc.max_code]
		inc eax
		push eax
		mov eax,[edi+deflate_state.l_desc.max_code]
		inc eax
		stdcall send_all_trees, edi, eax ;, ..., ...
		mov eax,edi
		add eax,deflate_state.dyn_dtree
		push eax
		add eax,deflate_state.dyn_ltree-deflate_state.dyn_dtree
		stdcall compress_block, edi, eax ;, ...
if DEBUG eq 1
		mov eax,[edi+deflate_state.opt_len]
		add eax,3
		add [edi+deflate_state.compressed_len],eax
end if
	.end4:
;    Assert (s->compressed_len == s->bits_sent, "bad compressed size");
	; The above check is made mod 2^32, for files larger than 512 MB
	; and uLong implemented on 32 bits.
 
	stdcall init_block,edi
 
	cmp dword[last],0
	je @f ;if (..)
		stdcall bi_windup,edi
if DEBUG eq 1
		add [edi+deflate_state.compressed_len],7 ;align on byte boundary
end if
	@@:
;    Tracev((stderr,"\ncomprlen %lu(%lu) ", s->compressed_len>>3,
;           s->compressed_len-7*last));
	ret
endp
 
; ===========================================================================
; Save the match info and tally the frequency counts. Return true if
; the current block must be flushed.
 
;int (s, dist, lc)
;    deflate_state* s
;    unsigned dist ;distance of matched string
;    unsigned lc   ;match length-MIN_MATCH or unmatched char (if dist==0)
align 4
proc _tr_tally uses ebx edi, s:dword, dist:dword, lc:dword
	mov edi,[s]
	zlib_debug '_tr_tally'
	mov eax,[edi+deflate_state.last_lit]
	shl eax,1
	add eax,[edi+deflate_state.d_buf]
	mov ebx,[dist]
	mov word[eax],bx
	mov eax,[edi+deflate_state.last_lit]
	add eax,[edi+deflate_state.l_buf]
	mov ebx,[lc]
	mov byte[eax],bl
	inc dword[edi+deflate_state.last_lit]
	cmp dword[dist],0
	jne @f ;if (..==0)
		; lc is the unmatched char
		mov eax,[lc]
		imul eax,sizeof.ct_data
		add eax,edi
		inc word[eax+deflate_state.dyn_ltree+Freq]
		jmp .end0
	@@: ;else
		inc dword[edi+deflate_state.matches]
		; Here, lc is the match length - MIN_MATCH
		dec dword[dist] ;dist = match distance - 1
		MAX_DIST edi
		cmp word[dist],ax
		jge @f
		cmp word[lc],MAX_MATCH-MIN_MATCH
		jg @f
		d_code [dist]
		cmp ax,D_CODES
		jl .end2
		@@:
			zlib_assert '_tr_tally: bad match' ;Assert(..<.. && ..<=.. && ..<..)
		.end2:
		mov eax,[lc]
		add eax,_length_code
		movzx eax,byte[eax]
		add eax,LITERALS+1
		imul eax,sizeof.ct_data
		add eax,edi
		inc word[eax+deflate_state.dyn_ltree+Freq]
		d_code [dist]
		imul eax,sizeof.ct_data
		add eax,edi
		inc word[eax+deflate_state.dyn_dtree+Freq]
	.end0:
 
if TRUNCATE_BLOCK eq 1
	; Try to guess if it is profitable to stop the current block here
	mov eax,[edi+deflate_state.last_lit]
	and eax,0x1fff
	cmp eax,0
	jne .end1
	cmp word[edi+deflate_state.level],2
	jle .end1 ;if (..==.. && ..>..)
	; Compute an upper bound for the compressed length
;        ulg out_length = (ulg)s->last_lit*8L;
;        ulg in_length = (ulg)((long)s->strstart - s->block_start);
;        int dcode;
;        for (dcode = 0; dcode < D_CODES; dcode++) {
;            out_length += (ulg)s->dyn_dtree[dcode].Freq *
;                (5L+extra_dbits[dcode]);
;        }
;        out_length >>= 3;
;        Tracev((stderr,"\nlast_lit %u, in %ld, out ~%ld(%ld%%) ",
;               s->last_lit, in_length, out_length,
;               100L - out_length*100L/in_length));
;        if (s->matches < s->last_lit/2 && out_length < in_length/2) return 1;
	.end1:
end if
	mov ebx,[edi+deflate_state.last_lit]
	mov edi,[edi+deflate_state.lit_bufsize]
	dec edi
	xor eax,eax
	cmp ebx,edi
	jne @f
		inc eax ;return (..==..)
	@@:
	; We avoid equality with lit_bufsize because of wraparound at 64K
	; on 16 bit machines and because stored blocks are restricted to
	; 64K-1 bytes.
	ret
endp
 
; ===========================================================================
; Send the block data compressed using the given Huffman trees
 
;void (s, ltree, dtree)
;    deflate_state* s
;    ct_data *ltree ;literal tree
;    ct_data *dtree ;distance tree
align 4
proc compress_block uses eax edi, s:dword, ltree:dword, dtree:dword
locals
	dist  dd ? ;unsigned ;distance of matched string
	lc    dd ? ;int      ;match length or unmatched char (if dist == 0)
	lx    dd 0 ;unsigned ;running index in l_buf
	u_code dd ? ;unsigned ;the code to send
	extra  dd ? ;int      ;number of extra bits to send
endl
	mov edi,[s]
	cmp dword[edi+deflate_state.last_lit],0
	je .end0 ;if (..!=0)
	.cycle0: ; do
		mov eax,[lx]
		shl eax,1
		add eax,[edi+deflate_state.d_buf]
		movzx eax,word[eax]
		mov [dist],eax
		mov eax,[lx]
		add eax,[edi+deflate_state.l_buf]
		movzx eax,byte[eax]
		mov [lc],eax
		inc dword[lx]
		cmp dword[dist],0
		jne @f ;if (..==0)
			send_code edi, [lc], [ltree] ;send a literal byte
;            Tracecv(isgraph(lc), (stderr," '%c' ", lc));
			jmp .end1
		@@: ;else
			; Here, lc is the match length - MIN_MATCH
			mov eax,[lc]
			add eax,_length_code
			movzx eax,byte[eax]
			mov [u_code],eax
			add eax,LITERALS+1
			send_code edi, eax, [ltree] ;send the length code
			mov eax,[u_code]
			shl eax,2
			add eax,extra_lbits
			mov eax,[eax]
			mov [extra],eax
			cmp eax,0
			je @f ;if (..!=0)
				mov eax,[u_code]
				shl eax,2
				add eax,base_length
				mov eax,[eax]
				sub [lc],eax
				stdcall send_bits, edi, [lc], [extra] ;send the extra length bits
			@@:
			dec dword[dist] ;dist is now the match distance - 1
			d_code [dist]
			mov [u_code],eax
			cmp eax,D_CODES
			jl @f
				zlib_assert 'bad d_code' ;Assert(..<..)
			@@:
			send_code edi, [u_code], [dtree] ;send the distance code
			mov eax,[u_code]
			shl eax,2
			add eax,extra_dbits
			mov eax,[eax]
			mov [extra],eax
			cmp eax,0
			je .end1 ;if (..!=0)
				mov eax,[u_code]
				shl eax,2
				add eax,base_dist
				mov eax,[eax]
				sub [dist],eax
				stdcall send_bits, edi, [dist], [extra] ;send the extra distance bits
		.end1: ;literal or match pair ?
 
		; Check that the overlay between pending_buf and d_buf+l_buf is ok:
		mov eax,[lx]
		shl eax,1
		add eax,[edi+deflate_state.lit_bufsize]
		cmp [edi+deflate_state.pending],eax
		jl @f
			zlib_assert 'pendingBuf overflow' ;Assert(..<..)
		@@:
		mov eax,[edi+deflate_state.last_lit]
		cmp [lx],eax
		jl .cycle0 ;while (..<..)
align 4
	.end0:
 
	send_code edi, END_BLOCK, [ltree]
	ret
endp
 
; ===========================================================================
; Check if the data type is TEXT or BINARY, using the following algorithm:
; - TEXT if the two conditions below are satisfied:
;    a) There are no non-portable control characters belonging to the
;       "black list" (0..6, 14..25, 28..31).
;    b) There is at least one printable character belonging to the
;       "white list" (9 {TAB}, 10 {LF}, 13 {CR}, 32..255).
; - BINARY otherwise.
; - The following partially-portable control characters form a
;   "gray list" that is ignored in this detection algorithm:
;   (7 {BEL}, 8 {BS}, 11 {VT}, 12 {FF}, 26 {SUB}, 27 {ESC}).
; IN assertion: the fields Freq of dyn_ltree are set.
 
;int (s)
;    deflate_state* s
align 4
proc detect_data_type uses ebx ecx edi, s:dword
	; black_mask is the bit mask of black-listed bytes
	; set bits 0..6, 14..25, and 28..31
	; 0xf3ffc07f = binary 11110011111111111100000001111111
locals
	black_mask dd 0xf3ffc07f
endl
	mov edi,[s]
	zlib_debug 'detect_data_type'
 
	; Check for non-textual ("black-listed") bytes.
	xor ecx,ecx
	mov ebx,edi
	add ebx,deflate_state.dyn_ltree+Freq
	.cycle0:
	cmp ecx,31
	jg .cycle0end ;for (..;..<=..;..,..)
		bt dword[black_mask],0
		jnc @f
		cmp word[ebx],0
		je @f ;if (..&.. && ..!=0)
			mov eax,Z_BINARY
			jmp .end_f
		@@:
		shr dword[black_mask],1
		add ebx,sizeof.ct_data
		inc ecx
		jmp .cycle0
	.cycle0end:
 
	; Check for textual ("white-listed") bytes.
	mov ebx,edi
	add ebx,deflate_state.dyn_ltree+Freq+9*sizeof.ct_data
	cmp word[ebx],0
	jne @f
	add ebx,sizeof.ct_data
	cmp word[ebx],0
	jne @f
	add ebx,3*sizeof.ct_data
	cmp word[ebx],0
	je .end0
	@@: ;if (..!=0 || ..!=0 || ..!= 0)
		mov eax,Z_TEXT
		jmp .end_f
	.end0:
	mov ecx,32
	mov ebx,edi
	add ebx,deflate_state.dyn_ltree+Freq+32*sizeof.ct_data
	.cycle1:
	cmp ecx,LITERALS
	jge .cycle1end ;for (..;..<..;..,..)
		cmp word[ebx],0
		je @f ;if (..!=0)
			mov eax,Z_TEXT
			jmp .end_f
		@@:
		add ebx,sizeof.ct_data
		inc ecx
		jmp .cycle1
	.cycle1end:
 
	; There are no "black-listed" or "white-listed" bytes:
	; this stream either is empty or has tolerated ("gray-listed") bytes only.
 
	mov eax,Z_BINARY
.end_f:
	ret
endp
 
; ===========================================================================
; Reverse the first len bits of a code, using straightforward code (a faster
; method would use a table)
; IN assertion: 1 <= len <= 15
 
;unsigned (code, len)
;    unsigned code ;the value to invert
;    int len       ;its bit length
align 4
proc bi_reverse uses ebx, p1code:dword, len:dword
	zlib_debug 'bi_reverse'
	xor eax,eax
	@@: ;do
		mov ebx,[p1code]
		and ebx,1
		or eax,ebx
		shr dword[p1code],1
		shl eax,1
		dec dword[len]
		cmp dword[len],0
		jg @b ;while (..>..)
	shr eax,1
	ret
endp
 
; ===========================================================================
; Flush the bit buffer, keeping at most 7 bits in it.
 
;void (s)
;    deflate_state* s
align 4
proc bi_flush uses eax ecx edi, s:dword
	mov edi,[s]
	cmp dword[edi+deflate_state.bi_valid],16
	jne @f ;if (..==..)
		mov cx,[edi+deflate_state.bi_buf]
		put_short edi,cx
		mov word[edi+deflate_state.bi_buf],0
		mov dword[edi+deflate_state.bi_valid],0
		jmp .end0
	@@: ;else if (..>=..)
		cmp dword[edi+deflate_state.bi_valid],8
		jl .end0
		mov cl,byte[edi+deflate_state.bi_buf]
		put_byte edi,cl
		shr word[edi+deflate_state.bi_buf],8
		sub dword[edi+deflate_state.bi_valid],8
	.end0:
	ret
endp
 
; ===========================================================================
; Flush the bit buffer and align the output on a byte boundary
 
;void (s)
;    deflate_state* s
align 4
proc bi_windup uses eax ecx edi, s:dword
	mov edi,[s]
	cmp dword[edi+deflate_state.bi_valid],8
	jle @f ;if (..>..)
		mov cx,[edi+deflate_state.bi_buf]
		put_short edi, cx
		jmp .end0
	@@: ;else if (..>0)
		cmp dword[edi+deflate_state.bi_valid],0
		jle .end0
		mov cl,byte[edi+deflate_state.bi_buf]
		put_byte edi, cl
	.end0:
	mov word[edi+deflate_state.bi_buf],0
	mov dword[edi+deflate_state.bi_valid],0
if DEBUG eq 1
	mov eax,[edi+deflate_state.bits_sent]
	add eax,7
	and eax,not 7
	mov [edi+deflate_state.bits_sent],eax
end if
	ret
endp
 
; ===========================================================================
; Copy a stored block, storing first the length and its
; one's complement if requested.
 
;void (s, buf, len, header)
;    deflate_state* s
;    charf    *buf   ;the input data
;    unsigned len    ;its length
;    int      header ;true if block header must be written
align 4
proc copy_block uses eax ebx ecx edi esi, s:dword, buf:dword, len:dword, p4header:dword
	mov edi,[s]
	stdcall bi_windup,edi ;align on byte boundary
 
	cmp dword[p4header],0
	je @f ;if (..)
		mov ecx,[len]
		put_short edi, cx
		not cx
		put_short edi, cx
if DEBUG eq 1
		add dword[edi+deflate_state.bits_sent],2*16
end if
	@@:
if DEBUG eq 1
	mov ecx,[len]
	shl ecx,3
	add [edi+deflate_state.bits_sent],ecx
end if
	mov ecx,[len]
	mov esi,[buf]
	@@: ;while (len--)
		lodsb
		mov bl,al
		put_byte edi, bl
		loop @b
	ret
endp
>
>
align>
>
>
>
>
;>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>

Rev 6799	Rev 6813
1	; trees.asm -- output deflated data using Huffman coding	1	; trees.asm -- output deflated data using Huffman coding
2	; Copyright (C) 1995-2012 Jean-loup Gailly	2	; Copyright (C) 1995-2012 Jean-loup Gailly
3	; detect_data_type() function provided freely by Cosmin Truta, 2006	3	; detect_data_type() function provided freely by Cosmin Truta, 2006
4	; For conditions of distribution and use, see copyright notice in zlib.h	4	; For conditions of distribution and use, see copyright notice in zlib.h
5		5
6	; ALGORITHM	6	; ALGORITHM
7		7
8	; The "deflation" process uses several Huffman trees. The more	8	; The "deflation" process uses several Huffman trees. The more
9	; common source values are represented by shorter bit sequences.	9	; common source values are represented by shorter bit sequences.
10		10
11	; Each code tree is stored in a compressed form which is itself	11	; Each code tree is stored in a compressed form which is itself
12	; a Huffman encoding of the lengths of all the code strings (in	12	; a Huffman encoding of the lengths of all the code strings (in
13	; ascending order by source values). The actual code strings are	13	; ascending order by source values). The actual code strings are
14	; reconstructed from the lengths in the inflate process, as described	14	; reconstructed from the lengths in the inflate process, as described
15	; in the deflate specification.	15	; in the deflate specification.
16		16
17	; REFERENCES	17	; REFERENCES
18		18
19	; Deutsch, L.P.,"'Deflate' Compressed Data Format Specification".	19	; Deutsch, L.P.,"'Deflate' Compressed Data Format Specification".
20	; Available in ftp.uu.net:/pub/archiving/zip/doc/deflate-1.1.doc	20	; Available in ftp.uu.net:/pub/archiving/zip/doc/deflate-1.1.doc
21		21
22	; Storer, James A.	22	; Storer, James A.
23	; Data Compression: Methods and Theory, pp. 49-50.	23	; Data Compression: Methods and Theory, pp. 49-50.
24	; Computer Science Press, 1988. ISBN 0-7167-8156-5.	24	; Computer Science Press, 1988. ISBN 0-7167-8156-5.
25		25
26	; Sedgewick, R.	26	; Sedgewick, R.
27	; Algorithms, p290.	27	; Algorithms, p290.
28	; Addison-Wesley, 1983. ISBN 0-201-06672-6.	28	; Addison-Wesley, 1983. ISBN 0-201-06672-6.
29		29
30	; ===========================================================================	30	; ===========================================================================
31	; Constants	31	; Constants
32		32
33		33
34	MAX_BL_BITS equ 7	34	MAX_BL_BITS equ 7
35	; Bit length codes must not exceed MAX_BL_BITS bits	35	; Bit length codes must not exceed MAX_BL_BITS bits
36		36
37	END_BLOCK equ 256	37	END_BLOCK equ 256
38	; end of block literal code	38	; end of block literal code
39		39
40	REP_3_6 equ 16	40	REP_3_6 equ 16
41	; repeat previous bit length 3-6 times (2 bits of repeat count)	41	; repeat previous bit length 3-6 times (2 bits of repeat count)
42		42
43	REPZ_3_10 equ 17	43	REPZ_3_10 equ 17
44	; repeat a zero length 3-10 times (3 bits of repeat count)	44	; repeat a zero length 3-10 times (3 bits of repeat count)
45		45
46	REPZ_11_138 equ 18	46	REPZ_11_138 equ 18
47	; repeat a zero length 11-138 times (7 bits of repeat count)	47	; repeat a zero length 11-138 times (7 bits of repeat count)
48		48
49	align 4	49	align 4
50	extra_lbits dd \ ;int [LENGTH_CODES] ;extra bits for each length code	50	extra_lbits dd \ ;int [LENGTH_CODES] ;extra bits for each length code
51	0,0,0,0,0,0,0,0,1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4,5,5,5,5,0	51	0,0,0,0,0,0,0,0,1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4,5,5,5,5,0
52		52
53	align 4	53	align 4
54	extra_dbits dd \ ;int [D_CODES] ;extra bits for each distance code	54	extra_dbits dd \ ;int [D_CODES] ;extra bits for each distance code
55	0,0,0,0,1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,11,11,12,12,13,13	55	0,0,0,0,1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,11,11,12,12,13,13
56		56
57	align 4	57	align 4
58	extra_blbits dd \ ;int [BL_CODES] ;extra bits for each bit length code	58	extra_blbits dd \ ;int [BL_CODES] ;extra bits for each bit length code
59	0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,3,7	59	0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,3,7
60		60
61	align 4	61	align 4
62	bl_order db 16,17,18,0,8,7,9,6,10,5,11,4,12,3,13,2,14,1,15	62	bl_order db 16,17,18,0,8,7,9,6,10,5,11,4,12,3,13,2,14,1,15
63	; The lengths of the bit length codes are sent in order of decreasing	63	; The lengths of the bit length codes are sent in order of decreasing
64	; probability, to avoid transmitting the lengths for unused bit length codes.	64	; probability, to avoid transmitting the lengths for unused bit length codes.
65		65
66		66
67	; ===========================================================================	67	; ===========================================================================
68	; Local data. These are initialized only once.	68	; Local data. These are initialized only once.
69		69
70		70
71	DIST_CODE_LEN equ 512 ;see definition of array dist_code below	71	DIST_CODE_LEN equ 512 ;see definition of array dist_code below
72		72
73	if GEN_TREES_H eq 1 ;\| !(STDC)	73	if GEN_TREES_H eq 1 ;\| !(STDC)
74	; non ANSI compilers may not accept trees.inc	74	; non ANSI compilers may not accept trees.inc
75		75
76	align 4	76	align 4
77	static_ltree rb sizeof.ct_data * (L_CODES+2)	77	static_ltree rb sizeof.ct_data * (L_CODES+2)
78	; The static literal tree. Since the bit lengths are imposed, there is no	78	; The static literal tree. Since the bit lengths are imposed, there is no
79	; need for the L_CODES extra codes used during heap construction. However	79	; need for the L_CODES extra codes used during heap construction. However
80	; The codes 286 and 287 are needed to build a canonical tree (see _tr_init	80	; The codes 286 and 287 are needed to build a canonical tree (see _tr_init
81	; below).	81	; below).
82		82
83	align 4	83	align 4
84	static_dtree rb sizeof.ct_data * D_CODES	84	static_dtree rb sizeof.ct_data * D_CODES
85	; The static distance tree. (Actually a trivial tree since all codes use	85	; The static distance tree. (Actually a trivial tree since all codes use
86	; 5 bits.)	86	; 5 bits.)
87		87
88	align 4	88	align 4
89	_dist_code rb DIST_CODE_LEN ;uch[]	89	_dist_code rb DIST_CODE_LEN ;uch[]
90	; Distance codes. The first 256 values correspond to the distances	90	; Distance codes. The first 256 values correspond to the distances
91	; 3 .. 258, the last 256 values correspond to the top 8 bits of	91	; 3 .. 258, the last 256 values correspond to the top 8 bits of
92	; the 15 bit distances.	92	; the 15 bit distances.
93		93
94	align 4	94	align 4
95	_length_code rb MAX_MATCH-MIN_MATCH+1 ;uch[]	95	_length_code rb MAX_MATCH-MIN_MATCH+1 ;uch[]
96	; length code for each normalized match length (0 == MIN_MATCH)	96	; length code for each normalized match length (0 == MIN_MATCH)
97		97
98	align 4	98	align 4
99	base_length rd LENGTH_CODES ;int[]	99	base_length rd LENGTH_CODES ;int[]
100	; First normalized length for each code (0 = MIN_MATCH)	100	; First normalized length for each code (0 = MIN_MATCH)
101		101
102	align 4	102	align 4
103	base_dist rd D_CODES ;int[]	103	base_dist rd D_CODES ;int[]
104	; First normalized distance for each code (0 = distance of 1)	104	; First normalized distance for each code (0 = distance of 1)
105		105
106	else	106	else
107	include 'trees.inc'	107	include 'trees.inc'
108	end if ;GEN_TREES_H	108	end if ;GEN_TREES_H
109		109
110	struct static_tree_desc ;_s	110	struct static_tree_desc ;_s
111	static_tree dd ? ;const ct_data * ;static tree or NULL	111	static_tree dd ? ;const ct_data * ;static tree or NULL
112	extra_bits dd ? ;const intf * ;extra bits for each code or NULL	112	extra_bits dd ? ;const intf * ;extra bits for each code or NULL
113	extra_base dd ? ;int ;base index for extra_bits	113	extra_base dd ? ;int ;base index for extra_bits
114	elems dd ? ;int ;max number of elements in the tree	114	elems dd ? ;int ;max number of elements in the tree
115	max_length dd ? ;int ;max bit length for the codes	115	max_length dd ? ;int ;max bit length for the codes
116	ends	116	ends
117		117
118	align 4	118	align 4
119	static_l_desc static_tree_desc static_ltree, extra_lbits, LITERALS+1, L_CODES, MAX_BITS	119	static_l_desc static_tree_desc static_ltree, extra_lbits, LITERALS+1, L_CODES, MAX_BITS
120		120
121	align 4	121	align 4
122	static_d_desc static_tree_desc static_dtree, extra_dbits, 0, D_CODES, MAX_BITS	122	static_d_desc static_tree_desc static_dtree, extra_dbits, 0, D_CODES, MAX_BITS
123		123
124	align 4	124	align 4
125	static_bl_desc static_tree_desc 0, extra_blbits, 0, BL_CODES, MAX_BL_BITS	125	static_bl_desc static_tree_desc 0, extra_blbits, 0, BL_CODES, MAX_BL_BITS
126		126
127	; ===========================================================================	127	; ===========================================================================
128	; Local (static) routines in this file.	128	; Local (static) routines in this file.
129		129
130		130
131	macro send_code s, c, tree	131	macro send_code s, c, tree
132	{	132	{
133	if DEBUG eq 1	133	if DEBUG eq 1
134	; if (z_verbose>2) fprintf(stderr,"\ncd %3d ",(c))	134	; if (z_verbose>2) fprintf(stderr,"\ncd %3d ",(c))
135	end if	135	end if
136	push eax ebx	136	push eax ebx
137	if c eq eax	137	if c eq eax
138	else	138	else
139	mov eax,c	139	mov eax,c
140	end if	140	end if
141	imul eax,sizeof.ct_data	141	imul eax,sizeof.ct_data
142	add eax,tree	142	add eax,tree
143	movzx ebx,word[eax+Len]	143	movzx ebx,word[eax+Len]
144	push ebx	144	push ebx
145	movzx ebx,word[eax+Code]	145	movzx ebx,word[eax+Code]
146	push ebx	146	push ebx
147	stdcall send_bits, s ;tree[c].Code, tree[c].Len	147	stdcall send_bits, s ;tree[c].Code, tree[c].Len
148	pop ebx eax	148	pop ebx eax
149	}	149	}
150	; Send a code of the given tree[c] and tree must not have side effects	150	; Send a code of the given tree[c] and tree must not have side effects
151		151
152	; ===========================================================================	152	; ===========================================================================
153	; Output a short LSB first on the stream.	153	; Output a short LSB first on the stream.
154	; IN assertion: there is enough room in pendingBuf.	154	; IN assertion: there is enough room in pendingBuf.
155		155
156	macro put_short s, w	156	macro put_short s, w
157	{	157	{
158	mov eax,[s+deflate_state.pending]	158	mov eax,[s+deflate_state.pending]
159	add eax,[s+deflate_state.pending_buf]	159	add eax,[s+deflate_state.pending_buf]
160	mov word[eax],w	160	mov word[eax],w
161	add dword[s+deflate_state.pending],2	161	add dword[s+deflate_state.pending],2
162	}	162	}
163		163
164	; ===========================================================================	164	; ===========================================================================
165	; Send a value on a given number of bits.	165	; Send a value on a given number of bits.
166	; IN assertion: length <= 16 and value fits in length bits.	166	; IN assertion: length <= 16 and value fits in length bits.
167		167
168	;void (s, value, length)	168	;void (s, value, length)
169	; deflate_state* s	169	; deflate_state* s
170	; int value ;value to send	170	; int value ;value to send
171	; int length ;number of bits	171	; int length ;number of bits
172	align 4	172	align 4
173	proc send_bits uses eax ecx edi, s:dword, value:dword, length:dword	173	proc send_bits uses eax ecx edi, s:dword, value:dword, length:dword
174	; Tracevv((stderr," l %2d v %4x ", length, value));	174	; Tracevv((stderr," l %2d v %4x ", length, value));
175	zlib_debug 'send_bits value = %d',[value]	175	zlib_debug 'send_bits value = %d',[value]
176	;if DEBUG eq 1	176	;if DEBUG eq 1
177	mov eax,[length]	177	mov eax,[length]
178	cmp eax,0	178	cmp eax,0
179	jle @f	179	jle @f
180	cmp eax,15	180	cmp eax,15
181	jle .end1	181	jle .end1
182	@@:	182	@@:
183	zlib_assert 'invalid length' ;Assert(..>0 && ..<=15)	183	zlib_assert 'invalid length' ;Assert(..>0 && ..<=15)
184	.end1:	184	.end1:
185	mov edi,[s]	185	mov edi,[s]
186	add [edi+deflate_state.bits_sent],eax	186	add [edi+deflate_state.bits_sent],eax
187		187
188	; If not enough room in bi_buf, use (valid) bits from bi_buf and	188	; If not enough room in bi_buf, use (valid) bits from bi_buf and
189	; (16 - bi_valid) bits from value, leaving (width - (16-bi_valid))	189	; (16 - bi_valid) bits from value, leaving (width - (16-bi_valid))
190	; unused bits in value.	190	; unused bits in value.
191		191
192	mov ecx,Buf_size	192	mov ecx,Buf_size
193	sub ecx,eax	193	sub ecx,eax
194	cmp [edi+deflate_state.bi_valid],ecx	194	cmp [edi+deflate_state.bi_valid],ecx
195	jle @f ;if (..>..)	195	jle @f ;if (..>..)
196	mov eax,[value]	196	mov eax,[value]
197	mov ecx,[edi+deflate_state.bi_valid]	197	mov ecx,[edi+deflate_state.bi_valid]
198	shl eax,cl	198	shl eax,cl
199	or [edi+deflate_state.bi_buf],ax	199	or [edi+deflate_state.bi_buf],ax
200	mov cx,[edi+deflate_state.bi_buf]	200	mov cx,[edi+deflate_state.bi_buf]
201	put_short edi, cx	201	put_short edi, cx
202	mov eax,[value]	202	mov eax,[value]
203	mov ecx,Buf_size	203	mov ecx,Buf_size
204	sub ecx,[edi+deflate_state.bi_valid]	204	sub ecx,[edi+deflate_state.bi_valid]
205	shr eax,cl	205	shr eax,cl
206	mov [edi+deflate_state.bi_buf],ax	206	mov [edi+deflate_state.bi_buf],ax
207	mov eax,[length]	207	mov eax,[length]
208	sub eax,Buf_size	208	sub eax,Buf_size
209	jmp .end0	209	jmp .end0
210	@@: ;else	210	@@: ;else
211	mov eax,[value]	211	mov eax,[value]
212	mov ecx,[edi+deflate_state.bi_valid]	212	mov ecx,[edi+deflate_state.bi_valid]
213	shl eax,cl	213	shl eax,cl
214	or [edi+deflate_state.bi_buf],ax	214	or [edi+deflate_state.bi_buf],ax
215	mov eax,[length]	215	mov eax,[length]
216	.end0:	216	.end0:
217	add [edi+deflate_state.bi_valid],eax	217	add [edi+deflate_state.bi_valid],eax
218	;else ;!DEBUG	218	;else ;!DEBUG
219		219
220	;{ int len = length;	220	;{ int len = length;
221	; if (s->bi_valid > (int)Buf_size - len) {	221	; if (s->bi_valid > (int)Buf_size - len) {
222	; int val = value;	222	; int val = value;
223	; s->bi_buf \|= (uint_16)val << s->bi_valid;	223	; s->bi_buf \|= (uint_16)val << s->bi_valid;
224	; put_short(s, s->bi_buf);	224	; put_short(s, s->bi_buf);
225	; s->bi_buf = (uint_16)val >> (Buf_size - s->bi_valid);	225	; s->bi_buf = (uint_16)val >> (Buf_size - s->bi_valid);
226	; s->bi_valid += len - Buf_size;	226	; s->bi_valid += len - Buf_size;
227	; } else {	227	; } else {
228	; s->bi_buf \|= (uint_16)(value) << s->bi_valid;	228	; s->bi_buf \|= (uint_16)(value) << s->bi_valid;
229	; s->bi_valid += len;	229	; s->bi_valid += len;
230	; }	230	; }
231	;}	231	;}
232	;end if ;DEBUG	232	;end if ;DEBUG
233	ret	233	ret
234	endp	234	endp
235		235
236	; the arguments must not have side effects	236	; the arguments must not have side effects
237		237
238	; ===========================================================================	238	; ===========================================================================
239	; Initialize the various 'constant' tables.	239	; Initialize the various 'constant' tables.
240		240
241	;int static_init_done = 0	241	;int static_init_done = 0
242		242
243	;void ()	243	;void ()
244	align 4	244	align 4
245	proc tr_static_init	245	proc tr_static_init
246	if GEN_TREES_H eq 1	246	if GEN_TREES_H eq 1
247		247
248	; int n ;iterates over tree elements	248	; int n ;iterates over tree elements
249	; int bits ;bit counter	249	; int bits ;bit counter
250	; int length ;length value	250	; int length ;length value
251	; int code ;code value	251	; int code ;code value
252	; int dist ;distance index	252	; int dist ;distance index
253	; uint_16 bl_count[MAX_BITS+1];	253	; uint_16 bl_count[MAX_BITS+1];
254	; number of codes at each bit length for an optimal tree	254	; number of codes at each bit length for an optimal tree
255		255
256	; if (static_init_done) return;	256	; if (static_init_done) return;
257		257
258	; For some embedded targets, global variables are not initialized:	258	; For some embedded targets, global variables are not initialized:
259	;if NO_INIT_GLOBAL_POINTERS	259	;if NO_INIT_GLOBAL_POINTERS
260	; static_l_desc.static_tree = static_ltree;	260	; static_l_desc.static_tree = static_ltree;
261	; static_l_desc.extra_bits = extra_lbits;	261	; static_l_desc.extra_bits = extra_lbits;
262	; static_d_desc.static_tree = static_dtree;	262	; static_d_desc.static_tree = static_dtree;
263	; static_d_desc.extra_bits = extra_dbits;	263	; static_d_desc.extra_bits = extra_dbits;
264	; static_bl_desc.extra_bits = extra_blbits;	264	; static_bl_desc.extra_bits = extra_blbits;
265	;end if	265	;end if
266		266
267	; Initialize the mapping length (0..255) -> length code (0..28)	267	; Initialize the mapping length (0..255) -> length code (0..28)
268	; length = 0;	268	; length = 0;
269	; for (code = 0; code < LENGTH_CODES-1; code++) {	269	; for (code = 0; code < LENGTH_CODES-1; code++) {
270	; base_length[code] = length;	270	; base_length[code] = length;
271	; for (n = 0; n < (1<	271	; for (n = 0; n < (1<
272	; _length_code[length++] = (uch)code;	272	; _length_code[length++] = (uch)code;
273	; }	273	; }
274	; }	274	; }
275	; Assert (length == 256, "tr_static_init: length != 256");	275	; Assert (length == 256, "tr_static_init: length != 256");
276	; Note that the length 255 (match length 258) can be represented	276	; Note that the length 255 (match length 258) can be represented
277	; in two different ways: code 284 + 5 bits or code 285, so we	277	; in two different ways: code 284 + 5 bits or code 285, so we
278	; overwrite length_code[255] to use the best encoding:	278	; overwrite length_code[255] to use the best encoding:
279		279
280	; _length_code[length-1] = (uch)code;	280	; _length_code[length-1] = (uch)code;
281		281
282	; Initialize the mapping dist (0..32K) -> dist code (0..29)	282	; Initialize the mapping dist (0..32K) -> dist code (0..29)
283	; dist = 0;	283	; dist = 0;
284	; for (code = 0 ; code < 16; code++) {	284	; for (code = 0 ; code < 16; code++) {
285	; base_dist[code] = dist;	285	; base_dist[code] = dist;
286	; for (n = 0; n < (1<	286	; for (n = 0; n < (1<
287	; _dist_code[dist++] = (uch)code;	287	; _dist_code[dist++] = (uch)code;
288	; }	288	; }
289	; }	289	; }
290	; Assert (dist == 256, "tr_static_init: dist != 256");	290	; Assert (dist == 256, "tr_static_init: dist != 256");
291	; dist >>= 7; /* from now on, all distances are divided by 128 */	291	; dist >>= 7; /* from now on, all distances are divided by 128 */
292	; for ( ; code < D_CODES; code++) {	292	; for ( ; code < D_CODES; code++) {
293	; base_dist[code] = dist << 7;	293	; base_dist[code] = dist << 7;
294	; for (n = 0; n < (1<<(extra_dbits[code]-7)); n++) {	294	; for (n = 0; n < (1<<(extra_dbits[code]-7)); n++) {
295	; _dist_code[256 + dist++] = (uch)code;	295	; _dist_code[256 + dist++] = (uch)code;
296	; }	296	; }
297	; }	297	; }
298	; Assert (dist == 256, "tr_static_init: 256+dist != 512");	298	; Assert (dist == 256, "tr_static_init: 256+dist != 512");
299		299
300	; Construct the codes of the static literal tree	300	; Construct the codes of the static literal tree
301	; for (bits = 0; bits <= MAX_BITS; bits++) bl_count[bits] = 0;	301	; for (bits = 0; bits <= MAX_BITS; bits++) bl_count[bits] = 0;
302	; n = 0;	302	; n = 0;
303	; while (n <= 143) static_ltree[n++].Len = 8, bl_count[8]++;	303	; while (n <= 143) static_ltree[n++].Len = 8, bl_count[8]++;
304	; while (n <= 255) static_ltree[n++].Len = 9, bl_count[9]++;	304	; while (n <= 255) static_ltree[n++].Len = 9, bl_count[9]++;
305	; while (n <= 279) static_ltree[n++].Len = 7, bl_count[7]++;	305	; while (n <= 279) static_ltree[n++].Len = 7, bl_count[7]++;
306	; while (n <= 287) static_ltree[n++].Len = 8, bl_count[8]++;	306	; while (n <= 287) static_ltree[n++].Len = 8, bl_count[8]++;
307	; Codes 286 and 287 do not exist, but we must include them in the	307	; Codes 286 and 287 do not exist, but we must include them in the
308	; tree construction to get a canonical Huffman tree (longest code	308	; tree construction to get a canonical Huffman tree (longest code
309	; all ones)	309	; all ones)
310		310
311	; gen_codes((ct_data *)static_ltree, L_CODES+1, bl_count);	311	; gen_codes((ct_data *)static_ltree, L_CODES+1, bl_count);
312		312
313	; The static distance tree is trivial:	313	; The static distance tree is trivial:
314	; for (n = 0; n < D_CODES; n++) {	314	; for (n = 0; n < D_CODES; n++) {
315	; static_dtree[n].Len = 5;	315	; static_dtree[n].Len = 5;
316	; static_dtree[n].Code = bi_reverse((unsigned)n, 5);	316	; static_dtree[n].Code = bi_reverse((unsigned)n, 5);
317	; }	317	; }
318	; static_init_done = 1;	318	; static_init_done = 1;
319		319
320	if GEN_TREES_H eq 1	320	if GEN_TREES_H eq 1
321	call gen_trees_header	321	call gen_trees_header
322	end if	322	end if
323	end if ;(GEN_TREES_H) \| !(STDC)	323	end if ;(GEN_TREES_H) \| !(STDC)
324	ret	324	ret
325	endp	325	endp
326		326
327	; ===========================================================================	327	; ===========================================================================
328	; Genererate the file trees.h describing the static trees.	328	; Genererate the file trees.h describing the static trees.
329		329
330	;# define SEPARATOR(i, last, width) \	330	;# define SEPARATOR(i, last, width) \
331	; ((i) == (last)? "\n};\n\n" : \	331	; ((i) == (last)? "\n};\n\n" : \
332	; ((i) % (width) == (width)-1 ? ",\n" : ", "))	332	; ((i) % (width) == (width)-1 ? ",\n" : ", "))
333		333
334	;void ()	334	;void ()
335	align 4	335	align 4
336	proc gen_trees_header	336	proc gen_trees_header
337	; FILE *header = fopen("trees.inc", "w");	337	; FILE *header = fopen("trees.inc", "w");
338	; int i;	338	; int i;
339		339
340	; Assert (header != NULL, "Can't open trees.inc");	340	; Assert (header != NULL, "Can't open trees.inc");
341	; fprintf(header,	341	; fprintf(header,
342	; "/* header created automatically with -DGEN_TREES_H */\n\n");	342	; "/* header created automatically with -DGEN_TREES_H */\n\n");
343		343
344	; fprintf(header, "local const ct_data static_ltree[L_CODES+2] = {\n");	344	; fprintf(header, "local const ct_data static_ltree[L_CODES+2] = {\n");
345	; for (i = 0; i < L_CODES+2; i++) {	345	; for (i = 0; i < L_CODES+2; i++) {
346	; fprintf(header, "{{%3u},{%3u}}%s", static_ltree[i].Code,	346	; fprintf(header, "{{%3u},{%3u}}%s", static_ltree[i].Code,
347	; static_ltree[i].Len, SEPARATOR(i, L_CODES+1, 5));	347	; static_ltree[i].Len, SEPARATOR(i, L_CODES+1, 5));
348	; }	348	; }
349		349
350	; fprintf(header, "local const ct_data static_dtree[D_CODES] = {\n");	350	; fprintf(header, "local const ct_data static_dtree[D_CODES] = {\n");
351	; for (i = 0; i < D_CODES; i++) {	351	; for (i = 0; i < D_CODES; i++) {
352	; fprintf(header, "{{%2u},{%2u}}%s", static_dtree[i].Code,	352	; fprintf(header, "{{%2u},{%2u}}%s", static_dtree[i].Code,
353	; static_dtree[i].Len, SEPARATOR(i, D_CODES-1, 5));	353	; static_dtree[i].Len, SEPARATOR(i, D_CODES-1, 5));
354	; }	354	; }
355		355
356	; fprintf(header, "const uch ZLIB_INTERNAL _dist_code[DIST_CODE_LEN] = {\n");	356	; fprintf(header, "const uch ZLIB_INTERNAL _dist_code[DIST_CODE_LEN] = {\n");
357	; for (i = 0; i < DIST_CODE_LEN; i++) {	357	; for (i = 0; i < DIST_CODE_LEN; i++) {
358	; fprintf(header, "%2u%s", _dist_code[i],	358	; fprintf(header, "%2u%s", _dist_code[i],
359	; SEPARATOR(i, DIST_CODE_LEN-1, 20));	359	; SEPARATOR(i, DIST_CODE_LEN-1, 20));
360	; }	360	; }
361		361
362	; fprintf(header,	362	; fprintf(header,
363	; "const uch ZLIB_INTERNAL _length_code[MAX_MATCH-MIN_MATCH+1]= {\n");	363	; "const uch ZLIB_INTERNAL _length_code[MAX_MATCH-MIN_MATCH+1]= {\n");
364	; for (i = 0; i < MAX_MATCH-MIN_MATCH+1; i++) {	364	; for (i = 0; i < MAX_MATCH-MIN_MATCH+1; i++) {
365	; fprintf(header, "%2u%s", _length_code[i],	365	; fprintf(header, "%2u%s", _length_code[i],
366	; SEPARATOR(i, MAX_MATCH-MIN_MATCH, 20));	366	; SEPARATOR(i, MAX_MATCH-MIN_MATCH, 20));
367	; }	367	; }
368		368
369	; fprintf(header, "local const int base_length[LENGTH_CODES] = {\n");	369	; fprintf(header, "local const int base_length[LENGTH_CODES] = {\n");
370	; for (i = 0; i < LENGTH_CODES; i++) {	370	; for (i = 0; i < LENGTH_CODES; i++) {
371	; fprintf(header, "%1u%s", base_length[i],	371	; fprintf(header, "%1u%s", base_length[i],
372	; SEPARATOR(i, LENGTH_CODES-1, 20));	372	; SEPARATOR(i, LENGTH_CODES-1, 20));
373	; }	373	; }
374		374
375	; fprintf(header, "local const int base_dist[D_CODES] = {\n");	375	; fprintf(header, "local const int base_dist[D_CODES] = {\n");
376	; for (i = 0; i < D_CODES; i++) {	376	; for (i = 0; i < D_CODES; i++) {
377	; fprintf(header, "%5u%s", base_dist[i],	377	; fprintf(header, "%5u%s", base_dist[i],
378	; SEPARATOR(i, D_CODES-1, 10));	378	; SEPARATOR(i, D_CODES-1, 10));
379	; }	379	; }
380		380
381	; fclose(header);	381	; fclose(header);
382	ret	382	ret
383	endp	383	endp
384		384
385	; ===========================================================================	385	; ===========================================================================
386	; Initialize the tree data structures for a new zlib stream.	386	; Initialize the tree data structures for a new zlib stream.
387		387
388	;void (s)	388	;void (s)
389	; deflate_state* s	389	; deflate_state* s
390	align 4	390	align 4
391	proc _tr_init uses eax edi, s:dword	391	proc _tr_init uses eax edi, s:dword
392	mov edi,[s]	392	mov edi,[s]
393	zlib_debug '_tr_init'	393	zlib_debug '_tr_init'
394	call tr_static_init	394	call tr_static_init
395		395
396	mov eax,edi	396	mov eax,edi
397	add eax,deflate_state.dyn_ltree	397	add eax,deflate_state.dyn_ltree
398	mov [edi+deflate_state.l_desc.dyn_tree],eax	398	mov [edi+deflate_state.l_desc.dyn_tree],eax
399	mov [edi+deflate_state.l_desc.stat_desc],static_l_desc	399	mov [edi+deflate_state.l_desc.stat_desc],static_l_desc
400		400
401	add eax,deflate_state.dyn_dtree-deflate_state.dyn_ltree	401	add eax,deflate_state.dyn_dtree-deflate_state.dyn_ltree
402	mov [edi+deflate_state.d_desc.dyn_tree],eax	402	mov [edi+deflate_state.d_desc.dyn_tree],eax
403	mov [edi+deflate_state.d_desc.stat_desc],static_d_desc	403	mov [edi+deflate_state.d_desc.stat_desc],static_d_desc
404		404
405	add eax,deflate_state.bl_tree-deflate_state.dyn_dtree	405	add eax,deflate_state.bl_tree-deflate_state.dyn_dtree
406	mov [edi+deflate_state.bl_desc.dyn_tree],eax	406	mov [edi+deflate_state.bl_desc.dyn_tree],eax
407	mov [edi+deflate_state.bl_desc.stat_desc],static_bl_desc;	407	mov [edi+deflate_state.bl_desc.stat_desc],static_bl_desc;
408		408
409	mov word[edi+deflate_state.bi_buf],0	409	mov word[edi+deflate_state.bi_buf],0
410	mov dword[edi+deflate_state.bi_valid],0	410	mov dword[edi+deflate_state.bi_valid],0
411	if DEBUG eq 1	411	if DEBUG eq 1
412	mov dword[edi+deflate_state.compressed_len],0	412	mov dword[edi+deflate_state.compressed_len],0
413	mov dword[edi+deflate_state.bits_sent],0	413	mov dword[edi+deflate_state.bits_sent],0
414	end if	414	end if
415		415
416	; Initialize the first block of the first file:	416	; Initialize the first block of the first file:
417	stdcall init_block,edi	417	stdcall init_block,edi
418	ret	418	ret
419	endp	419	endp
420		420
421	; ===========================================================================	421	; ===========================================================================
422	; Initialize a new block.	422	; Initialize a new block.
423		423
424	;void (s)	424	;void (s)
425	; deflate_state* s	425	; deflate_state* s
426	align 4	426	align 4
427	proc init_block uses eax ecx edi, s:dword	427	proc init_block uses eax ecx edi, s:dword
428	mov edi,[s]	428	mov edi,[s]
429		429
430	; Initialize the trees.	430	; Initialize the trees.
431	mov eax,edi	431	mov eax,edi
432	add eax,deflate_state.dyn_ltree+Freq	432	add eax,deflate_state.dyn_ltree+Freq
433	mov ecx,L_CODES	433	mov ecx,L_CODES
434	@@:	434	@@:
435	mov word[eax],0	435	mov word[eax],0
436	add eax,sizeof.ct_data	436	add eax,sizeof.ct_data
437	loop @b	437	loop @b
438	mov eax,edi	438	mov eax,edi
439	add eax,deflate_state.dyn_dtree+Freq	439	add eax,deflate_state.dyn_dtree+Freq
440	mov ecx,D_CODES	440	mov ecx,D_CODES
441	@@:	441	@@:
442	mov word[eax],0	442	mov word[eax],0
443	add eax,sizeof.ct_data	443	add eax,sizeof.ct_data
444	loop @b	444	loop @b
445	mov eax,edi	445	mov eax,edi
446	add eax,deflate_state.bl_tree+Freq	446	add eax,deflate_state.bl_tree+Freq
447	mov ecx,BL_CODES	447	mov ecx,BL_CODES
448	@@:	448	@@:
449	mov word[eax],0	449	mov word[eax],0
450	add eax,sizeof.ct_data	450	add eax,sizeof.ct_data
451	loop @b	451	loop @b
452		452
453	mov word[edi+sizeof.ct_data*END_BLOCK+deflate_state.dyn_ltree+Freq],1	453	mov word[edi+sizeof.ct_data*END_BLOCK+deflate_state.dyn_ltree+Freq],1
454	mov dword[edi+deflate_state.static_len],0	454	mov dword[edi+deflate_state.static_len],0
455	mov dword[edi+deflate_state.opt_len],0	455	mov dword[edi+deflate_state.opt_len],0
456	mov dword[edi+deflate_state.matches],0	456	mov dword[edi+deflate_state.matches],0
457	mov dword[edi+deflate_state.last_lit],0	457	mov dword[edi+deflate_state.last_lit],0
458	ret	458	ret
459	endp	459	endp
460		460
461	SMALLEST equ 1	461	SMALLEST equ 1
462	; Index within the heap array of least frequent node in the Huffman tree	462	; Index within the heap array of least frequent node in the Huffman tree
463		463
464		464
465	; ===========================================================================	465	; ===========================================================================
466	; Remove the smallest element from the heap and recreate the heap with	466	; Remove the smallest element from the heap and recreate the heap with
467	; one less element. Updates heap and heap_len.	467	; one less element. Updates heap and heap_len.
468		468
469	macro pqremove s, tree, top	469	macro pqremove s, tree, top
470	{	470	{
471	mov eax,s	471	mov eax,s
472	add eax,deflate_state.heap+2*SMALLEST	472	add eax,deflate_state.heap+2*SMALLEST
473	movzx top,word[eax]	473	movzx top,word[eax]
474	push ebx	474	push ebx
475	mov ebx,[s+deflate_state.heap_len]	475	mov ebx,[s+deflate_state.heap_len]
476	mov bx,[s+deflate_state.heap+2*ebx]	476	mov bx,[s+deflate_state.heap+2*ebx]
477	mov word[eax],bx	477	mov word[eax],bx
478	dec dword[s+deflate_state.heap_len]	478	dec dword[s+deflate_state.heap_len]
479	pop ebx	479	pop ebx
480	stdcall pqdownheap, s, tree, SMALLEST	480	stdcall pqdownheap, s, tree, SMALLEST
481	}	481	}
482		482
483	; ===========================================================================	483	; ===========================================================================
484	; Compares to subtrees, using the tree depth as tie breaker when	484	; Compares to subtrees, using the tree depth as tie breaker when
485	; the subtrees have equal frequency. This minimizes the worst case length.	485	; the subtrees have equal frequency. This minimizes the worst case length.
486		486
487	macro smaller tree, n, m, depth, m_end	487	macro smaller tree, n, m, depth, m_end
488	{	488	{
489	;if (..<.. \|\| (..==.. && depth[n] <= depth[m]))	489	;if (..<.. \|\| (..==.. && depth[n] <= depth[m]))
490	local .end0	490	local .end0
491	movzx eax,n	491	movzx eax,n
492	imul eax,sizeof.ct_data	492	imul eax,sizeof.ct_data
493	add eax,tree	493	add eax,tree
494	mov ax,word[eax+Freq]	494	mov ax,word[eax+Freq]
495	movzx ebx,m	495	movzx ebx,m
496	imul ebx,sizeof.ct_data	496	imul ebx,sizeof.ct_data
497	add ebx,tree	497	add ebx,tree
498	mov bx,word[ebx+Freq]	498	mov bx,word[ebx+Freq]
499	cmp ax,bx	499	cmp ax,bx
500	jl .end0	500	jl .end0

Subversion Repositories Kolibri OS

(root)/programs/fs/kfar/trunk/zlib/trees.asm – Rev 6799 → 6813