Mercurial > repos > public > sbplib
annotate +sbp/Upwind.m @ 163:51aaf67a7df5 feature/grids
Fixed naming and added functions for converting vectors to matrices for structured grid.
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Mon, 22 Feb 2016 13:55:19 +0100 |
parents | c56437d097de |
children |
rev | line source |
---|---|
52
4f5a65f49035
Fixed the upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
47
diff
changeset
|
1 classdef Upwind < sbp.OpSet |
47 | 2 properties |
3 norms % Struct containing norm matrices such as H,Q, M | |
4 boundary % Struct contanging vectors for boundry point approximations | |
5 derivatives % Struct containging differentiation operators | |
6 borrowing % Struct with borrowing limits for different norm matrices | |
7 m % Number of grid points. | |
8 h % Step size | |
9 end | |
10 | |
11 methods | |
52
4f5a65f49035
Fixed the upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
47
diff
changeset
|
12 function obj = Upwind(m,h,order) |
47 | 13 |
119
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
14 switch order |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
15 case 2 |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
16 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind2(m,h); |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
17 case 3 |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
18 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind3(m,h); |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
19 case 4 |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
20 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind4(m,h); |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
21 case 5 |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
22 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind5(m,h); |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
23 case 6 |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
24 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind6(m,h); |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
25 case 7 |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
26 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind7(m,h); |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
27 case 8 |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
28 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind8(m,h); |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
29 case 9 |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
30 [H, HI, Dp, Dm, e_1, e_m] = sbp.upwind9(m,h); |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
31 otherwise |
c56437d097de
I added remaining upwind operators.
Jonatan Werpers <jonatan@werpers.com>
parents:
52
diff
changeset
|
32 error('Invalid operator order %d.',order); |
47 | 33 end |
34 | |
35 obj.h = h; | |
36 obj.m = m; | |
37 | |
38 obj.norms.H = H; | |
39 obj.norms.HI = HI; | |
40 | |
41 obj.boundary.e_1 = e_1; | |
42 obj.boundary.e_m = e_m; | |
43 | |
44 obj.derivatives.Dp = Dp; | |
45 obj.derivatives.Dm = Dm; | |
46 end | |
47 end | |
48 | |
49 methods (Static) | |
50 function lambda = smallestGrid(obj) | |
51 error('Not implmented') | |
52 end | |
53 end | |
54 end | |
55 | |
56 | |
57 | |
58 | |
59 |