function    [Fd,Fc] = Fpmu(dLam,dX,Lam,X,V,H)
X=sparse(X);
V=sparse(V);
H=sparse(H);
dX=sparse(dX);
DLam=sparse(dLam);
Lam=sparse(Lam);
Fd= sparse(2*KWs(KW(dX,V,H),V,H)-dLam);
Fc = sparse(Lam*dX + dLam*X);
