Mercurial > repos > public > sbplib

comparison +scheme/Heat2dCurvilinear.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 5a9acf282b34
children 78db023a7fe3
comparison
equal deleted inserted replaced
741:5a9acf282b34 742:08f3ffe63f48
34 H_boundary_l, H_boundary_r 34 H_boundary_l, H_boundary_r
35 35
36 % Metric coefficients 36 % Metric coefficients
37 b % Cell matrix of size dim x dim 37 b % Cell matrix of size dim x dim
38 J, Ji 38 J, Ji
39 beta % Cell array of scale factors
39 40
40 % Numerical boundary flux operators 41 % Numerical boundary flux operators
41 flux_l, flux_r 42 flux_l, flux_r
42 43
43 end 44 end
137 end 138 end
138 end 139 end
139 obj.J = J; 140 obj.J = J;
140 obj.Ji = Ji; 141 obj.Ji = Ji;
141 obj.b = b; 142 obj.b = b;
143 obj.beta = beta;
142 %---------------------------- 144 %----------------------------
143 145
144 % Boundary operators 146 % Boundary operators
145 obj.e_l{1} = kron(e_l{1},I{2}); 147 obj.e_l{1} = kron(e_l{1},I{2});
146 obj.e_l{2} = kron(I{1},e_l{2}); 148 obj.e_l{2} = kron(I{1},e_l{2});