% Prop8MultiplierInversionOpGabor: 
% Inversion of Gabor Frame multiplier according to Prop. 8 for Gabor frames
% involves the program Prop8InvMultOpGabor(L,gPhi,a,M,gPsi,m,e)
% requires entry of the input-parameters
%function y= Prop8MultiplierInversionOpGabor;
ltfatstart;
L=input('Please enter the length L of the transform: ');
a=input('Please enter the time-shift a (should be divisor of L) for the Gabor frame Phi : ');
M=input('Please enter the number M of channels (should be divisor of L and bigger or equal to aPhi) for the Gabor frame Phi: ');
%so, we will have m (at frequency) from 0 to M-1 and n (at time) from 0 to N-1, where N=L/a,
gPhi=input('Please enter the window function of the Gabor frame Phi (e.g. gauss, hann, hamming, between text-brackets): ');
gG=input('Please enter the window function of the Gabor frame G(Psi-Phi): ');
m=input('Please enter the elements of m (ML/a positive numbers in a row): ');
mpos=1;
while(mpos==1)
am1=min(m);
    if am1<=0
       m=input('Your m contains non-positive elements. Please enter new m (c positive numbers in a row): ');
    else
        mpos=0;
    end
end
e=input('Please enter the error bound you aim: ');
v=1;
while(v==1)
   [TPhi,TPsi,M1,M2,M1inv,M2inv,n]=Prop8InvMultOpGabor(L,a,M,gPhi,gG,m,e)
  gG=input('Please enter the window function of the Gabor frame G or press zero if you want to stop: ');
   if gG==0 
       v=2;
       break
   end
 e=input('Please enter the error bound you aim: ');
end


   
%end