Mercurial > repos > public > sbplib
changeset 565:f4b0d0e84305 feature/grids/laplace_refactor
Add missing square
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Wed, 30 Aug 2017 17:42:22 +0200 |
| parents | 9bf49338f8e6 |
| children | 9c98a0526afc |
| files | +scheme/LaplaceCurvilinear.m |
| diffstat | 1 files changed, 4 insertions(+), 4 deletions(-) [+] |
line wrap: on
line diff
diff -r 9bf49338f8e6 -r f4b0d0e84305 +scheme/LaplaceCurvilinear.m --- a/+scheme/LaplaceCurvilinear.m Wed Aug 30 17:33:41 2017 +0200 +++ b/+scheme/LaplaceCurvilinear.m Wed Aug 30 17:42:22 2017 +0200 @@ -185,10 +185,10 @@ a22_n = spdiag(a22(I_n)); a12_n = spdiag(a12(I_n)); - s_w = sqrt((e_w'*x_v).^2 + (e_w'*y_v)'); - s_e = sqrt((e_e'*x_v).^2 + (e_e'*y_v)'); - s_s = sqrt((e_s'*x_u).^2 + (e_s'*y_u)'); - s_n = sqrt((e_n'*x_u).^2 + (e_n'*y_u)'); + s_w = sqrt((e_w'*x_v).^2 + (e_w'*y_v).^2); + s_e = sqrt((e_e'*x_v).^2 + (e_e'*y_v).^2); + s_s = sqrt((e_s'*x_u).^2 + (e_s'*y_u).^2); + s_n = sqrt((e_n'*x_u).^2 + (e_n'*y_u).^2); obj.d_w = -1*(a11_w*obj.du_w' + a12_w*obj.dv_w')'; obj.d_e = (a11_e*obj.du_e' + a12_e*obj.dv_e')';