function [selx,sely] = compute_quadrant_selection(j,q)

% compute_quadrant_selection - compute the indice
%   for selecting a given quadrant at a given scale
%   of a wavelet transform.
%
%   q=1 is for phi/psi (on X / on Y)
%   q=2 is for psi/phi (on X / on Y)
%   q=3 is for psi/psi (on X / on Y)
%
%   % Suppose that MW contains a wavelet transform 
%   j = 4; q = 2;
%   [selx,sely] = compute_quadrant_selection(j,q);
%   MWs = MW(selx,sely);    % MWs contains the quadrant.
%
%   Copyright (c) 2004 Gabriel Peyr