PowerIteration.m

classdef PowerIteration < ARE.GammaCalculation.GammaCalculatorInterface
    % Perform a power itertion for the calculation of gamma
    
    methods
        function data = gen_detailed_data(this, model, model_data)
        end
        function data = gen_reduced_data(this, model, model_data)
        end
        function g = calculate(this, model, model_data, dsim)
            g = -1;
            [E,A,B,C,Q,R] = model.assemble(model_data);
            Y = A - sparse(B*inv(R)*B'*dsim.Z*dsim.Z'*E);
            [L,U] = lu(Y);
            
        end
    end
end