Mercurial > repos > public > sbplib

view semiDefIneq.m @ 1115:0e5dc5dfa559 feature/laplace_curvilinear_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Temporary ugly fix for reversing direction of boundary in interface coupling
author Martin Almquist <malmquist@stanford.edu>
date Mon, 29 Apr 2019 12:55:51 -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