Mercurial > repos > public > sbplib
comparison +scheme/LaplaceCurvilinear.m @ 925:8f8f5ff23ead feature/utux2D
Bugfixes in LaplaceCurve
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Sun, 02 Dec 2018 18:04:33 -0800 |
| parents | 1d91c2a8aada |
| children | 4291731570bb |
comparison
equal
deleted
inserted
replaced
| 924:93a5c59dd2fb | 925:8f8f5ff23ead |
|---|---|
| 370 tau_v = tuning * spdiag(tau_v); | 370 tau_v = tuning * spdiag(tau_v); |
| 371 beta_u = tau_v; | 371 beta_u = tau_v; |
| 372 | 372 |
| 373 % Build interpolation operators | 373 % Build interpolation operators |
| 374 [I_u2v_good, I_u2v_bad, I_v2u_good, I_v2u_bad] = ... | 374 [I_u2v_good, I_u2v_bad, I_v2u_good, I_v2u_bad] = ... |
| 375 obj.InterpolationOperators(x_u, x_v, obj.order, opts.interpolation); | 375 obj.interpolationOperators(x_u, x_v, obj.order, opts.interpolation); |
| 376 | 376 |
| 377 closure = a*Hi*e_u*tau_u*H_b_u*e_u' + ... | 377 closure = a*Hi*e_u*tau_u*H_b_u*e_u' + ... |
| 378 a*Hi*e_u*H_b_u*I_v2u_bad*beta_u*I_u2v_good*e_u' + ... | 378 a*Hi*e_u*H_b_u*I_v2u_bad*beta_u*I_u2v_good*e_u' + ... |
| 379 a*1/2*Hi*d_u*H_b_u*e_u' + ... | 379 a*1/2*Hi*d_u*H_b_u*e_u' + ... |
| 380 -a*1/2*Hi*e_u*H_b_u*d_u'; | 380 -a*1/2*Hi*e_u*H_b_u*d_u'; |
| 423 otherwise | 423 otherwise |
| 424 error('No such boundary: boundary = %s',boundary); | 424 error('No such boundary: boundary = %s',boundary); |
| 425 end | 425 end |
| 426 | 426 |
| 427 X = obj.grid.getBoundary(boundary); | 427 X = obj.grid.getBoundary(boundary); |
| 428 if isa(obg.grid, 'grid.Curvilinear')) | 428 if isa(obj.grid, 'grid.Curvilinear') |
| 429 X_logic = obj.grid.logic.getBoundary(boundary); | 429 X_logic = obj.grid.logic.getBoundary(boundary); |
| 430 else | 430 else |
| 431 % Cartesian physical coordinates are also logical coordinates | 431 % Cartesian physical coordinates are also logical coordinates |
| 432 X_logic = X; | 432 X_logic = X; |
| 433 end | 433 end |