Ndependence_launch_solvers.m

% use a varitey of algorithms to find the maximum likelihood process 
% from various datasets 
% parpool(3)
%% check global variable set
if exist('ensemble')
   fprintf(['ensemble = ',ensemble])
else
    error('you must set the ensemble variable to either qp (quasi pure) or fr (full rank)')
end

if exist('drange')
    fprintf(['drange = ',drange])
else
    error('you must set the drange variable')
end

if exist('LIswitch')
    fprintf(['LIswitch = ',LIswitch])
else
    error('you must set the LIsiwtch variable (if 1 Linear Inversion is run, if 0 it is not)')
end

if exist('ensemble_size')
    fprintf(['ensemble_size = ',ensemble_size])
else
    error('you must set the ensemble_size variable')
end

if exist('Npows')
    fprintf(['Npows = ',ensemble_size])
else
    error('you must set the Npows variable')
end


%%

for d=drange
    for method={'mosek','gdapB','DIA','LinInversion'}
%     for method={'LinInversion'}
% for method={'DIA'}
% for method = {'DIA'}
        fprintf(char(10));
        fprintf(method{1});

        fprintf(char(10));
        fprintf('%d :', d);

        A = PM_minimal(d);
%          A = GGMall_IO(d);
        for Npow=Npows
%         for Npow=[Inf]
            N = 10^Npow;
            dir = sprintf('./Ndependence_benchmarking_results/d%i/Npow%i',d,Npow);
            fprintf(char(10));
            fprintf('[%d] ', N);
            for i = 1:ensemble_size % total number of simulated datasets
                fprintf('%d ', i); 
                load([dir,'/dataset',num2str(i)]);

                    switch char(method)
                        case'gdapB'
                            tic;
                            [choi_ml_vec, solution, costs] = gdapB(A,n);
                            elapsedTime = toc;
    %                         semilogy(costs)
    %                         hold on
                        case 'DIA'
                            tic;
                            [choi_ml_vec, solution, costs] = DIA(A,n);
                            elapsedTime = toc;
%                             semilogy(costs)
%                             hold on
                        case 'mosek'
                            cvx_solver mosek
                            tic;
                            [choi_ml_vec] = cvx_wrapper(A,n);
                            elapsedTime = toc;
                            solution=[]; costs = []; % cannot currently extract these
                        case 'sdpt3'
                            cvx_solver sdpt3
                            tic;
                            [choi_ml_vec] = cvx_wrapper(A,n);
                            elapsedTime = toc;
                            solution=[]; costs = []; % cannot currently extract these
                        case 'sedumi'
                            cvx_solver sedumi
                            tic;
                            [choi_ml_vec] = cvx_wrapper(A,n);
                            elapsedTime = toc;
                            solution=[]; costs = []; % cannot currently extract these
                        case 'LinInversion'
                            tic;
                            choi_ml_vec = LinInversion(A,n);
                            elapsedTime = toc;
                            solution = []; costs = [];
                    end
                    save([dir,'/',char(method),'_results',num2str(i)],'elapsedTime','choi_ml_vec','costs','solution')
            end
        end
    end
end