Mercurial > repos > public > sbplib

view semiDefIneq.m @ 795:1f6b2fb69225 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Revert bcSetup and update bc functions in elastic schemes to be compatible.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 25 Jul 2018 18:33:03 -0700
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