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;