Mercurial > repos > public > sbplib

annotate +scheme/+bc/closureSetup.m @ 873:dee5b5a57be6 bcSetupExperiment

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Clean up closureSetup
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 01 Nov 2018 12:05:53 +0100
parents fb91d12093f8
children ba10f24bf476
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
873
dee5b5a57be6 Clean up closureSetup
Jonatan Werpers <jonatan@werpers.com>
parents: 870
diff changeset
1 % Setup closure and penalty matrices for several boundary conditions at once.
870
fb91d12093f8 Change some naming of functions
Jonatan Werpers <jonatan@werpers.com>
parents: 869
diff changeset
2 function [closure, penalties] = closureSetup(diffOp, bcs)
873
dee5b5a57be6 Clean up closureSetup
Jonatan Werpers <jonatan@werpers.com>
parents: 870
diff changeset
3 scheme.bc.verifyFormat(bcs, diffOp);
869
d356f1a22d4f Start organizing the code
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4
d356f1a22d4f Start organizing the code
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5 % Setup storage arrays
d356f1a22d4f Start organizing the code
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6 closure = spzeros(size(diffOp));
d356f1a22d4f Start organizing the code
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7 penalties = cell(1, length(bcs));
d356f1a22d4f Start organizing the code
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
8
d356f1a22d4f Start organizing the code
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
9 % Collect closures and penalties
d356f1a22d4f Start organizing the code
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10 for i = 1:length(bcs)
d356f1a22d4f Start organizing the code
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
11 [localClosure, penalties{i}] = diffOp.boundary_condition(bcs{i}.boundary, bcs{i}.type);
d356f1a22d4f Start organizing the code
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
12 closure = closure + localClosure;
d356f1a22d4f Start organizing the code
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
13 end
d356f1a22d4f Start organizing the code
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
14 end