Mercurial > repos > public > sbplib

comparison +grid/evalOnTest.m @ 157:ea8103ad2cc5 feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added function eval on.
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 19 Feb 2016 14:28:27 +0100
parents
children 4c3f55a628c8
comparison
equal deleted inserted replaced
156:095cdd472688 157:ea8103ad2cc5
1 function tests = evalOnTest()
2 tests = functiontests(localfunctions);
3 end
4
5 function testInputConstant(testCase)
6 in = {
7 0,
8 47,
9 1,
10 [1; 2],
11 };
12
13 out = {
14 [0; 0; 0],
15 [47; 47; 47],
16 [1; 1; 1],
17 [1; 2; 1; 2; 1; 2],
18 };
19
20 g = getTestGrid('1d');
21
22 for i = 1:length(in)
23 gf = grid.evalOn(g,in{i});
24 testCase.verifyEqual(gf, out{i});
25 end
26 end
27
28 function testInputScalarFunction1d(testCase)
29 in = {
30 @(x)1+x*0,
31 @(x)x,
32 @(x)x.*x,
33 };
34
35 out = {
36 [1; 1; 1],
37 [0; 1; 2],
38 [0; 1; 4],
39 };
40
41 g = getTestGrid('1d');
42
43 for i = 1:length(in)
44 gf = grid.evalOn(g,in{i});
45 testCase.verifyEqual(gf, out{i});
46 end
47 end
48
49 function testInputScalarFunction2d(testCase)
50 in = {
51 @(x,y)1+x*0,
52 @(x,y)x-y,
53 @(x,y)x./(1+y),
54 };
55
56 out = {
57 [1; 1; 1; 1; 1; 1; 1; 1; 1],
58 [0; -1; -2; 1; 0; -1; 2; 1; 0],
59 [0; 0; 0; 1; 1/2; 1/3; 2; 1; 2/3],
60 };
61
62 g = getTestGrid('2d');
63
64 for i = 1:length(in)
65 gf = grid.evalOn(g, in{i});
66 testCase.verifyEqual(gf, out{i});
67 end
68 end
69
70
71 function testInputVectorFunction(testCase)
72 g = getTestGrid('1d');
73 in = @(x)[x; -2*x];
74 out = [0; 0; 1; -2; 2; -4];
75
76 gf = grid.evalOn(g,in);
77 testCase.verifyEqual(gf, out);
78
79 g = getTestGrid('2d');
80 in = @(x,y)[x.^2; -2*y];
81 out = [
82 0; 0;
83 0; -2;
84 0; -4;
85 1; 0;
86 1; -2;
87 1; -4;
88 4; 0;
89 4; -2;
90 4; -4;
91 ];
92
93 gf = grid.evalOn(g,in);
94 testCase.verifyEqual(gf, out);
95 end
96
97
98 function testInputErrorVectorValued(testCase)
99 in = {
100 [1,2,3],
101 @(x,y)[x,-y];
102 };
103
104 g = getTestGrid('2d');
105
106 for i = 1:length(in)
107 testCase.verifyError(@()grid.evalOn(g, in{i}),'grid:evalOn:VectorValuedWrongDim',sprintf('in(%d) = %s',i,toString(in{i})));
108 end
109 end
110
111 function g = getTestGrid(d)
112 switch d
113 case '1d'
114 g = grid.equidistant(3,{0,2});
115 case '2d'
116 g = grid.equidistant([3,3],{0,2},{0,2});
117 end
118 end