Mercurial > repos > public > sbplib

annotate +sbp/D1Staggered.m @ 638:17dffd56f2fe feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bug fix and improved comments in D1Staggered wrapper class.
author Martin Almquist <malmquist@stanford.edu>
date Fri, 10 Nov 2017 15:57:21 -0800
parents eead18a8964d
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
1 classdef D1Staggered < sbp.OpSet
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
2 properties
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
3 % x_primal: "primal" grid with m points. Equidistant. Called Plus grid in Ossian's paper.
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
4 % x_dual: "dual" grid with m+1 points. Called Minus grid in Ossian's paper.
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
5
638
17dffd56f2fe Bug fix and improved comments in D1Staggered wrapper class.
Martin Almquist <malmquist@stanford.edu>
parents: 637
diff changeset
6 % D1_primal takes FROM dual grid TO primal grid
17dffd56f2fe Bug fix and improved comments in D1Staggered wrapper class.
Martin Almquist <malmquist@stanford.edu>
parents: 637
diff changeset
7 % D1_dual takes FROM primal grid TO dual grid
17dffd56f2fe Bug fix and improved comments in D1Staggered wrapper class.
Martin Almquist <malmquist@stanford.edu>
parents: 637
diff changeset
8
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
9 D1_primal % SBP operator approximating first derivative
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
10 D1_dual % SBP operator approximating first derivative
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
11 H_primal % Norm matrix
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
12 H_dual % Norm matrix
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
13 H_primalI % H^-1
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
14 H_dualI % H^-1
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
15 e_primal_l % Left boundary operator
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
16 e_dual_l % Left boundary operator
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
17 e_primal_r % Right boundary operator
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
18 e_dual_r % Right boundary operator
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
19 m % Number of grid points.
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
20 m_primal % Number of grid points.
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
21 m_dual % Number of grid points.
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
22 h % Step size
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
23 x_primal % grid
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
24 x_dual % grid
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
25 x
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
26 borrowing % Struct with borrowing limits for different norm matrices
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
27 end
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
28
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
29 methods
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
30 function obj = D1Staggered(m,lim,order)
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
31
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
32 x_l = lim{1};
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
33 x_r = lim{2};
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
34 L = x_r-x_l;
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
35
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
36 m_primal = m;
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
37 m_dual = m+1;
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
38
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
39 switch order
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
40 case 2
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
41 [x_primal, x_dual, Pp, Pm, Qp, Qm] = sbp.implementations.d1_staggered_2(m, L);
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
42 case 4
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
43 [x_primal, x_dual, Pp, Pm, Qp, Qm] = sbp.implementations.d1_staggered_4(m, L);
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
44 case 6
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
45 [x_primal, x_dual, Pp, Pm, Qp, Qm] = sbp.implementations.d1_staggered_6(m, L);
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
46 otherwise
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
47 error('Invalid operator order %d.',order);
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
48 end
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
49
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
50 obj.m = m;
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
51 obj.m_primal = m_primal;
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
52 obj.m_dual = m_dual;
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
53 obj.x_primal = x_l + x_primal';
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
54 obj.x_dual = x_l + x_dual';
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
55
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
56 D1_primal = Pp\Qp;
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
57 D1_dual = Pm\Qm;
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
58
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
59 obj.D1_primal = D1_primal;
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
60 obj.D1_dual = D1_dual;
638
17dffd56f2fe Bug fix and improved comments in D1Staggered wrapper class.
Martin Almquist <malmquist@stanford.edu>
parents: 637
diff changeset
61 obj.H_primal = Pp;
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
62 obj.H_dual = Pm;
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
63
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
64 obj.e_primal_l = sparse(m_primal,1);
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
65 obj.e_primal_r = sparse(m_primal,1);
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
66 obj.e_primal_l(1) = 1;
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
67 obj.e_primal_r(m_primal) = 1;
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
68
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
69 obj.e_dual_l = sparse(m_dual,1);
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
70 obj.e_dual_r = sparse(m_dual,1);
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
71 obj.e_dual_l(1) = 1;
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
72 obj.e_dual_r(m_dual) = 1;
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
73
637
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
74 obj.H_primalI = inv(obj.H_primal);
eead18a8964d Rename +- grids primal and dual grids
Martin Almquist <malmquist@stanford.edu>
parents: 636
diff changeset
75 obj.H_dualI = inv(obj.H_dual);
636
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
76
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
77 obj.borrowing = [];
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
78 obj.x = [];
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
79
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
80 end
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
81
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
82 function str = string(obj)
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
83 str = [class(obj) '_' num2str(obj.order)];
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
84 end
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
85 end
476b7e411ce3 Add staggered operators
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
86 end