function B= param_2dconv_B(Bml1, Bml2)% function B= param_2dconv_b1(Bml1, Bml2)%   This file is part of the TFPM toolbox v1.0 (c)%   michael.jachan@tuwien.ac.at and underlies the GPL.% % Generates the TBT matrix of B.%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear;tfpm;MAR  =   7;LAR  =   5;MMA  =   5;LMA  =   3;N    = 128;re_im= 'i';mo_no= 'n';tfpm_file_gen;%-------------beta= -1/2;Bml1= Aml;Bml2= Bml;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Dimensions:[M1, L1]= param_dim(Bml1);[M2, L2]= param_dim(Bml2);M= M1+M2;L= L1+L2;B= zeros((M+1)*(2*L+1), (M2+1)*(2*L2+1));% Compute B1for m= 0:M   for l= -L:L      for mp= 0:M2         for lp= -L2:L2	    B((m*(2*L+1))+l+L+1, mp*(2*L2+1)+lp+L2+1)= param_get(Bml1, m-mp, l-lp);	 end;      end;   end;end;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Test B:norm(param_destack_ml(B*param_stack_ml(Bml2), M, L)-conv2(Bml1, Bml2))%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%