1 function infsup_beta = get_infsup(model,K,W)
2 %
function cont_gamma = get_continuity(model,A,W)
4 % This
function computes the inf-sup-constant beta
for a discrete
5 %
operator K with respect to the inner product matrix (mass matrix) W.
7 % A generalized eigenvalue problem is solved
for obtaining the continuity
11 % - K: Operator Matrix
12 % - W: inner product matrix
15 % - infsup_beta : continuity constant.
17 % This
function does exactly the same as get_coercivity.m
18 % but returns the smallest eigenvalue.
23 infsup_beta = sqrt(eigs(Y,W,1,'sm
'));