wolffd@0: function result = threshFull_mixed(A, lam)
wolffd@0: %return argmin_V lam*||V||_2,1 + 0.5 ||V - A||_F^2
wolffd@0: % no symmetry constraint on V
wolffd@0: 
wolffd@0: norm = sqrt(sum(A.^2));
wolffd@0: result = bsxfun(@times, A, max(1 - (lam ./ norm),0));
wolffd@0: 	
wolffd@0: end