Mercurial > repos > public > sbplib

view semiDefIneq.m @ 1295:cb053fabbedc feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add displacement BC on normal tangential components, but it is only provably stable if traction bc are used for the other component.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 02 Jul 2020 20:05:50 -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