Mercurial > repos > public > sbplib

view semiDefIneq.m @ 725:e9e15d64f9f9 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Implement symmetric Dirichlet BC with variable tuning parameter in Heat2D scheme.
author Martin Almquist <malmquist@stanford.edu>
date Fri, 30 Mar 2018 15:35:46 -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