Mercurial > repos > public > sbplib

annotate operator_def/assemble_opTest.m @ 826:b94bb6ffa38b feature/operator_files

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
rename d2_2 to D1_standard_2
author Ylva Rydin <ylva.rydin@telia.com>
date Mon, 10 Sep 2018 17:48:48 +0200
parents 32c360bb480e
children d1e5143d67ed
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
821
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
1
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
2 function tests = assemble_opTest()
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
3 tests = functiontests(localfunctions);
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
4 end
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
5
825
32c360bb480e Add test and fix typo
Ylva Rydin <ylva.rydin@telia.com>
parents: 821
diff changeset
6 function TestAssemble_op(testCase)
821
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
7 m = 10;
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
8 op = sbp.D2Standard(m,{0 1},4);
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
9 h = op.h;
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
10
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
11 boundary_block = op.D1(1:4,1:6)*h;
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
12 inner = op.D1(5,3:7)*h;
825
32c360bb480e Add test and fix typo
Ylva Rydin <ylva.rydin@telia.com>
parents: 821
diff changeset
13 D1_new = assemble_op(inner,boundary_block,m)/h;
821
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
14
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
15 verifyEqual(testCase,D1_new,op.D1,'AbsTol',1e-10)
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
16
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
17 end
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
18
825
32c360bb480e Add test and fix typo
Ylva Rydin <ylva.rydin@telia.com>
parents: 821
diff changeset
19 function TestAssembleD1(testCase)
32c360bb480e Add test and fix typo
Ylva Rydin <ylva.rydin@telia.com>
parents: 821
diff changeset
20 m = 10;
32c360bb480e Add test and fix typo
Ylva Rydin <ylva.rydin@telia.com>
parents: 821
diff changeset
21 op = sbp.D2Standard(m,{0 3},2);
32c360bb480e Add test and fix typo
Ylva Rydin <ylva.rydin@telia.com>
parents: 821
diff changeset
22 h = op.h;
821
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
23
826
b94bb6ffa38b rename d2_2 to D1_standard_2
Ylva Rydin <ylva.rydin@telia.com>
parents: 825
diff changeset
24 [D1,e_l,e_r] = assemble_D1('D1_standard_2',h,m);
821
95c26000c0ba Add files for parsing operator data stencil file
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
25
825
32c360bb480e Add test and fix typo
Ylva Rydin <ylva.rydin@telia.com>
parents: 821
diff changeset
26 verifyEqual(testCase,D1,op.D1,'AbsTol',1e-10)
32c360bb480e Add test and fix typo
Ylva Rydin <ylva.rydin@telia.com>
parents: 821
diff changeset
27 verifyEqual(testCase,e_l,op.e_l,'AbsTol',1e-10)
32c360bb480e Add test and fix typo
Ylva Rydin <ylva.rydin@telia.com>
parents: 821
diff changeset
28 verifyEqual(testCase,e_r,op.e_r,'AbsTol',1e-10)
32c360bb480e Add test and fix typo
Ylva Rydin <ylva.rydin@telia.com>
parents: 821
diff changeset
29 end