Mercurial > repos > public > sbplib

annotate semiDefIneq.m @ 1012:1e437c9e5132 feature/advectionRV

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Create residual viscosity package +rv and generalize the ResidualViscosity class - Generalize residual viscosity, by passing user-defined flux and calculating the time derivative outside of the update. - Create separate RungekuttaRV specifically using interior RV updates - Separate the artifical dissipation operator from the scheme AdvectionRV1D so that the same scheme can be reused for creating the diff op used by the ResidualViscosity class
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 05 Dec 2018 13:44:10 +0100
parents eaf557023fbe
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
234
eaf557023fbe Added verbose flag to semiDefIneq and minors.
Jonatan Werpers <jonatan@werpers.com>
parents: 129
diff changeset
1 function ineq = semiDefIneq(A, verbose)
eaf557023fbe Added verbose flag to semiDefIneq and minors.
Jonatan Werpers <jonatan@werpers.com>
parents: 129
diff changeset
2 default_arg('verbose', true);
eaf557023fbe Added verbose flag to semiDefIneq and minors.
Jonatan Werpers <jonatan@werpers.com>
parents: 129
diff changeset
3 [m, sub] = minors(A, verbose);
129
0a881a3dc9a9 Cleaned up the code for minors.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4
0a881a3dc9a9 Cleaned up the code for minors.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5 ineqsys = true;
0a881a3dc9a9 Cleaned up the code for minors.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6 for i = 1:length(m)
0a881a3dc9a9 Cleaned up the code for minors.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7 ineqsys = ineqsys & m(i) >= 0;
0a881a3dc9a9 Cleaned up the code for minors.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
8 end
0a881a3dc9a9 Cleaned up the code for minors.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
9
0a881a3dc9a9 Cleaned up the code for minors.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10 ineq = simplify(ineqsys);
0a881a3dc9a9 Cleaned up the code for minors.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
11
0a881a3dc9a9 Cleaned up the code for minors.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
12 str = toString(ineq);
0a881a3dc9a9 Cleaned up the code for minors.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
13 fprintf('%s\n',strjoin(strsplit(str,' & '), '\n'));
0a881a3dc9a9 Cleaned up the code for minors.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
14 end