Mercurial > repos > public > sbplib

comparison +scheme/Elastic2dCurvilinear.m @ 742:08f3ffe63f48 feature/poroelastic Heimisson2018 elastic1.0

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add metric scale factors to properties in heat and elastic curvilinear
author Martin Almquist <malmquist@stanford.edu>
date Fri, 11 May 2018 13:11:32 -0700
parents f4e2a6a2df08
children 1f6b2fb69225
comparison
equal deleted inserted replaced
741:5a9acf282b34 742:08f3ffe63f48
27 RHO, RHOi % Density, variable 27 RHO, RHOi % Density, variable
28 28
29 % Metric coefficients 29 % Metric coefficients
30 b % Cell matrix of size dim x dim 30 b % Cell matrix of size dim x dim
31 J, Ji 31 J, Ji
32 beta % Cell array of scale factors
32 33
33 D % Total operator 34 D % Total operator
34 D1 % First derivatives 35 D1 % First derivatives
35 36
36 % Second derivatives 37 % Second derivatives
173 end 174 end
174 end 175 end
175 obj.J = J; 176 obj.J = J;
176 obj.Ji = Ji; 177 obj.Ji = Ji;
177 obj.b = b; 178 obj.b = b;
179 obj.beta = beta;
178 %---------------------------- 180 %----------------------------
179 181
180 % Boundary operators 182 % Boundary operators
181 obj.e_l{1} = kron(e_l{1},I{2}); 183 obj.e_l{1} = kron(e_l{1},I{2});
182 obj.e_l{2} = kron(I{1},e_l{2}); 184 obj.e_l{2} = kron(I{1},e_l{2});