Mercurial > repos > public > sbplib

comparison +scheme/Elastic2dCurvilinear.m @ 1331:60c875c18de3 feature/D2_boundary_opt

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with feature/poroelastic for Elastic schemes
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 10 Mar 2022 16:54:26 +0100
parents cb4bfd0d81ea b8bd54332763
children
comparison
equal deleted inserted replaced
1330:855871e0b852 1331:60c875c18de3
154 154
155 J = x_xi.*y_eta - x_eta.*y_xi; 155 J = x_xi.*y_eta - x_eta.*y_xi;
156 156
157 b = cell(dim,dim); 157 b = cell(dim,dim);
158 b{1,1} = y_eta./J; 158 b{1,1} = y_eta./J;
159 b{1,2} = -x_eta./J; 159 b{2,1} = -x_eta./J;
160 b{2,1} = -y_xi./J; 160 b{1,2} = -y_xi./J;
161 b{2,2} = x_xi./J; 161 b{2,2} = x_xi./J;
162 162
163 % Scale factors for boundary integrals 163 % Scale factors for boundary integrals
164 beta = cell(dim,1); 164 beta = cell(dim,1);
165 beta{1} = sqrt(x_eta.^2 + y_eta.^2); 165 beta{1} = sqrt(x_eta.^2 + y_eta.^2);