Mercurial > repos > public > sbplib

comparison +scheme/LaplaceCurvilinear.m @ 559:e3f587bfd7f8 feature/grids/laplace_refactor

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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
comparison
equal deleted inserted replaced
557:2a856a589510 559:e3f587bfd7f8
183 a22_s = spdiag(a22(I_s)); 183 a22_s = spdiag(a22(I_s));
184 a12_s = spdiag(a12(I_s)); 184 a12_s = spdiag(a12(I_s));
185 a22_n = spdiag(a22(I_n)); 185 a22_n = spdiag(a22(I_n));
186 a12_n = spdiag(a12(I_n)); 186 a12_n = spdiag(a12(I_n));
187 187
188 d_w = -1*(a11_w*obj.du_w' + a12_w*obj.dv_w')'; 188 obj.d_w = -1*(a11_w*obj.du_w' + a12_w*obj.dv_w')';
189 d_e = (a11_e*obj.du_e' + a12_e*obj.dv_e')'; 189 obj.d_e = (a11_e*obj.du_e' + a12_e*obj.dv_e')';
190 d_s = -1*(a22_s*obj.dv_s' + a12_s*obj.du_s')'; 190 obj.d_s = -1*(a22_s*obj.dv_s' + a12_s*obj.du_s')';
191 d_n = (a22_n*obj.dv_n' + a12_n*obj.du_n')'; 191 obj.d_n = (a22_n*obj.dv_n' + a12_n*obj.du_n')';
192 192
193 obj.Dx = spdiag( y_v./J)*Du + spdiag(-y_u./J)*Dv; 193 obj.Dx = spdiag( y_v./J)*Du + spdiag(-y_u./J)*Dv;
194 obj.Dy = spdiag(-x_v./J)*Du + spdiag( x_u./J)*Dv; 194 obj.Dy = spdiag(-x_v./J)*Du + spdiag( x_u./J)*Dv;
195 195
196 196