Mercurial > repos > public > sbplib

comparison +scheme/LaplaceCurvilinear.m @ 564:9bf49338f8e6 feature/grids/laplace_refactor

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in creation of boundary inner products
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 30 Aug 2017 17:33:41 +0200
parents 8f780429dab2
children f4b0d0e84305
comparison
equal deleted inserted replaced
563:8f780429dab2 564:9bf49338f8e6
197 197
198 obj.Dx = spdiag( y_v./J)*Du + spdiag(-y_u./J)*Dv; 198 obj.Dx = spdiag( y_v./J)*Du + spdiag(-y_u./J)*Dv;
199 obj.Dy = spdiag(-x_v./J)*Du + spdiag( x_u./J)*Dv; 199 obj.Dy = spdiag(-x_v./J)*Du + spdiag( x_u./J)*Dv;
200 200
201 %% Boundary inner products 201 %% Boundary inner products
202 H_w = H_v*s_w; 202 obj.H_w = H_v*spdiag(s_w);
203 H_e = H_v*s_e; 203 obj.H_e = H_v*spdiag(s_e);
204 H_s = H_u*s_s; 204 obj.H_s = H_u*spdiag(s_s);
205 H_n = H_u*s_n; 205 obj.H_n = H_u*spdiag(s_n);
206 206
207 % Misc. 207 % Misc.
208 obj.m = m; 208 obj.m = m;
209 obj.h = [h_u h_v]; 209 obj.h = [h_u h_v];
210 obj.order = order; 210 obj.order = order;