Mercurial > repos > public > sbplib

comparison +scheme/LaplaceCurvilinear.m @ 1064:1341c6cea9c1 feature/laplace_curvilinear_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Set tuning to 1.0 in both old and new Dirichlet.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 17 Jan 2019 18:18:00 -0800
parents 78db023a7fe3
children 90c23fd08f59
comparison
equal deleted inserted replaced
1063:e913cdb34dcb 1064:1341c6cea9c1
242 H_b = obj.getBoundaryQuadrature(boundary); 242 H_b = obj.getBoundaryQuadrature(boundary);
243 gamm = obj.getBoundaryBorrowing(boundary); 243 gamm = obj.getBoundaryBorrowing(boundary);
244 switch type 244 switch type
245 % Dirichlet boundary condition 245 % Dirichlet boundary condition
246 case {'D','d','dirichlet'} 246 case {'D','d','dirichlet'}
247 tuning = 1.2; 247 tuning = 1.0;
248 % tuning = 20.2;
249 248
250 b1 = gamm*obj.lambda./obj.a11.^2; 249 b1 = gamm*obj.lambda./obj.a11.^2;
251 b2 = gamm*obj.lambda./obj.a22.^2; 250 b2 = gamm*obj.lambda./obj.a22.^2;
252 251
253 tau1 = tuning * spdiag(-1./b1 - 1./b2); 252 tau1 = tuning * spdiag(-1./b1 - 1./b2);