function    [tdX,tdLam] = FpsFp(dX,dLam,X,Lam,V,H)
X=sparse(X);
V=sparse(V);
H=sparse(H);
dX=sparse(dX);
DLam=sparse(dLam);
Lam=sparse(Lam);
[tFd,tFc]=Fpmu(dX,dLam,X,Lam,V,H);
[tdX,tdLam]=Fpmus(tFd,tFc,X,Lam,V,H);
