diff_rank_rho_pcp.m

clear;

config = diff_rank_rho_config();

method = 'pcp';

all_d = config.d;
all_rho = config.rho;
all_rep = config.repetitions;

p = config.p;
n = config.n;

%% add path
addpath('PROPACK/');
addpath('inexact_alm_rpca');

%% compute EV for PCP

for k=1:length(all_rep)
    rep = all_rep(k);
    
    for j=1:length(all_d)
        d = all_d(j);
        
        for i=1:length(all_rho)
            rho = all_rho(i);
            
            fprintf('PCP: rep = %d, d = %d, rho = %.2f\n', rep, d, rho);
            
            result_file = sprintf(config.result_file_format, method, method, d, rho, rep);            
            
            data_file = sprintf(config.data_file_format, d, rho, rep);
            load(data_file);
            
            U_all = cell2mat(U);
            UUt = U_all * U_all';
            traceUUt = trace(UUt);
            
            [X, ~, ~] = inexact_alm_rpca(Z);
            
            [UX, ~, ~] = svds(X, d);
            
            ev = js_compute_EV(UX, UUt, traceUUt);
            
            EV = ev * ones(1, n);
            
            save(result_file, 'EV');
            fprintf('save to %s\n', result_file);
            
        end
    end
end