Mercurial > repos > public > sbplib

view semiDefIneq.m @ 680:cd1a76c38565 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove schemes elasticDilation and elasticShear because they are subsets of elastic.
author Martin Almquist <malmquist@stanford.edu>
date Mon, 05 Feb 2018 14:58:48 -0800
parents eaf557023fbe
children
line wrap: on
line source

function ineq = semiDefIneq(A, verbose)
    default_arg('verbose', true);
    [m, sub] = minors(A, verbose);

    ineqsys = true;
    for i = 1:length(m)
        ineqsys = ineqsys & m(i) >= 0;
    end

    ineq = simplify(ineqsys);

    str = toString(ineq);
    fprintf('%s\n',strjoin(strsplit(str,' & '), '\n'));
end