function [y,nbr_bits] = perform_wavelet_arithmetic_coding(x, T, options)

% perform_wavelet_arithmetic_coding - arithmetic code a wavelet transformed image.
%
% Coding:
%   [stream,nbr_bits] = perform_wavelet_arithmetic_coding(MW, T, options);
% decoding:
%   MW = perform_wavelet_arithmetic_coding(stream, T, options);
%
%	This will calls the function perform_arithmetic_coding through the scales.
%   * options.coder_type: see perform_arithmetic_coding.
%   * options.ordering_mode: specifies the ordering of coefficients
%       1 : zig-zag ordering of each scale (respect both scale and space locality).
%       2 : Pack each band at a time (keep only scale locality).
%       3 : flush in line (destroy scale locality), should not be used.
%   * options.coding_mode:
%       1 : code independantly each scale.
%       2 : code at once all the scales.
%
%   You should provide options.Jmin
%
%   Copyright (c) 2005 Gabriel Peyr