% test out operators
n1=n-1;
tn1=(n)*(n-1)/2;   % for dim n-1
indn1=find(triu(ones(n1))==1);   % indices for upper triang. part
%  test upper part
Ku=Fpmuup(n1);
dX=randn(n1);
dX=dX+dX';
vdX=dX(indn1);
dLam=randn(n1);
dLam=dLam+dLam';
vdLam=dLam(indn1);
vds=[vdX
     vdLam];
resFmat=Ku*vds;
resFop=2*KWs(KW(dX))-dLam;
norm(resFmat-resFop(indn1))
%  test lower part
X=randn(n1);X=X*X';X=X+X';
Lam=randn(n1);Lam=Lam*Lam';Lam=Lam+Lam';
Kd=Fpmulow(Lam,X);
resFmat=Kd*vds;
resFop=Lam*dX + dLam*X;
norm(resFmat-reshape(resFop,n1*n1,1))
