Mercurial > repos > public > sbplib

comparison +scheme/LaplaceCurvilinear.m @ 565:f4b0d0e84305 feature/grids/laplace_refactor

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add missing square
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 30 Aug 2017 17:42:22 +0200
parents 9bf49338f8e6
children 9c98a0526afc
comparison
equal deleted inserted replaced
564:9bf49338f8e6 565:f4b0d0e84305
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)'); 188 s_w = sqrt((e_w'*x_v).^2 + (e_w'*y_v).^2);
189 s_e = sqrt((e_e'*x_v).^2 + (e_e'*y_v)'); 189 s_e = sqrt((e_e'*x_v).^2 + (e_e'*y_v).^2);
190 s_s = sqrt((e_s'*x_u).^2 + (e_s'*y_u)'); 190 s_s = sqrt((e_s'*x_u).^2 + (e_s'*y_u).^2);
191 s_n = sqrt((e_n'*x_u).^2 + (e_n'*y_u)'); 191 s_n = sqrt((e_n'*x_u).^2 + (e_n'*y_u).^2);
192 192
193 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')';
194 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')';
195 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')';
196 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')';