function [w,info,part] = perform_alpert_transform_nd(v,pos,k,dir, options)

% perform_alpert_transform_nd - transform a nD signal using a nD Alpert basis.
%
%   [w,info] = perform_alpert_transform_nd(v,pos,k, dir, part)
%
%   'pos' is a 2D vector, pos(:,i) is the ith point.
%   'k' is the number of vanishing moments (1=>Haar, 2=>linear basis ...).
%       * 'k' can be an integer, and then the algorithm will use the same
%         order for X and Y, etc. direction
%       * 'k' can be a vector of integers k=[kx,ky,...] and 'kx' will be the
%         order on the X direction, and 'ky' the order on the Y direction, etc.
%   'dir' is 1 for fwd transform and -1 for bwd.
%
%   'w' is the transformed data.
%   'info' is a struct containing the localisation information for each
%       basis Alpert vector.
%       'info.l' is the scale of the vector (0=coarse scale).
%       'info.n' is the space location of the vector.
%       'info.k' is the number of multiwavelet (in [1,...,k(1)*k(2)]).
%       'info.s' is the number of the slice.
%
%   Copyright (c) 2004 Gabriel Peyr