Mercurial > repos > public > sbplib

comparison +sbp/Upwind.m @ 119:c56437d097de

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
I added remaining upwind operators.
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 15 Dec 2015 16:38:17 +0100
parents 4f5a65f49035
children
comparison
equal deleted inserted replaced
118:5046ff7d13b8 119:c56437d097de
9 end 9 end
10 10
11 methods 11 methods
12 function obj = Upwind(m,h,order) 12 function obj = Upwind(m,h,order)
13 13
14 if order == 3 14 switch order
15 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind3_3bp(m,h); 15 case 2
16 elseif order == 4 16 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind2(m,h);
17 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind4(m,h); 17 case 3
18 elseif order == 6 18 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind3(m,h);
19 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind6(m,h); 19 case 4
20 elseif order == 8 20 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind4(m,h);
21 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind8(m,h); 21 case 5
22 else 22 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind5(m,h);
23 error('Invalid operator order %d.',order); 23 case 6
24 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind6(m,h);
25 case 7
26 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind7(m,h);
27 case 8
28 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind8(m,h);
29 case 9
30 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind9(m,h);
31 otherwise
32 error('Invalid operator order %d.',order);
24 end 33 end
25 34
26 obj.h = h; 35 obj.h = h;
27 obj.m = m; 36 obj.m = m;
28 37