Mercurial > repos > public > sbplib

annotate +grid/primalDual1DTest.m @ 1259:99f92bfc1157 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add staggered operators and 1D grids
author Martin Almquist <malmquist@stanford.edu>
date Tue, 14 Apr 2020 17:50:49 -0700
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1259
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
1 function tests = primalDual1DTest()
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
2 tests = functiontests(localfunctions);
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
3 end
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
4
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
5
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
6 function testErrorInvalidLimits(testCase)
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
7 in = {
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
8 {10,{1}},
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
9 {10,[0,1]},
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
10 {10,{1,0}},
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
11 };
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
12
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
13 for i = 1:length(in)
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
14 testCase.verifyError(@()grid.primalDual1D(in{i}{:}),'grid:primalDual1D:InvalidLimits',sprintf('in(%d) = %s',i,toString(in{i})));
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
15 end
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
16 end
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
17
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
18 function testCompiles(testCase)
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
19 in = {
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
20 {5, {0,1}},
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
21 };
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
22
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
23 out = {
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
24 {[0; 0.25; 0.5; 0.75; 1], [0; 0.125; 0.375; 0.625; 0.875; 1]},
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
25 };
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
26
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
27 for i = 1:length(in)
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
28 [gp, gd] = grid.primalDual1D(in{i}{:});
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
29 testCase.verifyEqual(gp.points(),out{i}{1});
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
30 testCase.verifyEqual(gd.points(),out{i}{2});
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
31 end
99f92bfc1157 Add staggered operators and 1D grids
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
32 end