% Prop8MultiplierInversionOp: 
% Inversion of Frame multiplier according to Prop. 8 for Gabor frames
% involves the program Prop8InvMultOpGabor(c,r,TPhi,TG,m,e)
% requires entry of the input-parameters
%function y= Prop8MultiplierInversionOp;
ltfatstart;
c=input('Please enter the number c of the frame vectors: ');
r=input('Please enter the number r of the coordinates of the frame vectors: ');
TPhi=input('Please enter the synthesis matrix (rxc) of a frame Phi: ');
TG=input('Please enter the synthesis matrix (rxc) of a frame G (Psi-Phi): ');
m=input('Please enter the elements of m (c numbers in a row): ');
mpos=1;
while(mpos==1)
a=min(m);
    if a<=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)
   [TPsi,M1,M2,M1inv,M2inv,n]=Prop8InvMultOp(c,r,TPhi,TG,m,e)
  TG=input('Please enter the synthesis matrix of the frame G or press zero if you want to stop: ');
   if TG==0 
       v=2;
    break 
   end
 e=input('Please enter the error bound you aim: ');
end


   
%end