thesis_ex5_nesdisred_vs_camd_connected_components.m

% Experiment 5: nesdisred - attempt to establish param values
testcases = dir(fullfile('test_suites', 'basic', '*.mat'));
testcases = testcases(20:end-5);
tc_count = length(testcases);

RESULTS_DIR = fullfile('results', 'ex5_nesdisred_param_sweep');
mkdir(RESULTS_DIR);

EXPORTS_DIR = fullfile('exports', 'ex5_nesdisred_param_sweep');
mkdir(EXPORTS_DIR);

REFERENCE_DIR = fullfile('results', 'reference');

% Test matrix definition - reduction params
testMatrix_costFuncs_empty = {[]};
testMatrix_costFuncs = {@cost_res2_nodes};  % @count_resistors,

testMatrix_nesDisRedDummy = {[5 0 1 0]};
testMatrix_nesDisRedCamd = {[1000 0 1 0]};
testMatrix_emptyNesDisOpts = {[]};

testMatrix = [
    allcomb({'nesdis_camd'}, testMatrix_costFuncs, testMatrix_nesDisRedCamd);
    allcomb({'camd'}, testMatrix_costFuncs, testMatrix_emptyNesDisOpts);
    ];

for c = testMatrix'
    elimScheme = c{1};
    costFunc = c{2};
    nesdisOpts = c{3};
    
    if isempty(costFunc) costFuncName = ''; else costFuncName = func2str(costFunc); end
    if isempty(nesdisOpts) nesDisSuffix = ''; else nesDisSuffix = num2str(nesdisOpts(1)); end
    results_group = strjoin({elimScheme costFuncName num2str(nesDisSuffix)}, '-');
    
    % build struct
    options = struct;
    options.verbose = 0;
    options.early_exit = 1;
    options.nodewise_algorithm = elimScheme;
    options.graph_algorithm = 'none';
    if ~isempty(costFunc)
        options.cost_function = costFunc;
    end
    if ~isempty(costFunc)
        options.nesdis_opts = nesdisOpts;
    end
    
    reducer_warmup(options);
    
    % Run testcases
    for f = testcases'
        fname = fullfile('test_suites', 'basic', f.name);
        [tc_dir, tc_name, ext] = fileparts(f.name);
        output_filename = @(suffix) fullfile(RESULTS_DIR, [tc_name suffix]);
        
        if exist(output_filename([results_group '.mat']), 'file') == 2
            continue
        end
        tc_name
        load(fname);
        input_circuit_info = circuit_info(G, is_ext_node);
        
        tic
        output = reducer(G, is_ext_node, options);
        output.total_time = toc;
        
        i = input_for_circuit_composite_info(output);
        output_circuit_info = circuit_composite_info(i{:});
        
        i = input_for_dump(output);
        netlists.dump_composite(output_filename([results_group '.cir']), i{:});
        ngspice.run(output_filename([results_group '.cir']));
        save(output_filename([results_group '.mat']), 'input_circuit_info', 'output', 'output_circuit_info', 'G', 'is_ext_node');
    end
end

selectMatrix = [
    allcomb({'nesdis_camd'}, testMatrix_costFuncs, testMatrix_nesDisRedCamd);
    allcomb({'camd'}, testMatrix_costFuncs, testMatrix_emptyNesDisOpts);
    ];

group_count = size(selectMatrix, 1);

time_solve_orig = nan(tc_count, group_count);
time_solve_red = nan(tc_count, group_count);
num_nodes_orig = nan(tc_count, 1);
num_nodes_red = nan(tc_count, group_count);
num_term_orig = nan(tc_count, 1);
num_term_red = nan(tc_count, group_count);
num_res_orig = nan(tc_count, 1);
num_res_red = nan(tc_count, group_count);
red_total_time = nan(tc_count, group_count);
ee = nan(tc_count, group_count);
sim_time = nan(tc_count, group_count);

tc_num = 0;
for f = testcases'
    [tc_dir, tc_name, ext] = fileparts(f.name);
    output_filename = @(suffix) fullfile(RESULTS_DIR, [tc_name suffix]);
    
    tc_num = tc_num+1;
    group_num = 0;
    for c=selectMatrix'
        group_num = group_num+1;
        
        elimScheme = c{1};
        costFunc = c{2};
        nesdisOpts = c{3};
        if isempty(costFunc) costFuncName = ''; else costFuncName = func2str(costFunc); end
        if isempty(nesdisOpts) nesDisSuffix = ''; else nesDisSuffix = num2str(nesdisOpts(1)); end
        results_group = strjoin({elimScheme costFuncName num2str(nesDisSuffix)}, '-');
        
        try
            load(output_filename([results_group '.mat']),  'input_circuit_info', 'output', 'output_circuit_info', 'G', 'is_ext_node');
        catch
            continue
        end
        
        num_nodes_orig(tc_num, 1) = input_circuit_info.num_nodes;
        num_nodes_red(tc_num, group_num) = output_circuit_info.num_nodes;
        num_term_orig(tc_num, 1) = input_circuit_info.num_external;
        num_term_red(tc_num, group_num) = output_circuit_info.num_external;
        num_res_orig(tc_num, 1) = input_circuit_info.num_resistors;
        num_res_red(tc_num, group_num) = output_circuit_info.num_resistors;
        red_total_time(tc_num, group_num) = output.c{1}.nodewise_processing_time;
        sim_time(tc_num, group_num) = ngspice.get_simulation_time(output_filename([results_group '.log']));
        ee(tc_num, group_num) = output.c{1}.early_exit;
    end
end

speed_improve_pct = 100 * (sim_time(:,2) ./ sim_time(:,1) - 1);

T = table(num_nodes_orig,num_term_orig,num_res_orig,...
    num_nodes_red,num_res_red, red_total_time, ...
    sim_time, speed_improve_pct, ee, ...
    'RowNames', arrayfun(@(s) s.name, testcases, 'UniformOutput', false))
writetable(T, fullfile(EXPORTS_DIR, 'nesdisred_param_sweep.csv'))