% Prop8MultiplierInversionf: Computation of M^{-1}f according to Prop. 8 
% involves the program Prop8InvMultf(c,r,TPhi,TG,m,f,e)
% requires entry of the input-parameters
%function y= Prop8MultiplierInversionf;
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 the frame Phi: ');
TG=input('Please enter the synthesis matrix (rxc) of the frame G (Psi-Phi): ');
m=input('Please enter the elements of the positive symbol m (c positive 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
f=input('Please enter the coordinates of f (r numbers as columns): ');
e=input('Please enter the error bound you aim: ');

v=1;
while(v==1)
   [TPsi,M1,M2,M1invf,M2invf,n]=Prop8InvMultf(c,r,TPhi,TG,m,f,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