Mercurial > repos > public > sbplib
changeset 559:e3f587bfd7f8 feature/grids/laplace_refactor
Fix bug where boundary operators were not stored in their properties
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Tue, 29 Aug 2017 13:03:45 +0200 |
| parents | 2a856a589510 |
| children | 6132c52bf923 |
| files | +scheme/LaplaceCurvilinear.m |
| diffstat | 1 files changed, 4 insertions(+), 4 deletions(-) [+] |
line wrap: on
line diff
diff -r 2a856a589510 -r e3f587bfd7f8 +scheme/LaplaceCurvilinear.m --- a/+scheme/LaplaceCurvilinear.m Tue Aug 29 12:32:57 2017 +0200 +++ b/+scheme/LaplaceCurvilinear.m Tue Aug 29 13:03:45 2017 +0200 @@ -185,10 +185,10 @@ a22_n = spdiag(a22(I_n)); a12_n = spdiag(a12(I_n)); - d_w = -1*(a11_w*obj.du_w' + a12_w*obj.dv_w')'; - d_e = (a11_e*obj.du_e' + a12_e*obj.dv_e')'; - d_s = -1*(a22_s*obj.dv_s' + a12_s*obj.du_s')'; - d_n = (a22_n*obj.dv_n' + a12_n*obj.du_n')'; + 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')'; + obj.d_s = -1*(a22_s*obj.dv_s' + a12_s*obj.du_s')'; + obj.d_n = (a22_n*obj.dv_n' + a12_n*obj.du_n')'; obj.Dx = spdiag( y_v./J)*Du + spdiag(-y_u./J)*Dv; obj.Dy = spdiag(-x_v./J)*Du + spdiag( x_u./J)*Dv;