Mercurial > repos > public > sbplib

comparison +scheme/LaplaceCurvilinear.m @ 562:11d8d6ccbcd7 feature/grids/laplace_refactor

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add boundary inner products
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 30 Aug 2017 15:43:51 +0200
parents 3a13916f8ff0
children 8f780429dab2
comparison
equal deleted inserted replaced
561:3a13916f8ff0 562:11d8d6ccbcd7
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 s_w = sqrt((e_w'*x_v).^2 + (e_w'*y_v)');
189 s_e = sqrt((e_e'*x_v).^2 + (e_e'*y_v)');
190 s_s = sqrt((e_s'*x_u).^2 + (e_s'*y_u)');
191 s_n = sqrt((e_n'*x_u).^2 + (e_n'*y_u)');
192
188 obj.d_w = -1*(a11_w*obj.du_w' + a12_w*obj.dv_w')'; 193 obj.d_w = -1*(a11_w*obj.du_w' + a12_w*obj.dv_w')';
189 obj.d_e = (a11_e*obj.du_e' + a12_e*obj.dv_e')'; 194 obj.d_e = (a11_e*obj.du_e' + a12_e*obj.dv_e')';
190 obj.d_s = -1*(a22_s*obj.dv_s' + a12_s*obj.du_s')'; 195 obj.d_s = -1*(a22_s*obj.dv_s' + a12_s*obj.du_s')';
191 obj.d_n = (a22_n*obj.dv_n' + a12_n*obj.du_n')'; 196 obj.d_n = (a22_n*obj.dv_n' + a12_n*obj.du_n')';
192 197
193 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;
194 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;
195 200
201 %% Boundary inner products
202 H_w = H_v*s_w;
203 H_e = H_v*s_e;
204 H_s = H_u*s_s;
205 H_n = H_u*s_n;
196 206
197 % Misc. 207 % Misc.
198 obj.m = m; 208 obj.m = m;
199 obj.h = [h_u h_v]; 209 obj.h = [h_u h_v];
200 obj.order = order; 210 obj.order = order;