Mercurial > repos > public > sbplib
changeset 1138:afd06a84b69c feature/laplace_curvilinear_test
Bugfixes in VirtaMin.interfaceStandard()
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Mon, 10 Jun 2019 14:39:14 +0200 |
| parents | 2ff1f366e64a |
| children | 6bc93c091682 |
| files | +scheme/LaplaceCurvilinearVirtaMin.m |
| diffstat | 1 files changed, 9 insertions(+), 9 deletions(-) [+] |
line wrap: on
line diff
diff -r 2ff1f366e64a -r afd06a84b69c +scheme/LaplaceCurvilinearVirtaMin.m --- a/+scheme/LaplaceCurvilinearVirtaMin.m Mon Jun 10 13:27:29 2019 +0200 +++ b/+scheme/LaplaceCurvilinearVirtaMin.m Mon Jun 10 14:39:14 2019 +0200 @@ -388,7 +388,7 @@ bp = v.bp; lambda_mat = reshape(lambda_v, my, mx); - switch boundary + switch neighbour_boundary case 'w' lambda_min = min(lambda_mat(:,1:bp), [], 2); lambda_mat = repmat(lambda_min, 1, mx); @@ -406,20 +406,20 @@ switch boundary case {'w', 'e'} - b1_u = gamm_u*u.lambda_min_u(I_u)./u.a11(I_u).^2; - b2_u = h11_u *u.lambda(I_u)./u.a22(I_u).^2; + b1_u = gamm_u*lambda_min_u(I_u)./u.a11(I_u).^2; + b2_u = h11_u *lambda_u(I_u)./u.a22(I_u).^2; case {'s', 'n'} - b1_u = h11_u *u.lambda(I_u)./u.a11(I_u).^2; - b2_u = gamm_u*u.lambda_min(I_u)./u.a22(I_u).^2; + b1_u = h11_u *lambda_u(I_u)./u.a11(I_u).^2; + b2_u = gamm_u*lambda_min_u(I_u)./u.a22(I_u).^2; end switch neighbour_boundary case {'w', 'e'} - b1_v = gamm_v*v.lambda_min_v(I_v)./v.a11(I_v).^2; - b2_v = h11_v *v.lambda(I_v)./v.a22(I_v).^2; + b1_v = gamm_v*lambda_min_v(I_v)./v.a11(I_v).^2; + b2_v = h11_v *lambda_v(I_v)./v.a22(I_v).^2; case {'s', 'n'} - b1_v = h11_v *v.lambda(I_v)./v.a11(I_v).^2; - b2_v = gamm_v*v.lambda_min(I_v)./v.a22(I_v).^2; + b1_v = h11_v *lambda_v(I_v)./v.a11(I_v).^2; + b2_v = gamm_v*lambda_min_v(I_v)./v.a22(I_v).^2; end tau1 = -1./(4*b1_u) -1./(4*b1_v) -1./(4*b2_u) -1./(4*b2_v);