Mercurial > repos > public > sbplib

comparison +scheme/Elastic2dVariable.m @ 973:f5a99cca4654 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove unused properties in Elastic2dVariable
author Martin Almquist <malmquist@stanford.edu>
date Tue, 25 Dec 2018 16:26:20 +0100
parents 104f0af001e0
children 1c334842bf23
comparison
equal deleted inserted replaced
972:104f0af001e0 973:f5a99cca4654
30 T_l, T_r 30 T_l, T_r
31 tau_l, tau_r 31 tau_l, tau_r
32 32
33 H, Hi, H_1D % Inner products 33 H, Hi, H_1D % Inner products
34 e_l, e_r 34 e_l, e_r
35 e_w, e_e, e_s, e_n
36 35
37 36
38 d1_l, d1_r % Normal derivatives at the boundary 37 d1_l, d1_r % Normal derivatives at the boundary
39 E % E{i}^T picks out component i 38 E % E{i}^T picks out component i
40 39
41 H_boundary % Boundary inner products 40 H_boundary % Boundary inner products
42 H_w, H_e, H_s, H_n
43 41
44 % Kroneckered norms and coefficients 42 % Kroneckered norms and coefficients
45 RHOi_kron 43 RHOi_kron
46 Hi_kron 44 Hi_kron
47 45
148 % Boundary operators 146 % Boundary operators
149 obj.e_l{1} = kron(e_l{1},I{2}); 147 obj.e_l{1} = kron(e_l{1},I{2});
150 obj.e_l{2} = kron(I{1},e_l{2}); 148 obj.e_l{2} = kron(I{1},e_l{2});
151 obj.e_r{1} = kron(e_r{1},I{2}); 149 obj.e_r{1} = kron(e_r{1},I{2});
152 obj.e_r{2} = kron(I{1},e_r{2}); 150 obj.e_r{2} = kron(I{1},e_r{2});
153 obj.e_w = obj.e_l{1};
154 obj.e_e = obj.e_r{1};
155 obj.e_s = obj.e_l{2};
156 obj.e_n = obj.e_r{2};
157 151
158 obj.d1_l{1} = kron(d1_l{1},I{2}); 152 obj.d1_l{1} = kron(d1_l{1},I{2});
159 obj.d1_l{2} = kron(I{1},d1_l{2}); 153 obj.d1_l{2} = kron(I{1},d1_l{2});
160 obj.d1_r{1} = kron(d1_r{1},I{2}); 154 obj.d1_r{1} = kron(d1_r{1},I{2});
161 obj.d1_r{2} = kron(I{1},d1_r{2}); 155 obj.d1_r{2} = kron(I{1},d1_r{2});
190 obj.Hi = inv(obj.H); 184 obj.Hi = inv(obj.H);
191 obj.H_boundary = cell(dim,1); 185 obj.H_boundary = cell(dim,1);
192 obj.H_boundary{1} = H{2}; 186 obj.H_boundary{1} = H{2};
193 obj.H_boundary{2} = H{1}; 187 obj.H_boundary{2} = H{1};
194 obj.H_1D = {H{1}, H{2}}; 188 obj.H_1D = {H{1}, H{2}};
195 obj.H_w = H{2};
196 obj.H_e = H{2};
197 obj.H_s = H{1};
198 obj.H_n = H{1};
199 189
200 % E{i}^T picks out component i. 190 % E{i}^T picks out component i.
201 E = cell(dim,1); 191 E = cell(dim,1);
202 I = speye(m_tot,m_tot); 192 I = speye(m_tot,m_tot);
203 for i = 1:dim 193 for i = 1:dim