function [yQ, totalBits] = performEncoding(y, Fs, bits, FFTLength, M, F, k, thres, totalBits, quantType);

%
% quantType = -1 : full-range quantization
%				   0 : narrow quantization
%              8 : 8-bit full-range quantization
%             16 : 16-bit full-range quantization
%             99 : ATH full-range quantization

% find low and high frequencies for this filter in the bank
center = (2*k - 1) * Fs * pi/(2*M);
low = center - Fs*pi/(2*M);
high = center + Fs*pi/(2*M);

% find the index into the frequency array that corresponds to low and high
lowIndex = max(find(F<low));
if (isempty(lowIndex))
   lowIndex = 1;
end

highIndex = max(find(F<high));

% range of amplitude levels (in our case 65536)
R = 2^bits;

% find the number of bits needed to quantize with by finding the minimum of the 
% masking threshold over the frequency range of this filter
if quantType ~= 99   
   b = ceil(log2(R./min(thres(lowIndex:highIndex)')) - 1);
else   
   ATH=(3.64.*(F./1000).^(-.8) - 6.5.*exp(-0.6.*(F./1000-3.3).^2) + 0.001.*(F./1000).^4);
   b = ceil(log2(R./min(ATH(lowIndex:highIndex)')) - 1);   
end

% perform narrow-range quantization
if quantType == 0
   [yQ, scaling] = quantization(y, b);
	yQ = invquantization(yQ, scaling);
else
   % perform 8-bit quantization
   if quantType == 8
      b = 8;
   else 
      % perform 16-bit quantization
      if quantType == 16
         b = 16;
      end
   end
   
   % perform full-range quantization
	yQ = q(y, b);
   yQ = iq(yQ, b)';   
end   
  
if (quantType == -1) | (quantType == 8) | (quantType == 16) | (quantType == 99)
   % calculate how many bits we need for full-range quantization  
   newBits = log2(bits) + b * FFTLength/M;
else
   % calculate how many bits we need for narrow quantization  
	newBits = 30 + b * FFTLength/M;
end

% add the bit count to the total
totalBits = totalBits + newBits;