Mercurial > repos > public > sbplib

comparison +scheme/Wave2dCurve.m @ 89:53fe4b64f65e

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Minor cleaning of Wave2dCurve.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 30 Nov 2015 12:20:22 +0100
parents 643bc513b8b8
children 19d0c9325a3e
comparison
equal deleted inserted replaced
85:643bc513b8b8 89:53fe4b64f65e
292 d_n = obj.du_w; 292 d_n = obj.du_w;
293 d_t = obj.dv_w; 293 d_t = obj.dv_w;
294 s = -1; 294 s = -1;
295 295
296 I = gridMatrix(:,1); 296 I = gridMatrix(:,1);
297 coeff_n = obj.a11(:,1); 297 coeff_n = obj.a11(I);
298 coeff_t = obj.a12(:,1); 298 coeff_t = obj.a12(I);
299 case 'e' 299 case 'e'
300 e = obj.e_e; 300 e = obj.e_e;
301 d_n = obj.du_e; 301 d_n = obj.du_e;
302 d_t = obj.dv_e; 302 d_t = obj.dv_e;
303 s = 1; 303 s = 1;
304 304
305 I = gridMatrix(:,end); 305 I = gridMatrix(:,end);
306 coeff_n = obj.a11(:,end); 306 coeff_n = obj.a11(I);
307 coeff_t = obj.a12(:,end); 307 coeff_t = obj.a12(I);
308 case 's' 308 case 's'
309 e = obj.e_s; 309 e = obj.e_s;
310 d_n = obj.dv_s; 310 d_n = obj.dv_s;
311 d_t = obj.du_s; 311 d_t = obj.du_s;
312 s = -1; 312 s = -1;
313 313
314 I = gridMatrix(1,:)'; 314 I = gridMatrix(1,:)';
315 coeff_n = obj.a22(1,:)'; 315 coeff_n = obj.a22(I);
316 coeff_t = obj.a12(1,:)'; 316 coeff_t = obj.a12(I);
317 case 'n' 317 case 'n'
318 e = obj.e_n; 318 e = obj.e_n;
319 d_n = obj.dv_n; 319 d_n = obj.dv_n;
320 d_t = obj.du_n; 320 d_t = obj.du_n;
321 s = 1; 321 s = 1;
322 322
323 I = gridMatrix(end,:)'; 323 I = gridMatrix(end,:)';
324 coeff_n = obj.a22(end,:)'; 324 coeff_n = obj.a22(I);
325 coeff_t = obj.a12(end,:)'; 325 coeff_t = obj.a12(I);
326 otherwise 326 otherwise
327 error('No such boundary: boundary = %s',boundary); 327 error('No such boundary: boundary = %s',boundary);
328 end 328 end
329 329
330 switch boundary 330 switch boundary