error_indicator.m

% Calculate the errors for the basis stored in the matrix RB
function [errs, avg_time, errs_abs, errs_true] = error_indicator(model, model_data, RB_W, RB_V, training_mus)
    errs = zeros(size(training_mus,1),1);
    errs_abs = zeros(size(training_mus,1),1);

    dd = ARE.detailed_data;
    dd.RB_V = RB_V;
    dd.RB_W = RB_W;

    rd = model.reduce(model_data, dd);
    t = tic;
    parfor i = 1:size(training_mus,1)
        tmpModel = model;
        tmpModel.set_mu(training_mus(i, :));
        tmpModel.RB_error_estimation = false;
        tmpModel.RB_calc_residual = true;
        try
            rsim = rb_simulation(tmpModel, rd);
            errs(i) = rsim.nresidual;
            errs_abs(i) = rsim.residual;
        catch e
            training_mus(i,:)
            warning(['Solution not calculable: ', e.message])
            errs(i) = -Inf;
        end
    end

    avg_time = toc(t);
end