Mercurial > repos > public > sbplib
diff +sbp/Upwind.m @ 119:c56437d097de
I added remaining upwind operators.
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Tue, 15 Dec 2015 16:38:17 +0100 |
| parents | 4f5a65f49035 |
| children |
line wrap: on
line diff
--- a/+sbp/Upwind.m Tue Dec 15 16:38:02 2015 +0100 +++ b/+sbp/Upwind.m Tue Dec 15 16:38:17 2015 +0100 @@ -11,16 +11,25 @@ methods function obj = Upwind(m,h,order) - if order == 3 - [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind3_3bp(m,h); - elseif order == 4 - [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind4(m,h); - elseif order == 6 - [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind6(m,h); - elseif order == 8 - [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind8(m,h); - else - error('Invalid operator order %d.',order); + switch order + case 2 + [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind2(m,h); + case 3 + [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind3(m,h); + case 4 + [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind4(m,h); + case 5 + [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind5(m,h); + case 6 + [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind6(m,h); + case 7 + [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind7(m,h); + case 8 + [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind8(m,h); + case 9 + [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind9(m,h); + otherwise + error('Invalid operator order %d.',order); end obj.h = h;