Mercurial > repos > public > sbplib
changeset 1115:0e5dc5dfa559 feature/laplace_curvilinear_test
Temporary ugly fix for reversing direction of boundary in interface coupling
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Mon, 29 Apr 2019 12:55:51 -0700 |
parents | 9ddd2e89d018 |
children | 0aed89043ad6 |
files | +scheme/LaplaceCurvilinear.m |
diffstat | 1 files changed, 7 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- a/+scheme/LaplaceCurvilinear.m Wed Apr 10 13:27:50 2019 -0700 +++ b/+scheme/LaplaceCurvilinear.m Mon Apr 29 12:55:51 2019 -0700 @@ -372,6 +372,13 @@ [th_H_v, ~, th_R_v] = v.getBoundaryBorrowing(neighbour_boundary); m_v = v.getBoundaryNumber(neighbour_boundary); + % BUGFIX?!?!? + if (strcmp(boundary,'s') && strcmp(neighbour_boundary,'e')) || (strcmp(boundary,'e') && strcmp(neighbour_boundary,'s')) + e_v = fliplr(e_v); + d_v = fliplr(d_v); + s_b_v = rot90(s_b_v,2); + end + % Coefficients, v K_v = v.K; J_v = v.J;