function U=ummp(L)
% U=UMMP(L)
%
% INPUT:
%
% L      Spherical harmonic bandwidth
%
% OUTPUT:
%
% U      Collection of transformation matrices
%
% The unitary block-transformation matrix that maps real to complex
% harmonics, for all degrees 0:L. 
% Note that to convert series of coefficients, you need the transpose!
% See also ULMMP.
%
% Last modified by fjsimons-at-alum.mit.edu, 09/25/2006

defval('L',3)

U=repmat(0,(L+1)^2,(L+1)^2);

e=0;
for degree=0:L
  b=e+1;
  e=b+(2*degree+1)-1;
  U(b:e,b:e)=ulmmp(degree);
end