Mercurial > repos > public > sbplib

comparison +scheme/Wave2dCurve.m @ 55:a8ed986fcf57

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Minor renaming and clean up in 2d wave schemes.
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 12 Nov 2015 09:47:04 -0800
parents 97a638f91fb8
children df642750e8f7
comparison
equal deleted inserted replaced
54:2194cd385419 55:a8ed986fcf57
232 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary) 232 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
233 % u denotes the solution in the own domain 233 % u denotes the solution in the own domain
234 % v denotes the solution in the neighbour domain 234 % v denotes the solution in the neighbour domain
235 tuning = 1.2; 235 tuning = 1.2;
236 % tuning = 20.2; 236 % tuning = 20.2;
237 [e_u, d_n_u, d_t_u, coeff_n_u, coeff_t_u, s_u, gamm_u, halfnorm_inv_u_n, halfnorm_inv_u_t, halfnorm_u_t, a_n_u, a_t_u] = obj.get_boundary_ops(boundary); 237 [e_u, d_n_u, d_t_u, coeff_n_u, coeff_t_u, s_u, gamm_u, halfnorm_inv_u_n, halfnorm_inv_u_t, halfnorm_u_t] = obj.get_boundary_ops(boundary);
238 [e_v, d_n_v, d_t_v, coeff_n_v, coeff_t_v, s_v, gamm_v, halfnorm_inv_v_n, halfnorm_inv_v_t, halfnorm_v_t, a_n_v, a_t_v] = neighbour_scheme.get_boundary_ops(boundary); 238 [e_v, d_n_v, d_t_v, coeff_n_v, coeff_t_v, s_v, gamm_v, halfnorm_inv_v_n, halfnorm_inv_v_t, halfnorm_v_t] = neighbour_scheme.get_boundary_ops(boundary);
239 239
240 a_n_u = spdiag(coeff_n_u); 240 a_n_u = spdiag(coeff_n_u);
241 a_t_u = spdiag(coeff_t_u); 241 a_t_u = spdiag(coeff_t_u);
242 a_n_v = spdiag(coeff_n_v); 242 a_n_v = spdiag(coeff_n_v);
243 a_t_v = spdiag(coeff_t_v); 243 a_t_v = spdiag(coeff_t_v);
312 case {'w','e'} 312 case {'w','e'}
313 halfnorm_inv_n = obj.Hiu; 313 halfnorm_inv_n = obj.Hiu;
314 halfnorm_inv_t = obj.Hiv; 314 halfnorm_inv_t = obj.Hiv;
315 halfnorm_t = obj.Hv; 315 halfnorm_t = obj.Hv;
316 gamm = obj.gamm_u; 316 gamm = obj.gamm_u;
317 a_n = obj.a11;
318 a_t = obj.a12;
319 case {'s','n'} 317 case {'s','n'}
320 halfnorm_inv_n = obj.Hiv; 318 halfnorm_inv_n = obj.Hiv;
321 halfnorm_inv_t = obj.Hiu; 319 halfnorm_inv_t = obj.Hiu;
322 halfnorm_t = obj.Hu; 320 halfnorm_t = obj.Hu;
323 gamm = obj.gamm_v; 321 gamm = obj.gamm_v;
324 a_n = obj.a22;
325 a_t = obj.a12;
326 end 322 end
327 end 323 end
328 324
329 function N = size(obj) 325 function N = size(obj)
330 N = prod(obj.m); 326 N = prod(obj.m);