function    value = Fceval(muu,Lam,X,B,V,H,sigmaa)
%  Call:    value = Fceval(muu,Lam,X,B,V,H,sigmaa)
%   calculate complementary (centered) slackness
%  sigmaa for centering for efficency, e.g. .5 rather than 1
X=sparse(X);
Lam=sparse(Lam);
n1=size(Lam,1);
value = Lam*X - sigmaa*muu * speye(n1);
value=sparse(value);
