Mercurial > repos > public > sbplib
view semiDefIneq.m @ 1227:02dfe3a56742 feature/poroelastic
Add Upwind ElasticAnisotropic schemes. Seem to work really well!
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Sat, 16 Nov 2019 14:26:06 -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