Mercurial > repos > public > sbplib

view semiDefIneq.m @ 718:71aa5828cbbf feature/utux2D

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add Schrödinger scheme for 2d single block. Will develop to multiblock with interpolation.
author Martin Almquist <malmquist@stanford.edu>
date Sat, 03 Mar 2018 16:18:33 -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