% This computes function value for IPM_p. 
%
% Sept 26, 2011

function fval = function_value(X,w,p,mu)

D = svd(X);

fval = sum(D.^p) - mu*sum(log(w)) - mu*log(1 - sum(w));