Mercurial > repos > public > sbplib
comparison +scheme/LaplaceCurvilinearNew.m @ 1064:1341c6cea9c1 feature/laplace_curvilinear_test
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 | e913cdb34dcb |
| children | 7a55a72729e6 |
comparison
equal
deleted
inserted
replaced
| 1063:e913cdb34dcb | 1064:1341c6cea9c1 |
|---|---|
| 274 a_b = e'*obj.a*e; | 274 a_b = e'*obj.a*e; |
| 275 | 275 |
| 276 switch type | 276 switch type |
| 277 % Dirichlet boundary condition | 277 % Dirichlet boundary condition |
| 278 case {'D','d','dirichlet'} | 278 case {'D','d','dirichlet'} |
| 279 tuning = 1.2; | 279 tuning = 1.0; |
| 280 % tuning = 20.2; | |
| 281 | |
| 282 % b1 = gamm*obj.lambda./obj.a11.^2; | |
| 283 % b2 = gamm*obj.lambda./obj.a22.^2; | |
| 284 | |
| 285 % tau1 = tuning * spdiag(-1./b1 - 1./b2); | |
| 286 % tau2 = 1; | |
| 287 | |
| 288 % tau = (tau1*e + tau2*d)*H_b; | |
| 289 | |
| 290 % closure = obj.a*Hi*tau*e'; | |
| 291 % penalty = -obj.a*Hi*tau; | |
| 292 | 280 |
| 293 sigma = 0; | 281 sigma = 0; |
| 294 for i = 1:obj.dim | 282 for i = 1:obj.dim |
| 295 sigma = sigma + e'*J*K{i,m}*K{i,m}*e; | 283 sigma = sigma + e'*J*K{i,m}*K{i,m}*e; |
| 296 end | 284 end |