Mercurial > repos > public > sbplib

annotate diracDiscrTest.m @ 757:179f234f6cbf feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with poroelastic to get blockEvalOn.
author Martin Almquist <malmquist@stanford.edu>
date Sun, 17 Jun 2018 14:44:05 -0700
parents 5264ce57b573
children c70131daaa6e
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
647
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
1 function tests = diracDiscrTest()
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
2 tests = functiontests(localfunctions);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
3 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
4
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
5 function testLeftGP(testCase)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
6
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
7 orders = [2, 4, 6];
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
8 mom_conds = orders;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
9
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
10 for o = 1:length(orders)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
11 order = orders(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
12 mom_cond = mom_conds(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
13 [xl, xr, m, h, x, H, fs] = setupStuff(order, mom_cond);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
14
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
15 % Test left boundary grid points
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
16 x0s = xl + [0, h, 2*h];
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
17
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
18 for j = 1:length(fs)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
19 f = fs{j};
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
20 fx = f(x);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
21 for i = 1:length(x0s)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
22 x0 = x0s(i);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
23 delta = diracDiscr(x0, x, mom_cond, 0, H);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
24 integral = delta'*H*fx;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
25 err = abs(integral - f(x0));
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
26 testCase.verifyLessThan(err, 1e-12);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
27 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
28 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
29 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
30 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
31
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
32 function testLeftRandom(testCase)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
33
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
34 orders = [2, 4, 6];
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
35 mom_conds = orders;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
36
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
37 for o = 1:length(orders)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
38 order = orders(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
39 mom_cond = mom_conds(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
40 [xl, xr, m, h, x, H, fs] = setupStuff(order, mom_cond);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
41
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
42 % Test random points near left boundary
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
43 x0s = xl + 2*h*rand(1,10);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
44
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
45 for j = 1:length(fs)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
46 f = fs{j};
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
47 fx = f(x);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
48 for i = 1:length(x0s)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
49 x0 = x0s(i);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
50 delta = diracDiscr(x0, x, mom_cond, 0, H);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
51 integral = delta'*H*fx;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
52 err = abs(integral - f(x0));
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
53 testCase.verifyLessThan(err, 1e-12);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
54 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
55 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
56 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
57 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
58
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
59 function testRightGP(testCase)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
60
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
61 orders = [2, 4, 6];
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
62 mom_conds = orders;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
63
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
64 for o = 1:length(orders)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
65 order = orders(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
66 mom_cond = mom_conds(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
67 [xl, xr, m, h, x, H, fs] = setupStuff(order, mom_cond);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
68
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
69 % Test right boundary grid points
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
70 x0s = xr-[0, h, 2*h];
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
71
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
72 for j = 1:length(fs)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
73 f = fs{j};
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
74 fx = f(x);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
75 for i = 1:length(x0s)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
76 x0 = x0s(i);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
77 delta = diracDiscr(x0, x, mom_cond, 0, H);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
78 integral = delta'*H*fx;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
79 err = abs(integral - f(x0));
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
80 testCase.verifyLessThan(err, 1e-12);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
81 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
82 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
83 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
84 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
85
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
86 function testRightRandom(testCase)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
87
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
88 orders = [2, 4, 6];
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
89 mom_conds = orders;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
90
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
91 for o = 1:length(orders)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
92 order = orders(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
93 mom_cond = mom_conds(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
94 [xl, xr, m, h, x, H, fs] = setupStuff(order, mom_cond);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
95
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
96 % Test random points near right boundary
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
97 x0s = xr - 2*h*rand(1,10);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
98
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
99 for j = 1:length(fs)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
100 f = fs{j};
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
101 fx = f(x);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
102 for i = 1:length(x0s)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
103 x0 = x0s(i);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
104 delta = diracDiscr(x0, x, mom_cond, 0, H);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
105 integral = delta'*H*fx;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
106 err = abs(integral - f(x0));
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
107 testCase.verifyLessThan(err, 1e-12);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
108 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
109 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
110 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
111 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
112
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
113 function testInteriorGP(testCase)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
114
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
115 orders = [2, 4, 6];
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
116 mom_conds = orders;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
117
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
118 for o = 1:length(orders)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
119 order = orders(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
120 mom_cond = mom_conds(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
121 [xl, xr, m, h, x, H, fs] = setupStuff(order, mom_cond);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
122
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
123 % Test interior grid points
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
124 m_half = round(m/2);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
125 x0s = xl + (m_half-1:m_half+1)*h;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
126
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
127 for j = 1:length(fs)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
128 f = fs{j};
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
129 fx = f(x);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
130 for i = 1:length(x0s)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
131 x0 = x0s(i);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
132 delta = diracDiscr(x0, x, mom_cond, 0, H);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
133 integral = delta'*H*fx;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
134 err = abs(integral - f(x0));
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
135 testCase.verifyLessThan(err, 1e-12);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
136 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
137 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
138 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
139 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
140
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
141 function testInteriorRandom(testCase)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
142
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
143 orders = [2, 4, 6];
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
144 mom_conds = orders;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
145
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
146 for o = 1:length(orders)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
147 order = orders(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
148 mom_cond = mom_conds(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
149 [xl, xr, m, h, x, H, fs] = setupStuff(order, mom_cond);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
150
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
151 % Test random points in interior
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
152 x0s = (xl+2*h) + (xr-xl-4*h)*rand(1,20);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
153
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
154 for j = 1:length(fs)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
155 f = fs{j};
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
156 fx = f(x);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
157 for i = 1:length(x0s)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
158 x0 = x0s(i);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
159 delta = diracDiscr(x0, x, mom_cond, 0, H);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
160 integral = delta'*H*fx;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
161 err = abs(integral - f(x0));
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
162 testCase.verifyLessThan(err, 1e-12);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
163 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
164 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
165 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
166 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
167
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
168 % x0 outside grid should yield 0 integral!
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
169 function testX0OutsideGrid(testCase)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
170
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
171 orders = [2, 4, 6];
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
172 mom_conds = orders;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
173
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
174 for o = 1:length(orders)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
175 order = orders(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
176 mom_cond = mom_conds(o);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
177 [xl, xr, m, h, x, H, fs] = setupStuff(order, mom_cond);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
178
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
179 % Test points outisde grid
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
180 x0s = [xl-1.1*h, xr+1.1*h];
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
181
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
182 for j = 1:length(fs)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
183 f = fs{j};
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
184 fx = f(x);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
185 for i = 1:length(x0s)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
186 x0 = x0s(i);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
187 delta = diracDiscr(x0, x, mom_cond, 0, H);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
188 integral = delta'*H*fx;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
189 err = abs(integral - 0);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
190 testCase.verifyLessThan(err, 1e-12);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
191 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
192 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
193 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
194 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
195
651
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
196 function testAllGP(testCase)
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
197
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
198 orders = [2, 4, 6];
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
199 mom_conds = orders;
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
200
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
201 for o = 1:length(orders)
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
202 order = orders(o);
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
203 mom_cond = mom_conds(o);
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
204 [xl, xr, m, h, x, H, fs] = setupStuff(order, mom_cond);
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
205
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
206 % Test all grid points
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
207 x0s = x;
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
208
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
209 for j = 1:length(fs)
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
210 f = fs{j};
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
211 fx = f(x);
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
212 for i = 1:length(x0s)
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
213 x0 = x0s(i);
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
214 delta = diracDiscr(x0, x, mom_cond, 0, H);
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
215 integral = delta'*H*fx;
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
216 err = abs(integral - f(x0));
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
217 testCase.verifyLessThan(err, 1e-12);
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
218 end
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
219 end
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
220 end
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
221 end
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
222
754
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
223 function testHalfGP(testCase)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
224
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
225 orders = [2, 4, 6];
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
226 mom_conds = orders;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
227
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
228 for o = 1:length(orders)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
229 order = orders(o);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
230 mom_cond = mom_conds(o);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
231 [xl, xr, m, h, x, H, fs] = setupStuff(order, mom_cond);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
232
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
233 % Test halfway between all grid points
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
234 x0s = 1/2*( x(2:end)+x(1:end-1) );
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
235
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
236 for j = 1:length(fs)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
237 f = fs{j};
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
238 fx = f(x);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
239 for i = 1:length(x0s)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
240 x0 = x0s(i);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
241 delta = diracDiscr(x0, x, mom_cond, 0, H);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
242 integral = delta'*H*fx;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
243 err = abs(integral - f(x0));
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
244 testCase.verifyLessThan(err, 1e-12);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
245 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
246 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
247 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
248 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
249
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
250 function testAllGPStaggered(testCase)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
251
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
252 orders = [2, 4, 6];
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
253 mom_conds = orders;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
254
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
255 for o = 1:length(orders)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
256 order = orders(o);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
257 mom_cond = mom_conds(o);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
258 [xl, xr, m, h, x, H, fs] = setupStaggered(order, mom_cond);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
259
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
260 % Test all grid points
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
261 x0s = x;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
262
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
263 for j = 1:length(fs)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
264 f = fs{j};
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
265 fx = f(x);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
266 for i = 1:length(x0s)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
267 x0 = x0s(i);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
268 delta = diracDiscr(x0, x, mom_cond, 0, H);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
269 integral = delta'*H*fx;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
270 err = abs(integral - f(x0));
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
271 testCase.verifyLessThan(err, 1e-12);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
272 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
273 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
274 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
275 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
276
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
277 function testHalfGPStaggered(testCase)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
278
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
279 orders = [2, 4, 6];
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
280 mom_conds = orders;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
281
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
282 for o = 1:length(orders)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
283 order = orders(o);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
284 mom_cond = mom_conds(o);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
285 [xl, xr, m, h, x, H, fs] = setupStaggered(order, mom_cond);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
286
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
287 % Test halfway between all grid points
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
288 x0s = 1/2*( x(2:end)+x(1:end-1) );
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
289
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
290 for j = 1:length(fs)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
291 f = fs{j};
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
292 fx = f(x);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
293 for i = 1:length(x0s)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
294 x0 = x0s(i);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
295 delta = diracDiscr(x0, x, mom_cond, 0, H);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
296 integral = delta'*H*fx;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
297 err = abs(integral - f(x0));
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
298 testCase.verifyLessThan(err, 1e-12);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
299 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
300 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
301 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
302 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
303
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
304 function testRandomStaggered(testCase)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
305
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
306 orders = [2, 4, 6];
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
307 mom_conds = orders;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
308
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
309 for o = 1:length(orders)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
310 order = orders(o);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
311 mom_cond = mom_conds(o);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
312 [xl, xr, m, h, x, H, fs] = setupStaggered(order, mom_cond);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
313
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
314 % Test random points within grid boundaries
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
315 x0s = xl + (xr-xl)*rand(1,300);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
316
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
317 for j = 1:length(fs)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
318 f = fs{j};
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
319 fx = f(x);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
320 for i = 1:length(x0s)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
321 x0 = x0s(i);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
322 delta = diracDiscr(x0, x, mom_cond, 0, H);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
323 integral = delta'*H*fx;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
324 err = abs(integral - f(x0));
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
325 testCase.verifyLessThan(err, 1e-12);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
326 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
327 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
328 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
329 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
330
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
331
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
332 % ============== Setup functions =======================
647
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
333 function [xl, xr, m, h, x, H, fs] = setupStuff(order, mom_cond)
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
334
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
335 % Grid
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
336 xl = -3;
754
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
337 xr = 900;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
338 L = xr-xl;
651
4ee7d15bd8e6 Fix roundoff bug in diracDiscr and improve test.
Martin Almquist <malmquist@stanford.edu>
parents: 647
diff changeset
339 m = 101;
647
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
340 h = (xr-xl)/(m-1);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
341 g = grid.equidistant(m, {xl, xr});
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
342 x = g.points();
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
343
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
344 % Quadrature
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
345 ops = sbp.D2Standard(m, {xl, xr}, order);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
346 H = ops.H;
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
347
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
348 % Moment conditions
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
349 fs = cell(mom_cond,1);
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
350 for p = 0:mom_cond-1
754
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
351 fs{p+1} = @(x) (x/L).^p;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
352 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
353
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
354 end
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
355
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
356 function [xl, xr, m, h, x, H, fs] = setupStaggered(order, mom_cond)
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
357
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
358 % Grid
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
359 xl = -3;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
360 xr = 900;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
361 L = xr-xl;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
362 m = 101;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
363 [~, g_dual] = grid.primalDual1D(m, {xl, xr});
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
364 x = g_dual.points();
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
365 h = g_dual.h;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
366
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
367 % Quadrature
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
368 ops = sbp.D1Staggered(m, {xl, xr}, order);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
369 H = ops.H_dual;
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
370
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
371 % Moment conditions
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
372 fs = cell(mom_cond,1);
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
373 for p = 0:mom_cond-1
5264ce57b573 Bugfix diracDiscr to make it work for grids that are non-equidistant near boundaries.
Martin Almquist <malmquist@stanford.edu>
parents: 651
diff changeset
374 fs{p+1} = @(x) (x/L).^p;
647
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
375 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
376
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
377 end
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
378
46c40711830f Add test for diracDiscr.m
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
379