riccati_stability_constant.m

Go to the documentation of this file.

function nu = riccati_stability_constant(model, model_data, sim)
%function nu = riccati_stability_constant
% This function can be used to calculate the stability constant for the
% model. 
%
% Note that this function uses high dimensional calculations!
%
%  Andreas Schmidt, 2015
%

[E,A,~,~] = riccati_assemble_system_matrices(model,model_data);
if ~issymmetric(E)
    error('Currently only symmetric E is supported')
end
nu = eigs(A+A', E, 1, 'LA')/2;

end